求助全息谱，瀑布图

qijunshuai

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([0 200 -40 40]);% 画出全息谱图
hold on
end
这是全息谱的程序，但在图中如何表示初振矢怎么在上面画出来？？
多谢多谢
还有，在启停车中，瀑布图如何画出？用plot3(x,y,z)没做出来。请高手指点一下。多谢

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

qijunshuai

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

qijunshuai

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

eaglezxw

你的瀑布图如何画出？

qijunshuai

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

miaomiao

这个程序有误
麻烦大师帮忙看看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([0 200 -40 40]);% 画出全息谱图
hold on
end

zhlong

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

miaomiao

运行不出图啊

zhangnan3509

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

miaomiao

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

zhangnan3509

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

miaomiao

Undefined function or variable 'x'.

zhlong

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

miaomiao

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

zhlong

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

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