帮我写一下7位m序列通过高斯信道的直接序列扩频matlab代码，并画出误码率

好的，以下是7位m序列通过高斯信道的直接序列扩频matlab代码：

clear all;
close all;
clc;

%% 参数设置
N = 1000000; % 发送比特数
Eb = 1; % 比特能量
SNR_dB = 0:2:20; % 信噪比范围
SNR = 10.^(SNR_dB/10); % 信噪比
m_seq = [1,1,1,0,1,0,0]; % 7位m序列
N_seq = length(m_seq); % 序列长度
L = 16; % 扩频因子
fc = 1; % 载波频率
T = 1/N_seq; % 符号周期
t = 0:T/L:T-T/L; % 扩频信号时间轴

%% 生成扩频信号
s = zeros(1,L*N_seq);
for i = 1:N_seq
    if m_seq(i) == 1
        s((i-1)*L+1:i*L) = ones(1,L);
    else
        s((i-1)*L+1:i*L) = -ones(1,L);
    end
end
c = exp(1j*2*pi*fc*t); % 载波信号
x = s.*c; % 发送信号

%% 高斯信道模拟
for k = 1:length(SNR)
    y = awgn(x,SNR_dB(k)+10*log10(log2(L)),'measured'); % 添加高斯白噪声
    r = y.*conj(c); % 接收信号
    r_seq = zeros(1,N_seq);
    for i = 1:N_seq
        r_seq(i) = sum(r((i-1)*L+1:i*L))/L;
    end
    m_hat = (r_seq > 0); % 判决
    err(k) = sum(m_hat ~= m_seq); % 统计误码数
end

%% 误码率计算与绘图
BER = err/N;
semilogy(SNR_dB,BER); % 绘制误码率曲线
xlabel('SNR(dB)');
ylabel('BER');
title('7位m序列通过高斯信道的直接序列扩频');
grid on;

运行上述代码，即可得到误码率曲线图像。