数字信号处理第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 与定点分析情况相同，因为前一级误差通过后一级蝶 形时其方差保持不变，所以浮点FFT总的输出误差与 从输入x(n)到输出X(k)经过的蝶形个数相关，与第k个 F(k)有关的蝶形： M(末级)级一个；M?1 (末前级)级 两个；M?2级四个；……，即在输出端总的输出误差 方差为 (8.5-30) E[|F (k)|2] 2 = 4?B E[|XM?1 (j)|2]+2E[|X M?2 (j)|2]+??? ??? + 2M?1E[|X (j)|2] 因为每级蝶形输出信号的方差是输入信号方差的2倍， 即 (8.5-31) 求式(8.5-31)的方差，可得 E[?Xm+1(i)?2]= E[?Xm (i)?2]+ E[?Xm(j)?2] (8.5-32) =2 E[?Xm (i)?2] =2 E[?Xm (j)?2] Xm+1(i)= Xm(i)+ WN Xm(j) r 利用这一关系类推，将所有Xm(j)的方差用X(k)的方差 表示为 E[?[XM?1(j)]?2]=0.5 E[?[X(k)]?2] ； E[?[XM?2(j)]?2]=(0.5)2 E[?[X(k)]?2] ； E[?[x(j)]?2]=(0.5)M E[?[X(k)]?2] ； … 将上面的关系代入(8.5-30)式，得到最后总的输出误差 方差为 输出噪信比 E[|F(k)|2] 2 = 4?B E[|X (k)|2] 2 1 1 22 2M 1 +2 +???+ 2M?1 = 4?B E[|X (k)|2] 2 2 M = 2M?B E[|X (k)|2] 2 (8.5-33) E[|F(k)|2] E[|X(k)|2] = 2M?B 2 M 2?2b 3 1 = (8.5-34) 变化，这也是所有浮点制运算的共同特点。 由上式可见，输出噪信比正比M，远远小于定点制的 N=2M。这一结果表明相同尾数字长情况下，浮点噪信 比比定点小(运算精度高)。当然，这是由浮点制数字表 达的复杂性换来的；另外浮点噪信比不随信号幅度大小 8.5.4、FFT的系数量化效应 无限精度表示的DFT为 化误差等效为随机噪声，用统计的方法分析系数量化 了解FFT系数量化效应影响很有必要，但要精确分析也 不易。所以分析时为方便要做许多假设，也采用输入、 输出方式分析，即输入为 x(n)，输出为X(k)。将系数量 效应。 X(k)= kn n=0 ? N?1 x(n)WN 均是WN的各次幂系数，其中一些可能是不需相乘的±1。 系数量化后，上式可表示为 式中F(k)是由系数量化引起的DFT计算误差由某个 x(n)