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（‘聲音信號(hào)的波形‘）;%y代表信號(hào)
%得到病人數(shù)據(jù)
pat=3;
[?hearloss?]?=?patientchar（pat）;
[IG40IG60IG95]?=insertgain（hearloss）;
%初始化
a1=zeros（164）;
a2=zeros（164）;
a3=zeros（164）;
a4=zeros（164）;
f2=zeros（1256）;
%信號(hào)處理
?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為取出的樣點(diǎn)長(zhǎng)度256
b1=f（1:128）;
b2=f（129:256）;
y1=b1+b2;
cir=mod（i*64128）;%循環(huán)移位
y2=cir