%qpsk各種仿真值
clear all
x1=0:1:7;
y1=[0.42 0.3 0.1665 0.045 0.0056 0.000987 0.000155 0.00001043];%高斯信道
x2=0:1:7;
y2=[0.46 0.335 0.1623 0.045 0.0056 0.001 0.000155 0.0000174];%單徑信道（視距外、近場開闊）
x3=0:1:7;
y3=[0.42 0.325 0.17 0.059 0.0128 0.00118 0.0001269 0.000018];%兩徑信道（視距內、近場開闊）
x4=0:1:7;
y4=[0.4388 0.3167 0.1643 0.054 0.00768 0.001586 0.000176 0.000005564];%四徑信道（遮擋、近場開闊）
x5=0:1:12;
y5=[0.457 0.4225 0.3148 0.2244 0.13 0.069 0.0288 0.01356 0.00486 0.00158 0.00038 0.00009165 0.0000458];%萊斯平坦信道K=3.8dB
x6=0:1:9;
y6=[0.435 0.37 0.258 0.1357 0.06 0.017 0.00368 0.00196 0.0001269 0.00004034];%萊斯平坦信道K=8.1dB
x7=0:5:20;
y7=[0.46 0.163 0.02 0.0011 0.00002];%瑞利平坦信道10Hz
x8=0:5:20;
y8=[0.5 0.36 0.0245 0.001232 0.00002115];%Suzuki信道
figure(1)
semilogy(x1,y1,'^k',x2,y2,'-k',x3,y3,'-.k',x4,y4,'-..k'),grid,hold on
semilogy(x5,y5,'-*k')
semilogy(x6,y6,'-sk')
semilogy(x7,y7,'-+k',x8,y8,'-ok')
legend('高斯信道','單徑信道（視距外、近場開闊）','兩徑信道（視距內、近場開闊）','四徑信道（遮擋、近場開闊）','萊斯平坦信道K=3.8dB（視距內、近場不開闊）','萊斯平坦信道K=8.1dB（視距內、近場不開闊）','瑞利平坦信道10Hz（視距外、近場不開闊）','Suzuki信道（遮擋、近場不開闊）')
xlabel('Eb/No(dB)')
ylabel('pi/4-QPSK超短波系統誤比特率')
x9=0:5:20;
y9=[0.4 0.13 0.008 0.00053 0.0000125];%瑞利平坦信道0.3Hzqpsk
x12=5:5:20;
y12=[0.064 0.0044 0.00034 0.000009];%qpsk
x10=0:5:35;
%x10=0:4:24;
y10=[0.4 0.2 0.05 0.014 0.0058 0.003 0.0025 0.001];%瑞利選擇信道
%y10=[0.47 0.284 0.055 0.00439 0.00034 0.0000847 0.0000139];
x11=0:5:30;
y11=[0.4 0.13 0.0356 0.01 0.0009 0.000445 0.00005];%瑞利平坦信道0.3Hz8psk
x13=[5 10 15 20 30];
y13=[0.086 0.035 0.01 0.001 0.00018];%8psk
figure(2)
semilogy(x9,y9,'-k',x10,y10,'-ok',x11,y11,'-.k',x12,y12,'sk',x13,y13,'^k'),grid
legend('qpsk瑞利平坦信道0.3Hz','qpsk瑞利選擇信道','8psk瑞利平坦信道0.3Hz','qpsk系統美軍標中的性能標準','8psk系統美軍標中的性能標準')
xlabel('Eb/No(dB)')
ylabel('短波系統誤比特率')