| 查看: 2610 | 回复: 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级博士申请 PRB*1，四级484
已经有1人回复
-
求助VISSIM破解版软件
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有216人回复
-
求2026年在台湾举行的物理和材料领域国际学术会议信息
已经有0人回复
-
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有14人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
萝花山人
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 303.2
- 沙发: 2
- 帖子: 169
- 在线: 16.7小时
- 虫号: 7186758
- 注册: 2017-09-15
- 性别: GG
- 专业: 机械动力学
2楼2019-06-01 08:41:25