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

源码开发语言/平台
当前位置：首页 > 源码/资料 > Windows编程 > 通讯编程 > 查看源码
PER.cc
资源名称：ns-allinone-2.33.tar.gz　[点击查看]
上传用户：rrhhcc
上传日期：2015-12-11
资源大小：54129k
文件大小：17k
源码类别：

通讯编程

开发平台：

Visual C++

PER.cc：源码内容
  1. /*
  2.  * Copyright (c) 2007 Regents of the SIGNET lab, University of Padova.
  3.  * All rights reserved.
  4.  *
  5.  * Redistribution and use in source and binary forms, with or without
  6.  * modification, are permitted provided that the following conditions
  7.  * are met:
  8.  * 1. Redistributions of source code must retain the above copyright
  9.  *    notice, this list of conditions and the following disclaimer.
  10.  * 2. Redistributions in binary form must reproduce the above copyright
  11.  *    notice, this list of conditions and the following disclaimer in the
  12.  *    documentation and/or other materials provided with the distribution.
  13.  * 3. Neither the name of the University of Padova (SIGNET lab) nor the 
  14.  *    names of its contributors may be used to endorse or promote products 
  15.  *    derived from this software without specific prior written permission.
  16.  *
  17.  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 
  18.  * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED 
  19.  * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 
  20.  * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR 
  21.  * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, 
  22.  * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, 
  23.  * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; 
  24.  * OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, 
  25.  * WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR 
  26.  * OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF 
  27.  * ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  28.  */
  30. /* -*- Mode:C++; c-basic-offset:3; tab-width:3; indent-tabs-mode:t -*- */
  32. #include "mac.h"
  33. #include "mac-802_11mr.h"
  34. #include "PER.h"
  35. #include <math.h>
  36. #include <float.h>
  37. #include <rng.h>
  39. #define ABS(x) (((x)>=0) ? (x) : (-x) )
  41. //Completion for the ADT
  44. /// Container for the length
  45. class leng : public Object {
  46. public:
  47. int len;
  48. DLList* snr_list;
  50. leng(int len_){
  51. len=len_;
  52. snr_list= new DLList();
  53. }
  55. ~leng() {
  56. delete snr_list;
  57. }
  59. };
  62. /// Container for SNR and PER
  63. class Snr : public Object {
  64. public:
  65. double snr;
  66. double per;
  68. Snr(double snr_, double per_){
  69. snr= snr_;
  70. per= per_;
  71. }
  73. };
  76. /**
  77.  * TCL wrapper for the PER class
  78.  * 
  79.  */
  80. static class PERClass : public TclClass {
  81. public:
  82. PERClass() : TclClass("PER") {}
  83. TclObject* create(int, const char*const*) {
  84. return (new PER());
  85. }
  86. } class_per;
  89. PER::PER()
  90. {
  91. for(int i = 0 ; i < NumPhyModes ; i++)
  92. PHYmod[i] = new DLList();
  94. #ifndef COMPILING_TEST_PROGRAM
  95. bind("debug_", &debug_);
  96. bind("noise_", &noise_);
  97. #endif
  99. }
  101. PER::~PER(){
  102.   for(int i = 0 ; i < NumPhyModes ; i++)
  103. delete PHYmod[i];
  104. }
  108. int PER::command(int argc, const char*const* argv){
  110. if(debug_>2) printf("Called PER::command argc=%in",argc);
  112. if (argc == 2) {
  113. if (strcmp(argv[1], "print") == 0) {
  114. print_data();
  115.  return TCL_OK;
  116. }
  117. } else if (argc == 6) {
  119. // Command to add PER entries to PER table
  120. if (strcmp(argv[1], "add") == 0) {
  121.                        
  122. PhyMode pm = str2PhyMode(argv[2]);
  124. if (pm==ModeUnknown) return(TCL_ERROR);
  126. int len =   atoi(argv[3]);
  127. double snr = atof(argv[4]);
  128. double  per = atof(argv[5]);
  130. if(debug_>2)  printf("%s=%f ------- %s=%fn",argv[4],snr,argv[5],per);
  132. set_per(pm,len,snr,per);
  133. return TCL_OK;    
  134. }
  135. }
  136. }
  140. int PER::get_err(PowerProfile* pp, Packet* pkt, double* snrptr, double* snirptr, double* interfptr){
  143. if(debug_>3)
  144. {
  146. printf("PER::get_err(): dumping everythingn");
  148. for (int i=0; i< 12; i++)
  149. {
  150. if(PHYmod[i]!=0)
  151. {
  152. NSNode* nl  = (NSNode*)(PHYmod[i]->first());
  153. while(nl !=0  )
  154. {
  155. cout << ((leng*)nl->element())->len << "n";
  156. DLList* sl = ((leng*)nl->element())->snr_list;
  157. NSNode* ns =(NSNode*)sl->first();
  158. while( ns !=0 )
  159. {
  160. cout << ((Snr*)ns->element())->snr << "  " ;
  161. cout << ((Snr*)ns->element())->per << "n" ;
  162. ns  = (NSNode*)(sl->after(ns));
  163. }
  164. nl  = (NSNode*)(PHYmod[i]->after(nl));
  165. }
  169. }
  170. }
  171. }
  173. /* Hack needed not to link the whole NS with the test program */ 
  174. #ifdef COMPILING_TEST_PROGRAM         
  175. double pkt_pow = 0 ;
  176. int len = 0;
  177. PhyMode pm=Mode1Mb;
  178. double duration = 1;
  179. #else        
  180. struct hdr_cmn *hdr = HDR_CMN(pkt);
  181. struct MultiRateHdr *hdm = HDR_MULTIRATE(pkt);
  182. double pkt_pow = pkt->txinfo_.RxPr;
  183. double duration = hdr->txtime(); 
  184. PhyMode pm = hdm->mode(); 
  185. int len = hdr->size(); 
  186. #endif
  188. double avgpower = avg(pp,duration);
  189. double interf = avgpower-pkt_pow;
  190. double snr1=DBL_MAX; 
  191. double snr2;
  192. double per1;
  193. double per2;
  196. if (ABS(noise_) > DBL_MIN)
  197. {
  198. snr1 = 10*log10(pkt_pow/(noise_));
  199. per1 = get_per( pm, snr1, len);
  200. }
  201. else
  202. per1 = 0;
  204. if (ABS(noise_ + interf) > DBL_MIN)
  205. {
  206. snr2 = 10*log10(pkt_pow/(noise_ + interf));
  207. per2 = get_per( pm, snr2, len);
  208. }
  209. else
  210. per2=0;
  212. /* Return SNR */
  213. if (snrptr!=NULL)
  214. *snrptr = snr1;
  216. /* Return SNIR */
  217. if (snirptr!=NULL)
  218. *snirptr = snr2;
  220. /* Return interfering power */
  221. if (interfptr!=NULL)
  222. *interfptr = interf;
  224. // double x = (double) rand()/(double)RAND_MAX;
  225. double x = RNG::defaultrng()->uniform_double();
  227. //if ((debug_)&&(interf>1e-15))
  228. if ((debug_))
  229. {
  230. cout << "PER::get_err() at " << Scheduler::instance().clock() << "sn" ;
  231. cout << "PER::get_err() Average received power      "<< avgpower << " Wn"; 
  232. cout << "PER::get_err() Signal power                "<< pkt_pow << " Wn";
  233. cout << "PER::get_err() Noise power                 "<< noise_ << " Wn";
  234. cout << "PER::get_err() Interference power          "<< interf<< " Wn";
  235. cout << "PER::get_err() SNR (noise only)            "<< snr1<< " dBn";
  236. cout << "PER::get_err() SNIR (noise + interference) "<< snr2<< " dBn";
  237. cout << "PER::get_err() PER due to noise only "<< per1 << "n"; 
  238. cout << "PER::get_err() PER due to noise and interference "<< per2<< "n";
  239. cout << "PER::get_err() Rand "<< x<< "n";
  240. }
  242.   if (x<per1) 
  243. return PKT_ERROR_NOISE; 
  244. else if (x<per2) 
  245. return PKT_ERROR_INTERFERENCE ; 
  246. else 
  247. return PKT_OK; 
  249. }
  254. int PER::get_per(PhyMode pm, double tsnr, double tlen, per_set* ps_){
  256. DLList* snrlist;
  257. NSNode* pl1 =0;
  258. NSNode* pl2 =0;
  259. NSNode* ps1 =0;
  260. NSNode* ps2 =0;
  261. NSNode* ps3 =0;
  262. NSNode* ps4 =0;
  264. //cerco la lunghezza
  265. pl1 = (NSNode*)PHYmod[pm]->first();
  266.  
  267. if(pl1 == 0) {
  268.   cerr << "No PER entry found for PHYMode " << PhyMode2str(pm) << endl;
  269. return 1;
  270. }
  272. while (pl1 != PHYmod[pm]->last() && ((leng*)pl1->element())->len < tlen ){
  273. pl1=(NSNode*)PHYmod[pm]->after(pl1);
  274. }
  275. //pl1 o e' il primo maggiore o uguale di tlen o e' l'ultimo
  277. //prendo quello precedente (potrebbe essere nullo)
  278. if (((leng*)pl1->element())->len > tlen)
  279. pl2 = (NSNode*)PHYmod[pm]->before(pl1);
  282. snrlist = ((leng*)pl1->element())->snr_list;
  283. //cerco snr nella prima lunghezza
  284. ps1 = (NSNode*)snrlist->first();
  286. if(ps1 == 0) {
  287. cerr<<"Tabella non inizializzata 2";
  288. return 1;
  289. }
  290. while ( ps1 != snrlist->last()  && ((Snr*)ps1->element())->snr < tsnr  ) {
  291. ps1=(NSNode*)PHYmod[pm]->after(ps1);
  292. }
  293. //ps1 o e' maggiore di tsnr o e' l'ultimo
  295. //prendo quello precedente
  296. if (((Snr*)ps1->element())->snr > tsnr)
  297. ps2 = (NSNode*)snrlist->before(ps1);
  300. //se le lunghezze erano due diverse cerco snr anche nella seconda
  301. if ((pl1 != pl2) & (pl2 != 0)){
  302. //cerr << "dicversi";
  303. snrlist = ((leng*)pl2->element())->snr_list;
  304. //cerco
  305. ps3 = (NSNode*)snrlist->first();
  307. if(ps3 == 0) {
  308. cerr<<"Tabella non inizializzata 3";
  309. return 1;
  310. }
  312. while (not (((Snr*)ps3->element())->snr >= tsnr or ps3 == snrlist->last() ) ){
  313. ps3=(NSNode*)PHYmod[pm]->after(ps3);
  314. }
  315. //prendo anche quello precedente
  316. if (((Snr*)ps3->element())->snr > tsnr)
  317. ps4 = (NSNode*)snrlist->before(ps3);
  319. }
  322. //settaggio delle variabili di uscita in base ai risultati della ricerca. 
  323. //I valori /1 indicano che no e
  324. if (ps1){
  325. //cerr << ps_;
  326. ((Snr*)ps1->element())->snr;
  328. ps_->snr[0] = ((Snr*)ps1->element())->snr;
  329. ps_->per[0] = ((Snr*)ps1->element())->per;
  330. ps_->len[0] = ((leng*)pl1->element())->len;
  331. }else{
  332. ps_->snr[0] = -1;
  333. ps_->per[0] = -1;
  334. ps_->len[0] = -1;
  335. }
  338. if (ps2){
  339. ps_->snr[1] = ((Snr*)ps2->element())->snr;
  340. ps_->per[1] = ((Snr*)ps2->element())->per;
  341. ps_->len[1] = ((leng*)pl1->element())->len;
  342. }else{
  343. ps_->snr[1] = -1;
  344. ps_->per[1] = -1;
  345. ps_->len[1] = -1;
  346. }
  348. if (ps3){
  349. ps_->snr[2] = ((Snr*)ps3->element())->snr;
  350. ps_->per[2] = ((Snr*)ps3->element())->per;
  351. ps_->len[2] = ((leng*)pl2->element())->len;
  352. }else{
  353. ps_->snr[2] = -1;
  354. ps_->per[2] = -1;
  355. ps_->len[2] = -1;
  356. }
  358. if (ps4){
  359. ps_->snr[3] = ((Snr*)ps4->element())->snr;
  360. ps_->per[3] = ((Snr*)ps4->element())->per;
  361. ps_->len[3] = ((leng*)pl2->element())->len;
  362. }else{
  363. ps_->snr[3] = -1;
  364. ps_->per[3] = -1;
  365. ps_->len[3] = -1;
  366. }
  368. if (debug_ > 2) {
  369. int i;
  370. printf("Pointers: %x %x %x %xn",ps1,ps2,ps3,ps4);
  371. fprintf(stderr,"4 PER table entries: n");
  372. for(i=0; i<4; i++)
  373. fprintf(stderr,"t %d %f %f  n", ps_->len[i], ps_->snr[i], ps_->per[i] );
  375. }
  377. }
  380. void PER::set_per(PhyMode phy, int l , double snr, double per)
  381. {
  383. NSNode* nl  = (NSNode*)(PHYmod[phy]->first());
  385. if (nl==0){
  386. leng* le = new leng(l);
  387. PHYmod[phy]->insertFirst(le);
  388. le->snr_list->insertFirst(new Snr(snr,per));
  391. }else{
  392. while(nl !=0 && l>((leng*)nl->element())->len)
  393. nl  = (NSNode*)(PHYmod[phy]->after(nl));
  395. if (nl==0){
  396. leng* le =  new leng(l);
  397. PHYmod[phy]->insertLast(le);
  398. le->snr_list->insertFirst(new Snr(snr,per));
  400. }else
  401. {
  403. if (((leng*)nl->element())->len==l){
  404. DLList* sl = ((leng*)nl->element())->snr_list;
  405. NSNode* ns =(NSNode*)sl->first();
  406. while( ns !=0 && ((Snr*)ns->element())->snr<snr )
  407. ns  = (NSNode*)(sl->after(ns));
  408. if (ns==0){
  409. sl->insertLast(new Snr(snr,per));
  410. }
  411. else{
  412. sl->insertBefore(ns,new Snr(snr,per));
  413. }
  415. }else{
  417. leng* le =  new leng(l);
  418. PHYmod[phy]->insertBefore(nl,le);
  419. le->snr_list->insertFirst(new Snr(snr,per));
  421. }
  422. }
  423. }
  424. }
  429. double PER::get_per(PhyMode pm, double tsnr, double tlen){
  431. struct per_set ps_;
  433. get_per( pm,  tsnr, tlen, &ps_);
  436. int len1 = -1;
  437. int len2 = -1;
  438. double per1 = -1;
  439. double per2 = -1;
  440. double per = -1;
  443. //not found entries are set equal  to the clesest one
  445. //cerr <<  ps_.len[0] << ps_.len[1] << ps_.len[2] << ps_.len[3];
  447. if(ps_.len[0]==-1){
  448. ps_.len[0]=ps_.len[1];
  449. ps_.snr[0]=ps_.snr[1];
  450. ps_.per[0]=ps_.per[1];
  451. }
  453. if(ps_.len[1]==-1){
  454. ps_.len[1]=ps_.len[0];
  455. ps_.snr[1]=ps_.snr[0];
  456. ps_.per[1]=ps_.per[0];
  457. }
  459. if(ps_.len[2]==-1){
  460. ps_.len[2]=ps_.len[3];
  461. ps_.snr[2]=ps_.snr[3];
  462. ps_.per[2]=ps_.per[3];
  463. }
  465. if(ps_.len[3]==-1){
  466. ps_.len[3]=ps_.len[2];
  467. ps_.snr[3]=ps_.snr[2];
  468. ps_.per[3]=ps_.per[2];
  469. }
  471. if (debug_ > 2) {
  472. int i;
  473. fprintf(stderr,"4 PER table entries: n");
  474. for(i=0; i<4; i++)
  475. fprintf(stderr,"t %d %f %f  n", ps_.len[i], ps_.snr[i], ps_.per[i] );
  477. }
  480. // intepolation: 1) linear interpolation at the varying of the SNR for each length
  481. // 2) linear interpolation among lenghts 
  483. if (ps_.snr[1] != ps_.snr[0])
  484. per1 = ps_.per[0] +  (ps_.per[1]- ps_.per[0])/(ps_.snr[1]- ps_.snr[0])*(tsnr- ps_.snr[0]);
  485. else 
  486. per1 = ps_.per[1];
  488. if (ps_.snr[3] != ps_.snr[2])
  489. per2 = ps_.per[2] +  (ps_.per[3]- ps_.per[2])/(ps_.snr[3]- ps_.snr[2])*(tsnr- ps_.snr[2]);
  490. else
  491. per2 = ps_.per[2];
  493. if (per1!=-1 && per2!=-1){
  494. if (ps_.len[3] != ps_.len[1])
  495. per = per1 +  (per2- per1)/(ps_.len[3]- ps_.len[1])*(tlen-ps_.len[1]);
  496. else
  497. per = per1;
  498. }
  499. else 
  500. if (per1!=-1){
  501. per = per1;
  502. }else 
  503. if (per2!=-1){
  504. per = per2;
  505. }else {cerr << "No PER available";}
  506. return per;
  507. }
  511. //compute the mean average power given by the time power profile.
  513. double PER::avg(PowerProfile* pp,double duration){
  515.   if (debug_>1)  printf("PER::avg()  starting calculationsn");
  518. double avg_pow = 0 ;
  519. for (int i=0; i< pp->getNumElements()-1; i++){
  520. double a =  pp->getElement(i).power;
  521. //double b = pp->getElement(i+1).power;
  522. //cerr << "pot ist" << a << " "<< b << "---n"; 
  523. double dt = (pp->getElement(i+1).time -pp->getElement(i).time);
  524. avg_pow = avg_pow + a * dt;
  526. if (debug_>1)  printf("PER::avg()  power=%2.10f t dt=%fn",a,dt);
  527. }
  529. double a = pp->getElement(pp->getNumElements()-1).power;
  530. double dt = (duration -pp->getElement(pp->getNumElements()-1).time);
  532. // avg_pow = avg_pow + pp->getElement(pp->getNumElements()-1).power * (duration -pp->getElement(pp->getNumElements()-1).time);
  533. avg_pow = avg_pow + a * dt;
  535. if (debug_>1)  printf("PER::avg()  power=%2.10f t dt=%f duration=%fn",a,dt,duration);
  537. //cerr << "average0" << pp->getElement(1).time - pp->getElement(0).time;
  539. avg_pow = avg_pow / duration;
  541. if(debug_>1) printf("PER::avg() finished calculations: at %f average power=%2.10fn", Scheduler::instance().clock(),avg_pow);
  543. return avg_pow;
  544. }
  554. int PER::get_snr(PhyMode pm, double tper, double tlen, per_set* ps_){
  556. DLList* snrlist;
  557. NSNode* pl1 =0;
  558. NSNode* pl2 =0;
  559. NSNode* ps1 =0;
  560. NSNode* ps2 =0;
  561. NSNode* ps3 =0;
  562. NSNode* ps4 =0;
  564. //cerco la lunghezza
  565. pl1 = (NSNode*)PHYmod[pm]->first();
  567. if(pl1 == 0) {
  568. cerr<<"Tabella non inizializzata 1";
  569. return 1;
  570. }
  572. while (pl1 != PHYmod[pm]->last() && ((leng*)pl1->element())->len < tlen ){
  573. pl1=(NSNode*)PHYmod[pm]->after(pl1);
  574. }
  575. //pl1 o e' il primo maggiore o uguale di tlen o e' l'ultimo
  577. //prendo quello precedente (potrebbe essere nullo)
  578. if (((leng*)pl1->element())->len > tlen)
  579. pl2 = (NSNode*)PHYmod[pm]->before(pl1);
  582. snrlist = ((leng*)pl1->element())->snr_list;
  583. //cerco snr nella prima lunghezza
  584. ps1 = (NSNode*)snrlist->first();
  586. if(ps1 == 0) {
  587. cerr<<"Tabella non inizializzata 2";
  588. return 1;
  589. }
  590. while ( ps1 != snrlist->last()  && ((Snr*)ps1->element())->per > tper  ) {
  591. ps1=(NSNode*)PHYmod[pm]->after(ps1);
  592. }
  593. //ps1 o e' minore di tper o e' l'ultimo
  595. //prendo quello precedente
  596. if (((Snr*)ps1->element())->per < tper)
  597. ps2 = (NSNode*)snrlist->before(ps1);
  600. //se le lunghezze erano due diverse cerco snr anche nella seconda
  601. if ((pl1 != pl2) & (pl2 != 0)){
  602. //cerr << "dicversi";
  603. snrlist = ((leng*)pl2->element())->snr_list;
  604. //cerco
  605. ps3 = (NSNode*)snrlist->first();
  607. if(ps3 == 0) {
  608. cerr<<"Tabella non inizializzata 3";
  609. return 1;
  610. }
  612. while (not (((Snr*)ps3->element())->per <= tper or ps3 == snrlist->last() ) ){
  613. ps3=(NSNode*)PHYmod[pm]->after(ps3);
  614. }
  615. //prendo anche quello precedente
  616. if (((Snr*)ps3->element())->per < tper)
  617. ps4 = (NSNode*)snrlist->before(ps3);
  619. }
  622. //settaggio delle variabili di uscita in base ai risultati della ricerca. 
  623. //I valori /1 indicano che no e
  624. if (ps1){
  625. //cerr << ps_;
  626. ((Snr*)ps1->element())->snr;
  628. ps_->snr[0] = ((Snr*)ps1->element())->snr;
  629. ps_->per[0] = ((Snr*)ps1->element())->per;
  630. ps_->len[0] = ((leng*)pl1->element())->len;
  631. }else{
  632. ps_->snr[0] = -123456789;
  633. ps_->per[0] = -1;
  634. ps_->len[0] = -1;
  635. }
  638. if (ps2){
  639. ps_->snr[1] = ((Snr*)ps2->element())->snr;
  640. ps_->per[1] = ((Snr*)ps2->element())->per;
  641. ps_->len[1] = ((leng*)pl1->element())->len;
  642. }else{
  643. ps_->snr[1] = -123456789;
  644. ps_->per[1] = -1;
  645. ps_->len[1] = -1;
  646. }
  648. if (ps3){
  649. ps_->snr[2] = ((Snr*)ps3->element())->snr;
  650. ps_->per[2] = ((Snr*)ps3->element())->per;
  651. ps_->len[2] = ((leng*)pl2->element())->len;
  652. }else{
  653. ps_->snr[2] = -123456789;
  654. ps_->per[2] = -1;
  655. ps_->len[2] = -1;
  656. }
  658. if (ps4){
  659. ps_->snr[3] = ((Snr*)ps4->element())->snr;
  660. ps_->per[3] = ((Snr*)ps4->element())->per;
  661. ps_->len[3] = ((leng*)pl2->element())->len;
  662. }else{
  663. ps_->snr[3] = -123456789;
  664. ps_->per[3] = -1;
  665. ps_->len[3] = -1;
  666. }
  668. }
  673. double PER::get_snr(PhyMode pm, double tper, double tlen){
  675. struct per_set ps_;
  677. get_snr( pm,  tper, tlen, &ps_);
  680. int len1 = -1;
  681. int len2 = -1;
  682. double snr1 = -123456789;
  683. double snr2 = -123456789;
  684. double snr = -123456789;
  687. //not found entries are set equal  to the clesest one
  689. //cerr <<  ps_.len[0] << ps_.len[1] << ps_.len[2] << ps_.len[3];
  691. if(ps_.len[0]==-1){
  692. ps_.len[0]=ps_.len[1];
  693. ps_.snr[0]=ps_.snr[1];
  694. ps_.per[0]=ps_.per[1];
  695. }
  697. if(ps_.len[1]==-1){
  698. ps_.len[1]=ps_.len[0];
  699. ps_.snr[1]=ps_.snr[0];
  700. ps_.per[1]=ps_.per[0];
  701. }
  703. if(ps_.len[2]==-1){
  704. ps_.len[2]=ps_.len[3];
  705. ps_.snr[2]=ps_.snr[3];
  706. ps_.per[2]=ps_.per[3];
  707. }
  709. if(ps_.len[3]==-1){
  710. ps_.len[3]=ps_.len[2];
  711. ps_.snr[3]=ps_.snr[2];
  712. ps_.per[3]=ps_.per[2];
  713. }
  716. // intepolation: 1) linear interpolation at the varying of the PER for each length
  717. // 2) linear interpolation among lenghts 
  719. if (ps_.snr[0] == ps_.snr[1])
  720. snr1 = ps_.snr[0];
  721. else
  722. snr1 = ps_.snr[0] +  (ps_.snr[1]- ps_.snr[0])/(ps_.per[1]- ps_.per[0])*(tper- ps_.per[0]);
  724. if (ps_.snr[2] == ps_.snr[3])
  725. snr2 = ps_.snr[2];
  726. else
  727. snr2 = ps_.snr[2] +  (ps_.snr[3]- ps_.snr[2])/(ps_.per[3]- ps_.per[2])*(tper- ps_.per[2]);
  729. if (snr1!=-123456789 && snr2!=-123456789){
  730. if (ps_.len[3] != ps_.len[1])
  731. snr = snr1 +  (snr2- snr1)/(ps_.len[3]- ps_.len[1])*(tlen-ps_.len[1]);
  732. else 
  733. snr = snr1;
  734. }
  735. else 
  736. if (snr1!=-123456789){
  737. snr = snr1;
  738. }else 
  739. if (snr2!=-123456789){
  740. snr = snr2;
  741. }else {cerr << "No SNR available " << snr << " " << pm << " " << tper <<" " << tlen << "/n";}
  742. return snr;
  743. }
  746. void PER::print_data(){
  748. for (int i = 0 ; i<= Mode54Mb; i++) {
  749. for (int l = 128 ; l<= 1024; l=l*2) {
  750. for (double per = 0 ; per<1; per=per+0.1) {
  752. printf("%d %d %f %f n",i,l,per,get_snr((PhyMode)i, per , l));
  754. }
  755. }
  756. }
  759. }