已经求出功率谱密度函数矩阵如何画频谱图
本帖最后由 gumengyinyin 于 2014-7-5 15:52 编辑功率谱矩阵中有傅里叶变换的w，不知道怎么画频谱图，求指教，用什么函数
我之前是想对w进行积分，但是算不出来
图片中的公式描述的是加速度频响和位移频响的关系，里面的所谓w，实际上是omega，等于2×PI×freq(i)
westrongmc 发表于 2014-7-5 18:53
图片中的公式描述的是加速度频响和位移频响的关系，里面的所谓w，实际上是omega，等于2×PI×freq(i)
恩恩，我用abs对功率谱密度求模，发现数字都好小，画出的曲线是直线，可不可以指导一下啊 gumengyinyin 发表于 2014-7-5 21:26
恩恩，我用abs对功率谱密度求模，发现数字都好小，画出的曲线是直线，可不可以指导一下啊
没看懂你的问题所在。数值小是一回事（数值小不一定相等），曲线是直线是另外一回事（直线代表数值都相等吧？）。
westrongmc 发表于 2014-7-5 22:42
没看懂你的问题所在。数值小是一回事（数值小不一定相等），曲线是直线是另外一回事（直线代表数值都相等 ...
恩，其实实不相等的，我后来修改了下，不过还在跑程序中，到时候如果还出错，求指导啊啊 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);
我求的功率谱图的纵坐标的值是这样求模得的，是正确的么 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楼】，可能有点启发意义。
关于频响函数估计H1,H2,H3,H4的问题
http://forum.vibunion.com/thread-110984-1-1.html
Gxx是矩阵，subs是我代入w的值所用的函数，因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过现在又有一个问题呢
x=;
x=x';
y=;
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();
xlabel('f/Hz');
ylabel('幅值/m/s^2');
这个是我simlink仿真的振动响应曲线，用上面的程序画频谱图，，发现采样频率的设置对结果影响很大，这个是什么原因。 gumengyinyin 发表于 2014-7-9 11:30
Gxx是矩阵，subs是我代入w的值所用的函数，因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过 ...
FFT谱分析中的能量泄漏 gumengyinyin 发表于 2014-7-9 11:30
Gxx是矩阵，subs是我代入w的值所用的函数，因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过 ...
个人水平有限, 实在看不清楚LZ的意思
若数据是simlink仿真的振动响应曲线,采样频率不就确定了吗? 为何还有设置的问题? ChaChing 发表于 2014-7-11 10:36
个人水平有限, 实在看不清楚LZ的意思
若数据是simlink仿真的振动响应曲线,采样频率不就确定了吗? 为何还 ...
恩恩，是的，现在弄明白了。之前都没好好思考采样频率的问题
页:
[1]