声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2194|回复: 10

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

[复制链接]
发表于 2014-7-5 15:51 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 gumengyinyin 于 2014-7-5 15:52 编辑

功率谱矩阵中有傅里叶变换的w,不知道怎么画频谱图,求指教,用什么函数
我之前是想对w进行积分,但是算不出来
360截图20140705154909015.jpg
回复
分享到:

使用道具 举报

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

 楼主| 发表于 2014-7-5 21:26 | 显示全部楼层
westrongmc 发表于 2014-7-5 18:53
图片中的公式描述的是加速度频响和位移频响的关系,里面的所谓w,实际上是omega,等于2×PI×freq(i)

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

点评

没看懂你的问题所在。数值小是一回事(数值小不一定相等),曲线是直线是另外一回事(直线代表数值都相等吧?)。  详情 回复 发表于 2014-7-5 22:42
发表于 2014-7-5 22:42 | 显示全部楼层
gumengyinyin 发表于 2014-7-5 21:26
恩恩,我用abs对功率谱密度求模,发现数字都好小,画出的曲线是直线,可不可以指导一下啊

没看懂你的问题所在。数值小是一回事(数值小不一定相等),曲线是直线是另外一回事(直线代表数值都相等吧?)。
 楼主| 发表于 2014-7-7 09:00 | 显示全部楼层
westrongmc 发表于 2014-7-5 22:42
没看懂你的问题所在。数值小是一回事(数值小不一定相等),曲线是直线是另外一回事(直线代表数值都相等 ...

恩,其实实不相等的,我后来修改了下,不过还在跑程序中,到时候如果还出错,求指导啊啊
 楼主| 发表于 2014-7-8 10:58 | 显示全部楼层
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);
我求的功率谱图的纵坐标的值是这样求模得的,是正确的么

点评

没看懂你到底要求什么?Gxx是什么?矩阵?向量?还是单值元素? subs函数用来干什么? Pxx是一个单值元素,plot有意义吗? 建议你step by step调试一下。看每一步中的变量值是否是你想要的。 另外,看一  详情 回复 发表于 2014-7-8 14:03
发表于 2014-7-8 14:03 | 显示全部楼层
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

评分

1

查看全部评分

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

点评

FFT谱分析中的能量泄漏  详情 回复 发表于 2014-7-9 22:58
发表于 2014-7-9 22:58 | 显示全部楼层
gumengyinyin 发表于 2014-7-9 11:30
Gxx是矩阵,subs是我代入w的值所用的函数,因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过 ...

FFT谱分析中的能量泄漏
发表于 2014-7-11 10:36 | 显示全部楼层
gumengyinyin 发表于 2014-7-9 11:30
Gxx是矩阵,subs是我代入w的值所用的函数,因为最开始我是 syms w;
那个功率谱曲线我已经画出来了。
不过 ...

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

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

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-10 23:04 , Processed in 0.066585 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表