新数字信号处理 第2版 教学课件 张小虹 8数字信号处理4.pptVIP

式中Re[Xm+1(i)]为无限精度实部运算结果，er为有限精度 = Re[Xm+1(i)]+er (8.5-26) 实部运算结果的误差。 由图8.5-7得到浮点实部运算的输出为 ^ Re[Xm+1(i)] =Re[WN]Re[Xm (j)] (1+?1) (1+?1)(1+?2) r ? Im[WN ]Im[Xm (j)](1+?2)(1+?1)(1+?2)+Re[Xm(i)](1+?2) r 从表达式中略去?、? 的所有高次项，可得： 为分析简便，作如下假设 (1)、所有误差源与信号x(n)不相关、独立；都是彼此独 立、白色、等概分布，具有相同方差，因此 er ? (?1 +?1 +?2) Re[WN]Re[Xm (j)] r ? (?2 +?1 +?2) Im[WN ]Im[Xm (j)] +?2 Re[Xm(i)] r E[|?1(n)|2] =E[|?2(n)|2]=E[|?1(n)|2]=E[|?2(n)|2] =?B 2 2?2b 3 = 1 (2)、x(n)是白色；实、虚部方差相同；且同一级各节点 方差相同，即 E[|Re[Xm(j)]|2]= E[|Im[Xm(j)]|2] E[|[Xm(i)]|2]= E[|[Xm(j)]|2] 2 1 = E[|[Xm(q)]|2] 由此得到实部误差的方差为 (8.5-27) = 2?B E[|Xm (j)|2] 2 E[er ] = 3?B [Re(WN)]2 E[|Xm (j)|2] r 2 1 2 ? 2 + 3?B [Im(WN)]2 E[|Xm (j)|2] r 1 2 ? +?B E[|Xm (i)|2] 1 2 ? 2 2 则输出节点噪声em+1的方差为 上式表明浮点制的输出节点噪声与其输入节点变量 相关。 (8.5-28) 同理可得到虚部误差的方差为 E[ei ] = 2?B E[|Xm (i)|2] 2 2 ?m+1 =E[em+1 ]=E[er ] +E[ei ] 2 ? 2 2 2 (8.5-29) = 4?B E[|Xm (i)|2] 2 * 在线教务辅导网： 教材其余课件及动画素材请查阅在线教务辅导网 QQ:349134187 或者直接输入下面地址： 即输入x(n)，输出为X(k)，分别讨论DFT与FFT运算及系 §8.5 DFT与FFT的有限字长影响 DFT与FFT在数字滤波器、频谱分析中应用广泛，了解 DFT有限字长影响很重要，不过要精确分析是很困难的， 只对选择所需要寄存器长度的简化分析就足以。分析时 为方便要做许多假设，也采用输入、输出方式进行分析， 数量化误差的影响。 8.5.1、DFT直接定点舍入计算误差分析 DFT定义为 k=0,1,2,…,N?1 (8.5-1) 的作用相当于单位脉冲响应， 的运算流图如图所示。 无限精度下第k个X(k)值 X(k) X(k)= kn n=0 ? N?1 x(n)WN x(0) WN 0 x(1) WN 1 x(2) WN 2 x(N?1) WN N?1 kn 若将上式中的X(k)当作输出，x(n)当作输入，而WN … 的分析类似，将每次舍入误差等效为一个噪声源，此时 第k个X(k)值的等效 统计模型如图所示。 ε(0,k) ε(1,k) ε(2,k) ε(N?1,k) kn 定点舍入运算x(n)WN 每次相乘都会产生误差，与前面 kn 不考虑WN 的量化误差，仅讨论运算引起的误差。因为 X(k) =X(k)+F(k) ^ x(0) WN 0 x(1) WN 1 x(2) WN 2 x(N?1) WN N?1 … 上图 X(k) 是无限精度下第 k个 X (k) 值的运算结果， k个X (k) 值运算结果的误差。由图8.5-2可见，各误差 直接加至输出端，因此总的误差输出为 有限精度运算下的输出为 (8.5-2) (8.5-3) X(k) =X(k)+F(k) ^ X (k)是有限精度下第k个X(k)值的运算结果，F(k)是第 ^ m=0 N?1 F(k) =? ? (n,m) 法是由四次实数乘法实现的，因此一次复数乘法会产 生四次误差，即 kn 通常x(n)、WN 均为复数，计算x(n)WN 的一次复数乘 kn kn Q[x(n)WN ] =Q {Re[x(n)]+jIm[x(n)]} N kn 2? N kn 2? cos ?jsin =Re[x(n)] cos +?1(n,k) +Im[x(n)]sin N kn 2? N kn 2? +?2(n,k) +jIm[x(n)] cos +?3(n,k) ?jRe[x(n)]sin N kn 2? N kn 2? +?4(n,k)