matlab编写2psk、qpsk、8psk的蒙特卡罗仿真

以下是MATLAB编写2PSK、QPSK、8PSK的蒙特卡罗仿真的示例代码：

2PSK：

N = 10^5; %仿真次数 
SNR_dB = 0:2:20; %信噪比范围 
SNR = 10.^(SNR_dB/10); %信噪比 
BER = zeros(1,length(SNR)); %误码率 
for ii=1:length(SNR) 
  for jj=1:N 
    b = randi([0 1],1,1); %随机生成二进制数据 
    s = 2*b-1; %2PSK调制 
    n = sqrt(1/(2*SNR(ii)))*randn(1,1); %高斯噪声 
    r = s + n; %接收信号 
    b_hat = (r>=0); %判决 
    BER(ii) = BER(ii) + (b~=b_hat); %计算误码率 
  end 
  BER(ii) = BER(ii)/N; %求平均 
end 
semilogy(SNR_dB,BER,'-o'); %绘制误码率曲线

QPSK：

N = 10^5; %仿真次数 
SNR_dB = 0:2:20; %信噪比范围 
SNR = 10.^(SNR_dB/10); %信噪比 
BER = zeros(1,length(SNR)); %误码率 
for ii=1:length(SNR) 
  for jj=1:N 
    b = randi([0 1],1,2); %随机生成二进制数据 
    s = 2*b-1; %2PSK调制 
    n = sqrt(1/(2*SNR(ii)))*randn(1,2); %高斯噪声 
    r = s + n; %接收信号 
    b_hat = (r>=0); %判决 
    BER(ii) = BER(ii) + sum(b~=b_hat); %计算误码率 
  end 
  BER(ii) = BER(ii)/(2*N); %求平均 
end 
semilogy(SNR_dB,BER,'-o'); %绘制误码率曲线

8PSK：

N = 10^5; %仿真次数 
SNR_dB = 0:2:20; %信噪比范围 
SNR = 10.^(SNR_dB/10); %信噪比 
BER = zeros(1,length(SNR)); %误码率 
for ii=1:length(SNR) 
  for jj=1:N 
    b = randi([0 1],1,3); %随机生成二进制数据 
    s = exp(1j*(2*pi/8)*(2*b(1)+b(2)+b(3))); %8PSK调制 
    n = sqrt(1/(2*SNR(ii)))*(randn(1,1)+1j*randn(1,1)); %高斯噪声 
    r = s + n; %接收信号 
    [~,index] = min(abs(r-exp(1j*(2*pi/8)*[0:7]))); %判决 
    b_hat = de2bi(index-1,3); %将索引转换为二进制数据 
    BER(ii) = BER(ii) + sum(b~=b_hat); %计算误码率 
  end 
  BER(ii) = BER(ii)/(3*N); %求平均 
end 
semilogy(SNR_dB,BER,'-o'); %绘制误码率曲线

以上代码中，N是仿真次数，SNR_dB是信噪比范围，SNR是信噪比，BER是误码率。对于2PSK和QPSK，随机生成的二进制数据b是一个长度为1或2的向量，对应于一个符号；对于8PSK，随机生成的二进制数据b是一个长度为3的向量，对应于一个符号。接收信号r是接收到的符号加上高斯噪声。判决是根据接收信号r是否大于0来判定，对于8PSK，通过最小距离判决。误码率是根据实际的二进制数据和判决得到的二进制数据计算得到的。注意，对于QPSK和8PSK，生成的符号是复数。