LMS算法学习曲线.docxVIP

%%Matlab 代码 本例中以通信原理中信号信道传输为例： clear all;close all;clc %%%%%%%%%%%%%%参数设 置%%%%%%%%%%%%%%%%%%%%%%% step_le n=0.05;%步长 varia nce_v=0.02;%方差 rep eat_times=100;%重 复次数 iteration_times=500;% 迭代次数 filt_len=2;%滤波器长度 %%%%%%%%%%%%%%%%%%%%变量及数组初始 化%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Square_err=zeros(repeat_times,iteration_times); Chan_factor1=[-0.195 0.95]; X=zeros(filt_len,1); W1=[]; W2=[]; Wm1=[repeat_times,iteration_times]; Wm2=[repeat_times,iteration_times]; %%%%%%%%%%%%%夕卜循环，重复做repeat_timws次实验，平均 值%%%%%%%%%%%%%%%% for loo p1=1:re peat_times%循环次数 %%%%%%%%%%%%%%内循环变量及数组初始 化%%%%%%%%%%%%%%%%%%%%%%%% Source=zeros(iteration_times,1); W=zeros(filt_len,1); %%%%%%%%%%%%% 内循环，做 iteration_times 次迭 代%%%%%%%%%%%%%%%%%%%%%%%%%%%% for i=1:iteration_times %%%%%%%%%%%%%%信源，产生等概率分布的正负 1%%%%%%%%%%%%%%%%%%%%%%%%%%%% if rand0.5 Source(i)=-1; else Source(i)=1; end end %%%%%%%%%%%%%%信道，滤波和加噪 声%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Revdata=filter(Chan_factor1,1,Source)+randn(iteration_times,1)*sqrt(variance_v); % filter the Source %%%%%%%%%%%%%%均衡器，基于 LMS 算 法%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% for loop2=1:iteration_times-filt_len+1 U=Source(loop2:loop2+filt_len-1); d=Revdata(loop2+1); e=d-W*U; W=W+step_len*U*e; Wm1(loop2,loop1)=W(1); Wm2(loop2,loop1)=W(2); Square_err(lo op 1,lo op 2)=e.^2; end end %%%%%%%%%%%%%%作图，显示结果 %%%%%%%%%%%%%%%%% Square_err_mean=sum(Square_err(:,:))/repeat_times; % Wm2 W1=Wm1(:,repeat_times); W2=Wm2(:,repeat_times); meanW1=mean(Wm1,2); meanW2=mean(Wm2,2); Square_err; %semilogy(Square_err_mean); figure(1); plot(Square_err_mean./max(Square_err_mean),LineWidth,2); hold on; suijishu=randperm(repeat_times); plot(Square_err(suijishu(1),:)./max(Square_err(suijishu(1),:)),r); legend([num2str(repeat_times)次重复试验平均误差],单次误差); xIabelC迭代次数);ylabel(均方误差);title([步长为 ,num2str(step_len),num2str(iteration_times)次迭代均方误差变化曲线]); grid on; 00 0 0 figure (2); plot(W1);hold on; plot(W2,r);hold on; plot(mea nW1,k);hold on; plot(mea nW2,--k); legend(W1,W2,[num2str(repeat_times),次重复试验 W1],[num2str(repeat_