帮我写一个通过高斯信道的直接序列扩频通信画出误码率的matlab代码

以下是一个简单的高斯信道的直接序列扩频通信的误码率计算的Matlab代码：

clear all;
close all;
clc;

% 系统参数
M = 2; % 2PSK调制
N = 10000; % 发送比特数
SNR_dB = 0:2:10; % 信噪比范围
SNR = 10.^(SNR_dB/10); % 信噪比
fc = 10^6; % 载波频率
fs = 20*fc; % 采样频率
T = 1/fc; % 周期
t = (0:N-1)*T; % 时间序列

% 发送端
msg = randi([0 1],1,N); % 随机生成发送比特流
c = ones(1,N); % 扩频码
s = 2*msg-1; % 2PSK调制
x = s.*c; % 发送序列
f = fc*cos(2*pi*fc*t); % 载波信号
tx_signal = x.*f; % 发送信号

% 信道
noise = randn(1,N); % 高斯噪声
rx_signal = tx_signal + sqrt(2*fs./SNR).*noise; % 加入高斯噪声

% 接收端
y = rx_signal.*f; % 解调
y = y.*c; % 移除扩频码
r = sum(y)/N; % 相干积分
demod_signal = sign(r); % 2PSK解调

% 误码率计算
BER = zeros(1,length(SNR));
for i=1:length(SNR)
    Nerr = sum(msg~=demod_signal);
    BER(i) = Nerr/N;
end

% 画图
semilogy(SNR_dB,BER,'-o');
grid on;
xlabel('信噪比 (dB)');
ylabel('误码率');
title('高斯信道的直接序列扩频通信误码率');

这个代码实现了一个基本的高斯信道的直接序列扩频通信系统，并且计算了不同信噪比下的误码率。你可以根据需要修改代码中的参数来进行实验和分析。