%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Sous-Programme calcul_pdf.m %%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function K_PDF=calcul_pdf(gain_K_PDF,phase_K_PDF,wc)

phi_d=phase_K_PDF;
wd=wc;
a=(1+sin(phi_d))/(1-sin(phi_d));
tau_d=1/(wd*sqrt(a));
k_d=gain_K_PDF/sqrt(a);
num_PDF=k_d*[a*tau_d 1];
den_PDF=[tau_d 1];
K_PDF=tf(num_PDF,den_PDF);