汽车系统动力学Matlab.docx

汽车系统动力学Matlab
作业报告
小组成员：
一、组内任务分配
二、Matlab程序与图形
1、不同转向特性车辆在不同车速下的系统特征根
m=1000;I=1500;a1=;b1=;Caf=53000/pi;D(i,2)=abs(Hw1(2));
C(i,1)=angle(Hw2(2))*180/pi;C(i,2)=abs(Hw2(2));
M(i,1)=angle(Hw3(2))*180/pi;M(i,2)=abs(Hw3(2));
i=i+1;
endsubplot(2,1,1)
f=::10;
semilogx(f,D(:,2),'k',f,C(:,2),'r-.',f,M(:,2))
gridonset(gca,'Xtick',[,,,3,10])set(gca,'FontName','Helvetica','FontSize',10)legend(缺乏转向','中性转向'，过多转向')
title]具有不同转向特性车辆的横摆角速度幅频响应],'FontSize',12);xlabel('频率
/Hz')
ylabel(横摆角速度增益/{(°/s)/(。)}')
subplot(2,1,2)
f=::10;
semilogx(f,D(:,1),'k',f,C(:,1),'r-.',f,M(:,1))m=1000;I=1500;
a1=;a2=;a3=;
b1=;b2=;b3=;
Caf=53000;Car=53000;
D=[];C=[];M=[];
uc=20;
i=1;
forf=::10
w=2*pi*f;
E=[10;01];
P=[m0;0I];
R1=[Caf;a1*Caf];
R2=[Caf;a2*Caf];
R3=[Caf;a3*Caf];
Q1=[(Caf+Car)/ucm*uc+(a1*Caf-b1*Car)/uc;(a1*Caf-b1*Car)/uc
(a1A2*Caf+b1A2*Car)/uc];
Q2=[(Caf+Car)/ucm*uc+(a2*Caf-b2*Car)/uc;(a2*Caf-b2*Car)/uc
(a2A2*Caf+b2A2*Car)/uc];
Q3=[(Caf+Car)/ucm*uc+(a3*Caf-b3*Car)/uc;(a3*Caf-b3*Car)/uc(a3A2*Caf+b3A2*Car)/uc];
A1=-inv(P)*Q1;
A2=-inv(P)*Q2;
A3=-inv(P)*Q3;
B1=inv(P)*R1;
B2=inv(P)*R2;
B3=inv(P)*R3;
Hw1=-inv(A1-1i*w*E)*B1;
Hw2=-inv(A2-1i*w*E)*B2;
Hw3=-inv(A3-1i*w*E)*B3;D(i,1)=angle(Hw1(2))*180/pi;D(i,2)=abs(Hw1(2));
C(i,1)=angle(Hw2(2))*180/pi;C(i,2)=abs(Hw2(2));M(i,1)=angle(Hw3(2))*180/pi;M(i,2)=abs(Hw3(2));
i=i+1;
end
subplot(2,1,1)f=::10;
semilogx(f,D(:,2),'k',f,C(:,2),'r-.',f,M(:,2))
gridonset(gca,'Xtick',[,,,3,10])set(gca,'FontName','Helvetica','FontSize',10)legend(缺乏转向','中性转向'，过多转向')
title]具有不同转向特性车辆的横摆角速度幅频响应],'FontSize',12);xlabel('频率/Hz')
ylabel(横摆角速度增益/{(°/s)/(。)}')
subplot(2,1,2)f=::10;
semilogx(f,D(:,1),'k',f,C(:,1),'r-.',f,M(:,1))gridon
set(gca,'Xtick',[,,,3,10])set(gca,'FontName','Helvetica','FontSize',10)legend(缺乏转向','中性转向'，过多转向')
title]具有不同转向特性车辆的横摆角速度相频响应],'FontSize',12);xlabel('频率/Hz')
ylabel(横摆角速度相位〔°〕’)
gridon
set(gca,'Xtick',