汽车系统动力学Matlab
作业报告
小组成员：
一、组内任务分配
二、Matlab程序与图形
m=1000;I=1500;a1=;b1=;Caf=53000;Car=53000;
i=1;R=;
for uc=10:5:100;
D=ICaf+Car+ma1^2Caf+b1^2Car/mIuc;
S=a1+b1^2CafCar/mIuc^2+b1Car-a1Caf/I;
P=1 D S;
r=rootsP;
Ri,1=r1,1;Ri,2=r2,1;i=i+1;
end
plotrealR:,1,imagR:,1,'bo';
hold
a2=;
b2=;
t=1;
S=;
for uc=10:5:100
P=m 0;0 I;
Q=Caf+Car/uc,muc+a2Caf-b2Car/uc;a2Caf-b2Car/uc,a2^2Caf+b2^2Car/uc;
R=Caf;a2Caf;
A=-P^-1Q;
d=eigA;
i=imagd;
r=reald;
St,1=r1;
St,2=i1;
t=t+1;
汽车动力系统end
plotS:,1,S:,2,''
a3=;
b3=;
for uc=10:5:100
P=m 0;0 I;
Q=Caf+Car/uc,muc+a3Caf-b3Car/uc;
a3Caf-b3Car/uc,a3^2Caf+b3^2Car/uc;
R=Caf;a3Caf;
A=-P^-1Q;
d=eigA;
i=imagd;
r=reald;
St,1=r1;
St,2=i1;
t=t+1;
end
grid on
plotS:,1,S:,2,'d';
axis-14 2 0 3;
xlabel'实轴Re';
ylabel'虚轴Im';
text-8,,'不足转向';
text0,,'过多转向';
text-3,,'中性转向'
setgca,'FontName','Helvetica','FontSize',10
title'不同转向特性车辆在不同车速下的系统特征根','FontSize',12;、具有不同转向特性车辆的横摆角速度幅频和相频响应
m=1000;I=1500;
a1=;a2=;a3=;
b1=;b2=;b3=;
Caf=53000;Car=53000;
D=;C=;M=;
uc=20;
i=1;
for f=::10
w=2pif;
E=1 0;0 1;
P=m 0;0 I;
R1=Caf;a1Caf;
R2=Caf;a2Caf;
R3=Caf;a3Caf;
Q1=Caf+Car/uc muc+a1Caf-b1Car/uc;a1Caf-b1Car/uc a1^2Caf+b1^2Car/uc;
Q2=Caf+Car/uc muc+a2Caf-b2Car/uc;a2Caf-b2Car/uc a2^2Caf+b2^2Car/uc;
Q3=Caf+Car/uc muc+a3Caf-b3Car/uc;a3Caf-b3Car/uc a3^2Caf+b3^2Car/uc;
A1=-invPQ1;
A2=-invPQ2;
A3=-invPQ3;
B1=invPR1;
B2=invPR2;
B3=invPR3;
Hw1=-invA1-1iwEB1;
Hw2=-invA2-1iwEB2;
Hw3=-invA3-1iwEB3;
Di,1=angleHw12180/pi;Di,2=absHw12;
Ci,1=angleHw22180/pi;Ci,2=absHw22;
Mi,1=angleHw32180/pi;Mi,2=absHw32;
i=i+1;
end
subplot2,1,1
f=::10;
semilogxf,D:,2,'k',f,C:,2,'r-.',f,M:,2
发布评论