请你用MATLAB写出一个bpsk调制的相干光通信系统

以下是一个简单的基于MATLAB的BPSK调制的相干光通信系统的代码示例：

% 定义系统参数
fc = 10e6; % 光载波频率
Rs = 1e6; % 符号速率
T = 1/Rs; % 符号时间
nBits = 1000; % 要发送的二进制比特数
SNRdB = 10; % 信噪比

% 生成随机比特流
data = randi([0, 1], nBits, 1);

% BPSK调制
s = 2*data-1;

% 生成QPSK调制的正弦余弦波形
t = linspace(0, T*nBits, nBits*100);
carrierI = cos(2*pi*fc*t);
carrierQ = sin(2*pi*fc*t);

% 将调制信号与正弦余弦波形相乘
sI = s.*carrierI(1:length(s));
sQ = s.*carrierQ(1:length(s));

% 将I、Q信号加起来
sTx = sI + 1j*sQ;

% 加入高斯白噪声
Es = mean(abs(sTx).^2);
SNR = 10^(SNRdB/10);
N0 = Es/SNR;
n = sqrt(N0/2)*randn(size(sTx)) + 1j*sqrt(N0/2)*randn(size(sTx));
rRx = sTx + n;

% 解调
rI = real(rRx).*carrierI(1:length(s));
rQ = imag(rRx).*carrierQ(1:length(s));
rI = lowpass(rI, fc/10, Rs);
rQ = lowpass(rQ, fc/10, Rs);
r = rI + rQ;
r = r(1:100:length(r));

% 比特误码率
numErrs = sum(abs((r+1)/2 - data));
BER = numErrs/nBits;

% 绘制调制信号、接收信号和解调信号
figure;
subplot(3,1,1);
plot(t, real(sTx));
title('BPSK调制信号');
xlabel('时间 (s)');
ylabel('幅度');
subplot(3,1,2);
plot(t, real(rRx));
title(sprintf('接收信号 (SNR = %d dB)', SNRdB));
xlabel('时间 (s)');
ylabel('幅度');
subplot(3,1,3);
plot(t, real(r));
title(sprintf('解调信号 (BER = %g)', BER));
xlabel('时间 (s)');
ylabel('幅度');

这个代码示例中生成了一个随机的二进制比特流，并使用BPSK调制将其调制到光载波上。然后，生成了正弦余弦波形，将调制信号分别与正弦余弦波形相乘，最后将信号加起来。接着，加入了高斯白噪声，解调接收信号，并计算了比特误码率。最后，绘制了调制信号、接收信号和解调信号。