giải phương trình bậc 2
ax^2+bx+c=0
th1: a=0
th2:b=0,c<0
th3:c=0
1/a
th4: a khác 0, b khác 0,c khác 0 deltal=b^2-4*a*c
+) deltal >0 có 2 nghiệm phân biệt:
27
_-b†v4 -7+1_
=— =——=
3
th2:
deltal =0 x1=x2=-b/2a
- ly
1⁄a rp x fp 3°
1/2 J th3: delta<0
pt vô nghiệm
code M-file: để tham số các hệ số a,b,c:
a=Input(nhap a);
b=input(‘nhap b');
c=mput(‘nhap c’');
simOut=sim('b | slx.slx');
b. Su dung khéi M-function:
Trong cách tiếp cận nảy, chúng ta sẽ sử dụng khối M-function đề viết mã MATL.AB tùy chỉnh thực hiện phép tính trong phương trình bậc 2. Khôi M-function cho phép chúng ta tích hợp mã MATLAB vào mô hình Simulink.
28
-Hệ số a,b,c nhập từ M-fie trên.
-Code M-function:
function [x1,x2,delta| = fen(a,b,e) delta=(b*2)-(4*a*c);
if (delta>0) x1=(-b-sqrt(delta))/(2*a);
x2=(-b+sqrt(delta))/(2*a);
elseif (delta==0) xI=-b(2*a);
x2=-b/(2*a);
else xl=0;
x2=0;
delta=0;
end
-Từ đó xuất ra kết quả xI,x2,deltal
® Ket qua:
Bang cach str dung ca hai cach tiép cận, chúng ta có thể mô phóng quá trình giải phương trình bậc 2 va kiểm tra tính chính xác của kết quả. So sánh hiệu suất vả tính linh hoạt của hai cách tiếp cận đề thấy rằng việc sử dụng khối M-function có thê là cách tốt hơn trong việc tích hợp mã MATLAB phic tap va tái sử dụng trong các mô hình mô phỏng phức tạp hơn.
Kết luận:
Bải toán mô phỏng giải phương trình bậc 2 băng Simulink đã cho thấy sự mạnh mẽ và linh hoạt của công cụ nảy trong việc mô phỏng và kiếm tra các khía cạnh của các hệ thông. Sự kết hợp giữa khối tính toán cơ bản và khối M-function đã giúp chúng ta thực hiện mô phỏng một cách hiệu quả và tùy chỉnh cho bải toán cụ thé.
Bài 2: Dùng simulink mô phỏng bải toán tính số sinh viên đậu/ rớt. Lưu ý: sử dụng khôi display đề hiện thị sô sv đậu và sô sv rớt, không vẽ biéu do tan suat.
Tiêu Đề: Mô phỏng tính số sinh viên đậu/rớt bằng Simulink
Giới thiệu:
Chung ta sé su dung Simulink dé tạo mô hình và sử dụng khối "Display" để hiện thị sô lượng sinh viên đậu và sô lượng sinh viên rớt.
29
Code M-fil
thiết lập các thông số:
while (true) N=input(‘nhap so hoc sinh la:');
a=Input(so hoe sinh de diem | la:');
b=imput('so hoc sinh de diem 2 la:');
c=mput('so hoc sinh de diem 3 la:');
d=input('so hoc sinh de diem 4 la:');
e=imput('so hoc sinh de diem 5 la:');
f=input('so hoc sinh de diem 6 la:');
g=input('so hoc sinh de diem 7 la:');
jJ=input('so hoc sinh de diem 8 la:');
q=input('so hoc sinh de diem 9 la:');
k=mput('so hoc sinh de diem 10 la:');
if ((atbt+c+d+et+fte+j+qtk)>N)
30
disp(‘moi nhap lai:');
end simOut=sim(‘untitled.slx');
svdau=erfrp+J+q+k;
svrot=a+b+c+d;
end
>> Untitled nhap al nhap b7 nhap c12
>> b2 nhap so hoc so so SƠ SO so so SƠ so so so
hoc hoc hoc hoc hoc hoc hoc hoc hoc hoc
sinh sinh sinh sinh sinh sinh sinh sinh sinh sinh inhap so hoc
sinh la:
dc sin
- code fen function function [dau,rot]
dau=etfte++qtk;
rof=a+b+c+d;
® Ket qua: end
diem dien dien diem dien diem 6 diem diem diem dien h la:
1 1 2 3 4 5 6 7 9 8 10
" fu BPRRPB RHE BBB
la:1
= fcn(a,b,c,d,e,fứ.J.q.k)
Bằng cách sử dụng Simulink và các khối tương ứng, chúng ta có thê thực hiện việc tính toán và hiện thị sô lượng sinh viên đậu và rớt một cách trực quan.
® Kết(luận:
Bải toán mô phỏng tính số sinh viên đậu/rớt bằng Simulink đã thê hiện tính đơn
giản và linh hoạt của công cụ này trong việc tạo các mô hình đơn giản.
Câu 3: Dùng Simulink mô phỏng lại quá trình săn mỗi của chim ưng. Biết phương trình chuyển động của quá trình này được mô tả bằng phương trình
A
=
sau:
- 1
mv 5? GC AY- m:
Trong đó,
A là điện tích cánh bằng 0.0123
Cd là hệ số cản bằng 0.1 ứ là gia tốc trọng trường bằng 9.81 (m/s2)
m là khối lượng bằng | (kg) tho băng 1.225
*Nhận xét kết quả trên khối scopeclc ;
®_ Tiêu Đề: Mô phỏng Quá trình Săn Môi của Chim Ung bang Simulink
® Giới thiệu:
Trong bài này, chúng ta sẽ mô phỏng quá trình săn mỗi của chim ưng bang Simulink. Quá trình này được mô tả bằng phương trình chuyển động cụ thê.
Chung ta sé str dung Simulink dé tao mé hinh va khéi "Scope" dé quan sat va nhan xét két qua.
Mô phỏng quá trình săn mỗi của chim ưng:
:)-[ E3]
ca
Li fa
® Kết Quả:
® Nhận xét kết quả trên khối Scope:
32
Kết quả mô phỏng sẽ hiển thị biểu đồ vị trí và vận tốc của chim ưng theo thời gian. Bằng cách xem biểu đồ trên khối Scope, chúng ta có thế nhận xét về hành vi của chim ưng trong quá trình săn môi. Cụ thế, chúng ta có thể quan sát làm thế nao chim ung tac d6ng vao van tốc và vị trí cla minh dé tiếp can va san mỗi.
® Kết luận:
Bải toán mô phỏng quá trình săn mỗi của chim ưng bằng Simulink đã thê hiện sức mạnh và linh hoạt của công cụ này trong việc mô hình hóa các quá trình vật lý phức tạp.
Câu 4: Thêm điều kiện đóng/ mở cánh của chim ưng vảo câu 3. Biết chim ưng mở rộng cánh khi cảm nhận được gân tiếp cận với mặt đât. Giả sử điêu này xảy ra khi chĩm ưng cách mặt đât 15m. Hành vị đó có thê miêu tả lại như sau:
lfh >= 15:
Cd = Cd1=0.1, A = A1 = 0.0123 Else
Cd = Cd2=0.4,A=A2=0
end
Gia str ban dau chim ung cach mat dat 100m, van téc ban dau -10 (m/s) Nhận xét kết quả trên khôi scope
-_ Mô hinh mô phỏng:
33
100 h 4 od fon AN T
code fen-function:
function [ed,A] = fen(h) if (h>=15)
cd=0.1;
A=0.0123;
else cd=0.4;
A=0;
End - _ Kết quả mô phỏng:
34