谱减法去噪matlab,最基本的谱减法去噪程序

各位大大好，新人刚学matlab，需要设计一个最基本的谱减法去噪程序；我按照课题要求写了，运行无错误，但是输出结果明显不对，还请高人解惑，感谢在先。我估计是哪里犯了严重错误……

clear

clc

[x,fs] = wavread('0db.wav');                           % read in the file, fs = frequency=48000

y = downsample(x,3);                                    % downsampling (fs/3)

figure;plot(y)

ind = vadsohn(y,16000,'n');                            %相当于vad检测silent或者voiced段，输出是0或者1

f=enframe(y,hamming(320),160);                   % enframing加窗，汉明窗，50%overlap

nmag = [];smag = [];r = [];                           %初始化

for i = 1:5                                                     %这里开始是获取初始噪音，按照课题要求直接拿开头5个frame

fi = f(i,:);

mag(i) = sum(abs(fft(fi,320)))/length(fi);

end

nmag = sum(mag)/length(mag);