%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）