已经求出功率谱密度函数矩阵如何画频谱图

gumengyinyin

功率谱矩阵中有傅里叶变换的w，不知道怎么画频谱图，求指教，用什么函数
我之前是想对w进行积分，但是算不出来

westrongmc

图片中的公式描述的是加速度频响和位移频响的关系，里面的所谓w，实际上是omega，等于2×PI×freq(i)

gumengyinyin

westrongmc 发表于 2014-7-5 18:53
图片中的公式描述的是加速度频响和位移频响的关系，里面的所谓w，实际上是omega，等于2×PI×freq(i)

恩恩，我用abs对功率谱密度求模，发现数字都好小，画出的曲线是直线，可不可以指导一下啊

westrongmc

gumengyinyin 发表于 2014-7-5 21:26
恩恩，我用abs对功率谱密度求模，发现数字都好小，画出的曲线是直线，可不可以指导一下啊

没看懂你的问题所在。数值小是一回事（数值小不一定相等），曲线是直线是另外一回事（直线代表数值都相等吧？）。

gumengyinyin

westrongmc 发表于 2014-7-5 22:42
没看懂你的问题所在。数值小是一回事（数值小不一定相等），曲线是直线是另外一回事（直线代表数值都相等 ...

恩，其实实不相等的，我后来修改了下，不过还在跑程序中，到时候如果还出错，求指导啊啊

gumengyinyin

westrongmc 发表于 2014-7-5 22:42
没看懂你的问题所在。数值小是一回事（数值小不一定相等），曲线是直线是另外一回事（直线代表数值都相等 ...

Gaa=(2*pi*f)^4*Gxx;
r=0:pi:100*pi;
Pxx=abs(Gaa(1,1));
Pss=subs(Gaa,w,r);
figure(1);
plot(r/(2*pi),Pxx);
我求的功率谱图的纵坐标的值是这样求模得的，是正确的么

westrongmc

gumengyinyin 发表于 2014-7-8 10:58
Gaa=(2*pi*f)^4*Gxx;
r=0:pi:100*pi;
Pxx=abs(Gaa(1,1));

没看懂你到底要求什么？Gxx是什么？矩阵？向量？还是单值元素？

subs函数用来干什么？

Pxx是一个单值元素，plot有意义吗？

建议你step by step调试一下。看每一步中的变量值是否是你想要的。

另外，看一下这个帖子中我的回复【见6楼】，可能有点启发意义。
[size=11.818181991577148px]关于频响函数估计H1,H2,H3,H4的问题
http://forum.vibunion.com/thread-110984-1-1.html

gumengyinyin

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仿真的振动响应曲线，用上面的程序画频谱图，，发现采样频率的设置对结果影响很大，这个是什么原因。

westrongmc

gumengyinyin 发表于 2014-7-9 11:30
Gxx是矩阵，subs是我代入w的值所用的函数，因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过 ...

FFT谱分析中的能量泄漏

ChaChing

gumengyinyin 发表于 2014-7-9 11:30
Gxx是矩阵，subs是我代入w的值所用的函数，因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过 ...

个人水平有限, 实在看不清楚LZ的意思
若数据是simlink仿真的振动响应曲线,采样频率不就确定了吗? 为何还有设置的问题?

gumengyinyin

ChaChing 发表于 2014-7-11 10:36
个人水平有限, 实在看不清楚LZ的意思
若数据是simlink仿真的振动响应曲线,采样频率不就确定了吗? 为何还 ...

恩恩，是的，现在弄明白了。之前都没好好思考采样频率的问题