| 查看: 2928 | 回复: 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 |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有287人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
萝花山人
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 303.2
- 沙发: 2
- 帖子: 169
- 在线: 16.7小时
- 虫号: 7186758
- 注册: 2017-09-15
- 性别: GG
- 专业: 机械动力学
2楼2019-06-01 08:41:25
50