08数字信号处理4.pptVIP

 Xm(j) Xm(i) ? (m,j) 8.5.3、FFT浮点舍入计算误差分析 与定点情况相同，对不同的FFT算法，相应的有限字长 效应不同。仍以基2、时选N=2M 为例，从一个蝶形运算 产生的误差开始讨论。 图中符号意义与定点运算统计模型相同。 单个蝶形浮点运算统计模型如图所示， ?1 蝶形上端输出为(下端输出的讨论相同) (8.5-24) 其中em+1具有实部和虚部，即 em+1 =er+j ei er、ei分别由Xm+1(i)的实部和虚部运算所决定。 浮点运算不论加法、乘法都产生误差，而在这个运算中 有两次乘法、两次加法，均引入舍入误差，均用相对误 差表示，即以?1、?2表示乘积引入的相对误差，以?1、 ?2表示加法引入的相对误差。 (8.5-25) 先分析实部误差er (虚部误差ei可以类推)，因为 Re[Xm (i)] 各噪声源与运算的关系如图所示。 1+?2 1+?1 1+?1 1+?2 由图8.5-7得到浮点实部运算的输出为 式中Re[Xm+1(i)]为无限精度实部运算结果，er为有限精度 = Re[Xm+1(i)]+er (8.5-26) 实部运算结果的误差。 从表达式中略去?、?的所有高次项，可得： 为分析简便，作如下假设 (1)、所有误差源与信号x(n)不相关、独立；都是彼此独 立、白色、等概分布，具有相同方差，因此 (2)、x(n)是白色；实、虚部