汽车理论习题Matlab程序68411 

XX理工大学

for i=1:M fs0=fs(i); lamtas=f/fs0;

deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/mu+1)*lamta.^2).^2;

p_z2=sqrt((1+(2*ypss*lamtas).^2)./((1-lamtas.^2).^2+(2*ypss*lamtas).^2)); z2_q=gama*sqrt((1+4*yps^2*lamta.^2)./deta); p_q=p_z2.*z2_q;

jfg_Gqddf=4*pi^2*sqrt(Gqn0*n0^2*ua)*f; jfg_Gaf=p_q.*jfg_Gqddf; kk=Wf.^2.*jfg_Gaf.^2; aw(i)=sqrt(trapz(f,kk)); end

Law=20*log10(aw/a0); figure(3)

plot(fs,aw);title('aw随fs的变化'),xlabel('“人体—座椅”系统的固有频率fs'),ylabel('aw/m*s^-2'); figure(4)

plot(fs,Law);title('Law随fs的变化'),xlabel('“人体—座椅”系统的固有频率fs'),ylabel('Law/dB');

46 / 56

XX理工大学

（3）问 程序1： clear

47 / 56

XX理工大学

figure(1)

fs=3;yps_s=0.25;g=9.8;

ua=20;Gqn0=2.56*10^(-8);n0=0.1;detaf=0.2;N=180; f0=1.5;yps=0.25;gama=9;mu=10;

ff0=[0.25:0.05:3];sigmaz2=0*ff0;sigmafd=0*ff0;sigmaFd_G=0*ff0; M=length(ff0); for i=1:M f0=ff0(i);

f=detaf*[0:N];lamta=f/f0;lamtas=f/fs;

deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/mu+1)*lamta.^2).^2;

z2_qdot=2*pi*f*gama.*sqrt((1+4*yps^2*lamta.^2)./deta); fd_qdot=gama*lamta.^2./(2*pi*f+eps)./sqrt(deta);

Fd_Gqdot=2*pi*f*gama/g.*sqrt(((lamta.^2/(mu+1)-1).^2+4*yps^2*lamta.^2)./deta); Gq_dotf=4*pi^2*Gqn0*n0^2*ua; Gz2f=(z2_qdot).^2*Gq_dotf; Gfd_qf=(fd_qdot).^2*Gq_dotf; GFd_Gf=(Fd_Gqdot).^2*Gq_dotf; sigmaz2(i)=sqrt(trapz(f,Gz2f)); sigmafd(i)=sqrt(trapz(f,Gfd_qf)); sigmaFd_G(i)=sqrt(trapz(f,GFd_Gf)); if f0==1.5

48 / 56

XX理工大学

sgmz2=sigmaz2(i); sgmfd=sigmafd(i); sgmFd_G=sigmaFd_G(i); end end

sz2=20*log10(sigmaz2/sgmz2); sfd=20*log10(sigmafd/sgmfd); sFd_G=20*log10(sigmaFd_G/sgmFd_G); plot(ff0,sz2,'r-',ff0,sfd,'b-.',ff0,sFd_G,'k--'); axis([0.25 3 -25 15]);

title('三个响应量均方根值随f0变化的曲线'),xlabel('车身部分固有频率f0/Hz'),ylabel('σz2/dB,σfd/dB,σFd/G/dB');

49 / 56

XX理工大学

程序2： clear figure(2)

fs=3;yps_s=0.25;g=9.8;

ua=20;Gqn0=2.56*10^(-8);n0=0.1;detaf=0.2;N=180; f0=1.5;yps=0.25;gama=9;mu=10; c=(0.5-0.125)/180;

yps0=[0.125:c:0.5];sigmaz2=0*yps0;sigmafd=0*yps0;sigmaFd_G=0*yps0; M=length(yps0); for i=1:M yps=yps0(i);

f=detaf*[0:N];lamta=f/f0;lamtas=f/fs;

deta=((1-lamta.^2).*(1+gama-1/mu*lamta.^2)-1).^2+4*yps^2*lamta.^2.*(gama-(1/mu+1)*lamta.^2).^2;

z2_qdot=2*pi*f*gama.*sqrt((1+4*yps^2*lamta.^2)./deta); fd_qdot=gama*lamta.^2./(2*pi*f+eps)./sqrt(deta);

Fd_Gqdot=2*pi*f*gama/g.*sqrt(((lamta.^2/(mu+1)-1).^2+4*yps^2*lamta.^2)./deta); Gq_dotf=4*pi^2*Gqn0*n0^2*ua; Gz2f=(z2_qdot).^2*Gq_dotf; Gfd_qf=(fd_qdot).^2*Gq_dotf; GFd_Gf=(Fd_Gqdot).^2*Gq_dotf; sigmaz2(i)=sqrt(trapz(f,Gz2f));

50 / 56