MATLAB程序

程序如下

t=0:1/pi:100*pi

x=sin(100*pi*t)+sin(200*pi*t)+sin(400*pi*t)

X=fft(x)

f = fftshift(X)

w = linspace(0, 250, length(t))%频率坐标

subplot(211)

plot(t,x)

xlabel('t/s')

ylabel('幅值')

subplot(212)

plot(w,abs(f))%频谱图

title('信号的频谱')

xlabel('频率(Hz)')

ylabel('幅值')

%Matlab仿真程序如下：

r1=36 %单位mm

r2=140

omiga1=60%单位d/sec

x11=1:720

for i=1:720

x1(i)=i*pi/180

%sin(x2(i))=-r1/r2*sin(x1(i))

x2(i)=asin(-r1/r2*sin(x1(i)))

x22(i)=x2(i)*180/pi

r3(i)=r1*cos(x1(i))+r2*cos(x2(i))

B=[-r1*omiga1*sin(x1(i))

r1*omiga1*cos(x1(i))]

A=[r2*sin(x2(i)) 1-r2*cos(x2(i)) 0] X=inv(A)*B

omiga2(i)=X(1,1)

v3(i)=X(2,1)

end

plot(x11/60,0.5*r1*sin(x1))

xlabel('时间 t/sec')

ylabel('连杆质心在Y轴上位置/mm')

figure(2)

plot(x11/60,r3)

xlabel('时间 t/sec')

ylabel('滑块位移r3/mm')

figure(3)

plot(x11/60,omiga2)

xlabel('时间 t/sec')

ylabel('连杆角速度omiga2/rad/sec')

figure(4)

plot(x11/60,v3*pi/180)

xlabel('时间 t/sec')

ylabel('滑块速度v3/mm/sec')

---