开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 游戏 > 模拟服务器 > 查看源码
vquality_g711.pl
资源名称：Server_win_src.rar　[点击查看]
上传用户：wzkunzhan
上传日期：2022-04-23
资源大小：2618k
文件大小：13k
源码类别：

模拟服务器

开发平台：

Visual C++

vquality_g711.pl：源码内容
  1. # File: vquality_g711.pl
  2. #
  3. # Modified OpenPhone based Voice Evaluation Tool (MOBVET) Version 0.1
  4. # --------------------------------------------------------------------
  5. # Voice Over IP Laboratory (http://www.voip.nce.ufrj.br)
  6. # Federal University of Rio de Janeiro
  8. # Copyright (c) 2002-2004 UFRJ (Federal University of Rio de Janeiro).
  9. # The contents of this file are subject to the Mozilla Public License
  10. # Version 1.0 (the "License"); you may not use this file except in
  11. # compliance with the License. You may obtain a copy of the License at
  12. # http://www.mozilla.org/MPL/
  13. # Versao de 10/11/2003 - Leandro
  15. # Versao de 11/11/2003 - Leandro
  16. # Modificacao no codigo para considerar perda de pacotes por "buffer full" 
  18. # Versao de 12/03/2004 - Leandro
  19. # Modificado para considerar tempo de dejitter_buffer como max_time/2
  21. use Emodel;
  23. #Definindo Codec como G711
  24. VoIP::Emodel->setIe(0);
  25. VoIP::Emodel->setBpl(25.1);
  26. $CODEC = "G711";
  29. $pasta = ".";
  30. my $dir = "side?.H323";
  31. #my $dir = "sideA.H323";
  33. my $linha;
  35. print "Processing the trace files...n";
  37. chdir $pasta;
  40. open DIR, "dir /b $dir |";
  41. #open DIR, "ls $dir |";
  42. my @arqs = <DIR>;
  43. foreach $linha (@arqs) {
  46. #print "Processing file: $linha";
  48. ($nome_arq, $ext_arq) = split("H",$linha);
  49. $log_vql = $nome_arq . "vql";
  50. $log_ppl = $nome_arq . "ppl";
  51. $log_mos = $nome_arq . "mos";
  53. #print "File PPL Nanme: $log_ppln";
  54. emodel($log_ppl);
  57. markov_chain($linha);
  60. }
  61. close DIR;
  64. sub emodel {
  66. $log_ppl = shift;
  68. open(LOG_PPL, $log_ppl) or die "CAnnot open file $log_ppln";
  70. foreach $line (<LOG_PPL>) {
  73. push(@emodel_array,$line);
  76. }
  78. close(LOG_PPL);
  83. }
  86. # Recebe string no formato "mm:ss.xxxx <lixo>" ou "hh:mm:ss.xxxx <lixo>" e converte para milisegundos
  87. sub converte_para_Milisegundos
  88. {
  89.         my $arg = shift;
  91.         if ($arg =~ /(d*):(d*):(d*).(d*)/o ) {
  92.                 return ( (($1*3600)+$2*60+$3)*1000 +$4 );
  93.         } elsif ($arg =~ /(d*):(d*).(d*)/o ) {
  94.                 return ( (($1*60)+$2)*1000 +$3 );
  95.         }
  96.         printf ("Parametro Invalido: %sn", $arg);
  97. }
  99. # Recebe valor em milisegundos e converte para string no formato "mm:ss.xxxx" ou "hh:mm:ss.xxxx" 
  100. sub ms
  101. {
  102. my $arg = shift;
  104. my ($h, $r, $m, $s, $ms);
  105. $h = int ($arg / 3600000);
  106. $r = $arg % 3600000;
  107. $m = int ($r / 60000);
  108. $r = $r % 60000;
  109. $s = int ($r / 1000); 
  110. $ms = $r % 1000;
  111. if ($h > 0) {
  112. return (sprintf("%d:%02d:%02d.%03d",$h,$m,$s,$ms));
  113. } else {
  114. return (sprintf("%d:%02d.%03d",$m,$s,$ms));
  115. }
  117. }
  119. sub markov_chain {
  123. my $arquivo = shift;
  125. my $status_old;
  126. my $time_old;
  128. $e = exp(1);
  131. $pkt = 0;   # contador de pacotes recebidos
  132. $lost = 0;  # contador de pacotes perdidos
  133. $gmin = 16; # tamanho do gap m�imo, normalmente 16
  134. $c5 = 0;    # utilizaco para calcular o tempo passado desde o ultimo burst significativo de perdas 
  135. $c14 = 0;
  136. $c13 = 0;
  137. $c11 = 0;
  138. $c33 = 0;
  139. $c23 = 0;
  140. $c22 = 0;
  141. $F = 0.03; # tamanho de frame e segundos do g711 - 1 frame per packet
  142.  
  144. # #Valores para G.723.1 6.3k + VAD
  145. # $d1 = 4.25;
  146. # $d2 = 4.8;
  147. # $d3 = 12;
  148. # $d4 = 1.75;
  149. $IeCODEC = 0;
  151. $t1 = 5;
  152. $t2 = 15;
  153. $i = 0; 
  155. ($time_emodel, $rtt, $psent, $plost, $too_late, $jitter, $buffer_size) = split(" ",shift(@emodel_array));
  157. open (LOG_VQL, ">$log_vql");
  159. open(LOG_H323, $arquivo) or die "lendo arquivo $arquivo";
  161. foreach $line (<LOG_H323>) {
  165. if ($line =~ /(.*)RTP Jitter:.*Receive statistics:.*packets=(d*).*octets=(d*).*tooLate=(d*)/o ) {
  167. $status = "received";
  168. $time = ms (converte_para_Milisegundos($1));
  169. if ($status_old eq "received") {
  171.  #print "n-> PACKET RECEIVED: $time_old -----------------------------------------n";
  172.  receive_event();
  173.  
  174.  
  175. }
  179. }elsif ($line =~ /(.*)RTP Jitter:.*Droppeds(d*)s/o ){
  181. $status = "dropped";
  182. $time = ms (converte_para_Milisegundos($1));
  183. $dropped = $2;
  185. #Considera pacotes multiframe 
  186. for ($i = 0; $i < $dropped ; $i++) {
  188. #print "n-> PACKET LOST: $time ---------------------------------------------n";
  189. loss_event();
  192. }
  195. }elsif ($line =~ /(.*)Jitter buffer oldest packet.*too late/o ){
  196.  
  197.   $status = "late";
  198. $time = ms (converte_para_Milisegundos($1));
  200. if ($status_old eq "received") {
  202. #print "n-> PACKET TOO LATE: $time_old ------------------------------------------n";
  203. loss_event();
  207. }
  208. # 2:33.566      RTP Jitter:80786e0 RTP Jitter buffer full, throwing away oldest frame (31920)
  209. }elsif ($line =~ /(.*)Jitter buffer full, throwing away oldest frame/o ){
  211. $status = "late";
  212. $time = ms (converte_para_Milisegundos($1));
  214. if ($status_old eq "received") {
  216. #print "n-> BUFFER FULL: $time_old ------------------------------------------n";
  217. loss_event();
  219. }
  220. }elsif($line =~ /(.*)SentReceiverReport:s*ssrc=(d*)s*fraction=(d*)s*lost=(d*)s*last_seq=(d*)s*jitter=(d*)s*lsr=(d*)s*dlsr=(d*)/o ) {
  222. $time = ms (converte_para_Milisegundos($1));
  223. if ($time == $time_emodel) {
  225. calcula_markov();
  226. $ta = onewaydelay($rtt,$CODEC,$buffer_size);
  228. VoIP::Emodel->setTa($ta);
  229. if ($psent != 0){
  231. VoIP::Emodel->setPpl(($plost + $too_late ) / $psent);
  232. }
  233. else {
  235. #Caso nenhum pacote tenha sido enviado
  236. VoIP::Emodel->setPpl(0);
  237. }
  239. #print ("Ppl: ",VoIP::Emodel->getPpl(),"  Bpl: ",VoIP::Emodel->getBpl(), " Ie: ",VoIP::Emodel->getIe(), "n");
  240. $r_emodel = VoIP::Emodel->r();
  241. $ie_emodel = VoIP::Emodel->getIeef();
  242. $id_emodel = VoIP::Emodel->getId();
  243. $mos_emodel = VoIP::Emodel->mos(VoIP::Emodel->r());
  245. $id_emodel_ext = $id_emodel; #ler o artigo outra vez para ajustar o valor de id
  246. $r_emodel_ext = $R - $id_emodel_ext;
  248. #print "$R $r_emodel_ext $id_emodel_extn";
  250. if ($r_emodel_ext < 0 ) {
  251. #nao deixa o r do emodel extendido ficar negativo
  252. $r_emodel_ext = 0;
  253. } 
  254. #print LOG_VQL ($time_emodel," ",$ie_emodel," ",$id_emodel," ",$r_emodel," ",$mos_emodel," ",$IeTotal," ",$id_emodel_ext," ",$r_emodel_ext," ",VoIP::Emodel->mos($r_emodel_ext)," ",$I1," ",$I2," ",$Ieav,"n");
  255. print LOG_VQL ($time_emodel," ",$ie_emodel," ",$id_emodel," ",$r_emodel," ",$mos_emodel," ",$IeTotal," ",$id_emodel_ext," ",$r_emodel_ext," ",VoIP::Emodel->mos($r_emodel_ext),"n");
  258. ($time_emodel, $lixo, $psent, $plost, $too_late, $jitter, $buffer_size) = split(" ",shift(@emodel_array));
  259. # a variavel lixo eh utilizada para nao se atualizar o rtt
  260. # o owd eh calculado sempre em relacao ao primeiro rtt calculado
  261. # o que varia eh o tamaho do buffer de compensacao de jitter
  263. }
  265. }
  267. $status_old = $status;
  268. $time_old = $time;
  272. }
  277. # calculo do MOS final da ligacao
  278. $Ie_end = $Ieav + (0.7 * ($I1 - $Ieav)) * $e ** (-$y/30);
  280. #Verificar como adicionar Id no R final da ligacao
  281. $R = 93.36 -  $Ie_end - $IeCODEC - $id_emodel_ext;
  282. #print "###########################################################n";
  283. #print "AVALIACAO FINAL DA LIGACAOn";
  284. #print "Ie(End of call) = $Ie_endn";
  285. #print "R-Factod(End of call) = $Rn";
  286. #print "MOS(End of call)  = ", mos($R),"n";
  287. open (LOG_MOS, ">$log_mos");
  288. print LOG_MOS (VoIP::Emodel->mos($R),"n");
  289. close (LOG_MOS);
  292. close(LOG_H323);
  293. close(LOG_VQL);
  295. }
  297. sub receive_event {
  298. # evento de recepcao de pacote
  300. $pkt++;
  302. }
  304. sub loss_event {
  305. # evento de perda de pacote
  307. $c5 = $c5 + $pkt;
  309. if ($pkt >= $gmin) {
  311.                 if ($lost == 1) {
  312.                
  313.                 $c14++;
  314.                 }
  315.                 else {
  316.                
  317.                 $c13++;
  318.                 }
  319.                
  320.                 $lost = 1;
  322.                 $c11 = $c11 + $pkt;
  324.         }
  325.         else {
  326.       
  327.       
  328.                 $lost++;
  329.                 if ($lost > 8) {
  330.                 
  331.                      $c5 = 0;
  332.                 }
  333.                 if( $pkt == 0) {
  334.                 $c33++;
  336. }
  337. else {
  339.                  $c23++;
  340.                  $c22 = $c22 + $pkt;
  341. }
  345. }
  347. $pkt = 0;
  349. }
  353. sub calcula_markov {
  355. $c31 = $c13;
  356. $c32 = $c23;
  357. if ( ($c11 + $c14 + $c13) == 0 ) {
  359. $p11 = 0;
  360. $p13 = 0;
  361. }
  362. else {
  363. $p11 = ( $c11 + $c14 ) / ( $c11 + $c14 + $c13);
  364. $p13 = $c13 / ( $c11 + $c14 + $c13);
  365. }
  367. #$p13 = 1 - $p11;
  370. if ( ( $c31 + $c32 + $c33) ==0 ) {
  372. $p31 = 0;
  373. $p32 = 0;
  374. $p33 = 0;
  376. }
  377. else {
  379. $p31 = $c31 / ( $c31 + $c32 + $c33);
  380. $p32 = $c32 / ( $c31 + $c32 + $c33);
  382. #$p33 = 1 - $p31 - $p32;
  383. $p33 = $c33 / ( $c31 + $c32 + $c33);
  384. }
  387. if ( ($c22 + $c23) == 0 ) {
  388. $p22 = 0;
  389. $p23 = 0;
  390. } else {
  391. $p22 = $c22 / ( $c22 + $c23 );
  392. $p23 = $c23 / ( $c22 + $c23 );
  393. }
  395. #$p23 = 1- $p22;
  398. $d = ($p23 * $p31 + $p13 * $p32 + $p13 * $p23);
  400. if ($d == 0) {
  401. $p1 = 0;
  402. $p2 = 0;
  403. $p3 = 0;
  404. }
  405. else {
  406. $p1 = ($p31 * $p23) / $d;
  407. $p2 = ($p13 * $p32) / $d;
  408. $p3 = ($p13 * $p23) / $d;
  409. }
  414. $L = 100 * $p3; # porcentagem media da taxa de pacotes perdidos
  416. $y = $F * $c5; # tempo passado desde o ultimo burst significativo de perdas em segundos
  418. if ( $p11 == 1) { 
  419. # significa que nao apareceu burst
  421. $g = $y;
  423. }else{
  425. $g = $F / ( 1 - $p11); # tamanho do gap em segundos
  426. }
  428. if ($c11 == 0) {
  429. $Dg = 0;
  430. }
  431. else {
  433. $Dg = 100 *  $c14 / ( $c11 + $c14); #densidade de perda no gap (porcentagem)
  434. }
  436. if ( ($p1 * $p13) == 0 ) {
  437. $b = 0;
  438. }
  439. else {
  440. $b = ( $F * (1 - $p1) ) / ( $p1 * $p13); # tamanho do burst em segundos
  441. }
  443. if ( ($p23 + $p32) == 0 ) {
  445. $Db = 0;
  446. }
  447. else {
  449. $Db = 100 * $p23 / ($p23 + $p32); # densidade de perda no burs (porcentagem)
  450. }
  451. # Ao meu ver, esta forma de calcular Db estah errada
  452. #$Db = 100 * $p23;
  457. #print "Average packet loss rate = $L %n";
  458. #print "Gap Length  = $g seg.n";
  459. #print "Gap loss density = $Dg %n";
  460. #print "Burst length = $b seg.n";
  461. #print "Burst loss density = $Db %n";
  462. #print "Delay since last burst = $y seg.n";
  464. # if ($Dg < 0.5) {
  465. #
  466. # $Ieg = 0;
  467. #
  468. # }
  469. # elsif ( (0.5 < $Dg) && ($Dg < $d2) ) {
  470. #
  471. # $Ieg = $d1 * $Dg;
  472. #
  473. # }
  474. # elsif ($d2 < $Dg) {
  475. #
  476. # $Ieg = $d3 + $d4 * $Dg;
  477. #
  478. # }
  479. #
  480. # if ($Db < 0.5) {
  481. #
  482. # $Ieb = 0;
  483. #
  484. # }
  485. # elsif ( (0.5 < $Db) && ($Db < $d2) ) {
  486. #
  487. #
  488. # $Ieb = $d1 * $Db;
  489. # }
  490. # elsif ($d2 < $Db) {
  491. #
  492. #
  493. # $Ieb = $d3 + $d4 * $Db;
  494. # }
  495. #
  496. VoIP::Emodel->setPpl($Dg/100);
  497. VoIP::Emodel->r();
  498. $Ieg = VoIP::Emodel->getIeef();
  500. VoIP::Emodel->setPpl($Db/100);
  501. VoIP::Emodel->r();
  502. $Ieb = VoIP::Emodel->getIeef();
  504. #print "Dg = $Dgn";
  505. #print "Db = $Dbn";
  508. #print "Ieg = $Iegn";
  509. #print "Ieb = $Iebn";
  511. if ( (1 - $e ** (-$b/$t1 - $g/$t2)) == 0) {
  512. $I2 = 0;
  514. }
  515. else {
  517. $I2 = ( $Ieg * ( 1 - $e ** (-$g /$t2)) + $Ieb * ( 1 - $e** (-$b/$t1)) * $e ** (-$g/$t2)) / (1 - $e ** (-$b/$t1 - $g/$t2) );
  519. }
  521. #$I1 = $Ieb - ($Ieg - $I2) * $e ** (-$b/$t1);
  522. $I1 = $Ieb - ($Ieb - $I2) * $e ** (-$b/$t1); # Our considerations
  524. $Ieav = ($b * $Ieb + $g * $Ieg - $t1 * ($Ieb - $I2) * (1 - $e ** (-$b/$t1)) + $t2 * ($I1 - $Ieg) * (1 - $e ** (- $g/$t2)) ) / ($b + $g);
  528. #print "I1 = $I1n";
  529. #print "I2 = $I2n";
  530. #print "g = $gn";
  531. #print "b = $bn";
  532. #print "Ie(av) = $Ieavn";
  533. $R = 94 -  $Ieav - $IeCODEC;
  534. $IeTotal = $Ieav + $IeCODEC;
  535. if ($IeTotal > 100) {
  536. #nao deixa o ietotal ser negativo 
  537. $IeTotal = 100;
  538. }
  540. #print "R-factior = $Rn";
  541. #print "MOS = ", mos($R),"  at  $time n";
  545. }
  548. sub onewaydelay
  549. #Parametros (rtt, codec, playoutbuffer)
  550. #Calcula o On Way Delay, tambe'm chamado de Absolute Delay (Ta).
  551. #Ta = RTT/2 + tempo de codificacao (dependente do codec utilizado) + tempo de playout buffer (fixo)
  553. {
  554. my $rtt = shift;
  555. my $codec = shift;
  556. my $plbuffer = shift;
  557. my $codec_time = 0;
  559. #confirmar tempo de codificacao
  561. $plbuffer = $plbuffer/2; #codecao do dejitter_buffer como max_time/2
  563. if ($codec eq "G711") { $codec_time = 0.25}
  564. elsif ($codec eq "GSM") { $codec_time = 40}
  565. elsif ($codec eq "G729") { $codec_time = 25}
  566. elsif ($codec eq "G723") { $codec_time = 157.5} #4 frames por pacote no caso do openam # ver g.114
  567. #considerando link de alta velocidade (melhor caso)
  568. #elsif ($codec eq "G723") { $codec_time = 277.5} #4 frames por pacote no caso do openam # ver g.114
  569. #considerando link de baixa velocidade (pior caso)
  570. #$a = ($rtt/2) + $plbuffer + $codec_time;
  571. #print "rtt= $rtt, buffer= $plbuffer, codec= $codec_time, total= $an";
  572. return ($rtt/2) + $plbuffer + $codec_time;
  573. }