BPSK的matlab仿真原理加源代码(12页).doc

上传人:1595****071 文档编号:34795557 上传时间:2022-08-18 格式:DOC 页数:12 大小:121.50KB
返回 下载 相关 举报
BPSK的matlab仿真原理加源代码(12页).doc_第1页
第1页 / 共12页
BPSK的matlab仿真原理加源代码(12页).doc_第2页
第2页 / 共12页
点击查看更多>>
资源描述

《BPSK的matlab仿真原理加源代码(12页).doc》由会员分享,可在线阅读,更多相关《BPSK的matlab仿真原理加源代码(12页).doc(12页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、-BPSK误码率的Matlab仿真实验作者:黄准南京航空航天大学,江苏省 南京市 430062摘要:对于刚刚从事通信领域研究性工作的研究生来说,很有必要通过Matlab强大的仿真功能来强化理解通信信号处理的过程。本文正是通过对BPSK的Matlab仿真来提供一个比较好的实例。关键词:BPSK,高斯白噪声,带通滤波器,载波,抽样判决1 BPSK的理论分析1.1 基本原理加法器s(t)相乘器BPF相乘器LPF抽样判决A(t)WGN定时脉冲A(t)图1.1 BPSK相干解调原理框图输入信号A(t)可表示为: (1.1-1)式中为矩形波,而表达式如下:对A(t)进行傅里叶变换,得到A(t)的功率谱密度

2、(式中Ts为码元时间长度):得到A(t)的频谱图如1.2所示(图中Tb为码元时间长度)图1.2 矩形基带信号频谱A(t)与载波相乘:其频谱图如图1.3所示:图1.3图中fc为载波频率,fs=1/Tb=Rb(码速率),故带通滤波器的中心频率为fc,带宽至少为2fs。经过带通滤波器的高斯白噪声变成了窄带高斯信号,可以表示为:经过带通滤波器后的信号为y(t)y(t)=s(t)+n(t)与载波相乘后再经过低通滤波器滤除高频分量得到x(t)由x(t)的功率谱密度图(图1.4所示)可知,最佳判决门限为0;误码率Pe=0.5*erfc(sqrt(r);式中信噪比r=A2/2*2,2为窄带高斯信号的方差,由于

3、其均值为0,故噪声功率就是其方差。图1.4附录:clear allclcnum = 5000; %码元数Rb = 1000; %码速率2*Rb2*(fp+Rb)fp = 100000; %载波频率t = 0:1/fs:num/Rb-1/fs; %时间t的长度与间隙。panjue1 = zeros(1,num); %基带传输的判决输出panjue2 = zeros(1,num); %频带传输的判决输出pe1 = zeros(1,15); %基带传输的误码率矩阵pe2 = zeros(1,15); %频带传输的误码率矩阵a = 2*(fp-2*Rb)/fs; %带通滤波器的第一截止频率b = 2*

4、(fp+2*Rb)/fs; %带通滤波器的第二截止频率s = 1; %循环累加数%a1,b1 = butter(2,a,b); %采用butter滤波器,得到butter滤波器的H(z)的分子分母参数noise = randn(1,num*tnum); %产生正态分布参数nwgnnoise = filter(a1,b1,noise); %对WGN进行带通滤波stdnwgn = std(nwgnnoise); %为了减少运算量。a2,b2 = butter(2,b); %低通滤波h1 w1 = freqz(a1,b1); %带通滤波器频谱h2 w2 = freqz(a2,b2); %低通滤波器频

5、谱%以上一段代码是为了防止加重计算机的计算量而从循环中提出来的。%for SNR_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*s

6、nr); %求出高斯白噪声标准差s2 = in2+noise*sigma; %基带传输s1 = inn+noise*sigma; %频带传输,暂未用上%带通滤波器v = filter(a1,b1,inn); %对载波信号进行带通滤波nwgnnoise = filter(a1,b1,noise); %对WGN进行带通滤波,此处分开滤波是为了对窄带高斯信号的功率进行处理nwgnnoise = nwgnnoise/stdnwgn; %使得滤波后的窄带高斯噪声任然为正态分布nwgnnoise = nwgnnoise*sigma; %改变噪声功率使其符合预设的信噪比。v = v+nwgnnoise; %

7、经过带通滤波器后的信号v2 = v.*sin(2*pi*fp*t); %与载波相乘v = filter(a2,b2,v2); %进行低通滤波i = 1; %进行抽样判决while i 0 %基带抽样判决 panjue1(i) = 1; else panjue1(i) = -1; end if v(i*tnum-tnum/2) 0 %频带抽样判决 panjue2(i) = 1; else panjue2(i) = -1; end i = i+1;end pe1(s) = length(find(panjue1=in1)/num; %基带误码率 pe2(s) = length(find(panju

8、e2=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(高斯白噪声信号频谱);fi

9、gure(1);subplot(3,2,6);pwelch(nwgnnoise,fs);title(窄带高斯信号频谱);figure(3);SNR_db = -5:10;pelilun = 0.5*erfc(sqrt(10.(SNR_db)./10);semilogy(SNR_db,pe1,-or);hold on;semilogy(SNR_db,pe2,-ob);hold on;semilogy(SNR_db,pelilun,-og);legend(基带误码率曲线,频带误码率曲线,理论误码率曲线);grid on;hold on;xlabel(SNR_db);ylabel(pe);第 12 页-

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 教育专区 > 单元课程

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