MSK、GMSK调制及相干解调Matlab仿真.pdfVIP

MSK、GMSK调制及相干解调Matlab仿真

实验目的：

1.掌握MSK调制、相干解调原理及特性；

2.了解MSK调制与GMSK调制的差别。

实验内容：

1.编写MATLAB程序仿真MSK调制及相干解调；

2.观察I、Q两路基带信号的特征及与输入NRZ码的关系；

3.观察I、Q调制解调过程中信号的变化；

4.对程序做修改，进行GMSK调制及解调仿真；

5.分析仿真中观察的数据，撰写实验报告。

仿真代码：

clearall

clc

globaldtdftfN

closeall

pi=3.1415926;

fc=5;

N=2^8;

L=8;

M=N/L;

Rb=2;

Tb=1/Rb;

dt=Tb/L;

df=1/(dt*N);

T=N*dt;

B=N*df/2;

t=[-T/2+dt/2:dt:T/2];

f=[-B+df/2:df:B];

EP=zeros(size(f));

EPg=zeros(size(f));

forii=1:10;

forj=1:50;

b=sign(randn(1,M));

fori=1:L,s(i+[0:M-1]*L)=b;end

P=t2f(s);

P=P.*conj(P)/T;

EP=(EP*(j-1)+P)/j;

end

Ps=10*log10(EP+eps);

Bb=Tb/0.3;

alpha=sqrt(logm(2)/2/Bb^2);

H=exp(-alpha^2*f.^2);

a(1)=b(1);

fori=M:-1:2,a(i)=b(i)*b(i-1);end

fori=1:L,sa(i+[0:M-1]*L)=a;end

send=real(f2t(t2f(s).*H));

It=zeros(size(t));

fork=0:2*L:N-1;

kk=1:2:2*L;

kkk=1:L;

It(k+kk)=send(k+kkk+L);

It(k+kk+1)=send(k+kkk+L);

end

fork=N:-1:L+1,It(k)=It(k-L);end

Qt=zeros(size(t));

fork=0:2*L:N-1;

kk=1:2:2*L;

kkk=1:L;

Qt(k+kk)=send(k+kkk);

Qt(k+kk+1)=send(k+kkk);

end

subplot(2,1,1);

stem(b);

title(原始);

%x=input(xxx);

%Itt

Itt=It.*cos(pi*t/2/Tb);

%Qtt

Qtt=Qt.*sin(pi*t/2/Tb);

%GMSK时域波形

gmsk=Itt.*cos(2*pi*fc*t)-Qtt.*sin(2*pi*fc*t);

%GMSK功率谱

PP=t2f(gmsk);

Pa=PP.*conj(PP)/T;

EPg=(EPg*(ii-1)+Pa)/ii;

end

Pgmsk=10*log10(EPg+eps);

%接收端

r=gmsk;

%接收端的低通滤波器，带宽为Rb

LPF=zeros(size(f));

ai=(B-Rb)/2/B*size(f);

aj=(B+Rb)/2/B*size(f);

fork=(ai(1,2):aj(1,2)),LPF(k)=1;end

%接收端上支路LPF的输出，与Itt相似（图九）

RI=r.*c