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

qijunshuai 发表于 2007-5-13 21:51

求助全息谱，瀑布图

z=x+y*i; fz=fft(z);
ZReal=real(fz)/512;ZImag=imag(fz)/512; ZR1=ZReal(1:256);ZI1=ZImag(1:256);
for i=1:256
ZR2(i)=ZReal(512-i+1); ZI2(i)=ZImag(512-i+1);
end;
XP=sqrt(ZR1.^2+ZI1.^2); XR=sqrt(ZR2.^2+ZI2.^2);
subplot(3,1,2);
for i=1:8;
t=1:128;
e=16*i;
w=(2*pi)/128;
q=pi/180;
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);
plot(x1,y1,'k'); axis();% 画出全息谱图
hold on
end
这是全息谱的程序，但在图中如何表示初振矢怎么在上面画出来？？
多谢多谢
还有，在启停车中，瀑布图如何画出？用plot3(x,y,z)没做出来。请高手指点一下。多谢

[ 本帖最后由 zhlong 于 2007-7-10 15:21 编辑 ]

qijunshuai 发表于 2007-5-14 10:55

怎么没人给俺回呀?
请高手指点,不吝赐教

qijunshuai 发表于 2007-5-15 09:14

问题已搞定!!看来没压力是不行的.:lol

eaglezxw 发表于 2007-5-15 12:19

你的瀑布图如何画出？

qijunshuai 发表于 2007-5-16 20:21

回复 #4 eaglezxw 的帖子

用函数plot3(x,y,z)
x为转速矩阵，y为频率轴，z为幅值

miaomiao 发表于 2007-7-8 11:41

全息谱程序问题

这个程序有误
麻烦大师帮忙看看z=x+y*i; fz=fft(z);
ZReal=real(fz)/512;ZImag=imag(fz)/512; ZR1=ZReal(1:256);ZI1=ZImag(1:256);
for i=1:256
ZR2(i)=ZReal(512-i+1); ZI2(i)=ZImag(512-i+1);
end;
XP=sqrt(ZR1.^2+ZI1.^2); XR=sqrt(ZR2.^2+ZI2.^2);
subplot(3,1,2);
for i=1:8;
t=1:128;
e=16*i;
w=(2*pi)/128;
q=pi/180;
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);
plot(x1,y1,'k'); axis();% 画出全息谱图
hold on
end

zhlong 发表于 2007-7-8 11:55

有什么样的误？
楼主先介绍一下你对全息谱的理解和你作全息谱的流程吧。

miaomiao 发表于 2007-7-8 16:17

运行不出图啊

zhangnan3509 发表于 2007-7-8 16:26

回复 #3 miaomiao 的帖子

这个程序是您自己编的吗？一开始的x,y指的是什么？他们所代表的数据是什么？这应该不是全部程序吧

miaomiao 发表于 2007-7-8 17:02

不是自己编的才问
自己编的太简单了
找不到全息谱的算法

zhangnan3509 发表于 2007-7-8 17:16

回复 #5 miaomiao 的帖子

也就是说这是您找得别人的程序，请把运行出错的提示语句粘贴上！

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

Undefined function or variable 'x'.

zhlong 发表于 2007-7-8 19:55

回复 #7 miaomiao 的帖子

你都没有给x和y赋值啊，就是说你都没有信号，怎么作谱呢？

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

我怎么给他赋值啊
我只知道全息谱是傅立叶变换的饿改进
只会一点点编程
其他好多不懂
请指教

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

你要作什么信号的全息谱？

实在不行先看点信号处理和matlab的基础书吧。