WOLA加权叠加滤波器组算法，matlab

clear all;
clc
bit=8;
fs=16000;
[yfsbit]=wavread（‘MaleVoice.wav‘）;
%wavplay（yfs）;
%y=[1;zeros（100001）];
[HW]=freqz（y1）;
mag=abs（H）;
fb=20*log10（mag/max（mag））;
figure（1）
subplot（211）
plot（W/pifb）; 
figure（2）
subplot（211）;
plot（y）;
title（‘声音信号的波形‘）;%y代表信号
%得到病人数据
pat=3;
[ hearloss ] = patientchar（pat）;
[IG40IG60IG95] =insertgain（hearloss）;
%初始化
a1=zeros（164）;
a2=zeros（164）;
a3=zeros（164）;
a4=zeros（164）;
f2=zeros（1256）;
%信号处理
 for i=1:fix（length（y）/64）
     a4=a3;
     a3=a2;
     a2=a1;
     a1=y（i*64-63:i*64）‘;
    x=[a1a2a3a4];
f=x‘.*hamming（length（x））;%x为取出的样点长度256
b1=f（1:128）;
b2=f（129:256）;
y1=b1+b2;
cir=mod（i*64128）;%循环移位
y2=cir