bpsk.rar_BPSK的LMS_LMS BPSK_LMS算法C_bpsk lms

clc num=10; %码元个数 tnum=200;%码元长度 Np=2000;%采样点数 M=11;%抽头数 q=1000;%仿真次数 u1=0.02;%步长 w1=zeros(1,M);%设置抽头加权值 pp=zeros(q,Np-M+1); N=num*tnum; a=randint(1,num,2); %产生1行num列的矩阵，矩阵内0和1随机出现 fc=0.1; %载波频率 t=0:0.05:9.99; s=[];c=[]; for i=1:num if(a(i)==0) A=zeros(1,tnum); else A=ones(1,tnum); end s=[s A]; cs=sin(2*pi*fc*t); c=[c cs]; %%%%%%%%%%%%%%%%%%%%%%%%产生载波信号 end figure(1);subplot(4,1,1);plot(s);grid on; xlabel('基带信号s(t)'); axis([0 N -2 2]); %%%%%%%%%%%%%%%%%%%%%产生BPSK已调信号。采用模拟调制方法 s_NRZ=[]; for i=1:num if(a(i)==0) A=ones(1,tnum); else A=-1*ones(1,tnum); end s_NRZ=[s_NRZ,A]; end d=s_NRZ.*c; figure(1);subplot(4,1,2);plot(d);grid on; xlabel('BPSK调制信号'); axis([0 N -2 2]); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%加白噪声 sh1=sqrt(10^(-20/10)); un1=sh1*randn(1,Np); un=un1+d; figure(1);subplot(4,1,3);plot(un);grid on; xlabel('加入白噪声的BPSK调制信号'); for qq=1:q %%%%%%%%%%%%%%%%用LMS算法迭代滤波 y1(1:M)=un(1:M); for m=M:Np; UN=[un(m) un(m-1) un(m-2) un(m-3) un(m-4) un(m-5) un(m-6) un(m-7) un(m-8) un(m-9) un(m-10)]; % UN=un((m-M+1):(m)); y1(m)=w1*UN'; e1(m)=d(m)-y1(m); w1=w1+2*u1*e1(m)*UN; end pp1(qq,:)=(e1(M:Np)).^2; end for c=1:Np-M+1%均方误差 b1(c)=sum(pp1(:,c))/q; end figure(1);subplot(4,1,4);plot(y1);grid on; xlabel('经过自适应滤波的输出信号'); figure(2);plot(b1,'r'); xlabel('迭代次数'); ylabel('均方误差'); axis([0,Np,0,1]); grid; hold on