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

[tonipe]

Bonjour

J'ai un soucis pour faire le calcul avec 3 coupures pour un HP.
Il s'agit, dans << //calculs grave complet >> des lignes $i et $j
Il y a 4 parties :
- Coupure électrique passe haut
- Coupure électrique passe bas
- Coupure acoustique passe haut
- Somme des trois coupures.

_fg_ concerne le grave
_fga_ concerne le grave, sur la coupure acoustique
_ph concerne le passe haut
_pb concerne le passe bas

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);

Merci pour votre aide.