关于小Matlab编程的小波变换

这个例子是书上的Morlet小波程序，运行时总是出错。包括
??? Error using ==> CWT_Morlet at 11
At least 1 parameter required!
这个问题是怎么回事？error的用法错在哪里了？还有，这个程序还有哪些错误呢？比如函数定义之类的问题？请各位大侠帮帮我吧。
function [WT,Freqbins,Scales]=CWT_Morlet(Sig,WinLen,nLevel);
%continuous Wavelet Transform using Morlet function
%Sig:信号
%WinLen:小波函数在尺度参数a=1时的长度（默认为10）
%nLevel:频率轴划分区间段（默认为1024）
%
%WT:返回的小波变换计算结果
%FreqBins:返回频率轴划分结果（归一化频率，最高频率为0.5）
%Scales:返回与频率轴划分值相对应的尺度划分（频率0.5对应的尺度为1）
if (nargin==0),
error('At least 1 parameter required!');
end;

if (nargin<4),
iShow=1;
elseif (nargin<3),
nLevel=1024;
elseif (nargin<2),
WinLen=10;
end
Sig=hilbert(real(Sig));%计算信号的解析信号
SigLen=length(Sig);%获取信号的长度
fmax=0.5;%设置最高分析频率
fmin=0.005;%设置最低分析频率
FreqBins=logspace(log10(fmin),log10(0.5),nLevel);%将频率轴在分析范围内等
%对数坐标划分
Scales=fmax*ones(size(FreqBins))./FreqBins;%计算响应的尺度参数
omg0=WinLen/6;
WT=zero(nLevel,SigLen);%分配计算结果的存储单元
wait=waitbar(0,'Under calculation,please wait');
for m=1:nLevel,
waitbar(m/nLevel,wait);
a=Scales(m);%提取尺度参数
t=-round(a*WinLen):1:round(a*WinLen);
Mor1=pi^(-1/4)*exp(i*2*pi*0.5*t/a).*exp(-t.^2/2/(2*omg0*a)^2);
%计算当前尺度下的小波函数
temp=conv(Sig,Mor1)/sqrt(a); %计算信号与小波函数的卷积
WT(m,:)=temp(round*(a*WinLen)+1:Length(temp)-round(a*WinLen));
end;
close(wait);
WT=WT/WinLen;

呵呵，这个程序是告诉你用morlet小波进行连续小波变换时的程序，相当于是matlab中的底层程序，function [WT,Freqbins,Scales]=CWT_Morlet(Sig,WinLen,nLevel);这个是定义了这样的一个函数，其中Sig,WinLen,nLevel这几个变量为输入变量，分别指信号，信号的长度以及分解的层数。WT,Freqbins,Scales这些为输出变量，分别指小波变换的系数，频率以及尺度。就是对信号Sig进行nLevel层小波变换。在实际使用的时候，只要把输入量给定，然后直接调用这个函数[WT,Freqbins,Scales]=CWT_Morlet(Sig,WinLen,nLevel)。就可以得到输出结果。

这是一个带参数的函数 根据下面的参数校验可以看出你应该是直接点运行了吧 所以报错了 你至少应该产生一个输入信号，比如一个正弦信号来做输入进行变换 。。。谁给你的程序啊 这bug也太多了 很多低级bug 我只是改的能运行了 function [wt,freqbins,scales]=cwt_morlet(sig,winlen,nlevel) %continuous wavelet transform using morlet function %sig:信号 %winlen:小波函数在尺度参数a=1时的长度（默认为10） %nlevel:频率轴划分区间段（默认为1024） % %wt:返回的小波变换计算结果 %freqbins:返回频率轴划分结果（归一化频率，最高频率为0.5） %scales:返回与频率轴划分值相对应的尺度划分（频率0.5对应的尺度为1） if (nargin==0), error('at least 1 parameter required!'); end; if (nargin<3), nlevel=1024; end if (nargin<2), winlen=10; end sig=hilbert(real(sig));%计算信号的解析信号 siglen=length(sig);%获取信号的长度 fmax=0.5;%设置最高分析频率 fmin=0.005;%设置最低分析频率 freqbins=logspace(log10(fmin),log10(0.5),nlevel);%将频率轴在分析范围内等 %对数坐标划分 scales=fmax*ones(size(freqbins))./freqbins;%计算响应的尺度参数 omg0=winlen/6; wt=zeros(nlevel,siglen);%分配计算结果的存储单元 wait=waitbar(0,'under calculation,please wait'); for m=1:nlevel, waitbar(m/nlevel,wait); a=scales(m);%提取尺度参数 t=-round(a*winlen):1:round(a*winlen); mor1=pi^(-1/4)*exp(1i*2*pi*0.5*t/a).*exp(-t.^2/2/(2*omg0*a)^2); %计算当前尺度下的小波函数 temp=conv(sig,mor1)/sqrt(a); %计算信号与小波函数的卷积 wt(m,:)=temp(round(a*winlen)+1:length(temp)-round(a*winlen)); end; close(wait); wt=wt/winlen; %%%%%%运行示例%%%%%%%%% fs=100; t=0:1/fs:2*pi; sig=sin(t); [wt,freqbins,scales]=cwt_morlet(sig);