MSK调制解调实现.docx

4.2最小频率键控(MSK)4.2.1 MSK基本原理MSK信号是一种相位连续、包络恒定并且占用带宽最小的二进制正交FSK信号它的 第k个码元可以表示为：s (t )= cos ( w t + a t + e) (k - 1)T < t < kT (2-1)k s 2 T k式中，w = 2吋为视在角载频；a二±1 ； T为第k个码元确定的初始相位s s k由上式可以看出，当a = +1时，码兀频率f等于f +1/4T ；当a = —1时，码兀频k 1 s k率f等于f — 1/4T故f和f的距离等于1/2T这是2FSK信号最小频率间隔0 s 1 0式(2-1)可以用频率为 f 的两个正交分量表示将式(2-1)进行三角公式变换，得 s到：s (t) = p cos cos w t 一 qk k 2T s k.Kt .sin sin w t2T s(k —1)T < t < kT3-1)式中，p = cos 申=±1kkq = a cos 申=±1k k k3-2)式(2-2)表示，此MSK信号可以分解为同相分量(I)和正交分量(Q)两部分 MSK信号的调制由式(2-2)可知，MSK信号可以用两个正交的分量表示。

根据该式构成 的MSK信 号的产生方框图如图 2-2所示Pk cos(兀 t / 2T)qkq^siniKt / 2T)q sin(兀 t / 2 T ) sin w tk sp cos(兀 t / 2 T ) cos w tk s图2-2 MSK调制原理图MSK信号的解调由于MSK信号是最小二进制FSK信号，所以它可以采用解调FSK信号的相干法和非 相干法解调图2-3是MSK信号的解调原理框图[2iT,(2i+l)T]MSK信号解 调 输 出[(2i-1)T,(2i+1)T]图2-3 MSK信号的解调原理图4.2.2 MSK仿真实现过程1.3设计步骤先定义MSK输入信号的参数然后用dmod函数做输入的调制，之后画相应的频谱图 再用dedmod做解调输出然后加入噪声，再进行解调输出1.定义载波频率Fs,输出信号频率Fd，采样频率Fs2.3.4.5.设置输入信号X,同时做信号的调制 画出相应的波形图做输出信号的解调波形，同时画解调波形图 加入高斯白噪声再进行信号的解调输出，画相应的波形6. 对结果进行比较分析MSK信号的产生因为MSK信号可以用两个正交的分量表示:/、 兀t .兀t .s (t) = p cos cos w t - q sin sin w tk k 2T c k 2T ss S(k—1) TsVtWkTs式中：右端第1项称作同相分量，其载波为cos w t ；c第2项称作正交分量，其载波为sin w t。

c在仿真时先设置输入信号的参数：载波频率Fc调制后，数字信号速率Fd,模拟信号 采样率Fs之后产生要调制的随机数字信号x;然后用dmod函数进行调制，产生调制信 号dmod函数是MATLAB中一个常用来进行信号调制的函数，它后面的参数包括被调制 信号，载波信号的频率，输出信号的速率以及采样速率和所进行调制的函数部分程序代码和仿真图如下所示：M=2;Fc=20;Fd=10;Fs=200;%定义输入随机信号的参数%载波信号的频率Fc为20；%输出信号的速率Fd为10；%采样速率Fs为200；x二randin t( 1,10000);y二dmod(x,Fc,Fd,Fs, 'msk');调制前的输入波形%产生随机的输入序列%对用dnod函数对输入的随机序列进行调制10.50212144 6 8 10调制后的输出液形图3-1 MSK调制信号的波形由调制波形图可以看出MSK的调制信号特性与2FSK调制信号的特性非常的相似，即: 当输入信号为“ 1”时，调制后的波形比输入信号为“0”时的波形要密同时MSK信号 的包络是恒定的，相位则是连续的带宽相对于一般的2FSK信号要小，而且正交MSK解调实现由于MSK信号是最小二进制FSK信号，所以它可以采用解调FSK信号的相干法和非相 干法解调。

