msk调制的MATLAB仿真程序

%E6_7_MskModem.m程序清单
 
ps=1*10^6;  %码速率为1MHz
Fs=16*10^6; %采样速率为16MHz
fc=3*10^6;  %载波频率为3MHz

N=100;        %数据码元个数
Len=N*Fs/ps;  %仿真数据的长度
x = randint（N12）‘; % 产生随机数据做为数据码元
%x=ones（1N）;
dx=ones（1N）;
for i=1:N
    if x（i）==0
        x（i）=-1;
    end
end
%求原码的相对码dx
for i=2:N
    if x（i）==1
        dx（i）=-dx（i-1）;
    else
        dx（i）=dx（i-1）;
    end
end

%将相对码按奇偶序号分成两路数据，形成Ik\Qk
di=ones（1N）;
dq=ones（1N）;
%取dx的偶数位，并列两位为di
for i=2:2:N
    di（i:i+1）=dx（i）;
end
%取dx的奇数位，并列两位为dq
for i=1:2:N-1
    dq（i:i+1）=dx（i）;
end

%对原始BIT数据进行Fs/ps倍重采样
udi=ones（1N*Fs/ps）;
udq=ones（1N*Fs/ps）;
for i=1:N
    udi（Fs/ps*（i-1）+1:Fs/ps*i）=di（i）;
    udq（Fs/ps*（i-1）+1:Fs/ps*i）=dq（i）;
end
% udi=rectpulse（diFs/ps）;
% udq=rectpulse（dqFs/ps）;

%产生MSK信号所需的载波信号
t=0:1/Fs:（Len-1）/Fs;
cf0c=cos（2*pi*fc.*t）;
sf0c=sin（2*pi*fc.*t）;
cfps=cos（pi*ps/2.*t）;
sfps=sin（pi*ps/2.*t）;

%正交调制法产生msk信号
msk=udi.*cfps.*cf0c-udq.*sfps.*sf0c;

 
%MSK 解调
%第一种相干解调，需要同时获取fcfb的载波频率
% demod_i=msk.*cf0c.*cfps;
% demod_q=msk.*sf0c.*sfps;
%平方环相干解调，只需获取fL、fH的载波频率
fL=cos（2*pi*fc.*t-2*pi*ps/4.*t）;%cos（2*pi*ps/4.*t）
fH=cos（2*pi*fc.*t+2*pi*ps/4.*t）;

demod_i=msk.*（fH+fL）;
demod_q=msk.*（fH-fL）;
%低通滤波后，获取I、Q支路基带波形
b=fir1（300.5*ps*2/Fs）;%设计低通滤波器
f_i=filter（b1demod_i）;
f_q=filter（b1demod_q）;

%求Tb/2的时钟信号波形
% fL2=cos（4*pi