| 查看: 2646 | 回复: 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 |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有234人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有29人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有52人回复
萝花山人
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 303.2
- 沙发: 2
- 帖子: 169
- 在线: 16.7小时
- 虫号: 7186758
- 注册: 2017-09-15
- 性别: GG
- 专业: 机械动力学
2楼2019-06-01 08:41:25
