CMA_Equalizer.doc_第1页
CMA_Equalizer.doc_第2页
CMA_Equalizer.doc_第3页
CMA_Equalizer.doc_第4页
CMA_Equalizer.doc_第5页
已阅读5页，还剩15页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

(6300:6800)(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:m a=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; % 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 onplot(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 offn=700 u=0.001 (100:680)n=700 u=0.01 (100:680)clear 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:m a=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; e1=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 的比较 clear 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:m a=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(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; 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);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 offclear 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:m a=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; % 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_av2+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 offu2=0.04 clear allM=2; 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:m a=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)=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 onscatterplot(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:m a=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; % 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 & imag(z(i)0; z1(i)=1+j; end;if real(z(i)0 & imag(z(i)0; z1(i)=1-j;end;if real(z(i)0 & imag(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*
人人文库> 全部分类> 应用文书 > 技术指导

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论