comp_exam10_2.m

来自「the source code of Ziemer s book Digital」· M 代码 · 共 18 行

clear all						% be safe
g = 0.2;							% gamma
zdB = 0:0.1:10;				% z in dB
z = 10.^(zdB/10);				% vector of z values
q1 = qfn((1-g)*sqrt(2*z));
q2 = qfn((1+g)*sqrt(2*z));
qt = qfn(sqrt(2*z));			% Gamma=0 case
p2 = q1-q2;						% P2
p1 = q2;							% P1
pe = p1./(1-p2);				% probability of error
semilogy(zdB,pe,'k',zdB,qt,'k' )
xlabel('z - dB')
ylabel('Probability of Error')
pause
N = 1./(1-p2);
plot(zdB,N,'k')
xlabel('z - dB')
ylabel('Expected Number of Transmissions')