| 查看: 2649 | 回复: 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 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有123人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
萝花山人
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 303.2
- 沙发: 2
- 帖子: 169
- 在线: 16.7小时
- 虫号: 7186758
- 注册: 2017-09-15
- 性别: GG
- 专业: 机械动力学
2楼2019-06-01 08:41:25
