MCRA噪声估计

function noise=MCRA（x）

count=256;
fs=8000;
%x=wavread（‘d:/0dB/sp01_airport_sn0.wav‘）;
%y=wavread（‘d:/clean/sp01.wav‘）;
%x2=x-y;
x1=zeros（1count）;
x6=zeros（1count）;
t=0:count/（2*fs）:ceil（length（x）-count）/fs;

powerspeech=zeros（count2*（ceil（length（x）/count）-1）-1）;
sf=zeros（count2*（ceil（length（x）/count）-1）-1）;
s=zeros（count2*（ceil（length（x）/count）-1）-1）;
stmp=zeros（count2*（ceil（length（x）/count）-1）-1）;
smin=zeros（count2*（ceil（length（x）/count）-1）-1）;
sr=zeros（count2*（ceil（length（x）/count）-1）-1）;
p=zeros（count2*（ceil（length（x）/count）-1）-1）;
psmooth=zeros（count2*（ceil（length（x）/count）-1）-1）;
alfad=zeros（count2*（ceil（length（x）/count）-1）-1）;
noise=zeros（count2*（ceil（length（x）/count）-1）-1）;
true=zeros（count2*（ceil（length（x）/count）-1）-1）;


d=80;alfas=0.8;deta=5;