... = 10.^(gamma_db/10); M = 16; Ps_exact=1-exp(2*log(( 1-( (2*(sqrt(M )-1 ))/(sqrt(M)))*Q(sqrt((3*gamma)/(M-1)))))); Ps_approx = ((4*(sqrt(M )-1 ))/sqrt(M))*Q(sqrt((3*gamma)/(M-1))); semilogy(gamma_db, ... 19.89µs Rb ≥ 50. 26 Kbps Combining the two, we have 50. 26 Kbps ≤ Rb ≤ 200 Kbps (or Mbps) 22 From figure 6. 5 with Pb = 10−3 , BPSK d = θTm /Ts , θTm = 3µs d = × 10−2 Ts = 60 µsec R = 1/Ts = 16. 7Kbps QPSK ... %Initializations avg_SNR = [0:0.1:20]; gamma_b_bar = 10.^(avg_SNR/10); m = [1 4]; line = [’-k’, ’-r’, ’-b’] for i = 1:size(m,2) for j = 1:size(gamma_b_bar, 2) Pb_bar(i,j) = (1/pi)*quad8(’nakag_MGF’,0,pi/2,[],[],gamma_b_bar(j),m(i),1);...