【扩频通信】直接序列扩频通信仿真【含Matlab源码 2650期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【扩频通信】基于matlab直接序列扩频通信仿真【含Matlab源码 2650期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏信号处理（Matlab）

备注：
点击上面蓝色字体付费专栏图像处理（Matlab），扫描上面二维码，付费299.9元订阅海神之光博客付费专栏，凭支付凭证，私信博主，可免费获得5份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：5份本博客上传CSDN资源代码

⛄二、简介

⛄三、部分源代码

 本程序用于多用户CDMA通信系统仿真,同时引入了对多径衰落的模拟，采用Walsh码作为扩频码

传输过程使用BPSK调制。
%}
close all;
user_num = 8; %用户数，必须少于或等于walsh码阶数
user_code = 1e4; %每个用户的码元数

walsh_order = 8; %walsh码阶数,也是扩频增益,必须是2的幂次方

multipath_fading_num = 5; %多径衰落通道数量

carrier_freq = 1e3; %载波频率
SampleRate = 25*1e3; %载波采样率
SamplePoint = 25; %载波采样点数

frame_size = 8; %汉明码一帧的二进制位数

mix_frame_size = 8; %交织编码一帧的二进制位数

%产生用户码元
userCode = genBipolar(user_num,user_code);

%产生多径衰落时滞码
mult_fading_code = multiPathCodeGen(multipath_fading_num);

%交织编码
foo = mixed(userCode,mix_frame_size,@mixed_encode);

%对每个用户分别进行汉明码编码
foo = hamming(foo,frame_size);

%产生Walsh矩阵
walshCode = walsh(walsh_order);

%扩频
foo = spreadSpectrum(foo,walshCode);

%扩频后的码元相加
foo = overlay2(foo);

%生成载波
carrier = carrierGen(carrier_freq,SampleRate,SamplePoint);
%调制
foo = modulate(foo,carrier);

%尝试不同的信噪比(dB)
snr_max = 30;
snr_min = -10;
snr_div = 1;
snr_array = floor((snr_max-snr_min)/snr_div);
%记录每个用户每个信噪比下的用户误码率
user_errorRate = zeros(snr_array,user_num);
%记录平均误码率
mean_errorRate = zeros(1,snr_array);
parfor snr_index = 1:snr_array

%当前信噪比
snr = snr_min + snr_div * snr_index;

%加噪
bar = awgn(foo,snr,powerCnt(carrier));

%添加时滞效应
length_bar = length(bar);
bar = conv(bar,mult_fading_code);
bar = arrayCut(bar,length_bar);

%解调(相干解调)
bar = demodulate(bar,carrier);

%解扩
bar = deSpreadSpectrum(bar,walshCode,user_num);

%码元判决
bar = codeJudge(bar,walsh_order);

%汉明码解码
bar = unhamming(bar,frame_size);

%逆交织编码
bar = mixed(bar,mix_frame_size,@mixed_decode);

%计算误码率
error_rate = 1 - compare(bar,userCode);

%保存误码率
user_errorRate(snr_index,:) = error_rate;
mean_errorRate(snr_index) = mean(error_rate);

fprintf('当前信噪比%f，平均误码率%f\n',snr,mean(error_rate));

end

%绘制误码率曲线
foo_x = linspace(snr_min,snr_max,snr_array);
%生成标签
legend_plot = cell(1,user_num);
figure;
for i = 1:user_num
semilogy(foo_x,user_errorRate(:,i));
legend_plot{i} = [‘用户’,num2str(i)];
hold on;
end
title(‘用户误码率曲线’);
legend(legend_plot);
xlabel(‘信噪比(dB)’);
ylabel(‘误码率’);

%绘制平均误码率曲线
figure;
%计算平均误码率
semilogy(foo_x,mean_errorRate);
title(‘平均误码率曲线’);
xlabel(‘信噪比(dB)’);
ylabel(‘误码率’);

⛄四、运行结果

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除