汽车系统动力学Matlab

汽车系统动力学Ｍａ tlab作业报告小组成员:一、组内任务分配二、Ma ｔ lab 程序与图形1、不同转向特性车辆在不同车速下得系统特征根m=１ 0 ０ 0;Ｉ=1500;a1＝１、15;b １=1、3 ５;Ｃ a ｆ＝５ 300 ０;Car=5 ３00 ０;i=1;R=［];for ｕ c＝1 ０:5:1 ００; Ｄ=(I＊(Caf＋C ａ r)+m＊(a1^２*Caf+b1^2*C ａｒ))/(m＊Ｉ＊ｕ c); S=(a1+b1)＾2＊Caf＊Ｃ ar/(m*I＊uc＾2)+(b1*Car-a1＊C ａ f)/I; P=［1 D Ｓ］; ｒ=r ｏ o ｔｓ(P); Ｒ(ｉ,1)=r(1,１);R(i,2)=ｒ(2,1);i=i+1;endplo ｔ(real(Ｒ(:,１)),im ａｇ(Ｒ(:,1)),＇bo＇);hol ｄａ 2＝1。25;ｂ 2=1.２ 5;ｔ=1;S=[］; f ｏ r uc=10:5:100 P＝［m 0;0 I］;Q＝[(Caf+Car)/ｕ c,ｍ*u ｃ+(a2＊Caf—ｂ 2＊C ａｒ)／ｕ c;(a2＊Ｃ af-b ２*Ｃａ r)/ｕｃ,(ａ２^2＊Caf+b2^２＊Ｃ ar)/ｕ c］; R＝［Caf;a ２*Caf］; Ａ=—Ｐ^(-１)*Q; d=eig(A); ｉ=imag(d); ｒ=real(d); S(ｔ,1)=r(1); S(t,2)=i(1); t＝t+１;e ｎ dｐ l ｏ t(S(:,1),S(:,2),'＊＇)ａ 3=１、３ 5;b3＝1。1 ５; fo ｒ uc=10:5:100 P=［m 0;０ Ｉ]; Q＝[(Caf＋Ｃ a ｒ)/uc,m＊u ｃ+(ａ 3*Ｃａ f—ｂ３＊C ａｒ)/uc;(ａ 3＊Ｃ af-b ３*Ca ｒ)/uc,(ａ 3^2*Caf+ｂ 3^2*C ａ r)/uc］; R＝[Caf;a3*Ｃ af］; A=-P^(－1)＊Q; d=ei ｇ(A); i=i ｍａｇ(d); ｒ＝r ｅａ l(ｄ); Ｓ(t,1)=r(1); S(t,2)=i(１); t=t+1;ｅ n ｄｇ r ｉｄ o ｎplo ｔ(Ｓ(:,１),Ｓ(:,２),'d＇);axis(［-14 2 0 3］);ｘ label(’实轴(Re)＇);ｙ label(＇虚轴(Ｉｍ)＇);ｔ ext(-８,２.8,'不足转向');tex ｔ(0,0.2,’过多转向');text(－3,0、2,'中性转向')s ｅ t(gca,'Ｆ o ｎ tName’,'Ｈ elve ｔ ic ａ’,’FontSize’,10)title(［’不同转向特性车辆在不同车速下得系统特征根'］,’FontS ｉze’,12); ２.1、具有不同转向特性车辆得横摆角速度幅频与相频响应m=1 ０ 00;I=1500;a1=1。１ 5;a2＝1。25;a ３=1.3 ５;ｂ１=1、35;ｂ２=1、２ 5;b3=1、15;Ｃ a ｆ=５３００ 0;C ａｒ=5 ３ 000;Ｄ=［］;C＝[］;M=[］;uc＝20;i=1;ｆｏｒ f=0。１:0。1:10 w=2＊p ｉ＊f; E=[1 0;0 1］; P=[ｍ 0;0 I］; Ｒ 1=[Caf;ａ１*Ca ｆ]; R2=[Caf;ａ２*Caf]; R3＝[Ca ｆ;ａ 3＊C ａ f]; Q １=[(Ｃ a ｆ＋C ａ r)/ｕ c m＊...