MATLAB读取gdas1气象数据,南京理工大学数字信号处理matlab上机完美版

1.已知3阶椭圆IIR数字低通滤波器的性能指标为：通带截止频率0.4π，通带波纹为0.6dB，最小阻带衰减为32dB。设计一个6阶全通滤波器对其通带的群延时进行均衡。绘制低通滤波器和级联滤波器的群延时。

%Q1_solution

%ellip(N,Ap,Ast,Wp)

%N--->The order of the filter

%Ap-->ripple in the passband

%Ast->a stopband Rs dB down from the peak value in the passband

%Wp-->the passband width

[be,ae]=ellip(3,0.6,32,0.4);

hellip=dfilt.df2(be,ae);

f=0:0.001:0.4;

g=grpdelay(hellip,f,2);

g1=max(g)-g;

[b,a,tau]=iirgrpdelay(6,f,[0 0.4],g1);

hallpass=dfilt.df2(b,a);

hoverall=cascade(hallpass,hellip);

hFVT=fvtool([hellip,hoverall]);

set(hFVT,'Filter',[hellip,hoverall]);

legend(hFVT,'Lowpass Elliptic filter','Compensated filter');

clear;

[num1,den1]=ellip(3,0.6,32,0.4);

[GdH,w]=grpdelay(num1,den1,512);

plot(w/pi,GdH); grid

xlabel('\omega/\pi'); ylabel('Group delay, samples');

F=0:0.001:0.4;

g=grpdelay(num1,den1,F,2); % Equalize the passband

Gd=max(g)-g;

% Design the allpass delay equalizer

[num2,den2]=iirgrpdelay(6,F,[0,0.4],Gd);

[GdA,w] = grpdelay(num2,den2,512);

hold on;

plot(w/pi,GdH+GdA,'r');