[idriss11]

السلام عليكم

clc
% close all
clear all
N = 64;
M= 16;
pfa = 10^-4;
Tc = 10^-6;
Td=N*Tc;
q0=1023;
K=1000;
k = 10000;
TCA =Fzero(*tempCA,0.5) ;
s = [-25:2:5];

for i = 1:16
snr1 = s(i);
snr = 10^(snr1/10);
S = N*sqrt(snr);
sigman = N/2;
S = S^2/sigman
cCA=0;
for j = 1:k
R = chi2rnd(2,1,M);
R1 = ncx2rnd(2,S,1,1);
XCA = sum (R);
q = R1;
if q >= XCA*TCA
cCA = cCA + 1;
end
end
pdCA(i)= cCA/k
TacqCA(i)=((2-pdCA(i))*(1+(K*pfa))*q0*Td)/(2*pdCA(i));
end
hold on
figure(1)

plot (s,pdCA)
hold on

figure(2)
plot (s,TacqCA)

[idriss11]

% calcule de TCA
function pfaCA = temp(TCA);
M =16;
pfaCA =((1+TCA)^-M)-10^-4;