【信号处理】构建直接序列扩频系统模型（Matlab代码实现）

news/2024/4/29 15:37:50/文章来源:https://blog.csdn.net/weixin_46039719/article/details/128086715

👨‍🎓个人主页：研学社的博客

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🌈3 Matlab代码实现

🎉4 参考文献

💥1 概述

编程产生周期长度为31的m序列；
产生随机发送信息，并分别进行BPSK调制；
采用m序列作为扩频码，产生扩频信号发送信号；
对扩频发送信号加入白噪声，得到接收扩频信号；
对接收信号进行解扩；
对接收信号进行BPSK解调；
通过蒙特卡洛方法，仿真系统误码率，画出误码率曲线。
产生周期长度为31的m序列，记为w；（实验6的内容，m序列周期T=2n-1）（需转换成1、-1）
产生发送信息0、1比特；（三种方法：循环判断、rand+round、randint/randi）（数据长度L自定，尽量大一些，比如1000）
对信息比特进行BPSK调制，记为d；（即码型变换：1→1, 0→-1）（注：出于简化，本实验省略载波调制部分）
对发送信息进行扩频，即d*w；（实现提示1：采用循环的方式，对每位数据进行扩频，运行速度较慢）（实现提示2：可将d的每位数据重复31次，w扩大L次(L为发送数据长度)，再点乘。可用函数：repmat函数、reshape函数）
信号经过信道，叠加白噪声，得到接收端接收信号；（噪声功率由信噪比决定，仿真结果横坐标为信噪比，可取-30:5，单位dB）（可使用awgn函数，也可自行实现）
对接收端信号进行恢复，包括解扩、解调、判决。方法如下：a. 首先与扩频码相乘进行解扩，；b. 对每位数据扩大的31次进行求和，以第一个数据b1为例，sum(b1)；c. 判决，sum(b1)>0，判为1，否则为-1。
误码率统计，并画图。（画图函数：semilogy）

📚2 运行结果

🌈3 Matlab代码实现

部分代码：

clc;clear all;
r=5;
N=2^r-1;%周期31
a=ones(1,r);
m=zeros(1,N);
for i=1:(2^r-1)
temp= mod((a(5)+a(2)),2);
for j=r:-1:2
a(j)=a(j-1);
end
a(1)=temp;
m(i)=a(r);
end
m=m*2-1;%双极性码
%产生随机发送信息，并进行BPSK/QPSK调制
n=1000;
Y_BIT_sum=zeros(1,36);
for jj=1:100%跑100次，得到平滑曲线
source=randi([0 1],1,n);
% s_imag=randi(1,n);
s_BPSK=2*source-1;
% s=source+j*s_imag;
% j=sqrt(-1);
% s_QPSK=s*2-(1+j);
%产生扩频信号：
% for z=1:1:100
% y_QPSK(31*(z-1)+(1:31))=m*s_QPSK(z);
% end
for z=1:1000
y_BPSK(31*(z-1)+(1:31))=m*s_BPSK(z);
end
%加入高斯白噪声
SNR=-30:5;
for i=1:length(SNR)
y(i,:)=awgn(y_BPSK,SNR(i));
end
%解扩
for z=1:1000
for i=1:length(SNR)
o_BPSK(i,z)=y(i,31*(z-1)+(1:31))/m;
end
end
Y_RE=o_BPSK;
%判决
Y_RE(o_BPSK>0)=1;
Y_RE(o_BPSK<0)=0;
Y_ERRO=zeros(length(SNR),1000);
Y_BIT=zeros(1,length(SNR));
for i=1:length(SNR)