Chương II : Lý thuyết điều tự động……… Chương III:Vẽ đồ thị……… Phần IV Giải tích mạch điện bằng Simulink………... Kết quả chưong trình trong Commanwindows:... Giải bài tập mạch điện dưới dạn
Trang 1Phục Lục
Chương I : Matlab cơ bản………
Chương II : Lý thuyết điều tự động……… Chương III:Vẽ đồ thị……… Phần IV Giải tích mạch điện bằng Simulink………
Trang 2Chương I : Matlab cơ bản
*2
Trang 3
Chương trình thực hiện trong comman Window:
Bài 3
[1 4]Hy lËp ch−¬ng tr×nh cho bµi to¸n néi suy Lagrange mét chiÒu trong kho¶ng c¸ch
Hy t×m kÕt qu¶ bÊt kú trong kho¶ng ®ã:
Trang 4
f=
0.6700 0.6200 0.5670 0.5120
n=4
nhap gia tri cua xA[1 4]=3
diem A can tim thuoc day ngoai suy
Có n công việc phải làm liên tiếp ,công việc phải làm xong trong thời gian t1
công việc phải hoàn thành trước thời gian c1.nếu công việc I đựoc hoàn thành
tại thời điểm d1 thì nó sẽ bị phạt với hệ số lương phạt như sau:pi(di-ci)nếu
hoàn thành tại di sau ci càng nhiều thì phạt càng nặng.Ngươcj lại nều hoàn
thành trước (di-ci)<0 thì đây chính là tiền thưởng Tìm cách sắp xếp thứ tự sao
cho tổng phạt là nhỏ nhất
Trang 5Kết quả chưong trình trong Commanwindows:
Trang 6KÕt qu¶ ch−¬ng tr×nh trong Commanwindows:
Trang 7Kết quả trong common windown
Function y= tinh luc (theta,T)
Trang 8Giải bài tập mạch điện dưới dạng script file: Cho mạch điện như hình vẽ
Trang 9Kết quả trong commonwindown :
Trang 11
dãy số đó theo thứ tự tăng dần ( kể cả số 0 nhập cuối cùng)
Trang 12
Bài 10:
Lập chương trình giải phương trình bậc 2:
ax2 + bx +c = 0 khi có nghiệm kép, vô nghiệm, hai nghiệm phân biệt
Giải:
Chương trình viết trong M-file như sau
Trang 13
Kết quả chương trình trong Commanwindows:
Bài 11:
Trang 14Vừa gà vừa chó, bó lại cho tròn 36 con, 100 chân chẵn, hỏi có bao nhiêu chó bao nhiêu gà
Bài 12:Chọn một số dương bất kì Nếu số đó là số chẵn thì chia hết cho 2, nếu
số đó là số lẻ thì nhân với 3 rồi cộng với 1 Lặp lại quá trình đó cho đến khi kết
quả là 1 (Dừng chương trình khi nhập số âm hoặc số 0)
Giải:
Chương trình viết trong M-file như sau:
Trang 16Kết quả trong Commanwindows:
*Bài 20:
Cho ma trận sau: A=[3 2/3;2 3/2] hy tÝnh gi¸ trÞ riªng,vÐct¬ riªng
ch−¬ng tr×nh ®−îc thùc hiÖn trªn cöa sæ comman Window
>> A=sym([3 2/3;2/3 2]);
>> [u,v]=eig(A)
Trang 17
Chương trình viết trong M-file
%tinh tich phan
Trang 18Cho sơ đồ mạch điện:
Bằng Việc áp dụng các định luật Kirchoff, Ohm, Faraday, Henry chúng ta có biểu thị sơ đồ mạch điện trên dưới dạng mô hình toán bởi hệ phương trình vi phân
L1 CHOKE HASH-IRON
E
R3 R2
R1
Trang 191 3 2
2 3
2
)(1)(
1
I R R
R t
q C t E R
R
dt
dq
++
=
Với điều kiện đầu: I1(0)=I10; q(0)=q0 (b) và E=sin(t)
Hãy lập chương trình giải phương trình vi phân trên (a),(b) bằng sử dụng phép biến đổi tích phân Laplace
Trang 20Cho một hệ thống có sơ đồ khối như sau:
Hàm Truyền của các khối như sau:
W1=1/(0.1s+1); W2=0.25/(s+0.3);W3=1/0.2s;W4=2;W5=0.5
1.Hãy xây dựng hàm truyền hệ kín
2.Kiểm tra độ ổn định của hệ thống
W4 W5
Trang 21title('bieu do nyquist va duong qua do');
xlabel('truc tung');
ylabel('truc ao');
xlabel('thoi gian [giay]');
ylabel('bieu do');
Trang 22bài 25:
Cho mạch điện như sau;
Với R=1.4 ohm, L=2H, C=0.32F, uc(0)=0.5V, iL(0)=0 Tại i=0 mạch điện được đóng vào điện áp u(t)=1(t) Hãy tính các đáp ứng quá độ uc(0), iL(0) và uRL(t), và
vẽ các đồ thị đó Sử dụng Matlab để tính quá trình quá độ trong mạch điện tuyến tính
ch-¬ng tr×nh ®-îc viÕt trong M-file nh- sau
Trang 23Chương trình thực hiện trong comman Window
Trang 24Với R=10ohm ,L=1h,C=10-3F, uc(0)=0 v,i1=0,Hãy tính các đáp ứng quá độ
uc(0),i1(0), uRL(t), và vẽ đồ thị trong các trương hợp sau:
Trang 26Chương III:Vẽ đồ thị
bài 27:
Vẽ hàm sin(x) và cos(x) và tan (x),cotan(x) trên cùng một đồ thị
Chương trình trong M-file như sau
title('do thi cac ham luong giac sinx,cosx,tanx');
set(gca,'xtick',-pi:pi/2:pi)
set(gca,'Xticklabel',{'-pi','-pi/2','0','pi/2','pi'})
grid on
Chương trình thực hiện trong comman Window
>> bai27
Trang 30set(gca,'Zscale','log');
Trang 32Đề bài 33:
Dùng câu lệnh Interp1 để tính nội suy cho hàm sau:
Y=sin(x), x thuộc khoảng [0:10] Hãy vẽ hàm đó bằng lệnh plot
Chương trình trong M-file nh- sau
x=0:10;
y=sin(x);
xi=0:0.25:10;
yi=interp1(x,y,xi);
plot(x,y,'o',xi,yi)
Chương trình trong common window
>> bai33
Trang 33Đề bài 34:
Cho mạch điên phi tuyến sau:
Với R1= 500 ohm, I(U) = 2*10-4U2A, U0=20V
Hãy tính I1, U1, U2=U
Chương trình viết trong M-file như sau
R=500;
I(U)
R
U1
Trang 36Phần IV Giải tích mạch điện bằng Simulink
Bài:35:
Minh họa chức năng của các khối Product (phép nhân), Gain (khuếch đại tín hiệu đầu vào), Math Function (trong khối này có một lượng khá lớn các hàm toán đã được chuẩn bị sẵn tại ô Function) Hãy mô phỏng phương trình sau:
F(t)=80exp(-1/80t)sin(0.25+pi/3)
Tham số của khối Sine wave
Amplitude:1; Frequency (rad/sec):0.25, Phase:pi/3, Sample Time:0 Bài làm
Trang 37Đề bài 36:
Hãy hiển thị độc lập đồng thời trên tọa độ XY hai tín hiệu sóng hình sin răng cưa
Bài làm:
Trang 39Đề bài 37:
Làm VD sơ đồ và đồ thị mô phỏng của hệ thống điều khiển như hình vẽ
Bài làm
Trang 41Đề bài 40:
Bài làm
Đề bài 41:
Trang 42Bài làm
Trang 43Bài 43:Cho sơ đồ như hình vẽ
Bài làm
Đề bài 44:
Tuyến tính hóa hệ thống:
Bài làm