设计作业
matlab和编程是懂一点
可这个实在不会啊
回复 #12 miaomiao 的帖子
还是那个问题,你要分析的信号在哪里? 我看不懂这个程序啊也运行不出来
老师让我编一个程序处理旋转机械的全息谱图
回复 #14 miaomiao 的帖子
你可以这样给x和y幅值:t=1:1000;
x=sin(2*pi*0.02*t);
y=cos(2*pi*0.02*t);
但是你程序中faiP()这个函数或是变量没有定义。
回复 #15 zhlong 的帖子
这个程序不是楼主的,那么里面要用到的一些函数自然也就没有了!所以我一直问这是不是完整的程序,是不是能画出来!看来是不行了回复 #16 zhangnan3509 的帖子
楼主在哪里获得的程序,应该把它弄全啊。强烈建议楼主尽快把faiP()给找出来。 楼主可以看一下这篇文章《基于MATLAB软件平台的全息谱图及在风机振动故障诊断中应用的研究》,里面有程序。 谢谢指点
这个程序就是不完整
faip这个函数我也没有找到
您说的那文章我也看了
太简单
而且没有源数据的时序 画不出图来,对这方面的了解太少,刚开始接触振动和信号处理这方面
知识太肤浅
望多多指教
[ 本帖最后由 zhlong 于 2007-7-9 07:48 编辑 ] 我的问题已经解决谢谢指点
回复 #26 miaomiao 的帖子
能不能具体说说怎么解决的? 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
输入数据就可以了
回复 #29 miaomiao 的帖子
能不能说得详细些 我都不太明白:loveliness:回复 #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 编辑 ]