对应的 MATLAB 代码
function mpma2
n=0:19;
h1=zeros(1,20);
h2=h1;
h1(1:5)=1/5;
h2(1:14)=1/14;
[H1,omega1]=freqz(h1);
[H2,omega2]=freqz(h2);
ph11=angle(H1);
ph21=angle(H2);
ph12=unwrap(ph11*180/pi);
ph22=unwrap(ph21*180/pi);
subplot(311)
plot(omega1/pi,abs(H1),omega2/pi,abs(H2));
title('amplitude')
legend('M=5','m=14')
subplot(312)
plot(omega1/pi,ph11,'g',omega2/pi,ph21,'b');
title('wrapped phase')
legend('M=5','M=14')
subplot(313)
plot(omega1/pi,ph12*pi/180,'--g',omega2/pi,ph22*pi/180,'--b')
title('unwrapped phase')
legend('M=5','M=14');
将一正弦信号加上高斯白噪声之后,投入 M=5 的 MPM 滤波器中,
得到的原信号和输出信号的图像如下
其中原信号的参数为
高斯噪声参数为