用MATLAB编程实现2DPSK信号的调制解调 实现模块：调制、信道（高斯白噪声）、解调（要求用相干检测＋码反变换） 选做：比较码反变换器前后的误码率性能曲线（要求用Nｅ／N总方式计算，信噪比范围至少－５～２０ｄＢ） 码元速率：学号的数字相加，单位KBaud，例如200021029，码元速率为2+0+0+0+2+1+0+2+2+9=18KBaud 载波频率：90。 传输信息内容：课号+个人学号+姓名 如：2222036+2000210229+MIXIANMIN

好的，这是一个关于2DPSK信号调制解调的MATLAB编程实现的问题，其中包括调制、信道（高斯白噪声）和解调（要求用相干检测＋码反变换），码元速率为18KBaud，载波频率为90。传输信息内容为课号+个人学号+姓名，例如2222036+2000210229+MIXIANMIN。选做任务是比较码反变换器前后的误码率性能曲线，要求使用Nｅ／N总方式计算，信噪比范围至少为-5～20dB。下面是一个可能的MATLAB实现，供参考：

1. 调制

fc = 90; % 载波频率
R = 18e3; % 码元速率
T = 1/R; % 码元时间
fs = 10*R; % 采样率
t = 0:1/fs:T-1/fs; % 采样时间序列
A = 1; % 信号幅度
M = 2; % 调制阶数
N = length(t); % 信号长度

% 生成调制信号
data = '2222036+2000210229+MIXIANMIN'; % 传输信息内容
bitStream = de2bi(double(data), 8, 'left-msb'); % 将字符转换为比特流
bitStream = reshape(bitStream', [], 1); % 将比特流转换为列向量
symbolStream = bi2de(reshape(bitStream, M, [])', 'left-msb'); % 将比特流转换为调制符号
phase = 2*pi/M*(symbolStream+1/2); % 计算相位
s = A*cos(2*pi*fc*t + phase); % 生成调制信号

2. 信道（高斯白噪声）

EbNo_dB = -5:1:20; % 信噪比范围
Ne_Ntotal = 10.^(EbNo_dB/10); % Ne/Ntotal
BER = zeros(size(EbNo_dB)); % 误码率

for i = 1:length(EbNo_dB)
    Ne = Ne_Ntotal(i)*R/2; % 噪声功率
    n = sqrt(Ne)*randn(1, N); % 生成高斯白噪声
    r = s + n; % 添加噪声
    BER_temp = zeros(1, M); % 记录每个符号的误码率
    for k = 0:M-1
        h = A*cos(2*pi*fc*t + 2*pi/M*k); % 信道频率响应
        y = r.*h; % 相乘
        z = conv(y, h, 'same'); % 相干检测
        z = z(N/2+1:end); % 去掉前导零
        z = z-mean(z); % 去除直流分量
        symbolHat = 2*(real(z)>0)-1; % 码元反变换
        BER_temp(k+1) = sum(symbolHat~=symbolStream(k+1:M:end))/length(symbolStream(k+1:M:end)); % 计算误码率
    end
    BER(i) = mean(BER_temp); % 计算平均误码率
end

3. 绘制误码率性能曲线

semilogy(EbNo_dB, BER, 'o-'); % 绘制误码率性能曲线
xlabel('Eb/No (dB)'); ylabel('BER');
grid on;

这是一个基本的MATLAB实现框架，您可以根据实际需要进行修改和完善。