帮我编写MANTLAB代码：4-PSK、16-PSK仿真，并与8-PSK比较

好的，下面是MANTLAB代码，用于4-PSK、16-PSK仿真，并与8-PSK比较。

clc;
clear all;
close all;

%定义常量
M1=4;   %4-PSK
M2=16;  %16-PSK
M3=8;   %8-PSK
N=10^6; %仿真次数
SNRdB=0:2:20;  %信噪比范围

%生成调制符号
data1=randi([0 M1-1],1,N);   %4-PSK
data2=randi([0 M2-1],1,N);   %16-PSK
data3=randi([0 M3-1],1,N);   %8-PSK

%调制
modData1=pskmod(data1,M1);   %4-PSK
modData2=pskmod(data2,M2);   %16-PSK
modData3=pskmod(data3,M3);   %8-PSK

%循环计算信噪比
for ii=1:length(SNRdB)
    %加噪声
    noise=(1/sqrt(2))*(randn(1,N/M1)+1i*randn(1,N/M1));  %4-PSK
    r1=modData1+10^(-SNRdB(ii)/20)*noise;
    noise=(1/sqrt(2))*(randn(1,N/M2)+1i*randn(1,N/M2));  %16-PSK
    r2=modData2+10^(-SNRdB(ii)/20)*noise;
    noise=(1/sqrt(2))*(randn(1,N/M3)+1i*randn(1,N/M3));  %8-PSK
    r3=modData3+10^(-SNRdB(ii)/20)*noise;
    
    %解调
    demodData1=pskdemod(r1,M1);   %4-PSK
    demodData2=pskdemod(r2,M2);   %16-PSK
    demodData3=pskdemod(r3,M3);   %8-PSK
    
    %计算误码率
    [num1,err1]=biterr(data1,demodData1); %4-PSK
    [num2,err2]=biterr(data2,demodData2); %16-PSK
    [num3,err3]=biterr(data3,demodData3); %8-PSK
    
    BER1(ii)=err1/N;  %4-PSK
    BER2(ii)=err2/N;  %16-PSK
    BER3(ii)=err3/N;  %8-PSK
end

%绘制误码率曲线
semilogy(SNRdB,BER1,'bo-','LineWidth',1.5);
hold on;
semilogy(SNRdB,BER2,'r*-','LineWidth',1.5);
semilogy(SNRdB,BER3,'g^-','LineWidth',1.5);
grid on;
axis([0 20 10^(-6) 1]);
legend('4-PSK','16-PSK','8-PSK');
xlabel('信噪比(dB)');
ylabel('误码率');
title('PSK调制三种方式比较');

这段代码可以生成三种PSK调制方式（4-PSK、16-PSK、8-PSK）的信号，并进行仿真比较，输出误码率曲线图。在代码中，我们设定了信噪比范围为0~20dB，仿真次数为10^6，可以根据需要进行修改。