完全重建QMF滤波器组实验

作业一

1.要求：

确定两个长度为100（或其他）的序列，一个脉冲在前一个脉冲在后，分别进行傅里叶变换并对比。

2.MATLAB代码：

clc;
clear;
close all;

N=100;
n=[0:1:N-1];

%脉冲在前的信号
x1=normpdf(n,30,1);
Xk1=fft(x1,N);
%画出该信号的波形
subplot(231);
stem(n,x1);
title('脉冲在前的原信号');
subplot(232);
stem(n,abs(Xk1));
title('脉冲在前的信号的幅度谱');
subplot(233);
stem(n,angle(Xk1)) 
title('脉冲在前的信号的相位谱');

%脉冲在后的信号
x2=normpdf(n,70,1);
Xk2=fft(x2,N);
%画出该信号的波形
subplot(234);
stem(n,x2);
title('脉冲在后的原信号'); 
subplot(235);
stem(n,abs(Xk2));
title('脉冲在后的信号的幅度谱');
subplot(236);
stem(n,angle(Xk2));
title('脉冲在后的信号的相位谱');

3.波形对比：

4.结论：

从得到的波形对比图可以看出，这两个脉冲位置不同的序列的幅度谱是一样的。结合傅里叶变换可知，时域上的移动对应到频域就是相位的变化，而这里确实两个序列的不同就体现在相位谱上，与理论是相符的。

作业二

1.引言：

一个两通道正交镜像滤波器组如图所示，在分析滤波器组一侧，输入信号（设为宽带信号）被分成K个子频带信号（窄带信号），通过抽取可降低采样率；在综合滤波器一侧，通过零值内插和带通滤波可以重建原来的信号。

2.MATLAB代码实现完全重建AMF滤波器组：

N=41;
w=0.43;
[h0,h1,g0,g1]=firpr2chfb(N,w);
[H1z