用M文件编写 4-PSK、16-PSK仿真，并与8-PSK比较；或者，完成2FSK/2ASK仿真。

好的，我可以为您提供一些相关的MATLAB代码，用于实现4-PSK、16-PSK、8-PSK、2FSK和2ASK的仿真。请注意，这些代码只是基础实现，您可以根据自己的需求进行修改和优化。

1. 4-PSK仿真

clc;
close all;
clear all;

M = 4; % 星座大小
k = log2(M); % 每个符号的比特数

% 创建4-PSK星座映射表
constellation = [1+1i, -1+1i, -1-1i, 1-1i];
constellation = constellation / sqrt(2);

% 随机生成M个符号
txBits = randi([0 1], 1, 10000*k);

% 将比特序列转换为符号序列
txSym = bi2de(reshape(txBits, k, []).', 'left-msb');
txSym = constellation(txSym+1);

% 加入高斯白噪声
snr = 10; % 信噪比
rxSym = awgn(txSym, snr, 'measured');

% 解调接收信号
rxBits = de2bi(qamdemod(rxSym, M, 'UnitAveragePower', true)-1, k, 'left-msb');
rxBits = rxBits(:).';

% 计算误比特率
ber = sum(txBits ~= rxBits) / length(txBits);
fprintf('4-PSK: SNR = %d dB, BER = %e\n', snr, ber);

2. 16-PSK仿真

clc;
close all;
clear all;

M = 16; % 星座大小
k = log2(M); % 每个符号的比特数

% 创建16-PSK星座映射表
constellation = exp(1i*(0:M-1)*2*pi/M);
constellation = constellation / sqrt(sum(abs(constellation).^2)/M);

% 随机生成M个符号
txBits = randi([0 1], 1, 10000*k);

% 将比特序列转换为符号序列
txSym = bi2de(reshape(txBits, k, []).', 'left-msb');
txSym = constellation(txSym+1);

% 加入高斯白噪声
snr = 10; % 信噪比
rxSym = awgn(txSym, snr, 'measured');

% 解调接收信号
rxBits = de2bi(qamdemod(rxSym, M, 'UnitAveragePower', true)-1, k, 'left-msb');
rxBits = rxBits(:).';

% 计算误比特率
ber = sum(txBits ~= rxBits) / length(txBits);
fprintf('16-PSK: SNR = %d dB, BER = %e\n', snr, ber);

3. 8-PSK仿真

clc;
close all;
clear all;

M = 8; % 星座大小
k = log2(M); % 每个符号的比特数

% 创建8-PSK星座映射表
constellation = exp(1i*(0:M-1)*2*pi/M + pi/8);
constellation = constellation / sqrt(sum(abs(constellation).^2)/M);

% 随机生成M个符号
txBits = randi([0 1], 1, 10000*k);

% 将比特序列转换为符号序列
txSym = bi2de(reshape(txBits, k, []).', 'left-msb');
txSym = constellation(txSym+1);

% 加入高斯白噪声
snr = 10; % 信噪比
rxSym = awgn(txSym, snr, 'measured');

% 解调接收信号
rxBits = de2bi(qamdemod(rxSym, M, 'UnitAveragePower', true)-1, k, 'left-msb');
rxBits = rxBits(:).';

% 计算误比特率
ber = sum(txBits ~= rxBits) / length(txBits);
fprintf('8-PSK: SNR = %d dB, BER = %e\n', snr, ber);

4. 2FSK仿真

clc;
close all;
clear all;

fs = 100; % 采样率
fc1 = 10; % 第一个载波频率
fc2 = 20; % 第二个载波频率
T = 1/fs;
t = 0:T:1-T;

% 随机生成0和1的比特序列
txBits = randi([0 1], 1, 10000);

% 将比特序列映射为调制信号
txSym = 2*txBits - 1;
txSig = txSym .* cos(2*pi*fc1*t + 2*pi*(fc2-fc1)*cumsum(t.*txBits));

% 加入高斯白噪声
snr = 10; % 信噪比
rxSig = awgn(txSig, snr, 'measured');

% 解调接收信号
rxBits = (sign(rxSig .* cos(2*pi*fc1*t + 2*pi*(fc2-fc1)*cumsum(t.*txBits))) + 1) / 2;

% 计算误比特率
ber = sum(txBits ~= rxBits) / length(txBits);
fprintf('2FSK: SNR = %d dB, BER = %e\n', snr, ber);

5. 2ASK仿真

clc;
close all;
clear all;

fs = 100; % 采样率
fc = 10; % 载波频率
T = 1/fs;
t = 0:T:1-T;

% 随机生成0和1的比特序列
txBits = randi([0 1], 1, 10000);

% 将比特序列映射为调制信号
txSym = 2*txBits - 1;
txSig = txSym .* cos(2*pi*fc*t);

% 加入高斯白噪声
snr = 10; % 信噪比
rxSig = awgn(txSig, snr, 'measured');

% 解调接收信号
rxBits = (sign(rxSig) + 1) / 2;

% 计算误比特率
ber = sum(txBits ~= rxBits) / length(txBits);
fprintf('2ASK: SNR = %d dB, BER = %e\n', snr, ber);