Definition at line 97 of file MakeGaussIIR1d.h. References Impala::Core::Array::complex::im, and Impala::Core::Array::complex::re. Referenced by MakeGaussIIR1d(). 00098 { 00099 00100 /* 00101 * Store order of recursive filters in pp[0] 00102 */ 00103 pp[0].re = (double) nn; pp[0].im = 0.0; 00104 00105 switch (nn) 00106 { 00107 case 5: 00108 switch (order) 00109 { 00110 case 2: 00111 pp[1].re = 0.70381; pp[1].im = 1.38271; 00112 pp[2].re = 0.70381; pp[2].im = -1.38271; 00113 pp[3].re = 1.42239; pp[3].im = 0.77978; 00114 pp[4].re = 1.42239; pp[4].im = -0.77978; 00115 pp[5].re = 1.69319; pp[5].im = 0.0; 00116 break; 00117 00118 case 1: 00119 pp[1].re = 0.70237; pp[1].im = 1.38717; 00120 pp[2].re = 0.70237; pp[2].im = -1.38717; 00121 pp[3].re = 1.43280; pp[3].im = 0.77903; 00122 pp[4].re = 1.43280; pp[4].im = -0.77903; 00123 pp[5].re = 1.70346; pp[5].im = 0.0; 00124 break; 00125 00126 case 0: 00127 default: 00128 pp[1].re = 0.85991; pp[1].im = 1.45235; 00129 pp[2].re = 0.85991; pp[2].im = -1.45235; 00130 pp[3].re = 1.60953; pp[3].im = 0.83009; 00131 pp[4].re = 1.60953; pp[4].im = -0.83009; 00132 pp[5].re = 1.87040; pp[5].im = 0.0; 00133 break; 00134 } 00135 break; 00136 00137 case 4: 00138 switch (order) 00139 { 00140 case 2: 00141 pp[1].re = 0.94576; pp[1].im = 1.21364; 00142 pp[2].re = 0.94576; pp[2].im = -1.21364; 00143 pp[3].re = 1.59892; pp[3].im = 0.42668; 00144 pp[4].re = 1.59892; pp[4].im = -0.42668; 00145 break; 00146 00147 case 1: 00148 pp[1].re = 1.04198; pp[1].im = 1.25046; 00149 pp[2].re = 1.04198; pp[2].im = -1.25046; 00150 pp[3].re = 1.69337; pp[3].im = 0.45006; 00151 pp[4].re = 1.69337; pp[4].im = -0.45006; 00152 break; 00153 00154 case 0: 00155 default: 00156 pp[1].re = 1.13231; pp[1].im = 1.28122; 00157 pp[2].re = 1.13231; pp[2].im = -1.28122; 00158 pp[3].re = 1.78532; pp[3].im = 0.46766; 00159 pp[4].re = 1.78532; pp[4].im = -0.46766; 00160 break; 00161 } 00162 break; 00163 00164 case 3: 00165 switch (order) 00166 { 00167 case 2: 00168 pp[1].re = 1.21969; pp[1].im = 0.91724; 00169 pp[2].re = 1.21969; pp[2].im = -0.91724; 00170 pp[3].re = 1.69485; pp[3].im = 0.0; 00171 break; 00172 00173 case 1: 00174 pp[1].re = 1.32094; pp[1].im = 0.97057; 00175 pp[2].re = 1.32094; pp[2].im = -0.97057; 00176 pp[3].re = 1.77635; pp[3].im = 0.0; 00177 break; 00178 00179 case 0: 00180 default: 00181 pp[1].re = 1.41650; pp[1].im = 1.00829; 00182 pp[2].re = 1.41650; pp[2].im = -1.00829; 00183 pp[3].re = 1.86131; pp[3].im = 0.0; 00184 break; 00185 } 00186 break; 00187 00188 case 2: 00189 pp[1].re = 1.69580; pp[1].im = 0.59955; 00190 pp[2].re = 1.69580; pp[2].im = -0.59955; 00191 break; 00192 00193 case 1: 00194 pp[1].re = 2.00000; pp[1].im = 0.00000; 00195 break; 00196 00197 default: 00198 break; 00199 } 00200 return; 00201 }
|