555.txt

可用于对信号的频谱进行分析,希望站长能整理好,以便大家互相学习

clear
s=wavread('C:\Documents and Settings\zhuwenhao\桌面\毕设\XJTU.wav');   %根据实际情况填入要被处理的文件的路径和文件名
%可在 COMMAD 窗口输入help wavread 查看其功能描述 
a=size(s);
u=0.002;      %收敛步长
M=64;          %滤波器阶数
N=a(1,1);
X=zeros(M,1);
W=zeros(M,1);
out=zeros(N,1);
i=1;
while i<=N
X(1,1)=s(i,1);
y=W'*X;
e=X(1,1)-y;
out(i,1)=y;
W=W+u*e*X;
i=i+1;
end
wavwrite(out,'C:\Documents and Settings\zhuwenhao\桌面\毕设\88882.wav');    
subplot(2,1,1),plot(s,'b'),axis([0,N,-2,2]),title('Input Signal')
subplot(2,1,2),plot(out,'r'),axis([0,N,-2,2]),title('Output Signal')