ch8example20main.m

来自「清华大学出版社 邵玉斌编写的《通信系统建模与仿真实例分析》一书的所有MATLAB」· M 代码 · 共 25 行

% ch8example20main.m
clear;clc;
data=[-5      0.29486 
    0       0.25006 
    5       0.16105 
    10      0.05192 
    15      0.00131 
    16      0.00036 
    17      0.00008 
    18      0.00009 
    19      0.00002];                  % 测试数据
SNR_dB=data(:,1); Pe_simu=data(:,2);
SNR=10.^(SNR_dB./10);                  % 测试信噪比范围
semilogy(SNR_dB, Pe_simu, 'o');hold on;% 测试误码率作图
Pe_0=0.99;                             % 信源符号概率猜测
[pp1,Jmin]=lsqcurvefit('ch8example20func',0.54,SNR,Pe_simu)% 拟合指令1
[pp2,r,J]=nlinfit(SNR,Pe_simu,'ch8example20func',Pe_0)     % 拟合指令2
ci = nlparci(pp2,r,J)                  % 拟合结果的置信区间
SNRThdB=-5:0.1:20;                     % 拟合曲线计算和作图
SNRTh=10.^(SNRThdB./10);
PeTh=ch8example20func(pp2,SNRTh);
semilogy(SNRThdB,PeTh,'k');
xlabel('A^2/\sigma^2 (dB)');
ylabel('错误率 P_e');
legend('仿真结果','拟合曲线');