Tải bản đầy đủ (.doc) (10 trang)

kiến thức CƠ BẢN để lập TRÌNH MATLAP

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (251.9 KB, 10 trang )



  !"
 #$%&'#$%&(%'#)%*+'+, /0
 123+(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
ABAC "
D BAC EF
% & !'('
) * !'#!+'
, $ /01)2
3 4#$ 5
6 7 $ 7 '8-92
: 7 %;$ 7 '8-9%;2
< 0-$ =>?@9'
A -'4$ B'4$
C -'4$ -'4$
%; '4$ DEFGH#!I;
%% 77$ DEFGH#!IJEGK
%) '$ DEFGH#!IL8 GK
%, 7>L4*$ !'M#!NL*0$
%3 E7LO0$ I!4*->@PQR.>+DEFS
%6 TT 0U
%: VT Q!(
%



%< WT X!8!7Y0U
%A Z [D
%C \ 7Y
); V ]!K
)% ^ !_!I!
))  RL 7 7``$
[ak=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
 AG HIF
Bbc
D ,%(J K%=
% d$ M#!N!e2d
) 'E d$ M#!N+72d
, 7fd$ M'S!"#2d
3 0-d$ ME7@>2d
6  d$ M >E2d
gh
% Ti%),j),
3k

l7E1)!D ,m
) Ti,j)j6k l7E1,!D %m
, '-#O0$ l7EmG7!D  nE%;;@'oE(!@H>!>O@'oE
@N>O@'oE>9'0
3 $ l7E1@8G=
6 EO$
l7E1@8G=Epm  T%O T;OR>'Q!(f$
**!SE\Eq\m*!D ;0S#!+'*LX'E1@8G=
)

: d7-EO$ l7E1;E!D m**->@P@'H(#*rGD7->
< 7-EO$ l7E1%E!D m
A L' G$ l7E1!s7GX'(#*rS@ !s7D(#*r
2G7G
C E '$
l7E1 t>!'SM# ,$O '(=(#*ru%
@R Ov !D 0U v m
wxy
% 'Of$ !E!'R>#!NrLF 'mf
) 'O$ !E!'R>LF '
, Of$ !E!'R>mf
3 'QO$ !E!'R>uLF '@RLF Q
6 OfQ$ !E!'R>umf@RmQ
z{hz|
% -'d$ ]I!}E1O+GHE1iEkGX'E!D m
) -'dO%$ B9LF E1O+GH%-9
, -'dO)$ B9mE1O+GH%-9
3 >E$ B9#!Nr2
6 -!#OEO$ !@v'QI!}E1!D!E!D mO@~!7
m**4E!#

: '-E#$ ]'oEP#!+'E1• Q!K
< Tik l7E1•
A 'O$Tik €PLF '2
C Of$Tik €Pmf2
%; $ !>o!D!E101!
%% i0kT$ TE101! uO0DG7!•!‚-9(m
80+
%) Q$ l 2E1
%, ƒ
R>!e$O'S!"#2R>#!•$
%3 `ƒ !>oG=2E1#!•
%6 $ [R2E1Tv (#!„rS@ !s7
%: … I!)E1
%< 5 E†GX'DE1G>K
%A
 
v !'?>!'E1
%C
… !J GX'u #*r2E1
,

);
‡ m  GX'u #*r2E1
)% L$ I!@=!!•2E1G>K 
)) 'G$ 1 !=!@+72E1
), & ˆ'+'!?#!8 ‰!4T
)3 *
€T* €…T!'#!+'$
)6 &
€T& …€T!'('$

): >$ 8-p2Q!K  ' !'?E!?4T;O>Oƒƒ$!7!?
 !'?ELl !Š>‚
)< '$ I!E1E '(LX'uE1
)A '>$ I!E1E '(SuE1
)C ' $ €>M='S 2E1
,; iakT' $
TaOR>@P''4• !eODE1e '7
 T%$
z{hz|[‹h
% 7EG$ mLD'2G7G!>„‹>'L2G$
)  !G$ B9#!Nr2G
, E4€$ +GH '(=X!M7 G7€
3 E'€$ +GH '(=!Œ!M7 G7€
6 L7>OG$ I!GK!X !I!•2>OG+GH%-9$
: 7 >OG$ I!Š>!X 2>OG+GH%G7$
< 7G>OG$ I! m#2>OG!'@!•!J!>$**Q!K L Q!'
P!E-9
A '#$G +7(#*r2Gu('- #!+'
C '#>LG$ +7(#*r2GuS4>9 LX'
ghzŽw•
% GLG$ l7E1[LE7LLeSG7G
) !LEL$ l7E1LELM#
, #-$ l7E1-M#
3 !'0 1'0
6 !7$ !JI!E1!D!I!)E1!7#!8 #!(#
!7-Q
: i•kT.$ !JI!E1!D!I!)E1•GD
< iykT>$ !JI!E1!D!I!)E1GDy
A 77- ‰E !'?E2@!•
C #7GE I! '(=@!•ED(0'RDE1

3

%; #7G I! '(=2@!•
%% #7' €M#4‚0U @!•
z{hz[‘wwyc’$wx-$
% -E- ]!'0(70'R**-E-0j
) >EOL$ !'R1##!J-9>E*L**QT%Oi%),k$DQT%*
45)‡)4‡,$
, >EL’$ I!rOEt>20'o>!•’+GHG7i>ELkGX'
>ErOLEt>$
3 'E'’O4O$
I!
6 77-$ ˆ'+'#!8 ‰!017GX'E1!?-K
: #7$ 'R !'?E‰El'#!8 ‰!GX'E1 !'?E
< 7-’O-$ l7G7I!(!?-9$-•#4s#!7-9E† '+ELN
A -E)#7’$ l7G7u0'o>!•I!(!?-9$
C #7)-EO-$ l70'o>!•uG7!?-9O0'R-
%; -7G’$ ˆ'+'#!8 ‰!’T;
%% 4#L$ ]!''o0'o>!•T…0$GX'O0D0'o>!•
%) 7’O-$ ]!''o0'o>!•’GD-•#4R#!7!'H> '+E2-9
E†
%, G#’O$ Mq%!Š-9!1##!JOP!oDE1!7Y0'o>!•
%3 #7'L$ I!@l7!DE’GX'E1!?-9
%6 B'E#'’$ 8 '+0'o>!•OE1`
 LDMNIOPQRS
 O+=+&43%8TJU
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%8TJU
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 -9wa

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)
%;

×