lmsjs.m

基于自适应滤波器的LMS算法仿真程序的计算部分

%写到附录中
close all;
fc=1e4;
fd=5000;
[fn,fp]=uigetfile('*.wav','pick a wav-file');%读取音频文件；
[y,fs,bit]=wavread([fp,fn]);
for i=1:fs;
    y1(i,1)=y(i,1);
end
figure(1);
plot(y1);
title('时域')
xlabel('t/s');
ylabel('y1(t)');
yy1=fft(y1);
py=abs(yy1);
figure(2);
plot(py);
title('频域')
xlabel('frequency/Hz');
ylabel('py(e^jw)');
n=1:fs;
opt=(fc/fs)*2*pi/(max(max(y1)));
x=cos(2*pi*fc*n'/fs+opt*cumsum(y1)); %信号调制;
x1=fft(x);
px=abs(x1);
carr=cos(2*pi*fd*n'/fs);
adata=x.*carr;  %信号搬移；
x2=fft(adata);
px2=abs(x2);
figure(3);
subplot(211);
plot(px);
axis([0 2.5e4 0 1000]);
title('信号调制')
xlabel('frequency/Hz');
ylabel('px(e^jw)');
subplot(212);
plot(px2);
axis([0 2.5e4 0 1000]);
title('信号搬移')
xlabel('frequency/Hz');
ylabel('px(e^jw)');