msk的调制解调MATLAB源代码.docx
《msk的调制解调MATLAB源代码.docx》由会员分享,可在线阅读,更多相关《msk的调制解调MATLAB源代码.docx(17页珍藏版)》请在冰豆网上搜索。
![msk的调制解调MATLAB源代码.docx](https://file1.bdocx.com/fileroot1/2023-1/23/35f08721-d468-416c-98d2-1d178f419d6b/35f08721-d468-416c-98d2-1d178f419d6b1.gif)
msk的调制解调MATLAB源代码
msk的调制解调MATLA源代码
functionout=delay(data,n,sample_number)
%data:
延迟的数据
%n:
延迟码元个数
%sample_number码元采样个数
out=zeros(1,length(data));
out(n*sample_numbe叶1:
length(data))=
data(1:
length(data)-n*sample_number);
function[data_diff]=difference(data)
%差分编码
%**********************************************************
****************
%data输入信号
%data_diff差分编码后信号
%**********************************************************
****************
data_diff=zeros(1,length(data));
data_diff
(1)=1*data
(1);%1为差分编码的初始参考值
fori=2:
length(data)
data_diff(i)=data_diff(i-1)*data(i);
end
%**********************************************************
****************
function[signal_out,l_out,Q_out]=mod_msk(data,data_len,sample_number,Rb)
%MS基带调制
%**********************************************************
****************
%data调制信号
%data_len码元个数
%sample_number每个码元采样点数
%Rb
码元速率
%signal_out
基带调制输出
%I_out
I路输出
%Q_out
Q路输出
%**********************************************************
****************
%data_len二10;%码元个数
%sample_number=8;%采样点数
%Rb=16000;%码元速率
%datal=randint(1,data_len);
%data=2*data1-1;%传输的序列
Tb=1/Rb;%码元时间
fs=Rb*sample_number;%采样速率
%
%差分编码
[data_diff]=difference(data);
%**********************************************************
****************
%并串转换,延时
1
(1)=1;%faiO=0,cos(fai0)=1
fori=1:
2:
data_len
Q(i)=data_diff(i);
Q(i+1)=data_diff(i);
end
fori=2:
2:
data_len
I(i+1)=data_diff(i);
I(i)=data_diff(i);
end
fori=1:
sample_number
I1(i:
sample_number:
data_len*sample_number)=
I(1:
data」en);
Q1(i:
sample_number:
data_len*sample_number)=
Q(1:
data_len);
end
%**********************************************************
****************
%乘加权函数
t=1/fs:
1/fs:
data_len*Tb;
l_out=I1.*cos(pi*t/2/Tb);
Q_out=Q1.*sin(pi*t/2/Tb);
%**********************************************************
****************
%
%调制信号产生
signal_out=I_out+j*Q_out;
%**********************************************************
****************
%%
%%|图
%subplot(221)
时间
%plot(data,'.-');title('MSK传输的数据');xlabel('
*3nH«
)y_abe_(®w・)
%subp_of(222)
%p-oadafaldiff:
-—.辽m®可3選箔・)x_abe-(兰凹
)y_abe_(®w・)
%subp_of(223)
%p_2-(-1・<)s_e(MOtf-^・)x_abe_(兰凹)y_abe三
®w)
%subp_of(224)
%p_of(Q1:
-—・)s_e(MOtfQ^)x_abe_(兰凹)y_abe三
®w)
%
%figu「e
(2)
%subpof(22」)
%p_2.(cos(pi*t/2Hb)「<辽f_e(8®^cos(mt/(2Tb)y)><一abe三鱼凹)y_abe_(®w・)
%subp_of(222)
%p_2.(sin(piG2Hb)・<)s_e(8®^
sin(mt/(2Tb)y)><一abe三鱼凹)y_abe_(s・)
%subp_of(223)
%p_of(_louu<辽f_e(岂商可-^・)x_abe-(昇「凹
)y_abe_(®w)
%subp_of(224)
6sfsyB
时间
%plot(Q_out,'.-');title('加权后Q路');xlabel('
');ylabel('幅度');
%%********************************************************
******************
function[signal_out,l_out,Q_out,phase]=mod_msk2(data,data_len,sample_number,Rb)
%MS基带调制
%**********************************************************
****************
%data调制信号
%data_len码元个数
%sample_number每个码元采样点数
%Rb码元速率
%signal_out基带调制输出
%I_outI路输出
%Q_outQ路输出
%**********************************************************
****************
%data_len二10;%码元个数
%sample_number=8;%采样点数
%Rb=16000;%码元速率
%datal=randint(1,data_len);
%data=2*data1-1;%传输的序列
Tb=1/Rb;%码元时间
fs=Rb*sample_number;%采样速率
%
%采样
fori=1:
sample_number
data_sample(i:
sample_number:
data_len*sample_number)=data;
end
%**********************************************************
****************
%
%十算相位
phase=zeros(1,data_len*sample_number);
phase
(1)=data_sample
(1)*pi/2/sample_number;
fori=2:
data_len*sample_number
phase(i)=phase(i-1)+data_sample(i-1)*
pi/2/sample_number;
end
%**********************************************************
****************
%
%I/Q
I_out=cos(phase);
Q_out=sin(phase);
%**********************************************************
****************
%
%调制信号产生signal_out=l_out+j*Q_out;
%**********************************************************
****************
%MS调制,差分解调方法一
clearall
closeall
%
%参数设置
data_len二10000;%码元个数
sample_number=8;%采样个数
Rb=24000;%码元速率
fc=96000;%载波频率
%**********************************************************
****************
%
%随机产生传输信号
data二rand_binary(data_len);
%**********************************************************
****************
%
%MS基带调制
[signal_out,l_out,Q_out]=
mod_msk(data,data_len,sample_number,Rb);
%**********************************************************
****************
%
帅频搬移
multi=fc/Rb;
I_temp=interp(I_out,multi);
Q_temp=interp(Q_out,multi);
Fs=fc*sample_number;
t=1/Fs:
1/Fs:
length(l_temp)*1/Fs;
signal」=l_temp.*cos(2*pi*fc*t);
signal_q二Q_temp.*sin(2*pi*fc*t);
signal_mod=l_temp.*cos(2*pi*fc*t)-Q_temp.*sin(2*pi*fc*t);
%**********************************************************
****************
%
%加噪声
forSNR=0:
8
signal_mod1=awgn(signal_mod,SNR);
%
%去载波
N=300;%滤波器的阶数为(N+1)
F二[0,fc-1000,fc+1000,Fs/2]*2/Fs;
A=[1,1,0,0];
lpf=firls(N,F,A);
[amp_lpf,w]=freqz(lpf);
l_dem二signal_mod1.*cos(2*pi*fc*t)*2;
I_dem=conv(l_dem,lpf);
l_dem=l_dem(N/2+1:
N/2+length(l_temp));
Q_dem=signal_mod1.*sin(2*pi*fc*t)*2;
Q_dem=conv(Q_dem,lpf);
Q_dem=-Q_dem(N/2+1:
N/2+length(l_temp));
%抽取
l_dem_out=zeros(1,length(l_dem)/multi);
Q_dem_out=zeros(1,length(Q_dem)/multi);
fori=1:
length(l_dem_out)
l_dem_out(i)=l_dem(multi*(i-1)+1);
Q_dem_out(i)=Q_dem(multi*(i-1)+1);
end;
%*******************************************************
*******************
%
%差分解调
demod_data=zeros(1,data_len);
demod_data
(1)=Q_dem_out(sample_number);
fori=2:
data_len
demod_data(i)=
Q_dem_out(i*sample_number)*l_dem_out((i-1)*sample_number)-l_dem_out(i*sample_number)*Q_dem_out((i-1)*sample_number);
end
%*******************************************************
*******************
%
%判决
demod_data=demod_data>0;
demod_data=2*demod_data-1;
%*******************************************************
*******************
%
%计算误码率
[num,ber(SNR+1)]=symerr(demod_data,data);
%*******************************************************
*******************
end
%**********************************************************
****************
%
頰码率曲线
semilogy([0:
8],ber,'r*-');
%**********************************************************
****************
%
頰码率理论值
snr=0:
0.1:
8;
fori=1:
length(snr)
snr1(1,i)=10A(snr(1,i)/10);
ps(1,i)=1/2*erfc(sqrt(snr1(1,i)));
pe(1,i)=2*ps(1,i);
end
holdon
semilogy([0:
.1:
8],pe);
%**********************************************************
****************
%MS调制,差分解调方法二
clearall
closeall
%
%参数设置
data_len二10000;
%码元个数
sample_number=8;
%采样个数
Rb=24000;
%码元速率
fc=96000;
%载波频率
%**********************************************************
****************
%
%随机产生传输信号
data二rand_binary(data_len);
%**********************************************************
****************
%
%MS基带调制
[signal_out,l_out,Q_out]=
mod_msk(data,data_len,sample_number,Rb);
%**********************************************************
****************
%
帅频搬移
multi=fc/Rb;
I_temp=interp(I_out,multi);
Q_temp=interp(Q_out,multi);
Fs二fc*sample_number;
t=1/Fs:
1/Fs:
length(l_temp)*1/Fs;
signal」=l_temp.*cos(2*pi*fc*t);
signal_q二Q_temp.*sin(2*pi*fc*t);
signal_mod=l_temp.*cos(2*pi*fc*t)-Q_temp.*sin(2*pi*fc*t);
%**********************************************************
****************
%
%加噪声
forSNR=0:
8
signal_mod1=awgn(signal_mod,SNR);
%
%去载波
N=300;%滤波器的阶数为(N+1)
F二[0,fc-1000,fc+1000,Fs/2]*2/Fs;
A=[1,1,0,0];
lpf=firls(N,F,A);
[amp_lpf,w]=freqz(lpf);
l_dem二signal_mod1.*cos(2*pi*fc*t)*2;
I_dem=conv(l_dem,lpf);
l_dem=l_dem(N/2+1:
N/2+length(l_temp));
Q_dem=signal_mod1.*sin(2*pi*fc*t)*2;
Q_dem=conv(Q_dem,lpf);
Q_dem=-Q_dem(N/2+1:
N/2+length(l_temp));
l_dem_out=zeros(1,length(l_dem)/multi);%抽取
Q_dem_out=zeros(1,length(Q_dem)/multi);
fori=1:
length(l_dem_out)
l_dem_out(i)=l_dem(multi*(i-1)+1);
Q_dem_out(i)=Q_dem(multi*(i-1)+1);
end;
%*******************************************************
*******************
%
%差分解调
demod_data=zeros(1,data_len);
demod_data
(1)=Q_dem_out(sample_number);
fori=2:
2:
data_len
demod_data(i)=
-l_dem_out(i*sample_number)*Q_dem_out((i-1)*sample_number);end
fori=3:
2:
data_len
demod_data(i)=
Q_dem_out(i*sample_number)*l_dem_out((i-1)*sample_number);end
%*******************************************************
*******************
%
%判决
demod_data=demod_data>0;
demod_data=2*demod_data-1;
%*******************************************************
*******************
%
%计算误码率
[num,ber(SNR+1)]=symerr(demod_data,data);
%*******************************************************
*******************
end
%**********************************************************
****************
%
頰码率曲线
semilogy([0:
8],ber,'r*');
%**********************************************************
****************
%
頰码率理论值
snr=0:
0.1:
8;
fori=1:
length(snr)
snr1(1,i)=10A(snr(1,i)/10);
ps(1,i)=1/2*erfc(sqrt(snr1(1,i)));
pe(1,i)=2*ps(1,i);
end
holdonsemilogy([0:
.1:
8],pe);
0/严*******************************************************
****************
°%MS^制,解调
clearall
closeall
%
%参数设置
data_len=30000;%码元个数
sample_number=8;%采样个数
Rb=24000;%码元速率
fc=96000;%载波频率
%**********************************************************
****************
%
%随机产生传输信号
data=rand_binary(data_len);
%**********************************************************
****************
%
%MS基带调制
[signal_out,l_out,Q_out]=
mod_msk(data,data_len,sample_number,Rb);
%**********************************************************