第5章离散时间系统的相位结构与逆系统

1、第5章 离散时间系统的相位、 结构与逆系统5.1 离散时间系统的相频响应；5.2 FIR 系统的线性相位；5.3 具有线性相位系统的零点分布；5.4 全通系统和最小相位系统；5.5 谱分解；5.6 FIR 系统的结构；5.7 离散时间系统的 Lattice 结构；5.8 逆系统（状态变量）5.1 离散时间系统的相频响应()()()jjjH eH ee 幅频响应相频响应() :( ):jH e 如果：( )k 我们称其为线性相位若( )k 也称为线性相位()()()()()xjjjjkjY eH eX eX ee()1jH e对输入 ，有( )x n假定：所以( )()y nx nk输出是输入的

2、简单移位，移位的大小正比于 因此不会发生失真。k( )k p.97 时移性例：令100( )cos()cos(2()y nnknk则没有发生相位失真00( )cos()cos(2)x nnn()jjkH ee具有线性相位例：令若200( )cos(/4)cos(2)y nnn则发生了相位失真00( )cos()cos(2)x nnn()()jjH ee 00/403/2( )3/2 -1001020-202-1001020-202-1001020-202( )x n1( )y n2( )y n 如果系统的相频响应不是线性的，而输入信号由多个正弦信号组成，那么系统的输出将不再是输入信号作线性移位

3、后的组合，因此，输出将发生相位失真。 ( )( )gdd 定义为系统的群延迟(Group Delay, GD)显然: 若系统具有线性相位，则其群延迟为常数00( )( )cos(),( ):Narrowband Signalacx nx nnx n若则0000( )()()cos()jagpy nH ex nn 即相位延迟 反映了载波信号的延迟， 而群延迟 反映了输出包络的延迟。 0()p0()g思考：如何实现对信号的零相位滤波？若 要保证系统是因果的，又如何实现？( )( )p )(nh5.2 FIR 系统的线性相位 在绝大部分信号处理的场合，人们都期盼系统具有线性相位，但是，如何实现线性相

4、位？对 FIR 系统，如果保证：( )(1)0,1,1h nh NnnN 则该系统具有线性相位。( )(1)h nh Nn :even:oddNN上述对称有四种情况：第一类 FIR 系统偶对称( )h n( )(1)h nh Nn :even:oddNN奇对称( )h n第二类 FIR 系统1. 为奇数N10()( )Njj nnH eh n e1(3) 2120(1)/21( )( )()2NNNjj nj nnnNNh n eh n ehe1mNn 令并利用 的对称性，有()jH e(3) 2(3) 2(1)00()( )(1)NNjj nj NmnmH eh n eh Nm e 1()2

5、1()2NjNhe第一类 FIR 系统(3) 2(1)/20112( )cos()()22Nj NmNNeh mmh()jH e(1) 2nNm令1()02( )12 ()1,2,(1) 22Nhna nNhnnN令(1) 2(1)/21112()cos()()22Nj NnNNehnnh()jH e实数(1) 2(1)/20()( )cos()Njj NnH eea nn最后有(1) 20( )(1)2( )( )cos()NgnNHa nn 相位增益 所以，只要保证滤波器的系数偶对称，该滤波器必然具有线性相位。 也呈偶对称，对呈偶对称，对20)(20)cos(gHn2(1)/211()2

6、()cos() 22Njj NnNH eehnn2. 为偶数N( )2 (),1,2,22NNb nhnn令2(1)/211()( )cos() 2Njj NnH eeb nn则呈偶对称呈偶对称或，对呈奇对称对奇对称对，）（为偶数时偶对称，)()2cos(0)21(cos,20)()21(cos)2(0) 1(, 0)(, 0)21(cos(1)(ggggHnnHnzHHnNnh( )(1)h nh Nn 第二类 FIR 系统3. 为奇数N1(1) 2221()( )sin()NNjjnH eec nn1( )2 ()1,2,(1) 22Nc nhnnN( )(1)/2/2N 呈奇对称，对处呈

7、奇对称，在，在）（为奇数时奇对称，20)(20)sin()2(0) 1(, 0)(, 0)sin(20)sin(1)(gggHnzHHnnNnh4. 为偶数N122211()( )sin ()2NNjjnH eed nn( )2 ()1,2,22Nd nhnnN( )(1)/2/2N 呈偶对称对呈奇对称，对处呈偶对称在处呈奇对称，在，在）（为偶数时奇对称，,)(20)(20)21sin()2(, 0) 1(, 0)(, 0)21sin(20)21sin(1)(ggggHHnzHHnnNnh的线性组合，在 时， 易取得最大值，因此这一类滤波器易体现低通特性，且是偶函数。通过频率移位，又可体现高通

8、、带通、带阻特性。所以，经典的低通、高通、带通和带阻滤波器的 都是偶对称的。说明： 第一类 FIR 系统是cos()n0()jH e( )h n的线性组合，在 时， 的值为零，且是奇函数。这一类滤波器都是作为特殊形式的滤波器，如 Hilbert变换器、差分器等。第二类 FIR 系统是sin()n0()jH e 最好取为奇数，以便以中心点为对称。N思考：四类滤波器的对称点在何处请掌握四种情况下线性相位表达式的推导方法所以， 的零点也是 的零点, 反之亦然( )H z1()H z5.3 具有线性相位的FIR系统的零点分布1100( )( )(1)NNnnnnH zh n zh Nn z 1mNn

9、令1(1)0( )NNnnzh n z 则1(1)0( )( )NNmmH zh m z (1)1()NzH z ( )H z 的零点: 零点分布可能有四种情况(单位圆上及其内）：1. 不在实轴也不在圆上，应是一对共轭零点，模1;2. 不在实轴，但在圆上，也是一对共轭零点；模1；3. 在实轴但不在圆上，无共轭，角度0， ，模1;4. 在实轴，且在圆上，无共轭，角度0， ，模1;( )H zkjkkerz .0, ,1kkar1111( )(1)(1)11(1)(1)kkkkjjkkkjjkkHzz r ez r ezezerr四个零点同时存在, 构成四阶系统.kz在单位圆内把该式展开，其系数也

10、是对称的，是具有线性相位的子系统。（不在实轴也不在圆上，应是一对共轭零点，模1）有镜象零点.0, ,1,.kkkbrz在实轴上)11)(1 ()(11kkmrzrzzH无共轭零点, 有镜象零点)1)(1 ()(11kkjjlezezzH.0, ,1,.kkkcrz在单位圆上无镜象对称零点, 有共轭零点.(在实轴但不在圆上，无共轭，角度0， ，模1)( 不在实轴，但在圆上，是一对共轭零点；模1 )kmlnnlmkzHzHzHzHzH)()()()()(一个具有线性相位的FIR数字滤波器的转移函数可表示为上述四类 FIR 子系统的级联，即：很容易证明，每一个子系统的系数都是对称的，因此它们都具有线

11、性相位。.0, ,11,1kkkdrz)1 ()(1zzHn无镜象零点, 也无共轭零点. ( 在实轴，且在圆上，无共轭，角度0， ，模1 )5.4 全通系统和最小相位系统 如果一个系统的幅频响应对所有的频率都等于1 (或一个常数), 即1| )(|japeH|0则称系统 为全通系统。)(zHap( )kapHzz最简单的全通系统，纯延迟全通系统1111( )11apzHzz一阶全通系统：Pole:,Zero:1zz镜像对称212( )( )()apapapHzHz Hz22()()()jjjapapapHeHeHe二阶全通系统： 111111(1)(1 ()( )(1)(1)apzzHzzz一

12、对位于单位圆内的共轭极点，一对共轭零点和极点以单位圆为镜像对称(共轭并求倒数）111( )11NkapkkkzHzz 高阶全通系统：1高阶全通系统的另一种表示形式1(1)111212( )1NNNNapNNaaza zzHza za za z 1()( )( )NapzA zHzA z 即2*()()()1jjjapapapHeHeHe对该全通系统，请自己证明：式中系数a1, a2,aN 均为实数 1 . 是IIR系统(不考虑纯延迟形式）； 2. 极点数和零点数相等； 3. 极点和零点是以单位圆镜像对称的； 4. 极点都在单位圆内,零点都在单位圆外； 5. 每一对极零点都是镜像对称的，且零点在

13、单位圆外， 由零变到 时，相频响应 是单调递减的(全通系统的群延迟始终为正值)全通系统的特点：)(IIR系统的 无限长，无法对称，即无法作到线性相位。在实际中，可以用一个全通系统和IIR系统相级联，在不改变幅频响应的情况下对相频响应做矫正，使其接近线性相位。( )h n全通系统的应用：全通系统还广泛应用在系统分析及一些特殊滤波器的设计方面-101-1-0.500.51(a)00.20.400.511.5(b)00.20.4-4-3-2-10(c)05101520-0.500.51(d) 一阶全通系统极零图幅频相频抽样响应)cos()(2)sin()(arctan)(arg11jeHrad)(2

14、/-101-1-0.500.51(a)00.20.400.51(b)00.20.4-8-6-4-20(c)0102030-0.500.51(d) 三阶全通系统相频响应rad)(2/极零图幅频 一个离散系统，其极点必须在单位圆内，但对零点没有限制，如果：1. 所有的零点都在单位圆内： 最小相位系统；2. 所有的零点都在单位圆外： 最大相位系统；3. 单位圆内、外都有零点 ： 混合相位系统。最小相位系统1. 在具有相同幅频响应的因果的稳定的滤波器 集合中, 最小相位滤波器具有最小的相位偏移最小相位系统的性质：例：作为练习，请证明如下两个系统具有相同 的幅频响应：1111( ),1,11zbbzH

15、zabzaaz1211( ),1,11bzbzHzabzaaz那一个是最小相位系统1111( ),1,11zbbzH zabzaaz1211( ),1,11bzbzHzabzaaz00.20.4024(a)00.20.4-1.501.53(b)0510152000.20.40.60.81(c)05101520-1-0.500.51(d)幅频相频21)(1nh)(2nh单位抽样响应更集中在n=0附近2. 在所有具有相同幅频响应的离散系统中, 最 小相位系统的 具有最小的延迟；( )h n令：20()( )0MnE Mh nM 累积能量有：22min00( )( )MMnnhnh n所以，最小相位

16、系统的单位抽样响应又称最小延迟序列。思考： 具有线性相位的FIR系统是否是最小相位系统？/31 2/31 212/31/31/31/3122/31 2/31 232(1 0.5) (1 0.5)( )1 0.81(1 0.5)(1 0.5)(0.5)(0.5)( )1 0.81(0.5) (0.5)( )1 0.81jjjjjjjjezezH zzezezezezHzzezezHzz例. 三个系统它们具有相同的幅频响应，试判断，那一个是最小相位系统？最大相位系统？混合相位系统？请注意：为保证系统具有相同的幅频响应（相同的定标）， 的表达式。123( ),( ),( )H zHzHz-101-1

17、01222-202-1012-202-10122200.250.50102000.250.5-20-1001001020-20201020-20201020-20201020051015121E3E)(1nh)(2nh幅频3.对于稳定因果系统，当且仅当其是最小相位 系统时, 该系统才有逆系统 （Inverse System) 令( )( )( )N zH zD z记INV1( )( )( )( )D zHzH zN z 的逆系统( )H z( )H zINV( )Hz( )y n( )x n( )x nDeconvolution(反卷积）System identification(系统辨识）1

18、)(),()()()()()( (z) 729. 011)(/1)( 729. 01)( )()(3 , 0, 0)(;729. 0)3(, 1)0( 33zYnnxzXzXzHzHYzzHzHzzHzYzHnnhhhINVINVINV则若则所以因为的输出求其逆滤波器例如4. 任一非最小相位的因果系统的转移函数均可由一个最小相位系统和一个全通系统 级联而成, 即：)()()(minzHzHzHap由于最小相位系统有着以上特殊的性质，因此有着广泛的应用，特别是在信号的建模与系统辨识方面。要理解，具有相同幅频响应的系统，它们所对应的转移函数可以是不相同的，区别就在于相位（或零点的位置）。那么，如何

19、由一个最小相位系统得到具有相同幅频响应的最大相位、混合相位系统？设 有一个零点在单位圆之外,即 其余的极零点均在单位圆内,则 可表示为, 是最小相位的, 又可表示为 )(zH, 1,/100zzz)(zH)()(011zzzHzH（)(1zH)(zH 例题：1*001min1*0011*011*01*00111)(1)1)(11)()(zzzzzHzzzzzzzHzzzzzzzHzH圆内的零点位于单位全通系统最小相位具有相同的幅频响应和)()(minzHzH5.5 谱分解（Spectral factorization)1,0 ,)1)(1 ()(1,0 ,)1)(1 ()(1,0 ,)1)(1

20、 ()(1,0 ,)1)(1 ()(12min11max11minbabzazzHbabzazzHbabzazzHbabzazzHmix系统还可构建两个混合相位统可构建一个最大相位系统若已知一个最小相位系具有相同的幅频特性)()()1)(1)(1)(1 ()()(max1max111minminzHzHbzazbzazzHzH显然：幅频特性一样1( )( )()P zH z H z若显然， 具有线性相位 或零相位( )P z零点以单位圆镜像对称将一个转移函数的极零点重新分配，得到两个转移函数， 这一过程（或方法）就称为“谱分解”。最常用的是将具有线性相位系统的转移函数作分解，并且往往是分解成两

21、个具有相同幅频响应的子系统。谱分解定义：1( )( )()P zH z H z)()()(10zHzHzP或( )p n=1.0000 ，4.0500，8.1000 ，14.9956，27.7248，43.2996，51.1831，43.2996，27.7248，14.9956，8.1000，4.0500，1.0000例. 令显然，该系统具有线性相位，共有12个零点：0.8,1/0.8,0.6,1/0.6,2/32/3/3/3,0.6,/0.6jjjjeeee 下图是对 作谱分解的结果，可以看出，分解后的两个系统具有相同的幅频响应。( )P z 谱分解的目的是想得到因果的、符合某种要求的系统，

22、这在信号建模、特殊滤波器的设计中经常要用到。分解的一般方法是： 令一个系统是最小相位系统； 则另一个系统必然是最大相位系统。这样，两个系统都有着相同的幅频响应。 另外一种分解方法是得到两个混合系统，目的是保证它们都具有线性相位。0510020406000.250.500.51-2-101-101221200.250.500.510510012-101-101600.250.500.510510051015-2-101-1016( )P z( )P z( )H z1()H z)(np)(0nh)(1nh幅频5.6 FIR 系统的结构直接实现:12012( )( )MMY zX z bb zb z

23、b z0( )MnnnH zb z一、 直接实现和级联实现级联实现:LkkkkzzzH122110)()( )(1)h nh Nn :oddM:evenM乘法量减少一半二、 具有线性相位的FIR系统的结构1 NM)1()(nNhnh1 NM )()0( ),1()0( )25()25( , )22()32( )23()23( ),12()22( )21()21( ),2() 12( MhhNhhMhMhNhNhMhMhNhNhMhMhNhNhN即即即即为偶数时 )()( )()()( )(2/ )1(00nMMnnMnnzznhzXznhzXzY所以)( )()1()( )()()()()1(

24、12/0)1(12/012/012/12/010nNnNnnNNnnNnnNNnnNnNnnzznhznNhznhznhznhznhzH:oddM)()()()()()()(2/ )1(0)(2/ )1(0nMMnnnMMnnzznhnxnyzznhzXzY:evenM)1()(nNhnh乘法量减少一半 所以: )2()( )( )()( )2()(12010MnMnMnNnnzMhzznhznhzHN为奇数时同理，当1/0( )Nh n 其它1, 1 , 0 Nn11011 1( )1NNnnzH zzNNzFIR 系统该系统实际上是一个N点平均器。11001( )() ( )()NNkky

25、 nx nk h kx nkNIIR系统三、 FIR系统的递归实现及梳状滤波器(p.289)该系统可由一FIR系统和一个一阶IIR系统级联而成，极零点抵消后，仍是一FIR系统。12111( ),( )1NzH zHzNz令令IIR 实现11 1( )1NzH zNz1,.,1 , 0 z 11)( )(2k11NkeNzzNzzHzHkNjNNN令其分子为零而得的零点由211101( )(1),NjkNkH zezN/21sin(2)()2jj NNH ej eN梳状滤波器频率响应 1)(1NeeHjNj1, 1 , 0,2,2NkNkkN得令211101( )(1),NjkNkH zezN/

26、21sin(2)()2jj NNH ej eN梳状滤波器N点平均器)()()(21jjjeHeHeH8N系统是一个一阶的IIRzzH1211)()2/sin()2/sin(1)()()(1)(, 12/)1(211NeNeHeHeHzzHzNjjjj保证了整个系统的稳定的零点相抵消，级联与其与不稳定极点在Sinc函数1100( )( ),( )( )NNnknNnnH zh n zH kh n W211001101( )( )11( )1NNNnknNnkNNjkkH zH k WzNzH kNez10121)(1NkkjNzekHNzN思路：用DFT系数 表示系统函数( )H k( )H z

27、四. 频率抽样实现 (p.202)令令11( )NzH zN1, 221)()(zekHzHkjkN梳状滤波器N个一阶IIR系统)()()(10, 21zHzHzHNkk则则21101( )( )1NNNjkkzH kH zNez可按上述级联方式得到系统的信号流图：该结构一方面反映了 Z 变换、DTFT、DFT之间的关系，另一方面，给出了FIR 滤波器设计的一种有效方法。取样型结构的特点。此为频率处的响应即为在频率点个零极点相抵消，好与梳状滤波器级联后正的无耗谐振器。是一个谐振频率为的响应是对。极点在单位圆上有一个每一个说明),(N2)(N2,N2)(z 1)()( :10, 210, 221

28、, 22kHkNzHkkzHezekHzHNkkNkkkNjkkjkN1011011)()1 ()( )()()()( )(, 1,)(1)()1 ()( ,11 1. :NkkNNNkNkNkNrrNkkNrNNzrWkHNzrzHkHWHrWHrWzzHkHrkHzrWkHNzrzHrFIRr故所以因是修正点上的取样值此时圆上相应移到即梳状滤波器的零点也系统而子的圆周上又接近设置在半径器的极点实践中只好将所有谐振，不能用；为零移动，系统稳定裕度当系数量化时，极点会问题讨论频率取样型结构的两个03 4.5 TN911.72 p.1741997, , )(4 . 6 . 5)( . 2北京理工

29、大学出版社，王世一处理信号数字可得到局部改善。见时序列为实复数运算。若系统的结构实现需进行大量的图都是复数，因此按书上nhkHWkN5.7 离散时间系统的 Lattice 结构Lattice 结构在基于模型的功率谱估计、语音信号处理、自适应滤波方面有着重要的应用。对下式所示的FIR系统，其Lattice 结构如下图：( )01( )( )( )1MMiiiMiiH zB zb i zbz )(1npm)(1nqm1zmkmk)(npm)(nqm) 1()()(11nqknpnpmmmm11( )( )(1)mmmmqnk pnqn mk反射系数Lattice 结构的基本单元1. 全零点系统(F

30、IR)的Lattice结构111111( )( )( )( )( )( )mmmmmmmmP zPzk z QzQzk Pzz Qz 如何实现滤波器系数和 的相互转换mk)()(1)()(1111zQzPzkzkzQzPmmmmmm( )010( )( ) /( )1( )( ) /( )miimmmimmBzPzP zbzBzQzQz定义：MmMm, 2 , 1, 2 , 1 00( )( )P zQzLattice结构中的基本关系( ),( )mmBzBz：是Lattice 结构中第 m 个上、下结点相对输入端的转移函数。( )01( )( )( )1MMiiiMiiH zB zb i z

31、bz 11111( )( )( )( )mmmmmmk zBzBzBzBzkz 1211( )( )/(1)( )( )mmmmmmmkBzBzkBzBzzkz得到由低阶倒高阶，或由高到低的递推关系。由高阶至低阶 或由低阶至高 阶转移函数的递推关系)1/()()()()()()(211111mmmmmmmmmmmkzBzkzBzBzBzkzBzB( )01( )( )( )1MMiiiMiiH zB zb i zbz 将下式关于 的定义代入分别上述两式)(),(1zBzBmm()( )( )()11()( )( )()21/(1)mmmiim immmmmmmiim immmmmbkbbk b

32、kbbbk bk 利用待定系数法得到时域递推关系：低到高阶高到低阶MATLAB中有相应的 m 文件。( )01( )( )( )1MMiiiMiiH zB zb i zbz 例：123( )( )1 1.71.50.648H zB zzzz (1)(2)(3)3331.7,1.5,0.648bbb (3)330.648kb (1)(1)(2)2233 33(2)(2)(1)2233 33/(1)1.221453/(1)0.738498bbk bkbbk bk 122( )1 1.2214530.738498B zzz (2)220.738498kb (1)(1)(1)212222/(1)0.7

33、0259bbk bk 11( )1 0.70259B zz (1)110.70259kb 111( )( )1MkkkH zA za z看作是FIR系统的逆形式。2. 全极点系统(IIR)的Lattice结构11( )( )(1)mmmmpnpnk qn11( )( )(1)mmmmqnk pnqn mkmk)(npm)(nqm1z)(1npm)(1nqm( )1( )1,( )( )( )( )mmmmY zY zP zAzQzAz( )1( )11( )( )( )1MiiMMMiY zH zPzAza z的求解方式同FIR系统Lattice结构的计算方法, 只是将多项式的系数 换成 .

34、Mkkk,21 Mmmiaim, 2 , 1, 2 , 1,)( )(ima系数系数及及( ) imb注意：在递推求解的过程中，反射系数1,1,2,mkkM有关反射系数的更多讨论见第12章信号建模。01( )( )( )1NkkrNkkkb zB zH zA za z3. 极零系统的Lattice结构上半部对应全极系统上半部对应全极系统下半部对应全零系统下半部对应全零系统12,Nk kkNk, 1 , 0 两组Lattice系数1201,NNk kkc cc1/( )( )A zB z求出同全极系统；()1Nm kkkmmm kcbc a 递推求解5.8 离散系统的状态变量描述1. LSI系统

35、的状态变量与状态方程)()()(01rnxbknyanyMrrNkk差分方程：01( )( )( )1NrrrNkkkb zB zH zA za z01( )( )1MrrrNkkkb zY zX za z11( )( )1NkkkV zX za z0( )( )MrrrY zV zb z1( )()( )Nkkv na v nkx n 0( )()Mkky nb v nk转移函数、差分方程、中间变量的关系1. “状态”指系统内一组变量, 它包含了系统全部 过去的信息, 由这一组变量和现在与将来的 输入，可求出现系统现在和将来的全部输出；2. 可用于分析多输入、多输出系统；如何选择状态变量？有

36、着不同的方法。方法之一是选择( ), (1), ()v n v nv nN作为系统的状态。12( )(1)( )(2)( )()Nw nv nw nv nwnv nN定义一组新的变量相互关系21321(1)( )(1)( )(1)( )NNw nw nw nw nwnwn111(1)( )( )( )( )NNw nv na w na wnx n 121112233(1)( )1(1)( )10000(1)0100( )0( )00100(1)( )NNNNaaaaw nw nw nw nw nw nx nwnwn 状态方程01011( )( )(1)()( )( )( )MMMy nb v

37、nbv nb v nMb v nbw nb wn0( )( )MrrrY zV zb z121,230( )( )( ),( )( )( )NNw nw ny nc ccw nb x nwn1101220201010,MMMMMNNcbb acbb acbb acb acb a 输出方程(1)( )( )( )( )( )nnnnnnwAwBxyCwDx上述内容讨论了如何由差分方程转换为状态方程。当然，反过来也可以。)(nxBD) 1( nw1z)(nwC)(nyA(1)( )( )( )( )( )nnnnnx nwAwBxyCwD两边取两边取Z变换：变换：( )( )( )( )( )(

38、)zzzzzzzWAWBXYCWDX)()()()(1zBXAzIzWzBXzWAzI)()()(1zDXzBXAzICzYDBAzICzXzYzH1)(/ )()(2.由状态方程求系统的转移函数)()()()()() 1(nDxnCwnynBxnAwnw) 3() 3()2()2()2() 1() 1() 1()(nBxnAwnwnBxnAwnwnBxnAwnw2( )(2)(2)(1)w nA w nABx nBx n) 1()2() 1()()() 1()2()3()3(2123nBxnABxinBxAinBxAinwAnBxnABxnBxAnwAiii状态方程输出方程3.由状态方程求输

39、出及单位抽样响应011000)()()(nnlnBxAnwAnwnnllnn,000nnininnn，则即，前的某时刻为设00110)()()()(nnllnnnDxlnBxACnwCAny反映系统输入序列在 n n0 后的输入零输入解)()(00nwCAnynnoi零状态解011)()()(nnllosnDxlnBxACny00110)()()()(nnllnnnDxlnBxACnwCAny00110)()()()(nnllnnnDxlnBxACnwCAny若系统矩阵A的所有特征值之模都小于1，whenn无穷大，An0)(,)(,0000nnnnnnnn为一常数，为一常数，11)()()(llnDxlnBxCAny0)()