abc.m

汽车振动的频谱分析。很有意思

[fz1,fz4,fz5,fz6,fz7,fz8,fy1,fy4,fy5,fy6,fy7,fy8,fx5]=free12(60,725.63);

%gragh the mounts force(z)
f=0.1:0.1:200;
h1=figure;
loglog(f,fz1(1:2000),'k-',f,fz5(1:2000),'k',f,fz6(1:2000),'k',f,fz7(1:2000),'k',f,fz8(1:2000),'k',f,fz4(1:2000),'k');
grid
xlabel('frequency(Hz)')
ylabel('power density')
title('economy speed 900r/min')
text(1,fz1(10),'---------------fz1')
text(2,fz5(20),'---------------fz5')
text(3,fz6(30),'---------------fz6')
text(4,fz7(40),'---------------fz7')
text(60,fz8(600),'---------------fz8')
text(50,fz4(500),'---------------fz4');

f=0.1:0.1:200;
h2=figure;
loglog(f,fy1(1:2000),'k-',f,fy5(1:2000),'k',f,fy6(1:2000),'k',f,fy7(1:2000),'k',f,fy8(1:2000),'k',f,fy4(1:2000),'k');
grid
xlabel('frequency(Hz)')
ylabel('power density')
title('economy speed 900r/min')
text(1,fy1(10),'---------------fy1')
text(20,fy5(200),'---------------fy5')
text(0.3,fy6(3),'---------------fy6')
text(4,fy7(40),'---------------fy7')
text(60,fy8(600),'---------------fy8')
text(50,fy4(500),'---------------fy4');