MSK调制MATLAB仿真程序 MSK为相位连续的FSK调制，其基本原理如图所示： 方案一：按照上图的流程编写代码：
根据上图a和d两组码元的关系，不难看出这一部分应属于差分编码，在调制系统中，有两类调制需要用到，一个是DPSK，另一个就是MSK，仔细观察两种调制的关系，可以发现，不是同一种差分编码，而且好多资料里也写的是MSK差分编码是a(i)和d(i-1)异或。其实不是这样，正确关系应该是d(1)=a(1);d(i)是a(i)和d(i-1)同或结果。这里可以去书上查看相应a和d组码元，自己计算一下就明白了。
%方法一：计算相位变量 c=zeros(1,length(a)); c(1)=0; for i=2:length(a) if a(i)==a(i-1) c(i)=c(i-1); else a(i)~=a(i-1); c(i)=c(i-1)+(a(i-1)-a(i))*((i-1)pi/2); end end I=cos©; %I(k) Q=a.cos©; %-Q(k) %扩展码元，每个码元采样点数为number1，一共number个码元 I1=zeros(1,numbernumber1); Q1=zeros(1,numbernumber1); for i=1:number; I1((i-1)number1+1:inumber1)=I(i); Q1((i-1)number1+1:inumber1)=Q(i); end
%加权系数 t1=1/fs:1/fs:length(a)tb; I=I1.cos(pit1/2/tb); Q=Q1.sin(pit1/2/tb); I和Q的波形如图所示： 然后进行中频搬移即可： multi=4; I=interp(I,multi); Q=interp(Q,multi); t=1/fs:1/fs:length(I)1/fs; I=I.cos(2pifct); Q=Q.sin(2pifct); out=I-Q;
方案二：利用ak和相位的联系，如下图： for i=1:number1 data_sample(i:number1:numbernumber1)=a; end %计算相位 phase=zeros(1,numbernumber1); phase(1)=0;%phase(1)=data_sample(1)pi/2/number1; for i=2:numbernumber1 phase(i)=phase(i-1)+data_sample(i-1)pi/2/number1; end %phase=phase+(0-phase(1)); I=cos(phase); Q=sin(phase); [r1,r2]=fen(a); r(1)=-r2pi; for i=2:(r1+r2)2+1 r(i)=-r2pi+(pi/2)(i-1); end figure; subplot(2,1,1);stem(a);axis([0,length(a)+1,-2,2]);grid; subplot(2,1,2);plot(phase);grid; set(gca,‘YTick’,-r2pi:pi/2:r1*pi); set(gca,‘YTickLabel’,{r/pi}); title(‘相位图’);ylabel(‘pi’); 结果如图所示： 然后再进行调制。 关于这一部分，有的资料上计算相位phase时，初始化令phase(1)=data_sample(1)*pi/2/number1; 我仿真之后发现是有误差的，I和Q路的图像也是有误差，经过仔细检查发现，这里phase(1)=0，也就是第一个相位值应该为1。具体可以自己仿真一下试试。
