帮忙解释下这个小波时频程序的问题吧

tonghuayuxing

fs=250E6;                    %采样频率
     f1=4E6;
     f2=5E6;                 
     t=0:1/fs:1;
     wavename='db4';
     totalscal=256;                    %尺度序列的长度，即scal的长度
     wcf=centfrq(wavename);            %小波的中心频率
     cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
     a=totalscal:-1:1;  
     scal=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列
     coefs=cwt(s,scal,wavename);       %得到小波系数
     f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率
     imagesc(t,f,abs(coefs));          %绘制色谱图
     colorbar;
     xlabel('时间 t/s');
     ylabel('频率 f/Hz');
??? Error using ==> :
Out of memory. Type HELP MEMORY for your options.
而且我不知道加上什么语言可以讲自己的原始信号加进去进行分析

tonghuayuxing

大家帮个忙

zhlong

fs=250E6;                    %采样频率]]
f1=4E6;
f2=5E6;                 
t=0:1/fs:1;

fs太大，导致t太长了。 将fs改小点。

coefs=cwt(s,scal,wavename);  

s就是你要分析的信号。

tonghuayuxing

谢谢了，真热心

对于已经采集到得信号，采样频率是确定的值，进行分析时是不是只能带确定的呀？

zhlong

分析数据少点，不要一次分析那么多

tonghuayuxing

load a;
     s=a;
     plot=(s);
     fs=250E6;                       %采样频率
     t=0:1/fs:0.00000001;
     wavename='db4';
     totalscal=128;                    %尺度序列的长度，即scal的长度
     wcf=centfrq(wavename);            %小波的中心频率
     cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
     a=totalscal:-1:1;  
     scal=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列
     coefs=cwt(s,scal,wavename);       %得到小波系数
     f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率
     imagesc(t,f,abs(coefs));          %绘制色谱图
     colorbar;
     xlabel('时间 t/s');
     ylabel('频率 f/Hz');
我最后按这个程序绘制出了时频图，但是图形有些问题，能不能指导一下？我的图形在附件中，我处理的信号的有效频带在4.15E6到5.85E6之间，最后我想分析在此频带内的信号的特征， 谢谢，麻烦了

tonghuayuxing

希望懂这个程序的给点意见吧