gửi bởi tvphuc » 27 Tháng 11 2009 01:02
@van persi: code Matlab de giai bai tap lon
Code Matlab tính nghiệm của bài toán ổn định dọc
disp(' Nhap cac thong so dau vao ');
disp(' ');
W=input('W (lbs)= ');
Ix=input('Ix (slug.ft^2) = ');
Iy=input('Iy (slug.ft^2) = ');
Iz=input('Iz (slug.ft^2) = ');
Ixz=input('Ixz (slug.ft^2) = ');
S=input('S (ft^2) = ');
b=input('b (ft) = ');
c=input('c (ft) = ');
klr=input('Khoi luong rieng (kg/m^3) = ');
u0=input('u0 (m/s) = ');
Cxo=input('Cxo = ');
Cx1=input('Cxalpha = ');
Cxu=input('Cxu = ');
Czo=input('Czo = ');
Cz1=input('Czalpha = ');
Cm1=input('Cmalpha = ');
Cm2=input('Cmalpha.
= ');
Cmu=input('Cmu = ');
Czu=input('Czu = ');
Cz2=input('Czalpha.
= ');
Czq=input('Czq = ');
Cmq=input('Cmq = ');
k1=0.454;k2=4812/3549;k3=0.334;
W=W*k1;
Ix=k2*Ix;Iy=k2*Iy;Iz=k2*Iz;Ixz=k2*Ixz;
S=S*k3^2;b=b*k3;c=c*k3;
L=c/2;
t1=L/u0;
m1=W/(klr*S*L);
ix=Ix/(klr*S*L^3);iy=Iy/(klr*S*L^3);iz=Iz/(klr*S*L^3);ixz=Ixz/(klr*S*L^3);
syms y
A(1,1)=2*m1*y-Cxu;A(1,2)=-Cx1;A(1,3)=-Czo;
A(2,1)=-2*Czo+Czu;A(2,2)=(2*m1-Cz2)*y-Cz1;A(2,3)=-(2*m1+Czq)*y;
A(3,1)=-Cmu;A(3,2)=-(Cm2*y+Cm1);A(3,3)=(iy*y-Cmq)*y;
disp(' Ma tran on dinh doc truc : ');
disp(A);
disp('Dinh thuc cua ma tran:');
dt=det(A)
disp('Cac nghiem :');
nghiem=solve(dt,y);
double(nghiem)
Tiền thì a không thiếu, chủ yếu các chú phải ngoan!