请根据上述这段话写出相应的代码

根据上面的描述，这里提供一个基本的MATLAB代码示例来实现AM调制、相干解调和非相干解调，以及输出信噪比的计算。假设我们已经定义了变量`t`,`m`, `carrier`, 和 `noise` (噪声)：

% 定义信号和相关参数
fs = 1e6;           % 采样率
t = 0:1/fs:3;       % 时间向量
x = sin(2*pi*3e3*t); % 信源信号
carrier_freq = 30 * 3e3;  % 载波频率

% AM调制
m = 2/3;
y_AM = carrier + m * x;

% 信噪比设定
SNR_1 = 1;         % 输入信噪比1
SNR_20 = 20;        % 输入信噪比20
noise = sqrt(SNR_1/2) * randn(size(y_AM));  % 白噪声

% 相关解调
y_coherent = y_AM;  %相干解调直接使用调制信号
y_coherent_noise = y_coherent + noise;  %相干解调后加噪声

% 非相干解调
y_noncoherent = y_coherent - mean(carrier);

% 输出信噪比计算
SNR_coherent = 10*log10(mean(abs(y_coherent_noise).^2) ./ mean(abs(noise).^2));
SNR_noncoherent = 10*log10(mean(abs(y_noncoherent).^2) ./ mean(abs(noise).^2));

% 波形绘制
figure;
subplot(2,1,1);
plot(t, x, 'b', t, y_AM, 'r'); % 调制前后波形
title('AM Modulation');
xlabel('Time (s)');
ylabel('Amplitude');

subplot(2,1,2);
hold on;
plot(t, y_coherent_noise, 'g', t, y_noncoherent, 'k');  % 解调波形
legend('Original Signal', 'Coherent Demodulation', 'Non-Coherent Demodulation');
xlabel('Time (s)');
ylabel('Amplitude');

% 输出信噪比关系图
figure;
semilogx([SNR_1 SNR_20], [SNR_coherent SNR_noncoherent], 'o-');
xlabel('Input SNR (dB)');
ylabel('Output SNR (dB)');
grid on;
title('Output SNR vs Input SNR');

这个代码演示了如何在MATLAB中完成所需的操作，但请注意实际运行时可能需要根据你的特定需求进行调整。如果你有其他疑问或需要更详细的帮助，请随时提问。