close?all;
clear?all;

N=400000;??????????%輸入數(shù)據(jù)長(zhǎng)度
[dnFsNb]=wavread（‘yuanshi.wav‘N）;??%讀入原始信號(hào)，即語(yǔ)音加噪聲
[noisenFsNb]=wavread（‘zao.wav‘N）;??%讀入噪聲
%---------------使原始信號(hào)和噪聲同步————————%
m=xcorr（dnnoisen）;
[Mnum]=max（m）;
plot（m）
grid?on;
title（‘原始信號(hào)和噪聲的互相關(guān)‘）;
flag=num-N;
if?flag>0
????shiftn=flag;
????d=dn（shiftn+1:N）;
????noise=noisen（1:N-shiftn）;
elseif?flag<0
????shiftn=-flag;
????d=dn（1:N-shiftn）;
????noise=noisen（shiftn+1:N）;
else
????d