|
楼主 |
发表于 2014-7-9 11:30
|
显示全部楼层
Gxx是矩阵,subs是我代入w的值所用的函数,因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过现在又有一个问题呢
x=[ScopeData2.time];
x=x';
y=[ScopeData2.signals.values(:,7)];
y=y';
subplot(2,1,1);
plot(x,y);
xlabel('t/s');
ylabel('加速度/(m/s^2)');
title('Zd')
axis auto;
%去除直流分量
Fs=100; %采样频率
N=length(y);
z=fft(y);
f=(0:N-1)*Fs/N;
mag=2*abs(z)/N; %幅值
Pyy=mag.^2; %功率
subplot(2,1,2);
plot(f(1:N/2),mag(1:N/2)) %频率-幅值图
% plot(f(1:N/2),Pyy(1:N/2)) %功率谱图
xlim([0 30]);
xlabel('f/Hz');
ylabel('幅值/m/s^2');
这个是我simlink仿真的振动响应曲线,用上面的程序画频谱图,,发现采样频率的设置对结果影响很大,这个是什么原因。 |
|