1、0 1000 2000 3000 4000 5000 6000 70000.40.60.811.21.41.61.822.22.4-1 -0.5 0 0.5 1-1-0.8-0.6-0.4-0.200.20.40.60.81QuadratureIn-PhaseScatter plot-2 -1 0 1 2-2.5-2-1.5-1-0.500.511.522.5QuadratureIn-PhaseScatter plot-2 -1 0 1 2-2-1.5-1-0.500.511.52QuadratureIn-PhaseScatter plot(6300:6800)-2 -1 0 1 2-2-1.
2、5-1-0.500.511.52QuadratureIn-PhaseScatter plot(1000:6800)clear allM=4; k=log2(M); n=7000; %u=0.05;u=0.001;m=400;%h=0.05 -0.063 0.088;%-0.126; -0.25;h=1 0.3 -0.3 0.1 -0.1;L=7; mse_av=zeros(1,n-L+1);for j=1:ma=randint(1,n,M); a1=qammod(a,M);m1=abs(a1).4;m2=abs(a1).2; r1=mean(m1);r2=mean(m2); R2=r1/r2;
3、% R2=sqrt(2);s=filter(h,1,a1); snr=15; x=awgn(s,snr,measured);c=0 0 0 1 0 0 0;for i=1:n-L+1 y=x(i+L-1:-1:i);z(i)=c*y;e=R2-(abs(z(i)2);%u=0.0005*(1-(1/sqrt(2*pi)*(exp(50*(-e2/2);%u=0.0005*(1-(exp(-50*e2);%u=e*(1-exp(-5*e2);c=c+u*e*y*z(i); mse(i)=e2;end;mse_av=mse_av+mse;end;mse_av=mse_av/m;hold onplo
4、t(1:n-L+1,mse_av,r)hold onscatterplot(a1,1,0,r*);hold onscatterplot(x,1,0,g*);hold onscatterplot(z(1000:6800),1,0,r*);hold off0 100 200 300 400 500 600 7000.811.21.41.61.822.22.42.6-1.5 -1 -0.5 0 0.5 1 1.5-1.5-1-0.500.511.5QuadratureIn-PhaseScatter plotn=700 u=0.001 (100:680)0 100 200 300 400 500 60
5、0 7000.811.21.41.61.822.22.42.6-2 -1 0 1 2-2-1.5-1-0.500.511.52QuadratureIn-PhaseScatter plotn=700 u=0.01 (100:680)0 500 1000 1500 2000 2500 3000 3500 4000 4500 500000.511.522.5clear allM=4; k=log2(M); n=5000; %u=0.05;u1=0.001;u2=0.0001;m=500;%h=0.05 -0.063 0.088;%-0.126; -0.25;h=1 0.3 -0.3 0.1 -0.1
6、;L=7; mse1_av=zeros(1,n-L+1);mse2_av=mse1_av;for j=1:ma=randint(1,n,M); a1=pskmod(a,M);m1=abs(a1).4;m2=abs(a1).2; r1=mean(m1);r2=mean(m2); R2=r1/r2;%R2=sqrt(2%);s=filter(h,1,a1); snr=15; x=awgn(s,snr,measured);c1=0 0 0 1 0 0 0;c2=c1;y=zeros(n-L+1,2);for i=1:n-L+1 y=x(i+L-1:-1:i);z1(i)=c1*y;z2=c2*y;e
7、1=R2-(abs(z1(i)2);e2=a1(i)-z2;c1=c1+u1*e1*y*z1(i);c2=c2+u2*e2*y;mse1(i)=e12;%u(i)=0.2*(1-exp(-(0.3*abs(e(i); mse2(i)=abs(e2)2; end;mse1_av=mse1_av+mse1;mse2_av=mse2_av+mse2;end;mse1_av=mse1_av/m;mse2_av=mse2_av/m;figureplot(1:n-L+1,mse1_av,r,1:n-L+1,mse2_av,b)axis(0,5100,0 2.8);lms 和 cma 的比较0 500 10
8、00 1500 2000 2500 3000 3500 4000 4500 500000.511.522.5-1 -0.5 0 0.5 1-1-0.8-0.6-0.4-0.200.20.40.60.81QuadratureIn-PhaseScatter plot-2 -1 0 1 2-2-1.5-1-0.500.511.52QuadratureIn-PhaseScatter plot-1.5 -1 -0.5 0 0.5 1 1.5-1.5-1-0.500.511.5QuadratureIn-PhaseScatter plot-1.5 -1 -0.5 0 0.5 1 1.5-1.5-1-0.50
9、0.511.5QuadratureIn-PhaseScatter plotclear allM=4; k=log2(M); n=5000; %u=0.05;u1=0.001;u2=0.0001;m=500;%h=0.05 -0.063 0.088;%-0.126; -0.25;h=1 0.3 -0.3 0.1 -0.1;L=7; mse1_av=zeros(1,n-L+1);mse2_av=mse1_av;for j=1:ma=randint(1,n,M); a1=pskmod(a,M);m1=abs(a1).4;m2=abs(a1).2; r1=mean(m1);r2=mean(m2); R
10、2=r1/r2;%R2=sqrt(2%);s=filter(h,1,a1); snr=15; x=awgn(s,snr,measured);c1=0 0 0 1 0 0 0 ;c2=c1;y=zeros(n-L+1,2);for i=1:n-L+1 y=x(i+L-1:-1:i);z1(i)=c1*y;z2(i)=c2*y;e1=R2-(abs(z1(i)2);e2=a1(i)-z2(i);c1=c1+u1*e1*y*z1(i);c2=c2+u2*e2*y;mse1(i)=e12;%u(i)=0.2*(1-exp(-(0.3*abs(e(i); mse2(i)=abs(e2)2; end;ms
11、e1_av=mse1_av+mse1;mse2_av=mse2_av+mse2;end;mse1_av=mse1_av/m;mse2_av=mse2_av/m;figureplot(1:n-L+1,mse1_av,r,1:n-L+1,mse2_av,b)axis(0,5100,0 2.8);scatterplot(a1,1,0,r*);hold onscatterplot(x,1,0,g*);hold onscatterplot(z1(2300:4800),1,0,r*);hold offscatterplot(z2(2300:4800),1,0,r*);hold off0 100 200 3
12、00 400 500 600 7000.811.21.41.61.822.22.42.6clear allM=4; k=log2(M); n=700; %u=0.05;u1=0.001;u2=0.01;m=400;%h=0.05 -0.063 0.088;%-0.126; -0.25;h=1 0.3 -0.3 0.1 -0.1;L=7; mse_av1=zeros(1,n-L+1);mse_av2=zeros(1,n-L+1);for j=1:ma=randint(1,n,M); a1=qammod(a,M);m1=abs(a1).4;m2=abs(a1).2; r1=mean(m1);r2=
13、mean(m2); R2=r1/r2;% R2=sqrt(2);s=filter(h,1,a1); snr=15; x=awgn(s,snr,measured);c1=0 0 0 1 0 0 0;c2=c1;for i=1:n-L+1 y=x(i+L-1:-1:i);z1(i)=c1*y; z2(i)=c2*y;e1=R2-(abs(z1(i)2);e2=R2-(abs(z2(i)2);c1=c1+u1*e1*y*z1(i); c2=c2+u2*e2*y*z2(i); mse1(i)=e12; mse2(i)=e22;end;mse_av1=mse_av1+mse1;mse_av2=mse_a
14、v2+mse2;end;mse_av1=mse_av1/m;mse_av2=mse_av2/m;hold onplot(1:n-L+1,mse_av1,r,1:n-L+1,mse_av2,g)hold onscatterplot(x,1,0,g*);hold off0 100 200 300 400 500 600 700024681012x 10158u2=0.040 1000 2000 3000 4000 5000 6000 70000.20.40.60.811.21.41.6-1 -0.5 0 0.5 1-1-0.8-0.6-0.4-0.200.20.40.60.81Quadrature
15、In-PhaseScatter plot-2 -1 0 1 2-2.5-2-1.5-1-0.500.511.522.5QuadratureIn-PhaseScatter plot-1.5 -1 -0.5 0 0.5 1 1.5-1.5-1-0.500.511.5QuadratureIn-PhaseScatter plot-1 -0.5 0 0.5 1-1-0.8-0.6-0.4-0.200.20.40.60.81QuadratureIn-PhaseScatter plotclear allM=2; n=7000; %u=0.05;u=0.001;m=400;%h=0.05 -0.063 0.0
16、88;%-0.126; -0.25;h=1 0.3 -0.3 0.1 -0.1;L=7; mse_av=zeros(1,n-L+1);for j=1:ma=randint(1,n,M); a1=pammod(a,M);m1=abs(a1).4;m2=abs(a1).2; r1=mean(m1);r2=mean(m2); R2=r1/r2;% R2=sqrt(2);s=filter(h,1,a1); snr=15; x=awgn(s,snr,measured);c=0 0 0 1 0 0 0;for i=1:n-L+1 y=x(i+L-1:-1:i);z(i)=c*y;if z(i)0;z1(i
17、)=1;else z1(i)=-1;end ;e=R2-(abs(z(i)2);%u=0.0005*(1-(1/sqrt(2*pi)*(exp(50*(-e2/2);%u=0.0005*(1-(exp(-50*e2);%u=e*(1-exp(-5*e2);c=c+u*e*y*z(i); mse(i)=e2;end;mse_av=mse_av+mse;end;mse_av=mse_av/m;hold onplot(1:n-L+1,mse_av,r)hold onscatterplot(a1,1,0,r*);hold onscatterplot(x,1,0,g*);hold onscatterpl
18、ot(z(1000:6800),1,0,r*);hold onscatterplot(z1(1000:6800),1,0,r*);hold offclear allM=4; k=log2(M); n=700; %u=0.05;u=0.001;m=400;%h=0.05 -0.063 0.088;%-0.126; -0.25;h=1 0.3 -0.3 0.1 -0.1;L=7; mse_av=zeros(1,n-L+1);for j=1:ma=randint(1,n,M); a1=qammod(a,M);m1=abs(a1).4;m2=abs(a1).2; r1=mean(m1);r2=mean
19、(m2); R2=r1/r2;% R2=sqrt(2);s=filter(h,1,a1); snr=15; x=awgn(s,snr,measured);c=0 0 0 1 0 0 0;for i=1:n-L+1 y=x(i+L-1:-1:i);z(i)=c*y;if real(z(i)0 z1(i)=1+j;end;if real(z(i)0 z1(i)=-1-j;end;e=R2-(abs(z(i)2);%u=0.0005*(1-(1/sqrt(2*pi)*(exp(50*(-e2/2);%u=0.0005*(1-(exp(-50*e2);%u=e*(1-exp(-5*e2);c=c+u*e*y*z(i); mse(i)=e2;end;mse_av=mse_av+mse;end;mse_av=mse_av/m;hold onplot(1:n-L+1,mse_av,r)hold onscatterplot(a1,1,0,r*);hold onscatterplot(x,1,0,g*);hold onscatterplot(z(100:680),1,0,r*);hold onscatterplot(z1(100:680),1,0,r*);hold off
