[讨论]关于正交镜像滤波器组的问题

[wanjuzdy]

--------------------------------------------------------------------------------
我要做个最基本正交镜像滤波器组（时域）的，我大概的弄了一个，发现是频域的，而调试出来有不对，能帮我改一下吗？谢谢了！
N=64; %滤波器的长度
div=0.1; %截止频率的初始步长
eval=exp(-j*pi*(0:N-1)/2); %在ω=n*pi/2处的复指数值
freq=0.5; %初始化截止频率
target=1/sqrt(2); %在ω=n*pi/2处,希望得到的幅值
for i=1:1000 %绘制两个滤波器频率响应模的平方和
w(i)=i*pi/1000;
hw1(i)=h*exp(-j*w(i)*[0:255]');
hw2(i)=h*exp(-j*(w(i)-pi)*[0:255]');
yy(i)=abs(hw1(i))^.2+abs(hw2(i))^.2;
end
plot(w,yy);