exe0705.m

Matlab程序为清华大学胡广书版信号处理的上机习题答案

%----------------------------------------------------
% exe0705.m
%----------------------------------------------------
clear all;

N=55;
M=N-1;
% 用切比雪夫最佳一致逼近设计线性相位多带FIR滤波器
f=[0 .05 .10 .15 .18 .25 .30 .36 .41 .50];
f=2.*f;
A=[0 0 1 1 0 0 1 1 0 0];
weigh=[1 10 1 10 1];	
b=remez(M,f,A,weigh);
%
[h,w]=freqz(b,1,256,1);
hr=abs(h);
hp=unwrap(angle(h));
h=20*log10(hr);

figure(1)
stem(b,'.');grid;
ylabel('h(n)')
xlabel('n');

figure(2)
plot(w,h);grid;
ylabel('20*log(Hr)')
xlabel('omg/2pi');

figure(3)
plot(w,hp);grid;
ylabel('hpase')
xlabel('omg/2pi');