BPSK的matlab仿真原理加源代码

BPSK的matlab仿真原理加源代码BPSK误码率的Matlab仿真实验作者：黄准南京航空航天大学，江苏省南京市430062摘要：对于刚刚从事通信领域研究性工作的研究生来说，很有必要通过Matlab强大的仿真功能来强化理解通信信号处理的过程。本文正是通过对BPSK的Matlab仿真来提供一个比较好的实例。关键词：BPSK，高斯白噪声，带通滤波器，载波，抽样判决1BPSK的理论分析1.1基本原理图1.1BPSK相干解调原理框图输入信号A(t)可表示为：A(t)，€a*g(t-NT)nbN(1.1-1)式中g(t-NT)为矩形波，而a表达式如下：bn+I•槪率为P-1,概率为(1-n对A(t)进行傅里叶变换，得到A(t)的功率谱密度P(f)(式中Ts为码元时间长度):sP3胡監［瞬］『=s(t)€A(t)*coswtcn(z)=nc(t)cosw^t-nH(f)sinwct经过带通滤波器后的信号为y(t)y(t)=s(t)+n(t)与载波相乘后再经过低通滤波器滤除高频分量得到x(t)『+%(门 发送“广符号x(t)-jl"+吆⑴发送符号

附录：clearallclcnum=5000;Rb=1000;tnum=500;fs=Rb*tnum;fp=100000;t=0:1/fs:num/Rb-1/fs;panjue1=zeros(1,num);panjue2=zeros(1,num);pe1=zeros(1,15);pe2=zeros(1,15);a=2*(fp-2*Rb)/fs;b=2*(fp+2*Rb)/fs;s=1;%[a1,b1]=butter(2,[a,b]);母参数noise=randn(1,num*tnum);nwgnnoise=filter(a1,b1,noise);stdnwgn=std(nwgnnoise);[a2,b2]=butter(2,b);[h1w1]=freqz(a1,b1);[h2w2]=freqz(a2,b2);%码元数%码速率2*Rb<=fs%码元数%码速率2*Rb<=fs%没一个码元用tnum个点代替%采样频率fs>2*(fp+Rb)%载波频率%时间t的长度与间隙。%基带传输的判决输出%频带传输的判决输出%基带传输的误码率矩阵%频带传输的误码率矩阵%带通滤波器的第一截止频率%带通滤波器的第二截止频率%循环累加数%采用butter滤波器，得到butter滤波器%产生正态分布参数%对WGN进行带通滤波%为了减少运算量。%低通滤波%带通滤波器频谱%低通滤波器频谱%forSNR_db=-5:10;in=randint(1,num);in1=2*in-1;g=ones(1,tnum);inn=upsample(in,tnum);inn=filter(g,1,inn);in2=inn*2-1;inn=in2.*sin(2*pi*fp*t);noise=randn(1,num*tnum);snr=10^(SNR_db/10);sigma=sqrt(1/(2*snr));s2=in2+noise*sigma;s1=inn+noise*sigma;%以上一段代码产生基带信号%产生频带信号%产生标准正态分布的高斯白噪声信号%信噪比换算%求出高斯白噪声标准差%基带传输%频带传输，暂未用上%带通滤波器v=filter(a1,b1,inn);nwgnnoise=filter(a1,b1,noise);高斯信号的功率进行处理nwgnnoise=nwgnnoise/stdnwgn;nwgnnoise=nwgnnoise*sigma;v=v+nwgnnoise;v2=v.*sin(2*pi*fp*t);%对载波信号进行带通滤波%对WGN进行带通滤波,此处分开滤波是为了对窄带%使得滤波后的窄带高斯噪声任然为正态分布%改变噪声功率使其符合预设的信噪比。%经过带通滤波器后的信号%与载波相乘v=filter(a2,b2,v2);%进行低通滤波i=1;whilei<=numifs2(i*tnum-tnum/2)>0panjue1(i)=1;%进行抽样判决%基带抽样判决elsepanjue1(i)=-1;endifv(i*tnum-tnum/2)>0panjue2(i)=1;%频带抽样判决elsepanjue2(i)=-1;endi=i+1;endpe1(s)=length(find(panjue1~=in1))/num; %基带误码率pe2(s)=length(find(panjue2~=in1))/num; %频带误码率s=s+1;endfigure(1);subplot(3,2,1);pwelch(in2,[],[],[],fs);title('基带信号频谱')；figure(1);subplot(3,2,2)；pwelch(inn,[],[],[],fs)；title('频带信号')；figure(1)；subplot(3,2,4)；plot(w1,h1)；title('带通信号频谱')；figure(1)；subplot(3,2,3)；plot(w2,h2)；title('低通信号频谱')；figure(1)；subplot(3,2,5)；pwelch(noise,[],[],[],fs)；title('高斯白噪声信号频谱')；figure(1)；subplot(3,2,6)；pwelch(nwgnnoise,[],[],[],fs)；title('窄带高斯信号频谱')；figure(3)；SNR_db=-5:10；pelilun=0.5*erfc(sqrt(10.A((SNR_db)./10)));semilogy(SNR_db,pe1,'-or')；holdon；semilogy(SNR_db,pe2,'-ob')；holdon；