matlab 谱减法降噪程序

function output=simplesubspec（signalwlenincNISab）
wnd=hamming（wlen）;                      % 设置窗函数
N=length（signal）;                       % 计算信号长度

y=enframe（signalwndinc）‘;             % 分帧
fn=size（y2）;                           % 求帧数

y_fft = fft（y）;                         % FFT
y_a = abs（y_fft）;                       % 求取幅值
y_phase=angle（y_fft）;                   % 求取相位角
y_a2=y_a.^2;                            % 求能量
Nt=mean（y_a2（:1:NIS）2）;               % 计算噪声段平均能量
nl2=wlen/2+1;                           % 求出正频率的区间