【原创】：利用Pmcalgorithm算法设计低通滤波器(LPF)

ljelly

归一化设计，代码如下：
% ****************************************
% Design a LPF with equal ripple method  *
% ****************************************
wp=0.2*pi; //通带频率
ws=0.3*pi; //阻带频率
Rp=0.15; // 通带纹波
As=80;  // 阻带衰减
delta1=(10^(Rp/20)-1)/(10^(Rp/20)+1)
delta2=(1+delta1)*(10^(-As/20))

f = [wp/pi ws/pi];
m=[1 0];
dev=[delta1 delta2];
[N,fo,ao,weights] = remezord(f,m,dev)
h=remez(N,fo,ao,weights);
[db,mag,pha,grd,w]=freqz_m(h,[1]);
Asd=-max(db(ws:1:501));
N=N+4
h=remez(N,fo,ao,weights);
[db,mag,pha,grd,w]=freqz_m(h,[1]);
delta_w=2*pi/1000;
wsi=ws/delta_w+1;                                                                                                                                                                                                                                                                                                                                                                  
wpi=wp/delta_w;
Asd=-max(db(wsi:1:501))
Rpd=-min(db(1:1:wp/delta_w+1))

figure(1);
subplot(2,1,1);stem([0:1:N],h);title('Actual Impulse Response');
axis([0 N -0.1 0.3]);xlabel('n');ylabel('h(n)');
set(gca,'XTickMode','manual','XTick',[0,N]);
set(gca,'YTickMode','manual','YTick',[-0.1:0.1:0.3]);
grid on;
subplot(2,1,2);plot(w/pi,db);title('Magnitude Response in dB');
axis([0 1 -100 10]);xlabel('frequency in pi units');ylabel('DECIBELS');
set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]);
set(gca,'YTickMode','manual','YTick',[-70,0]);
grid on;

不知道为什么，我上传不了附件，图只能大家自己运行看一下了!

[ 本帖最后由 ljelly 于 2008-11-4 11:12 编辑 ]