tonghuayuxing 发表于 2009-9-5 12:11

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

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 发表于 2009-9-8 16:43

大家帮个忙

zhlong 发表于 2009-9-8 17:03

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

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

coefs=cwt(s,scal,wavename);

s就是你要分析的信号。

tonghuayuxing 发表于 2009-9-9 15:47

回复 板凳 zhlong 的帖子

本帖最后由 wdhd 于 2016-9-19 13:13 编辑

谢谢了,真热心

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

zhlong 发表于 2009-9-10 00:32

回复 5楼 tonghuayuxing 的帖子

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

tonghuayuxing 发表于 2009-9-17 11:08

视频程序有点问题了,帮个忙吧

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 发表于 2009-9-18 15:06

回复 6楼 tonghuayuxing 的帖子

希望懂这个程序的给点意见吧
页: [1]
查看完整版本: 帮忙解释下这个小波时频程序的问题吧