qpsk_ber_sevelchannel.m

这是一个仿真qpsk信号通过不同信道模型（包括单径

%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('短波系统误比特率')