武汉理工大学《信息群处理课程设计》报告书
a=0.5; %a取小于等于1
Bz1=[a,0,0,0,0,0,0,0,0,0,1]; %分子的系数 Az1=[1,0,0,0,0,0,0,0,0,0,a]; %分母的系数
yy4=filter(Bz1,Az1,x); figure(8); %滤波器进行滤波
subplot(2,1,1);plot(yy4); title('全通结构混响器的时域图')%全通结构的混响器时域波形
YY4=fft(yy4,2001); %经全通结构的混响器后的信号做2001点的FFT变换
subplot(2,1,2);plot(n1(1:1000),YY4(1:1000)); %全通结构的混响器频谱图
title('全通结构混响器的频谱图')sound(5*yy4,fs,bits);
%回放经全通结构的混响器后的语音信号
22