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

源码开发语言/平台
当前位置：首页 > 源码/资料 > 通讯/手机编程 > 3G开发 > 查看源码
bilin_transf.cpp
资源名称：simulation_wireless_c++.rar　[点击查看]
上传用户：jtjnyq9001
上传日期：2014-11-21
资源大小：3974k
文件大小：3k
源码类别：

3G开发

开发平台：

Visual C++

bilin_transf.cpp：源码内容
  1. //
  2. //  File = bilin_transf.cpp
  3. //                     
  5.  #include <stdlib.h> 
  6.  #include <iostream>
  7.  #include <fstream> 
  8.  #include "misdefs.h"
  9.  #include "typedefs.h" 
  10.  #include <complex>
  11.  #include "bilin_transf.h"
  12.  extern ofstream *DebugFile;
  13.  
  14. //====================================================
  15. //
  16. //----------------------------------------------------
  17. void BilinearTransf( DenormalizedPrototype* filt_proto,
  18.                      double big_t,
  19.                      double **a_ret,
  20.                      double **b_ret)
  22. {
  23.   int max_poles;
  24.   int j, m, n;
  25.   int max_coeff;
  26.   double little_h;
  27.   complex<double> *mu;
  28.   complex<double> alpha;
  29.   complex<double> beta;
  30.   complex<double> gamma;
  31.   complex<double> delta;
  32.   complex<double> eta;
  33.   complex<double> work;
  34.   complex<double> c2;
  35.   complex<double> *pole;
  36.   complex<double> *zero;
  38.   double *a, *b;
  39.   int num_poles, num_zeros;
  40.   double h0;
  41.   double adjust;
  43.   pole = filt_proto->GetPoles();
  44.   num_poles = filt_proto->GetNumPoles();
  45.   zero = filt_proto->GetZeros();
  46.   num_zeros = filt_proto->GetNumZeros();
  47.   h0 = filt_proto->GetHZero();
  49.   if(num_poles > num_zeros){
  50.     max_poles = num_poles;
  51.     }
  52.   else{
  53.     max_poles = num_zeros;
  54.     }
  55.   //-------------------------------------------
  56.   a = new double[max_poles+1];
  57.   *a_ret = a;
  58.   b = new double[max_poles+1];
  59.   *b_ret = b;
  60.   
  61.   mu = new std::complex<double>[max_poles+1];
  62.   
  63.   for(j=0; j<=max_poles; j++){
  64.     mu[j] = 0.0;
  65.     a[j] = 0.0;
  66.     b[j] = 0.0;
  67.     }
  68.   //------------------------------------------
  69.   little_h = 1.0;
  70.   work = 1.0;
  71.   c2 = 2.0;
  72.   adjust = 1.0;
  74.   for(n=0; n<num_poles; n++){
  75.     work = work * (c2 - (big_t * pole[n]));
  76.     little_h = little_h * big_t;
  77.     adjust *= 4.0;
  78.     }
  79.   adjust *= num_poles;
  80.   little_h = h0 * little_h / std::real(work);
  82.   //-------------------------------------------------
  83.   mu[0] = 1.0;
  85.   for( n=1; n<=max_poles; n++){
  86.     mu[n] = 0.0;
  87.     }
  88.   for( n=1; n<=num_poles; n++){
  89.     gamma = complex<double>( 2.0/big_t, 0.0) - pole[n-1];
  90.     delta = complex<double>( -2.0/big_t, 0.0) - pole[n-1];
  91.     eta = delta/gamma;
  93.     for( j=n; j>=1; j--){
  94.       mu[j] = gamma * mu[j] + (delta * mu[j-1]);
  95.       #ifdef DEBUG
  96.         *DebugFile << "mu[" << j << "] = " << mu[j] << endl;
  97.       #endif
  98.       }
  99.     mu[0] = gamma * mu[0];
  100.     #ifdef DEBUG
  101.       *DebugFile << "mu[0] = " << mu[0] << endl;
  102.     #endif
  103.     }
  104.   for( j=1; j<=num_poles; j++){
  105.     a[j] = -little_h * mu[j].real() / h0;
  106.     }
  107.   //--------------------------------------------
  108.   mu[0] = 1.0;
  109.   for( n=1; n<=max_poles; n++){
  110.     mu[n] = 0.0;
  111.     }
  112.   max_coeff = 0;
  113.   //--------------------------------------------
  114.   for(m=1; m<=(num_poles-num_zeros); m++){
  115.     max_coeff++;
  116.     for( j=max_coeff; j>=1; j--){
  117.       mu[j] = mu[j] + mu[j-1];
  118.       }
  119.     }
  120.   for( m=1; m<=num_zeros; m++){
  121.     max_coeff++;
  122.     alpha = complex<double>( 2.0/big_t, 0.0) - zero[m-1];
  123.     beta = complex<double>( -2.0/big_t, 0.0) - zero[m-1];
  125.     for( j=max_coeff; j>=1; j--){
  126.       mu[j] = alpha * mu[j] + (beta * mu[j-1]);
  127.       }
  128.     mu[0] = alpha * mu[0];
  129.     }
  130.   for(j=0; j<=max_coeff; j++){
  131.     b[j] = little_h * mu[j].real();
  132.     }
  133.   return;
  134. }