基于MATLAB的PSK解调Word下载.doc
《基于MATLAB的PSK解调Word下载.doc》由会员分享,可在线阅读,更多相关《基于MATLAB的PSK解调Word下载.doc(19页珍藏版)》请在冰豆网上搜索。
本课程设计主要是实现2PSK信号的解调,在Matlab6.5环境下用M文件来编制程序并且分别绘制出有无噪声时解调前后在时域和频域波形并进行对比以及进行频谱分析,加深对2PSK解调原理的掌握。
1.1课程设计目的
本课程设计是基于Matlab6.5中用M文件来实现PSK信号的解调,通过该课程设计可以加深我对PSK解调原理的理解以及通过分析解调前后的频谱图更好的掌握PSK解调的原理和特性,对以后学习其它信号的解调过程奠定了基础,提高了自己的实践动手能力和独立思考问题能力。
1.2课程设计要求
1)熟悉Matlab6.5集成环境下M文件的使用方法,通过理解和掌握PSK信号解调原理,编写程序来实现PSK信号的解调。
2)分别绘制出PSK信号解调前后在时域和频域中的波形图,并观察解调前后频谱有何变化并分析频谱变化原因,从而可以加深对PSK信号解调原理的理解和掌握。
3)对已调信号叠加噪声,并进行解调,分别绘制出解调前后信号的时频波形图,
从而比较叠加噪声(高斯白噪声)和未叠加噪声时解调后信号的波形和频谱有何区别,进而分析噪声对信号解调所造成的影响。
4)独立的完成课程设计的全部内容,能正确阐述和分析设计和实验结果。
1.3课程设计步骤
首先是产生随机数字基带信号,接着用调制函数dmod实现PSK调制产生PSK信号(简单介绍),再对PSK信号用解调函数ddemod实现PSK解调,且绘制出解调后时频域的波形和频谱,并与调制前信号(即原始数字基带信号)进行比较,给出相应的分析,最后再在解调前叠加上噪声(高斯白噪声),然后再进行解调,绘制出解调后在时频域的波形和频谱,并与无噪声时进行比较分析,并且再叠加大噪声(信噪比小)和叠加小噪声(信噪比大)时比较时频域的波形和频谱,分析噪声对信号的影响。
2PSK解调原理
2PSK信号可以看成是双极性数字基带信号作用下的DSB调幅信号[2]。
2PSK的信号形式一般表示为:
(2-1)
2PSK的典型波形如图2-1所示:
图2-12PSK信号的典型波形
因为2PSK信号是属于DSB信号,所以对于它的解调不能够采用包络检测的方法,而只能够采用相干解调和鉴相法解调方法进行解调,其解调的框图如图2-2所示:
鉴相法解调
相干解调
解调器
y(t)z(t)x(t)
图2-22PSK信号的解调框图
而主要是采用相干解调的方法对PSK信号进行解调,工作原理简要分析如下:
当不考虑噪声时,带通滤波器输出可以表示为:
y(t)=cos(wct+n)(2-2)
式中n为2PSK信号的某一个码元的初相。
当n=0时,代表数字“0”;
而当n=时,代表数字“1”。
再通过乘法器与同步载波coswct相乘后,输出为:
z(t)=cos(wct+n)coswct=cosn+cos(2wct+n)(2-3)
接着经过低通滤波器滤除高频分量,可以得到解调器的输出为:
x(t)=cosn(2-4)
当n=0时,x(t)=;
当n=1时,x(t)=-。
根据输入端产生2PSK信号时n(0或)代表数字信息(“1”或“0”)的规定,以及接收端x(t)与n关系的特性,所以抽样判决器的判决准则应该为:
当x0时,判为“0”;
当x<
0时,判为“1”。
其中x为x(t)在抽样时刻的值。
由此可见,2PSK信号的相干解调的过程实际上就是对输入的已调信号与本地载波信号进行极性比较的过程,故2PSK信号的相干解调常称为极性比较法解调,也称为同步检测法解调[5]。
3PSK解调的实现
3.1PSK信号的产生
由原理上说2PSK信号的产生方法有模拟调制法和键控法。
本课程设计是在Matlab6.5集成环境下用M文件在产生基带信号的基础上调用数字调制函数dmod设置参数Fc,Fd,Fs,在Matlab6.5下具体产生PSK信号的实现为:
y=dmod(x,Fc,Fd,Fs,’psk’,M)(M=2表示二进制)
Fc表示载波频率,Fd表示数字基带信号与已调信号取样,Fs表示采样频率。
3.2PSK解调实现
本课程设计由原理上说是对数字基带信号在无噪声干扰的情况下采用相干解调方法对2PSK信号进行解调,具体是通过Matlab6.5集成环境下用M文件调用数字解调函数ddemod编制程序来实现PSK信号的解调。
首先是随机产生长度为32的二进制数字基带信号,在Matlab6.5下具体实现为:
x=randint(32,1,M)(M=2表示二进制)
然后是用调制函数dmod产生PSK信号,其中的参数Fc,Fd,Fs可以设置为:
Fc=10;
%载波频率
Fd=20;
%数字基带信号与已调信号取样
Fs=500;
%采样频率
Fc,Fd,Fs之间的关系为:
Fc/Fd=1/2(表示基带信号的一个码元对应的PSK信号波形有两个正弦载波),Fs/Fc=50(表示在一个正弦载波中采样50个点)。
再用解调函数ddemod对已调信号进行PSK解调,具体实现为:
z=ddemod(y,Fc,Fd,Fs,’psk’,M)(M=2表示二进制)
用图形窗口分割函数subplot[3]和绘制离散序列图函数stem[3]对比绘制出调制前数字基带信号时域波形和解调后信号的时域波形,如图3-1所示:
图3-1数字基带信号PSK调制前与解调后的时域波形
由图3-1对比波形图进行分析:
输入的数字基带信号为:
10011001011100011001000010101010
解调后的输出信号为:
10011001011100011001000010101010
由PSK解调原理可知在对已调信号不叠加噪声再进行解调所得到的输出信号与输
入的数字基带信号相一致,不会产生误码。
接着用fft函数对数字基带信号和解调后的信号进行傅里叶变换,再用图形窗口分割函数subplot[3]和绘制二维图形函数plot[3]对比绘制出调制前数字基带信号频谱和解调后信号的频谱,如图3-2所示:
图3-2数字基带信号调制前与解调后的频谱图
由图3-2对比频谱图进行分析:
同样与时域波形图分析一样,在无噪声的情况下,理论上进行解调后不会产生误码,所以数字基带信号调制前与经PSK解调后的频谱图是相一致的,而且该频谱图是关于x轴的正负半轴是对称的。
由上述的时频域图对比分析可以知道当无噪声时对PSK信号进行解调所得到的实际结果与理论上的分析是相一致的。
3.3叠加噪声时的PSK解调
因为在信道中噪声会对信号的传输有影响,主要是高斯白噪声,所以需要叠加噪声(主要是高斯白噪声)来分析噪声对信号影响。
为了能够更好的研究噪声对PSK解调信号的影响,可以用函数randint[3]产生长度为64的二进制随机数字基带信号,因为输入的数字基带信号的码元越多就可以更好的研究误码率。
现在先对已调信号叠加高斯白噪声,再进行PSK解调。
因为高斯白噪声是一组高
斯分布的随机数,所以可以用函数awgn[4]来实现叠加高斯白噪声,具体实现为:
SRNperBit=6[4];
(信噪比为6)
adjSRN=SRNperBit-10*log10(Fs/Fd)+10*log10(log2(M))[4];
yniose=awgn(y,adjSRN,'
measured'
,[],'
dB'
)[4];
然后再用解调函数ddemod对加入高斯白噪声的已调信号进行PSK解调。
同样用图形窗口分割函数subplot[3]和绘制离散序列图函数stem[3]对比绘制出调制前数字基带信号时域波形,解调后信号的时域波形以及加入高斯白噪声解调后信号的时域波形,如图3-3所示:
图3-3数字基带信号PSK调制前与有无噪声解调后的时域波形
由图3-3对比波形图进行分析:
1000110110110000010011000110101011111010111000100110000001000010
无噪声解调后的输出信号为:
加入高斯白噪声解调后的输出信号为:
1000110110110000010011000110101011111110011000100110000001000010
由此可知对已调信号叠加高斯白噪声再进行PSK解调,得到PSK信号解调波形与原来的基带数字信号(原始信号)以及无噪声的PSK信号解调波形相对比产生了误码。
其误码Ne=2;
N=64。
由公式Pe=Ne/N,可得出实际误码率为0.03125。
根据PSK信号加入噪声后的解调原理由理论上分析:
当发送“1”时,由于叠加了高斯白噪声,使抽样值在抽样判决时刻变为小于0时,将发生“1”错判为“0”的误码,其概率为:
(x<
0,发送“1”时)[1](3-1)
当发送“0”时,同理由于叠加了高斯白噪声,使抽样值在抽样判决时刻变为大于0时,将发生“0”错判为“1”的误码,其概率为:
(x>
0,发送“0”时)[1](3-2)
因为Pe1=Pe2,故通过计算可得2PSK信号采用同步检测法(相干解调法)时的系统总误码率为:
[1](3-3)
在大信噪比下,上式成为:
[1](3-4)
式中r表示信噪比。
由于加入的高斯白噪声的信噪比为6(噪声较大),所以产生误码的概率也比较大。
与理论上的分析是相符合的,r越小,则理论上的误码率越大。
接着同样用f