汽车系统动力学Ma tlab作业报告小组成员:一、组内任务分配二、Ma t lab 程序与图形1、不同转向特性车辆在不同车速下得系统特征根m=1 0 0 0;I=1500;a1=1、15;b 1=1、3 5;C a f=5 300 0;Car=5 300 0;i=1;R=[];for u c=1 0:5:1 00; D=(I*(Caf+C a r)+m*(a1^2*Caf+b1^2*C ar))/(m*I*u c); S=(a1+b1)^2*Caf*C ar/(m*I*uc^2)+(b1*Car-a1*C a f)/I; P=[1 D S]; r=r o o ts(P); R(i,1)=r(1,1);R(i,2)=r(2,1);i=i+1;endplo t(real(R(:,1)),im ag(R(:,1)),'bo');hol da 2=1。25;b 2=1.2 5;t=1;S=[]; f o r uc=10:5:100 P=[m 0;0 I];Q=[(Caf+Car)/u c,m*u c+(a2*Caf—b 2*C ar)/u c;(a2*C af-b 2*Ca r)/uc,(a2^2*Caf+b2^2*C ar)/u c]; R=[Caf;a 2*Caf]; A=—P^(-1)*Q; d=eig(A); i=imag(d); r=real(d); S(t,1)=r(1); S(t,2)=i(1); t=t+1;e n dp l o t(S(:,1),S(:,2),'*')a 3=1、3 5;b3=1。1 5; fo r uc=10:5:100 P=[m 0;0 I]; Q=[(Caf+C a r)/uc,m*u c+(a 3*Ca f—b3*C ar)/uc;(a 3*C af-b 3*Ca r)/uc,(a 3^2*Caf+b 3^2*C a r)/uc]; R=[Caf;a3*C af]; A=-P^(-1)*Q; d=ei g(A); i=i mag(d); r=r ea l(d); S(t,1)=r(1); S(t,2)=i(1); t=t+1;e n dg r id o nplo t(S(:,1),S(:,2),'d');axis([-14 2 0 3]);x label(’实轴(Re)');y label('虚轴(Im)');t ext(-8,2.8,'不足转向');tex t(0,0.2,’过多转向');text(-3,0、2,'中性转向')s e t(gca,'F o n tName’,'H elve t ic a’,’FontSize’,10)title([’不同转向特性车辆在不同车速下得系统特征根'],’FontS ize’,12); 2.1、具有不同转向特性车辆得横摆角速度幅频与相频响应m=1 0 00;I=1500;a1=1。1 5;a2=1。25;a 3=1.3 5;b1=1、35;b2=1、2 5;b3=1、15;C a f=5300 0;C ar=5 3 000;D=[];C=[];M=[];uc=20;i=1;for f=0。1:0。1:10 w=2*p i*f; E=[1 0;0 1]; P=[m 0;0 I]; R 1=[Caf;a1*Ca f]; R2=[Caf;a2*Caf]; R3=[Ca f;a 3*C a f]; Q 1=[(C a f+C a r)/u c m*...