cdma的MATLAB仿真源程序

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

%*************************************************************************************%ThisfunctionpertainstotheadditionofAWGNwithmeanzeroand%parameter'variance'toaninputsignal.%%AUTHOR:WenbinLuo%DATE:04/12/01%%SYNOPSIS:y=awgn(x,var)%x---inputsignal%var---variance%y---y=x+AWGN%***********************************************************************************functiony=awgn(x,var)w=randn(1,length(x));w=w-mean(w)*ones(size(w));w=sqrt(var)*(w/std(w));x=x(:);w=w(:);y=x+w;%*************************************************************************************%ThisfunctiondoestheDS-SSmodulation%%AUTHOR:WenbinLuo%DATE:04/28/01%%SYNOPSIS:y=ds_mod(c,x)%c---usercode(columnvector)%x---inputsignal(rowvector)%y---tmp=c*x,y=tmp(:)(ds-ssmodulatedsignal,columnvector)%***********************************************************************************functiony=ds_mod(c,x)tmp=c*x;y=tmp(:);%*************************************************************************************%Thisfunctiongeneratesrandom+1/-1sequencewithindependentidentically%distributedsymbols%%AUTHOR:WenbinLuo%DATE:04/28/01%%SYNOPSIS:x=bingen(L)%L---numberofrandomsymbols%***********************************************************************************functionx=bingen(L)%generateLsymbolsrandomlywithvalue+1or-1x=rand(1,L);x(find(x)=-1;x(find(x=)=1;%*************************************************************************************%ThisfunctiondoestheDS-SSmodulation%%AUTHOR:WenbinLuo%DATE:04/28/01%%SYNOPSIS:x=ds_demod(c,y)%c---usercode(columnvector)%y---tmp=c*x,y=tmp(:)(ds-ssmodulatedsignal,columnvector)%x---inputsignal(rowvector)%***********************************************************************************functionx=ds_demod(c,y)tmp=reshape(y,length(c),length(y)/length(c));tmp=tmp';%xisacolumnvectorx=tmp*c;%converttorowvectorx=x';%*************************************************************************************%ThisfunctiondoestheDS-SSmodulation%%AUTHOR:WenbinLuo%DATE:04/28/01%%SYNOPSIS:y=ds_mod(c,x)%c---usercode(columnvector)%x---inputsignal(rowvector)%y---tmp=c*x,y=tmp(:)(ds-ssmodulatedsignal,columnvector)%***********************************************************************************functiony=ds_mod(c,x)tmp=c*x;y=tmp(:);%***********************************************************%Thismfunctiongeneratesfadedenvelopeandphase%correspondingtoRayleighfading%%AUTHOR:WenbinLuo%DATE:04/27/01%%FUNCTIONSYNOPSIS:%[env,phi]=fade(L,para)%%ParameterDescription:%L:numberofsamplesneeded%variance:variance%**********************************************************function[env,phi]=fade(L,variance)%Errorcheckifvariance=0error('Positivevarianceneeded')elseifnargin~=2error('Insufficientinputparameters')end%GeneratebivariateGaussianuncorrelated%randomvariablesmu=zeros(1,2);C=variance*eye(2,2);r=mvnrnd(mu,C,L);%Converttopolarcoordinatesandcompute%magnitudeandphasez=r(:,1)+j*r(:,2);env=abs(z);phi=angle(z);%**********************************************************%***********************************************************%Thismfunctiongeneratestwochannelsoffaded%envelopeandphasecorrespondingto%Rayleighfading%%AUTHOR:WenbinLuo%DATE:04/27/01%%FUNCTIONSYNOPSIS:%[env,phi]=fade_diversity(L,para)%%ParameterDescription:%L:numberofsamplesneeded%variance:variance%**********************************************************function[env1,env2]=fade_diversity(L,variance)%Errorcheckifvariance=0error('Positivevarianceneeded')elseifnargin~=2error('Insufficientinputparameters')end%GeneratebivariateGaussianuncorrelated%randomvariablesmu=zeros(1,4);C=variance*eye(4,4);r=mvnrnd(mu,C,L);%Converttopolarcoordinatesandcompute%magnitudeandphasez1=r(:,1)+j*r(:,2);z2=r(:,3)+j*r(:,4);env1=abs(z1);env2=abs(z2);%**********************************************************%***********************************************************%Thismfunctiongeneratesfrequencyselective%Rayleighfading%%AUTHOR:WenbinLuo%DATE:05/02/01%%FUNCTIONSYNOPSIS:%y=fade_fs(x,L)%%ParameterDescription:%y:outputsignal%x:inputsignal%L:numberofindependentRayleigh%fadingprocess%**********************************************************functiony=fade_fs(x,L)%GeneratebivariateGaussianuncorrelated%randomvariablestmp1=0:1:(L-1);tmp1=exp(-tmp1);tmp(1:2:2*L-1)=tmp1;tmp(2:2:2*L)=tmp1;mu=zeros(1,2*L);C=*diag(tmp);x_len=length(x);r=mvnrnd(mu,C,x_len);%Converttopolarcoordinatesandcomputemagnitudex=x(:);y=zeros(x_len,1);fori=1:L,z=r(:,2*i-1)+j*r(:,2*i);env=abs(z);%phi=angle(z);tmp_y=env.*x;tmp_y=[zeros(i-1,1);tmp_y(1:x_len-i+1)];y=y+tmp_y;end%**********************************************************%**********************************************************************%ThisprogramcomputestheaverageBERofaDS-SS/BPSK%communicationsystemwithbinaryBCHcodeintheAWGNchannel%%AUTHOR:WenbinLuo%DATE:04/28/01%%%%**********************************************************************%functionPlot_Pe=final11_extra()clearall;%closeall;formatlong;%setupthethresholdVtVt=0;Plot_Pe=[];N=16;x_num=2500;plot_EbNo=-20:2:10;forEbNo=-20:2:10,%convertb

1 / 50
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功