我是个新手,这是我做的一个谐波小波的程序请大家帮我看一下:
function A=hmapdn(f)
% Copyright(c) D.E.Newland 1992
% All right reserved
%
N=length(f);
n=round(log(N)/log(2));
a=hwtdn(f);
%
for j=1:2^(n-2)%为什么进行谐波变换后数据量减少为1/4
A(1,j)=(abs(a(1))^2);%log
A(n+1,j)=(abs(a(2^(n-1)+1))^2);%log
end
for j=2:n
for k=1:2^(j-2)
for m=1:2^(n-j)
A(j,(k-1)*2^(n-j)+m)=(abs(a(2^(j-2)+k))^2+abs(a(N+2-2^(j-2)-k))^2);%log
end
end
end
错误:
??? Input argument 'f' is undefined.
Error in ==> E:\cc\data\复小波与谐波小波\谐波小波\hmapdn.m
On line 5 ==> N=length(f);