信号处理方法 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

zhangnan3509 发表于 2007-7-8 20:12

如果仅仅是引用数据的问题，那么不用从头做起。这就是全部程序吗？是不是以前有人用这个程序画出过全息谱图

miaomiao 发表于 2007-7-8 20:19

就是要做一般转子信号的全息谱分析
设计作业
matlab和编程是懂一点
可这个实在不会啊

zhlong 发表于 2007-7-8 20:21

回复 #12 miaomiao 的帖子

还是那个问题，你要分析的信号在哪里？

miaomiao 发表于 2007-7-8 20:26

我看不懂这个程序啊
也运行不出来
老师让我编一个程序处理旋转机械的全息谱图

zhlong 发表于 2007-7-8 20:35

回复 #14 miaomiao 的帖子

你可以这样给x和y幅值：
t=1:1000;
x=sin(2*pi*0.02*t);
y=cos(2*pi*0.02*t);

但是你程序中faiP()这个函数或是变量没有定义。

zhangnan3509 发表于 2007-7-8 20:40

回复 #15 zhlong 的帖子

这个程序不是楼主的，那么里面要用到的一些函数自然也就没有了！所以我一直问这是不是完整的程序，是不是能画出来！看来是不行了

zhlong 发表于 2007-7-8 20:44

回复 #16 zhangnan3509 的帖子

楼主在哪里获得的程序，应该把它弄全啊。
强烈建议楼主尽快把faiP()给找出来。

zhlong 发表于 2007-7-8 21:02

楼主可以看一下这篇文章《基于MATLAB软件平台的全息谱图及在风机振动故障诊断中应用的研究》，里面有程序。

miaomiao 发表于 2007-7-8 22:18

谢谢指点
这个程序就是不完整
faip这个函数我也没有找到
您说的那文章我也看了
太简单
而且没有源数据的时序

miaomiao 发表于 2007-7-8 22:20

画不出图来，对这方面的了解太少，刚开始接触振动和信号处理这方面
知识太肤浅
望多多指教

[ 本帖最后由 zhlong 于 2007-7-9 07:48 编辑 ]

miaomiao 发表于 2007-7-14 22:01

我的问题已经解决谢谢指点

zhangnan3509 发表于 2007-7-14 22:02

回复 #26 miaomiao 的帖子

能不能具体说说怎么解决的？

miaomiao 发表于 2007-7-15 16:01

for i=1:4;
a=abs(fft(x));
b=abs(fft(y));
px=angle(fft(x));
py=angle(fft(y));
end

for i=1:4
for j=1:128
   vx(j)=a(j)*cos(2*pi*(i-1)/32+px(j));
   vy(j)=b(j)*cos(2*pi*(i-1)/32+py(j));
   z(j)=j*2;
end
plot3(z,vx,vy,'*');
hold on
grid on
输入数据就可以了

zhangnan3509 发表于 2007-7-15 17:09

回复 #29 miaomiao 的帖子

能不能说得详细些我都不太明白:loveliness:

qijunshuai 发表于 2007-7-16 11:22

回复 #23 zhlong 的帖子

fai()这个函数是自己定义的。就是为了求相位角，你直接用Matlab里的angle（）就可以了。
而且这里面的算法利用的是全矢谱的快速算法，其中相位精度没有像全息谱所说的那样进行改进。
而且一开始给出的程序是不完整的。x，y是两个数组，就是你要分析的截面的两通道的数据。

x1=XP(e)*cos(w*t+faiP(e)*q)+XR(e)*cos(w*t+faiR(e)*q)+e; y1=XP(e)*sin(w*t+faiP(e)*q)-XR(e)*sin(w*t+faiR(e)*q);
这一步就是根据钟一鄂编《转子动力学》中图1.4几何关系得到的。

[ 本帖最后由 qijunshuai 于 2007-7-16 11:26 编辑 ]

页: 1 [2] 3