Đang tải... (xem toàn văn)
CƠ BẢN ĐỂ LẬP TRÌNH MATLAP A. THIẾT LẬP 1. Desktop –desktop layout –defaut : giao diện mặc định 2. Chỉnh cùng 1 nơi lưu file 3. File – new – script : mở cửa sổ soạn thảo 4. Cấu trúc mặc định Clc : xóa dòng lệnh Clear all : xóa giá trị biến, hàm Clf ( tùy trường hợp): xóa hình, đồ thị hiện tại B. CÁC TOÁN TỬ MATLAP STT TOÁN TỬ MÔ TẢ 1 Chia trái 2 Chia phải 3 sqrt(a) Căn bậc 2 của a 4 exp(a) ea 5 log(a) Logarit cơ số e của a 6 log10(a) Logarit cơ số 10 của a 7 abs(a) Trị tuyệt đối a 8 sin(x) Sin(x) 9 asin(x) Arcsin(x) 10 fix(a) Làm tròn về phía 0 11 floor(a) Làm tròn về phía âm vô cùng 12 ceil(a) Làm tròn về phía dương vô cùng 13 round(xn) Chia lấy phần dư (ab) 14 mod(a,b) Tính xn sau đó kết quả làm tròn lên 15 == bằng 16 ~= khác
!" #$%&'#$%&(%'#)%*+'+, /0 123+(4)+(# 5 6+(# #7$8+&%*-9:$; $<%= > ?%8@ /0 xóa dòng lệnh : xóa giá trị biến, hàm !"#$xóa hình, đồ thị hiện tại ABAC " D BAC EF % & !'(' ) * !'#!+' , $ /01)2 3 4#$ 5 6 7 $ 7 '8-92 : 7 %;$ 7 '8-9%;2 < 0-$ =>?@9' A -'4$ B'4$ C -'4$ -'4$ %; '4$ DEFGH#!I; %% 77$ DEFGH#!IJEGK %) '$ DEFGH#!IL8 GK %, 7>L4*$ !'M#!NL*0$ %3 E7LO0$ I!4*->@PQR.>+DEFS %6 TT 0U %: VT Q!( % %< WT X!8!7Y0U %A Z [D %C \ 7Y ); V ]!K )% ^ !_!I! )) RL 7 7``$ [ak=zeros(3,3); for i=1:3 for j=1:3 if i==j k(i,j)=1; elseif i>j k(i,j)=2; else k(i,j)=3; end end end AG HIF Bbc D ,%(J K%= % d$ M#!N!e2d ) 'E d$ M#!N+72d , 7fd$ M'S!"#2d 3 0-d$ ME7@>2d 6 d$ M >E2d gh % Ti%),j), 3k l7E1)!D ,m ) Ti,j)j6k l7E1,!D %m , '-#O0$ l7EmG7!D nE%;;@'oE(!@H>!>O@'oE @N>O@'oE>9'0 3 $ l7E1@8G= 6 EO$ l7E1@8G=Epm T%O T;OR>'Q!(f$ **!SE\Eq\m*!D ;0S#!+'*LX'E1@8G= ) : d7-EO$ l7E1;E!D m**->@P@'H(#*rGD7-> < 7-EO$ l7E1%E!D m A L' G$ l7E1!s7GX'(#*rS@ !s7D(#*r 2G7G C E '$ l7E1 t>!'SM# ,$O '(=(#*ru% @R Ov !D 0U v m wxy % 'Of$ !E!'R>#!NrLF 'mf ) 'O$ !E!'R>LF ' , Of$ !E!'R>mf 3 'QO$ !E!'R>uLF '@RLF Q 6 OfQ$ !E!'R>umf@RmQ z{hz| % -'d$ ]I!}E1O+GHE1iEkGX'E!D m ) -'dO%$ B9LF E1O+GH%-9 , -'dO)$ B9mE1O+GH%-9 3 >E$ B9#!Nr2 6 -!#OEO$ !@v'QI!}E1!D!E!D mO@~!7 m**4E!# : '-E#$ ]'oEP#!+'E1• Q!K < Tik l7E1• A 'O$Tik €PLF '2 C Of$Tik €Pmf2 %; $ !>o!D!E101! %% i0kT$ TE101! uO0DG7!•!‚-9(m 80+ %) Q$ l 2E1 %, ƒ R>!e$O'S!"#2R>#!•$ %3 `ƒ !>oG=2E1#!• %6 $ [R2E1Tv (#!„rS@ !s7 %: … I!)E1 %< 5 E†GX'DE1G>K %A v !'?>!'E1 %C … !J GX'u #*r2E1 , ); ‡ m GX'u #*r2E1 )% L$ I!@=!!•2E1G>K )) 'G$ 1 !=!@+72E1 ), & ˆ'+'!?#!8 ‰!4T )3 * €T* €…T!'#!+'$ )6 & €T& …€T!'('$ ): >$ 8-p2Q!K ' !'?E!?4T;O>Oƒƒ$!7!? !'?ELl !Š>‚ )< '$ I!E1E '(LX'uE1 )A '>$ I!E1E '(SuE1 )C ' $ €>M='S 2E1 ,; iakT' $ TaOR>@P''4• !eODE1e '7 T%$ z{hz|[‹h % 7EG$ mLD'2G7G!>„‹>'L2G$ ) !G$ B9#!Nr2G , E4€$ +GH '(=X!M7 G7€ 3 E'€$ +GH '(=!Œ!M7 G7€ 6 L7>OG$ I!GK!X !I!•2>OG+GH%-9$ : 7 >OG$ I!Š>!X 2>OG+GH%G7$ < 7G>OG$ I! m#2>OG!'@!•!J!>$**Q!K L Q!' P!E-9 A '#$G +7(#*r2Gu('- #!+' C '#>LG$ +7(#*r2GuS4>9 LX' ghzŽw• % GLG$ l7E1[LE7LLeSG7G ) !LEL$ l7E1LELM# , #-$ l7E1-M# 3 !'0 1'0 6 !7$ !JI!E1!D!I!)E1!7#!8 #!(# !7-Q : i•kT.$ !JI!E1!D!I!)E1•GD < iykT>$ !JI!E1!D!I!)E1GDy A 77- ‰E !'?E2@!• C #7GE I! '(=@!•ED(0'RDE1 3 %; #7G I! '(=2@!• %% #7' €M#4‚0U @!• z{hz[‘wwyc’$wx-$ % -E- ]!'0(70'R**-E-0j ) >EOL$ !'R1##!J-9>E*L**QT%Oi%),k$DQT%* 45)‡)4‡,$ , >EL’$ I!rOEt>20'o>!•’+GHG7i>ELkGX' >ErOLEt>$ 3 'E'’O4O$ I! 6 77-$ ˆ'+'#!8 ‰!017GX'E1!?-K : #7$ 'R !'?E‰El'#!8 ‰!GX'E1 !'?E < 7-’O-$ l7G7I!(!?-9$-•#4s#!7-9E† '+ELN A -E)#7’$ l7G7u0'o>!•I!(!?-9$ C #7)-EO-$ l70'o>!•uG7!?-9O0'R- %; -7G’$ ˆ'+'#!8 ‰!’T; %% 4#L$ ]!''o0'o>!•T…0$GX'O0D0'o>!• %) 7’O-$ ]!''o0'o>!•’GD-•#4R#!7!'H> '+E2-9 E† %, G#’O$ Mq%!Š-9!1##!JOP!oDE1!7Y0'o>!• %3 #7'L$ I!@l7!DE’GX'E1!?-9 %6 B'E#'’$ 8 '+0'o>!•OE1` LDMNIOPQRS O+=+&43%8TJU clc; clear all; a=input('nhap so a: '); b=input('nhap so b: '); c=input('nhap so c: '); if a==0 disp('ptrinh co nghiem don la:'); x=-c/b; disp('x = ') ; disp(x); else a~=0 delta=b^2-4*a*c; 6 if delta==0 disp('ptrinh co nghiem kep la:'); x=-b/(2*a); disp('x1 = x2 = '); disp(x); elseif delta>0 disp('ptrinh co 2 nghiem phan biet la:'); x1=(-b+sqrt(delta))/(2*a); x2=(-b-sqrt(delta))/(2*a); disp('x1 = '); disp(x1); disp('x2 = '); disp(x2); else disp('ptrinh co 2 nghiem ao la:'); x1=(-b+sqrt(delta))/(2*a); x2=(-b-sqrt(delta))/(2*a); disp('x1 = '); disp(x1); disp('x2 = '); disp(x2); end end VW/%0&43%8TJU clc; clear all; % do thi y=ax^2+bx+c a=input('nhap a:'); b=input('nhap b:'); c=input('nhap c:'); if a==0 disp('ham bac nhat'); x=-10:0.1:10; y=b*x+c; plot(x,y); else disp('pt bac hai'); x=-10:0.1:10; y=a*x.^2+b*x+c; plot(x,y); end ,` ~e clc; clear all; : d = input('nhap diem vao: '); if d>=0 & d<=10 if d<5 disp('hoc sinh yeu') elseif d>=5 & d<6.5 disp('hoc sinh trung binh') elseif d>=6.5 & d<8 disp('hoc sinh kha') elseif d>=8 & d<9 disp('hoc sinh gioi') elseif d>=9 disp('xuat sac') end end 3` !K -9wa clc; clear all; < syms s Kp Kd Ki a; disp('nhap ham truyen Gs'); disp('vi du:Gs=50/(s*(s+5)) ta nhap:'); disp('num2=50*s^0 ; den=s*(s+5)'); num=input('nhap phuong trinh tu='); den=input('nhap phuong trinh mau='); num1=sym2poly(num); den1=sym2poly(den); Kv=input('nhap he so van toc Kv='); e =input('nhap e='); Wn=input('nhap wn='); c=length(den1); h=num1/den1(1,c); Ki=Kv/h; den2 =[1 0]; num3 =conv(den1,den2); num2= [0 Kd Kp Ki]; num4=[num2*num1]; pt1= [num3 + num4]; num5=[1 a]; num6=[1 2*e*Wn Wn*Wn]; num7=[num6 0]; num8=[0 a*num6]; pt2=[num7+num8]; [Kd,Kp,a]=solve(pt1-pt2); num10= [Kd Kp Ki]; den10 = [1 0]; disp('Bo dieu khien PID:'); disp('Gc=Kd*s+Ki/s+Kp'); Gc = vpa((Kp + Kd*s +Ki/s),4) 6` [7 A clc; clear all; syms s k1 k2 k3; % Buoc 1:-Nhap diem cuc a vao,nhap thong so:e,Wn vao % -Nhap ba matran A,B,C vao và vector K =[k1 k2 k3] a=input('nhap diem cuc a = '); Wn=input('nhap gia tri Wn = '); e=input('nhap gia tri e = '); disp('Nhap ma tran A,B,C : '); A=input('nhap ma tran A = '); B=input('nhap ma tran B = '); C=input('nhap ma tran C = '); K=[k1 k2 k3]; C % Buoc 2:-Xet phuong trinh tinh det(s*I-A+B*K)=0 suy ra duoc pt1 % -Lap pt2=(s-a)*(s^2+2*e*Wn*s+Wn^2) tinhdet=det(s*eye(3)-A+B*K); pt1=collect(tinhdet,s); pt2=collect(((s-a)*(s^2+2*e*Wn*s+Wn^2)),s); % Buoc 3:Giai he phuong trình cua cac he so tuong ung cua pt1 và pt2 % tim duoc các an: k1, k2, k3 ham1=collect((pt1-pt2),s); ham2=coeffs(ham1,s); [k1,k2,k3]=solve(ham2); % Buoc 4:Xuat vector hoi tiep trang thai K disp('Vector hoi tiep trang thai : '); K=vpa([k1 k2 k3],4) %;