请给出matlab仿真代码

以下是一个简单的MATLAB仿真代码，用于模拟4发1收STBC发射分集多天线通信链路，并比较STBC BPSK信号和单发单收BPSK信号的误比特率性能：

clear all;
clc;

% 设置参数
M = 2;  % 调制阶数，这里使用BPSK
nT = 4; % 发射天线数
nR = 1; % 接收天线数
N = 100000; % 发送符号数
EbN0dB = 0:2:16; % 信噪比范围
EbN0 = 10.^(EbN0dB/10);

% 生成发送符号
s = randi([0 M-1], nT, N);
x = pskmod(s, M);

% 计算STBC矩阵
C = [1 0 0 0; 0 1 0 0; 0 0 1 0; 0 0 0 -1]; % STBC矩阵
X = reshape(x, nT/2, 2*N); % 将发送符号按列合并成2N个2x1的矩阵
Y = C * X; % 通过STBC矩阵计算4个发射天线的输出

% 计算噪声
for i = 1:length(EbN0)
    % 生成瑞利衰落信道
    h = (randn(nR, nT) + 1i*randn(nR, nT)) / sqrt(2);
    H = repmat(h, 1, 2*N); % 将信道复制2N次，用于发射2N个信号
    % 计算噪声
    n = sqrt(0.5*1/EbN0(i))*(randn(nR, 2*N) + 1i*randn(nR, 2*N));
    % 接收到的信号
    r = sum(H.*Y, 1) + n;
    % 最小误码率译码
    r = reshape(r, 2, N);
    r1 = r(1,:);
    r2 = r(2,:);
    s_hat = zeros(1,N);
    s_hat(r1>0) = s_hat(r1>0) + 1;
    s_hat(r2>0) = s_hat(r2>0) + 2;
    % 计算误比特率
    BER(i) = sum(s_hat~=s(:))/N;
end

% 单发单收BPSK信号性能比较
for i = 1:length(EbN0)
    % 生成瑞利衰落信道
    h = (randn(nR, 1) + 1i*randn(nR, 1)) / sqrt(2);
    % 计算噪声
    n = sqrt(0.5*1/EbN0(i))*(randn(nR, N) + 1i*randn(nR, N));
    % 发送和接收信号
    x1 = pskmod(randi([0 M-1], 1, N), M);
    y1 = h*x1 + n;
    % 最小误码率译码
    s_hat1 = pskdemod(y1, M);
    % 计算误比特率
    BER1(i) = sum(s_hat1~=randi([0 M-1], 1, N))/N;
end

% 画图
semilogy(EbN0dB, BER, 'bo-', EbN0dB, BER1, 'r^-');
grid on;
legend('4发1收STBC BPSK', '单发单收BPSK');
xlabel('Eb/N0 (dB)');
ylabel('BER');
title('误比特率性能比较');

需要注意的是，以上代码仅供参考，实际应用中需要根据具体情况进行调整和优化。