testwave.m

利用matlab读入wav音频

filename='Windows XP 电话拨入声.wav';
% Fs=12000;       % 12KHz
[y,Fs] = wavread(filename);
gain = 24;
y = gain * y;   %  提高音量
wavplay(y,Fs );          % gong signal begins playing.

subplot(2,1,1),plot(y(1 : fix(end/8)));title('原声')
x=y;
b=[ 
    0.09393423303564  
   -0.16284746686676
    0.09332757238627
    0.00324603545919
    0.00519151173625
   -0.04572863261770
    0.03284820495831
   -0.03020744855625
    0.07211329425296
   -0.07532762770276
    0.02501463773315
   -0.00022433000271
    0.01859168007976
   -0.03693965320228
    0.03418269061323
   -0.03635041949239
    0.03574972598228
   -0.02828369831337
    0.02405308778505
   -0.01437328116435
    0.00035267234096
    0.00150611385749
    0.00750470750932
   -0.01364292812061
    0.00636513782032
];
a=[
    1.0
   -0.60321120349729
   -0.43642949988157
    0.57020912710538
    0.05095591874974
   -0.38485083647156
   -0.00071267967606
   -0.09847674399986
    0.27096900161101
    0.14497007369623
   -0.18209318108990
    0.25538642287687
   -0.12176700888090
   -0.21574310691803
    0.16446450584643
    0.00363859185489
   -0.10186630200717
   -0.09872384968285
    0.11675452664540
    0.11820736487318
   -0.19973332096141
    0.15334567116746
   -0.09076585614211
    0.02006147771053
   -0.00680573138016
];
y=24.0*filter(b,a,x);    % 音量提高点
wavplay(y,Fs );          % gong signal begins playing.
subplot(2,1,2),plot(y(1:fix(end/8)));title('均衡后')