完全重建QMF滤波器组

任务一 完全重建QMF滤波器组的设计

在QMFB的设计仿真中，发现1<N<127时，N越大，低通滤波器的幅度响应越理想，并且达到通带尽量平，过渡带尽量窄，且阻带尽可能快速衰减的要求，具有实际意义。这里用matlab实现完全重建QMFB的设计，具有直观、简便、快捷的优点，极大的减轻了工作量，得到的均方误差也较小，具有相对较高的精度，最优参数一般可以选取N=41，=0.43。

结果及分析:

图1为滤波器h0和h1的幅度响应,h0与h1镜像对称，h0的通带截止频率为w=0.43：

图2为幅度响应关系误差，误差很小：

图3为滤波器h0的幅度：
图4为输入信号也就是原信号波形：

图5为重建输出信号：

图6为理想输出信号与重建输出信号的偏差。观察可以发现偏差的绝对值在4x10 ^-4之间，差值很小，因此可以实现QMFB的近似完全重建。

任务二：仿真分析两种波形的傅里叶变换，一个在序列前端有值，一个在后端有值

信号y1为长20的序列，n=3处有一单位脉冲；信号y2为长20的序列，n=18处有一单位脉冲.

观察比较两信号的相位谱与幅度谱发现，时域移动，幅度谱几乎不变，相位谱发生改变。

代码：

clear;
clc;
close all;

x1=[zeros(1,2),1,zeros(1,17)];
x2=[zeros(1,17),1,zeros(1,2)];

%绘制原始时域信号波形
figure(1)
subplot(2,1,1)
stem(x1)
title('信号1的时域波形')
subplot(2,1,2)
stem(x2)
title('信号2的时域波形')

Y1=fft(x1);
Y2=fft(x2);
figure(2)
subplot(3,1,1)
stem(abs(Y1));
title('Y1的幅度谱');
subplot(3,1,2)
stem(abs(Y2));
title('Y2的幅度谱');
subplot(3,1,3)
stem(abs(Y1)-abs(Y2));
title('Y1-Y2');

figure(3)
subplot(2,1,1)
stem(angle(Y1));
title('Y1的相位谱');
subplot(2,1,2)
stem(angle(Y2));
title('Y2的相位谱');