• 大小: 3KB
    文件类型: .m
    金币: 1
    下载: 0 次
    发布日期: 2021-06-08
  • 语言: Matlab
  • 标签: MSK:matlab  

资源简介

MSK 的调制与平方环解调方法,有详细的注释ps=1*10^6; %码速率为1MHz Fs=16*10^6; %采样速率为16MHz fc=3*10^6; %载波频率为3MHz N=100; %数据码元个数

资源截图

代码片段和文件信息

%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 = (randi(2N1)-1)‘; % 产生随机数据做为数据码元
%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);
figure(4);
subplot(211); plot(udi);
subplot(212); plot(udq);

%产生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)

评论

共有 条评论

相关资源