振动论坛 Matlab讨论区 紧急求助这段循环错那里.doc

《振动论坛 Matlab讨论区 紧急求助这段循环错那里.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、振动论坛»Matlab讨论区»紧急求助：这段循环错那里了？？查看完整版本:紧急求助：这段循环错那里了？？糖葫芦2006-9-301:34紧急求助：这段循环错那里了？？根据图中二式我编写的程序如下：clearall;clc;T  =1/2400;matchout=[1:2:100];%数据L      =length(matchout);%数据长度Corr_matchout    =xcorr(matchout);%自相关Cut_Corr_matchout=Corr_matchout(L:1:2*L-1);circle_

2、range=[0:100:400];%循环频率a    freq_range  =[0:500:4000];%频率fforff=1:length(freq_range)  forcc=1:length(circle_range)      fortao=0:L-1        forii=1:L          R(tao+1,ii)=Cut_Corr_matchout(tao+1)*exp(-j*2*pi*circle_range(cc)*ii*T);        end;        Rx(cc,tao+1

3、)=sum(R(tao+1,:))/(L*T);%循环子相关函数        SS(cc,tao+1)=Rx(cc,tao+1)*exp(-j*2*pi*freq_range(ff)*tao);      end;      Sx(ff,cc)=abs(sum(SS(cc,:)));%循环谱幅度  end;        end;surf(circle_range,freq_range,Sx);xlabel('循环频率/Hz');ylabel('频率/Hz');zlabel('幅度');最后得到的是随着循环频率a和频

4、率f变化的循环幅度谱3维图....问题是最后得到的循环幅度谱幅度并没有随着f和a的变化而变化，在每一个频率f下得到的都是同一组数值....看了2天了没有找到问题，求高人指点，不胜感激！:'(AaronSpark2006-9-308:07画三位谱图有这么麻烦吗？怎么那么多循环？用FFT不就行了吗？另外一般三位谱图也没有用surf画的页:[1]查看完整版本:紧急求助：这段循环错那里了？？PoweredbyDiscuz!Archiver6.0.0 ©2001-2006ComsenzInc.