| 查看: 2628 | 回复: 1 | |||
[交流]
关于matlab求加速度功率谱
|
|
最近碰到一个频谱分析的问题，要将如图所示的一个a-w振动信号，转换为加速度功率谱密度，以求得与频率相对应的PSD值。 我查了一些资料找到了一个大致的思路，a（w）通过逆傅里叶变换得到a（t），再求得a（t）的自相关函数R(s)，最后对R（s）做傅里叶变换，得到加速度的功率谱密度S(w)。我用Matlab编写了这些程序，最终得到的功率谱密度函数是一个关于w,t,s的函数。因此再求频率对应的功率谱密度值时，得到的是一个函数，而不是数值解，有谁知道问题出在哪里吗？ 下面是matlab的代码： clear all syms w t T s a1=10*exp(i*w*t) a2=(22-2/3)*w*exp(i*w*t) a3=2*exp(i*w*t) a(t)=(int(a1,w,7,18)+int(a2,w,18,30)+int(a3,w,30,50))/(2*pi) t=1:1:450 plot(t,abs(a(t))) syms w t T s m(s)=int(a(t)*a(t+s),t,0,T) R(s)=limit(m(s)/T,T,inf) S(w)=int(R(s)*exp(-i*w*s),s,-inf,inf) z1=S(7) z2=S(18) z3=S(30) z4=S(50)  PSD.PNG |
回复此楼» 猜你喜欢
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有86人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有1人回复
-
投稿chemical physical letters不送审？
已经有2人回复
萝花山人
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 303.2
- 沙发: 2
- 帖子: 169
- 在线: 16.7小时
- 虫号: 7186758
- 注册: 2017-09-15
- 性别: GG
- 专业: 机械动力学
2楼2019-06-01 08:41:25