在进行程序仿真时，用ddemod函数进行调制信号的解调同时画出解调前后的时域 与频域的波形图在画频域的波形图时先对已调信号与解调信号进行DTFT变换，之后画 出相应的波形ddemod函数是与dmod函数相对应，用来对已调信号进行解调的其后面 的参数与dmod函数后的一模一样部分程序与仿真波形图如下所示：z二ddemod(y,Fc,Fd,Fs, 'msk' ); %对调制后的MSK信号进行解调解调前的波形%定义频域参数nl为y的长度，n2为 输入原始信号的长度%对已调信号进行DTFT变换%对解调信号进行DTFT变换图3-2 MSK信号的解调波形由解调信号的时域波形可以看出，解调后的波形与原始输入信号的波形完全一致 同时不难发现解调后的信号很稳定 n1=1:1:length(y);w二[0:1:1000]*pi/500; n2=1:1:10000;w二[0:1:1000]*pi/500;Y二y* exp(-j *n1，* w); Z=z *exp(_j *n2，* w); magY二abs(Y); magZ二abs(Z);解调前的频域波形300200 - -100 1 -口 , ， J0 0.5 1 1.5 2图3-3 MSK信号解调前后的频域波形由解调后的频域波形可以看出MSK信号的稳定性很好，说明MSK信号的抗噪声性能很 强。

100 0.5 1 1.5 2图3-4调制前解调后的频域波形通过调制前和解调后的频域波形比较，发现解调非常的成功，波形被完整无误的输 出来叠加噪声的MSK解调由于信号在传输的过程中难免要受到外来噪声的影响，所以在进行通信仿真时务必 要在理想的模拟通信系统中加入高斯白噪声对系统进行影响，以此来判断一个通信系统 抗噪声性能的好坏利用awagn函数对已调信号加入信噪比为0.001的高斯白噪声然后对加入了噪声的已调信号进行解调，同时画出其时域与频域的波形，将其与没有加噪声的调制信号进行 比较，不难发现：加入噪声后对已调信号的影响很大，但对解调信号而言，噪声对解调 后的影响还是有一定限度的这说明MSK解调系统的抗噪声性能很强仿真的部分程序与相应的仿真波形图如下所示：yl二awgn(y,0.001); %对已调信号加入高斯白噪声zl二ddemod(y,Fc,Fd,Fs, 'msk'); %对加入高斯白噪声的调制信号进行解调[num1,ra ti o1]=symerr(x,z2) %对加入噪声后解调的信号进行误码率的测定通过symerr函数对解调信号的误码率进行测定，得到的误码率为0.0024.因为输入 的随机序列长度为10000,错了24个。

说明MSK信号的抗噪声性能很强当将信噪比变为 1时，输出的误码率很快发生了变化，由原来的0.0024变为了6.0000e-004.图形也同时发 生了变化通过图形的比较，可以发现噪声对调制信号的影响很大，但通过解调后影响 相对减小了很多加噪声后解调前的波形加嗥芦后解调后的波形图3-5信噪比为0.001的MSK信号加噪声解调前后的时域波形加噪声后解调前的液形加噪声后解调后的液形Y1=y1*exp(-j*n1' *w);Z1=z1 *exp(-j*n2' *w);图3-6信噪比为1的MSK信号加噪声解调前后的时域波形%对加噪声后的已调信号进行DTFT变换%对加噪声后的解调信号进行DTFT变换magYl二abs(Yl); magZl二abs(Zl);图3-7信噪比为2的MSK信号加噪声解调前后的频域波形在多次改变信噪比后，将信噪比与误码率联合做了如图3-8所式的波形图图3-8信噪比与产生误码率曲线图通过图形可知，随着系统信噪比的增加，解调后输出的误码率随之减小，当信噪比 增加到一定程度时，误码率变为0说明噪声对信号的影响在一定程度上是可以减小， 但很难消除，因为每一个系统都不可能是理想系统，总会受到外界的干扰。