matlabBPSK调制与解调clc;closeall;clear;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 假定:% 2倍载波频率采样的bpsk信号% 调制速率为在波频率的N/2m% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%m=128;N=512;n=1:1:N;N0=0.5*randn(1,N)%噪声h0=zeros(1,N);% 30阶低通滤波器h0f=[00.30.31];w0=[1100];b=fir2(30,f,w0);[h,w]=freqz(b,1,N/2);h0(1,1:N/2)=abs(h');fori=1:N/2h0(1,N-i1)=h0(1,i);end;%%%%%%%%% 随机序列a=rand(1,m);fori=1:mif(a(1,i)>0.5)a(1,i)=1; elsea(1,i)=-1;end;end;a%%%生成BPSK信号bpsk_m=zeros(1,N);j=1;k=1;fori=1:Nif(j==(N/m1))j=1;k=k1;end; %0.05*pi为初始相位,可以任意改变bpsk_m(1,i)=a(1,k)*sin(2*pi*0.5*i0.05*pi)a(1,k)*cos(2*pi*0.5*i0.05*pi);j=j1;end;bpsk_m=bpsk_mN0;%信号加噪声,模拟过信道%接收处理 用正交本振与信号相乘,变频bpsk_m1=bpsk_m.*sin(2*pi*0.5*n);bpsk_m2=bpsk_m.*cos(2*pi*0.5*n);%滤波tempx=fft(bpsk_m1);tempx=tempx.*h0; %低通滤波tempx=ifft(tempx);real_x=real(tempx);tempx=h0.*fft(bpsk_m2);tempx=tempx.*h0; %低通滤波tempx=ifft(tempx);real_x1=real(tempx);subplot(2,1,1);plot(real_x1real_x,'b'); axis([1 N-2.52.5]);gridon;holdon;In=real_x1real_x; %可只取一路,这里取了两路之和fori=1:N %滤波后整形if(In(1,i)>0) %判决,得到解调结果 In(1,i)=1;elseIn(1,i)=-1;end;end;plot(In,'r');an=zeros(1,m); fori=1:man(1,i)=In(1,(i-1)*N/mN/(2*m));end;subplot(2,1,2); % 比较误码plot(an,'r*');holdon;axis([1 m-22]);plot(a,'b^');
本文档为【matlab BPSK 调制与解调】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。