參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  盼╃字JUN團 2017-11-24 17:50:22

  clear all close all N=10; %濾波器階數(shù) sample_N=500; %采樣點數(shù) A=1; %信號幅度 snr=10; %信噪比 t=1:sample_N; length_t=100; %期望信號序列長度 d=A*sin(2*pi*t/length_t); %期望信號 M=length(d); %M為接收數(shù)據(jù)長度 x=awgn(d，snr); %經(jīng)過信道（加噪聲） delta=1/(10*N*(A^2)); %計算能夠使LMS算法收斂的delta y=zeros(1，M); h=zeros(1，N); %LMS濾波器系數(shù) h_normalized=zeros(1，N); %歸一化LMS濾波器系數(shù) y1=zeros(1，N); for n=N:M %系數(shù)調(diào)整LMS算法 x1=x(n:-1:n-N+1); %LMS算法 y(n)=h*x1'; e(n)=d(n)-y(n); h=h+delta*e(n)*x1; %NLMS算法 y_normalized(n)=h_normalized*x1'; e_normalized(n)=d(n)-y_normalized(n); h_normalized=h_normalized+e_normalized(n)*x1/(x1*x1'); end error=e.^2; %LMS算法每一步迭代的均方誤差 error_normalized=e_normalized.^2; %NLMS算法每一步迭代的均方誤差 for n=N:M %利用求解得到的h，與輸入信號x做卷積，得到濾波后結(jié)果 x2=x(n:-1:n-N+1); y1(n)=h*x2'; y2(n)=h_normalized*x2'; end subplot(411) plot(t，d); axis([1，sample_N，-2，2]); subplot(412) plot(t，x); subplot(413) plot(t，y); subplot(414) plot(t，y_normalized); figure(2) plot(t，error，'r'，t，error_normalized，'b');

  贊(10)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論