2022年直接序列扩频通信系统matlab仿真程序 .pdf

《2022年直接序列扩频通信系统matlab仿真程序 .pdf》由会员分享，可在线阅读，更多相关《2022年直接序列扩频通信系统matlab仿真程序 .pdf（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、直接序列扩频通信系统仿真程序 code_length=20; %信息码元个数 N=1:code_length; rand(seed,0); x=sign(rand(1,code_length)-0.5); %信息码 for i=1:20 s(1+(i-1)*800):i*800)=x(i); %每个信息码元内含fs/f=800 个采样点end 生成的信息码的波形图如图1 所示。图 1 信源信息码 %产生伪随机码，调用的mgen 函数见附录 length=100*20; %伪码频率5MHz, 每个信息码内含5MHz/50kHz=100个伪码 x_code=sign(mgen(19,8,lengt

2、h)-0.5); %把 0,1 序列码变换为-1,1 调制码 for i=1:2000 w_code(1+(i-1)*8):i*8)=x_code(i); %每个伪码码元内含8 个采样点end 生成的 PN 码波形如图2 所示。 %扩频 k_code=s.*w_code; %k_code 为扩频码扩频码如图3 所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 5 页 - - - - - - - - - 图 2 PN 码图 3 扩频码 %调制 fs=20e6; f0=3

3、0e6; for i=1:2000 AI=2; dt=fs/f0; n=0:dt/7:dt; %一个载波周期内采样八个点cI=AI*cos(2*pi*f0*n/fs); signal(1+(i-1)*8):i*8)=k_code(1+(i-1)*8):i*8).*cI; end 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 5 页 - - - - - - - - - PSK 调制后的波形如图4 所示。图 4 PSK 调制后的波形 %解调 AI=1; dt=fs/f0;

4、n=0:dt/7:dt; %一个载波周期内采样八个点 cI=AI*cos(2*pi*f0*n/fs); for i=1:2000 signal_h(1+(i-1)*8):i*8)=signal(1+(i-1)*8):i*8).*cI; end 解调后的波形如图5 所示。图 5 解调后的波形名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 5 页 - - - - - - - - - %解扩jk_code=signal_h.*w_code; %低通滤波 wn=5/1000000

5、0; %截止频率wn=fn/(fs/2), 这里的 fn 为信息码 (扩频码 )的带宽 5M b=fir1(16,wn); H=freqz(b,1,16000); signal_d=filter(b,1, jk_code); 解扩并滤波后的波形如图6 所示。图 6 解扩并滤波后的波形从图形整体看， 解扩出来的信息码与信源信息码基本相同。输入与输出对比如图7 所示。图 7 输入与输出对比名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 5 页 - - - - - - - -

6、- 附录：%mgen.m functionout=mgen(g,state,N)gen=dec2bin(g)-48;M=length(gen);curState=dec2bin(state,M-1)-48;for k=1:N out(k)=curState(M-1); a=rem(sum(gen(2:end).*curState),2); curState=a curState(1:M-2);end名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 5 页 - - - - - - - - -