版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、基于matlab的2PSK的系统仿真 摘要:Simulink是Mathworks公司推出的基于Matlab平台的著名仿真环境Simulin作为一种专业和功能强大且操作简单的仿真工具,目前已被越来越多的工程技术人员所青睐,它搭建积木式的建模仿真方式既简单又直观,而且已经在各个领域得到了广泛的应用。本次课程设计是基于MATLAB的2PSK和2DPSK仿真,通过系统分析,步骤来完成本次设计任务。通过课程设计从理论学习的轨道逐步引向实际应用,把理论上熟悉的定性分析、定量计算逐步和工程估算、实验调整等手段结合起来,掌握工程设计的步骤和方法,了解科学实验的程序和实施方法,为以后毕业设计和从事信息处理技术的
2、实际工作打下基础。关键词:MATLAB;2PSK,2DPSK;仿真1.二进制相移键控设计流程图如图2-1所示。图 2-1 2PSK调制解调框图1.1 PSK调制原理在二进制数字调制中,当正弦载波的相位随二进制数字基带信号离散变化时,则产生二进制移相键控(2PSK)信号。2PSK信号调制有两种方法,即模拟调制法和键控法。通常用已调信号载波的 0°和 180°分别表示二进制数字基带信号的 1 和 0,模拟调制法用两个反相的载波信号进行调制。2PSK以载波的相位变化作为参考基准的,当基带信号为0时相位相对于初始相位为0°,当基带信号为1时相对于初始相位为180°
3、;。键控法,是用载波的相位来携带二进制信息的调制方式。通常用0°和180°来分别代表0和1。其时域表达式为:其中,2PSK的调制中an必须为双极性码。本次设计中采用模拟调制法。两种方法原理图分别如图2-2和图2-3所示。图 2-2 模拟调制法原理图图 2-3 键控法原理图2PSK信号的时间波形1.2 PSK解调原理由于2PSK的幅度是恒定的,必须进行相干解调。经过带通滤波的信号在相乘器中与本地载波相乘,然后用低通滤波器滤除高频分量,在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1,负抽样值判为0。2PSK信号的相干解调原理图如图2-4所示,各点的波形如图2-5所示。
4、由于2PSK信号的载波回复过程中存在着180°的相位模糊,即恢复的本地载波与所需相干载波可能相同,也可能相反,这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反,即“1”变成“0”吗“0”变成“1”,判决器输出数字信号全部出错。这种现象称为2PSK方式的“倒”现象或“反相工作”。但在本次仿真中是直接给其同频同相的载波信号,所以不存在此问题。图 2-4 2PSK的相干解调原理图图 2-5 相干解调中各点波形图图中,假设相干载波的基准相位与2PSK信号的基准一致(通常默认为0相位)。但是由于2PSK信号的载波回复过程中存在着180°的相位模糊,即恢复的本
5、地载波与所需相干载波可能相同,也可能相反,这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反,即“1”变成“0”吗“0”变成“1”,判决器输出数字信号全部出错。这种现象称为2PSK方式的“倒”现象或“反相工作”。 2PSK信号在一个码元的持续时间Ts内可以表示为 u1T(t) 发送“1”时 ST(t)= uoT(t)=- u1T(t) 发送“0”时期中 Acosct 0< t < Ts u1T(t)= 0 其他 设发送端发出的信号如上式所示,则接收端带通滤波器输出波形y(t)为 a+nc(t)cosct-ns(t)sinct 发送“1”时 y(t)= -a+
6、nc(t)cosct-ns(t)sinct 发送“0”时y(t)经过想干解调(相乘低通)后,送入抽样判决器的输入波形为 a+nc(t) 发送“1”时 x(t)= -a+nc(t) 发送“0”时由于nc(t)是均值为0,方差为2的高斯噪声,所以x(t)的一维概率密度函数为由最佳判决门限分析可知,在发送“1”和“0”概率相等时,即P(1)=P(0)时,最佳门限b*=0.此时,发“1”而错判为“0”的概率为 P(0/1)=P(x0)=0-f1(x)dx=1/2erfc() 式中:r=a2/22n 同理,发“0而错判为“1”的概率为P(1/0)=P(x0)=0-f0(x)dx=1/2erfc() 故2
7、PSK信号相干解调系统的总误码率为 Pe=P(1)P(0/1)+P(0)P(0/1)= 1/2erfc() 在大信噪比(r>>1)的条件下,上式可近似为 Pee-r/22.基于matlab的2PSK的系统仿真通过编写M文件程序:2PSK调制解调程序及注释clear allclose alli=10;j=5000;fc=4; %载波频率fm=i/5; %码元速率B=2*fm;t=linspace(0,5,j);a=round(rand(1,i); %随机序列,基带信号figure(3);stem(a);st1=t;for n=1:10 if a(n)<1; for m=j/i*
8、(n-1)+1:j/i*n st1(m)=0; end else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end endendfigure(1);subplot(411);plot(t,st1);title('基带信号st1');axis(0,5,-1,2);%由于PSK中的是双极性信号,因此对上面所求单极性信号取反来与之一起构成双极性码st2=t; for k=1:j; if st1(k)>=1; st2(k)=0; else st2(k)=1; endend;subplot(412);plot(t,st2);title('基带信号反
9、码st2');axis(0,5,-1,2);st3=st1-st2;subplot(413);plot(t,st3);title('双极性基带信号st3');axis(0,5,-2,2);s1=sin(2*pi*fc*t);subplot(414);plot(s1);title('载波信号s1');e_psk=st3.*s1;figure(2);subplot(511);plot(t,e_psk);title('e_2psk');noise=rand(1,j);psk=e_psk+noise; %加入噪声subplot(512);plot
10、(t,psk);title('加噪后波形');psk=psk.*s1; %与载波相乘subplot(513);plot(t,psk);title('与载波s1相乘后波形');f,af = T2F(t,psk); %通过低通滤波器t,psk = lpf(f,af,B);subplot(514);plot(t,psk);title('低通滤波后波形');for m=0:i-1; if psk(1,m*500+250)<0; for j=m*500+1:(m+1)*500; psk(1,j)=0; end else for j=m*500+1:(
11、m+1)*500; psk(1,j)=1; end endendsubplot(515);plot(t,psk);axis(0,5,-1,2);title('抽样判决后波形')2DPSK调制解调程序及注释clear allclose alli=10;j=5000;fc=4; %载波频率fm=i/5; %码元速率B=2*fm;t=linspace(0,5,j);a=round(rand(1,i);figure(4);stem(a);st1=t;for n=1:10 if a(n)<1; for m=j/i*(n-1)+1:j/i*n st1(m)=0; end else f
12、or m=j/i*(n-1)+1:j/i*n st1(m)=1; end endendfigure(1);subplot(321);plot(t,st1);title('绝对码');axis(0,5,-1,2);b=zeros(1,i); %全零矩阵b(1)=a(1);for n=2:10 if a(n)>=1; if b(n-1)>=1 b(n)=0; else b(n)=1; end else b(n)=b(n-1); endendst1=t;for n=1:10 if b(n)<1; for m=j/i*(n-1)+1:j/i*n st1(m)=0; e
13、nd else for m=j/i*(n-1)+1:j/i*n st1(m)=1; end endendsubplot(323);plot(t,st1);title('相对码st1');axis(0,5,-1,2);st2=t; for k=1:j; if st1(k)>=1; st2(k)=0; else st2(k)=1; endend;subplot(324);plot(t,st2);title('相对码反码st2');axis(0,5,-1,2);s1=sin(2*pi*fc*t);subplot(325);plot(s1);title('
14、载波信号s1');s2=sin(2*pi*fc*t+pi);subplot(326);plot(s2);title('载波信号s2');d1=st1.*s1;d2=st2.*s2;figure(2);subplot(411);plot(t,d1);title('st1*s1');subplot(412);plot(t,d2);title('st2*s2');e_dpsk=d1+d2;subplot(413);plot(t,e_dpsk);title('调制后波形');noise=rand(1,j);dpsk=e_dpsk+
15、noise; %加入噪声subplot(414);plot(t,dpsk);title('加噪声后信号');dpsk=dpsk.*s1; %与载波s1相乘figure(3);subplot(411);plot(t,dpsk);title('与载波相乘后波形');f,af=T2F(t,dpsk); %通过低通滤波器t,dpsk=lpf(f,af,B);subplot(412);plot(t,dpsk);title('低通滤波后波形');st=zeros(1,i); %全零矩阵for m=0:i-1; if dpsk(1,m*500+250)<
16、0; st(m+1)=0; for j=m*500+1:(m+1)*500; dpsk(1,j)=0; end else for j=m*500+1:(m+1)*500; st(m+1)=1; dpsk(1,j)=1; end endendsubplot(413);plot(t,dpsk);axis(0,5,-1,2);title('抽样判决后波形')dt=zeros(1,i); %全零矩阵dt(1)=st(1);for n=2:10; if (st(n)-st(n-1)<=0&&(st(n)-st(n-1)>-1; dt(n)=0; else dt
17、(n)=1; endendst=t;for n=1:10 if dt(n)<1; for m=j/i*(n-1)+1:j/i*n st(m)=0; end else for m=j/i*(n-1)+1:j/i*n st(m)=1; end endendsubplot(414);plot(t,st);axis(0,5,-1,2);title('码反变换后波形');%利用FFT计算信号的频谱并与信号的真实频谱的抽样比较。%脚本文件T2F.m定义了函数T2F,计算信号的傅立叶变换。function f,sf= T2F(t,st)%This is a function using
18、 the FFT function to calculate a signal's Fourier%Translation%Input is the time and the signal vectors,the length of time must greater%than 2%Output is the frequency and the signal spectrumdt = t(2)-t(1);T=t(end);df = 1/T;N = length(st);f=-N/2*df:df:N/2*df-df; sf = fft(st);sf = T/N*fftshift(s);用
19、到的低通滤波器函数function t,st=lpf(f,sf,B)%This function filter an input data using a lowpass filter%Inputs: f: frequency samples% sf: input data spectrum samples% B: lowpass's bandwidth with a rectangle lowpass%Outputs: t: time samples% st: output data's time samplesdf = f(2)-f(1);T = 1/df;hf = zeros(1,length(f); %全零矩阵bf = -floor( B/df ): floor( B/df ) + floor( l
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 石河子大学《园林植物栽培养护》2023-2024学年第一学期期末试卷
- 大学生个人实习总结集合3篇
- 石河子大学《饲料学》2022-2023学年第一学期期末试卷
- 石河子大学《律师实务》2023-2024学年第一学期期末试卷
- 石河子大学《程序设计》2022-2023学年期末试卷
- 沈阳理工大学《模拟电路基础》2021-2022学年期末试卷
- 沈阳理工大学《机械设计》2022-2023学年第一学期期末试卷
- 银屑病的辩证施护
- 沈阳理工大学《复变函数与积分变换》2021-2022学年第一学期期末试卷
- 骨灰安放合同
- 起重机械自查报告
- 2021年至2023年广东省公务员遴选笔试真题、面试真题及答案解析(各地市、省直共12套)
- ZJ40J钻机技术参数
- 提高冠脉介入手术术前准备的合格率
- 《学习态度与习惯》主题班会课件
- 创建国家级旅游度假区自评报告
- 英语1-基础模块-unit3-Shopping-教案
- 水池防腐涂层施工方案范本
- 路面水稳层施工方案(完整版)
- 沉井下沉监测方案
- 2023年医养结合工作情况汇报范文(4篇)
评论
0/150
提交评论