feuille calcul simulation enceinte 2 voies d'après mesures

Code :
// Calculs grave passe haut

$puls_fg_ph = $f/$f_fg_ph;

$ct_fg_ph = $a0_fg_ph-($a2_fg_ph*pow($puls_fg_ph,2))+($a4_fg_ph*pow($puls_fg_ph,4))-($a6_fg_ph*pow($puls_fg_ph,6))+($a8_fg_ph*pow($puls_fg_ph,8));

$dt_fg_ph = +($a1_fg_ph*$puls_fg_ph)-($a3_fg_ph*pow($puls_fg_ph,3))+($a5_fg_ph*pow($puls_fg_ph,5))-($a7_fg_ph*pow($puls_fg_ph,7));

$et_fg_ph = +$b0_fg_ph-($b2_fg_ph*pow($puls_fg_ph,2))+($b4_fg_ph*pow($puls_fg_ph,4))-($b6_fg_ph*pow($puls_fg_ph,6))+($b8_fg_ph*pow($puls_fg_ph,8));

$ft_fg_ph = +($b1_fg_ph*$puls_fg_ph)-($b3_fg_ph*pow($puls_fg_ph,3))+($b5_fg_ph*pow($puls_fg_ph,5))-($b7_fg_ph*pow($puls_fg_ph,7));

$gt_fg_ph = +($ct_fg_ph*$et_fg_ph+$dt_fg_ph*$ft_fg_ph)/(pow($et_fg_ph,2)+pow($ft_fg_ph,2));

$ht_fg_ph = +($dt_fg_ph*$et_fg_ph-$ft_fg_ph*$ct_fg_ph)/(pow($et_fg_ph,2)+pow($ft_fg_ph,2));

// Calculs grave passe bas

$puls_fg_pb = $f/$f_fg_pb;

$ct_fg_pb = $a0_fg_pb-($a2_fg_pb*pow($puls_fg_pb,2))+($a4_fg_pb*pow($puls_fg_pb,4))-($a6_fg_pb*pow($puls_fg_pb,6))+($a8_fg_pb*pow($puls_fg_pb,8));

$dt_fg_pb = +($a1_fg_pb*$puls_fg_pb)-($a3_fg_pb*pow($puls_fg_pb,3))+($a5_fg_pb*pow($puls_fg_pb,5))-($a7_fg_pb*pow($puls_fg_pb,7));

$et_fg_pb = +$b0_fg_pb-($b2_fg_pb*pow($puls_fg_pb,2))+($b4_fg_pb*pow($puls_fg_pb,4))-($b6_fg_pb*pow($puls_fg_pb,6))+($b8_fg_pb*pow($puls_fg_pb,8));

$ft_fg_pb = +($b1_fg_pb*$puls_fg_pb)-($b3_fg_pb*pow($puls_fg_pb,3))+($b5_fg_pb*pow($puls_fg_pb,5))-($b7_fg_pb*pow($puls_fg_pb,7));

$gt_fg_pb = +($ct_fg_pb*$et_fg_pb+$dt_fg_pb*$ft_fg_pb)/(pow($et_fg_pb,2)+pow($ft_fg_pb,2));

$ht_fg_pb = +($dt_fg_pb*$et_fg_pb-$ft_fg_pb*$ct_fg_pb)/(pow($et_fg_pb,2)+pow($ft_fg_pb,2));

// calculs grave acoustique passe haut

$puls_fga_ph = $f/$f_fga_ph;

$ct_fga_ph = $a0_fga_ph-($a2_fga_ph*pow($puls_fga_ph,2))+($a4_fga_ph*pow($puls_fga_ph,4))-($a6_fga_ph*pow($puls_fga_ph,6))+($a8_fga_ph*pow($puls_fga_ph,8));

$dt_fga_ph = +($a1_fga_ph*$puls_fgas_ph)-($a3_fga_ph*pow($puls_fga_ph,3))+($a5_fga_ph*pow($puls_fga_ph,5))-($a7_fga_ph*pow($puls_fga_ph,7));

$et_fga_ph = +$b0_fga_ph-($b2_fga_ph*pow($puls_fga_ph,2))+($b4_fga_ph*pow($puls_fga_ph,4))-($b6_fga_ph*pow($puls_fga_ph,6))+($b8_fga_ph*pow($puls_fga_ph,8));

$ft_fga_ph = +($b1_fga_ph*$puls_fga_ph)-($b3_fga_ph*pow($puls_fga_ph,3))+($b5_fga_ph*pow($puls_fga_ph,5))-($b7_fga_ph*pow($puls_fga_ph,7));

$gt_fga_ph = +($ct_fga_ph*$et_fga_ph+$dt_fga_ph*$ft_fga_ph)/(pow($et_fga_ph,2)+pow($ft_fga_ph,2));

$ht_fga_ph = +($dt_fga_ph*$et_fga_ph-$ft_fga_ph*$ct_fga_ph)/(pow($et_fga_ph,2)+pow($ft_fga_ph,2));

// calculs grave complet

$it_fg_pbh = $gt_fg_ph*$gt_fg_pb-$ht_fg_ph*$ht_fg_pb;

$jt_fg_pbh = $ht_fg_ph*$gt_fg_pb+$gt_fg_ph*$ht_fg_pb;

$kt_fg_pbh = sqrt(pow($it_fg_pbh,2)+pow($jt_fg_pbh,2));

$lt_fg_pbh = +atan2($jt_fg_pbh,$it_fg_pbh);

$kt_fg_pb = $g_fg_pbh+20*log10($kt_fg_pbh);

$lt_fg_pb = 360*$lt_fg_pbh/2/M_PI+$bph_fg_pbh+(360/$cjmlc)*$f*$d_fg_pbh;

$mt_fg_pb = pow(10,$kt_fg_pb/20);

$nt_fg_pb = 2*M_PI*$lt_fg_pb/360;

$ot_fg_pb = $mt_fg_pb*cos($nt_fg_pb);

$pt_fg_pb = $mt_fg_pb*sin($nt_fg_pb);