专利名称：：基于双电阻校准和麦夸尔特法的emi内阻抗测量方法
技术领域：
：本发明涉及一种EMI内阻抗测量方法，具体说是一种基于双电阻校准和麦夸尔特法的EMI内阻抗测量方法，属于电磁兼容
技术领域：
。
背景技术：
：目前，光伏发电、风能发电、航空电源、开关磁阻电机、电动汽车、高速动车组、城市轨道交通等产生的电磁兼容问题越来越复杂，从而严重的影响社会经济发展。因此，为了尽量降低电磁干扰噪声，节省开关时间，节约开发成本，加快新能源发展进程，同时也为产品通过质量检测部门的检验做好前期准备，必须设计合适的EMI滤波器。然而，由于噪声源内阻抗的不确定性，因此在设计EMI滤波器前，对噪声源进行内阻抗测定是必不可少的。噪声源内阻抗测定方法包括插入损耗法、单电流探头法和双电流探头法，如ZL200610098038.2、CN1996030、CN101082655、CN101093235所保护的噪声源内阻抗测定方法。然而，插入损耗法、单电流探头法和双电流探头法精度较低且无法获取阻抗的相位，另一方面，上述三种方法对被测产品的噪声源内阻抗依赖性较大，即三种方法提出的简化计算模型均与被测产品的噪声源内阻抗有关，但在实际中，被测产品的噪声源内阻抗基本特性无法获知。可见，为了提高噪声源内阻抗测定精度，同时获取内阻抗(复阻抗)的幅值与相位，并解决简化计算模型中存在的不确定性，有必要设计一个简单、高效的噪声源内阻抗测定方法。
发明内容本发明所要解决的技术问题，在于克服现有技术存在的缺陷，提出了一种基于双电阻校准和麦夸尔特法的EMI内阻抗测量方法。通过使用两个电流探头，从而获取被测产品的EMI内阻抗，并采用麦夸尔特法对测得的EMI内阻抗幅频特性进行处理，最终得到被测产品的EMI内阻抗(复阻抗)，包括幅值与相位，从而为EMI滤波器的设计提供充分的理论依据。本发明中所有阻抗均为复阻抗，包括幅值与相位，缺一不可。本发明为一种基于双电阻校准和麦夸尔特法的EMI内阻抗测量方法，其步骤是第一步将被测产品(其EMI内阻抗为Zx)、标准电阻(ZSTD)用线缆串联成闭合线路，连接频谱分析仪的检测探头和连接信号源的注入探头套在线缆上；第二步根据第一步，在闭合线路中，移去被测产品和标准电阻(即线缆自身闭合)，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vp。；第三步根据第一步，在闭合线路中，移去被测产品(其EMI内阻抗为Z》，并采用标准电阻1(ZSTD1)作为标准电阻，打开信号源使其输出电压信号V^至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpl;第四步根据第一步，在闭合线路中，移去被测产品(其EMI内阻抗为Z》，并采用标准电阻2(ZSTD2)作为标准电阻，打开信号源使其输出电压信号VsiK至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vp2;第五步根据第一步，在闭合线路中，保留被测产品(其EMI内阻抗为Z》，并移去标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为V一；第六步根据第一步，在闭合线路中，保留被测产品(其EMI内阻抗为Z》，并采用标准电阻1(ZSTD1)作为标准电阻，打开信号源使其输出电压信号V^至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpxl;第七步根据第一步，在闭合线路中，保留被测产品(其EMI内阻抗为Z》，并采用标准电阻2(ZSTD2)作为标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpx2;第八步采用如下方法处理第二、三、四步中的测量结果，计算线缆等效阻抗Z^和被测产品的噪声源内阻抗被测产品EMI内阻抗为<formula>formulaseeoriginaldocumentpage7</formula>式中，Zy为被测产f比例因子(复常数)。由式(1)可得到(1)的EMI内阻抗，Vsig为注入探头电压，VD为检测探头电压，K为<formula>formulaseeoriginaldocumentpage7</formula>根据式(5)、(6)，结合第二、三、四步中测量结果，可得(RSTDi+RiN)2+(ImSTD1+ImIN):(RSTD2+RIN)+(ImSTD2+ImIN)'KK2=K(7)VsigApl，KK2式中，KK。二K'V一/Vp。，KK丄由式(7)可得2R,Rw+2Im画=|KKf-|KK0「-RV./V。。，sigZ，p2oSTDlIm:STDl(s、2Rstd2^in+2ImSTD2Imn由式(8)可得=KK，KKn_RSTD2Im;STD2RIN=ImSTDi1P2ImSTD22RsTDlImSTDiRsTD2RsTDl1RsTD2P22RSTDlImSTDiRSTD2ImSTD2(9)式(9)其中P^lKK^-|KK。|-R〖TD1即可得到线缆等效阻抗。根据式(5)、(6)，结合第五、六、七步中》(Rx+RIN)2+(Imx+ImIN)2=|KKX0|2(RSTD1+RX+RIN)2+(ImSTD1+Imx+ImIN)2=(Rstd2+Rx+R!n)2+(ImSTD2+Imx+ImIN)2=式中，KKx。K'V一/Vpx。，KI^:K由式(10)可得ImTD1'P^lKK^-lKK,'RsTD2—ImSTD2。根据KKKKV,:结果，可得(10)sig/vpxl，KKYVsig/Vpx2oImSTDl1P2ImSTD22RsTDlImSTD1RSTD2RsTDl1"STD2P22RSTDlImSTDi"STD2ImSTD2in(11)其中P'KK|-|KK0|—R測■Im:'STDl，P2=|KK2|2-|KKRSTD2ImSTD2根据式(9)、(11)即可得到被测产品的噪声源内阻抗Zp本发明与双电流探头相比，采用了两个标准电阻，在原理中，考虑到实际中。在100KHz-30MHz频段内，标准电阻很难实现，因此采用标准阻抗进行分析，提高了测量精度，同时获取了EMI内阻抗的幅值和相位。第九步计算被测产品的全频段(100KHz-30MHz)的噪声源内阻抗；由于信号源输出的是特性频率的信号，因此要获取全频段(100KHz-30MHz)的噪声源内阻抗，需要采用麦夸尔特法进行推算，其方法为在实际中，噪声源内阻多为纯电阻(R)、纯电容(C)、电感(L)中的两者或三者串联构成，如boost电路中的内阻为纯电阻与纯电感的串联。因此，其阻抗方程为Z=R+jL+1(12)由(12)可得<formula>formulaseeoriginaldocumentpage9</formula>令<formula>formulaseeoriginaldocumentpage9</formula>式(13)可表示为<formula>formulaseeoriginaldocumentpage9</formula>式中<formula>formulaseeoriginaldocumentpage9</formula>对其进行N个频点测量，即得到N组数据(x—先给定Z的一组初始值，记为Z(°)，即<formula>formulaseeoriginaldocumentpage9</formula>将f(xT，Z)在Z(0)处按泰勒级数展开，并略去二次及二次以上的项，可得<formula>formulaseeoriginaldocumentpage9</formula>式中，Z=[R，L，C]。这是关于R，L，C的线性函数，式中除R，L，C之外均为已知数，对此用最小二乘法原则，可得<formula>formulaseeoriginaldocumentpage9</formula>其中，d>0称为阻尼因子，当取d=0时，就是高斯_牛顿法，高斯_牛顿法是，麦夸尔特法的特殊形式。欲使Q值达到最小，令Q分别对R，L，C的一阶偏导数等于零，可得<formula>formulaseeoriginaldocumentpage9</formula>式中，k=1，2，3;可化为以下形式<formula>formulaseeoriginaldocumentpage9</formula>其中，<formula>formulaseeoriginaldocumentpage9</formula>1，2，3;j=1，2，3，可得<formula>formulaseeoriginaldocumentpage9</formula>据此，可以得到被测产品的全频段噪声源内阻抗，包括幅值和相位。本发明简单、高效，测量精度较高且解决了插入损耗法、单电流探头、双电流探头法的相位缺损问题。通过采用本发明方法能够精确的获取被测产品的噪声源内阻抗，从而为EMI滤波器的设计提供理论依据。图1是本发明基于双电阻校准和麦夸尔特法的EMI内阻抗测量方法噪声源内阻抗测定实验装置电路示意图。图2是标称值为200Q的标准电阻(ZSTD1)特性曲线，其中，(a)为幅频特性曲线，(b)为相频特性曲线。图3是标称值为220Q的标准电阻(ZSTD2)特性曲线，其中，(a)为幅频特性曲线，(b)为相频特性曲线。图4是分别采用本发明方法、阻抗仪测试法、双电流探头法测量被测电容(2200pF)的对比结果。具体实施例方式下面结合附图和具体实施例，对本发明作进一步详细说明。实施例、第一步如附图1所示，将一个标称2200pF的电容(被测阻抗，Zx)、标准电阻(ZSTD)用线缆连接成闭合线路，连接固纬GSP-827频谱分析仪的检测探头和连接射频信号源的注入探头套在闭合线路上；第二步在图1的闭合线路中，移去被测电容，并采用导线代替标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为vp。；第三步在图1的闭合线路中，移去被测电容，并采用一个阻抗为200Q的标准电阻(ZSTD1)作为标准电阻，打开信号源使其输出电压信号V一至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpl;第四步在图1的闭合线路中，移去被测电容，并采用一个阻抗为220Q的标准电阻(ZSTD2)作为标准电阻，打开信号源使其输出电压信号V一至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vp2;第五步在图1的闭合线路中，保留被测电容，并采用导线代替标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpx。；第六步在图1的闭合线路中，保留被测电容，并采用一个阻抗为200Q的标准电阻(ZSTD1)作为标准电阻，打开信号源使其输出电压信号V一至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpxl;第七步在图1的闭合线路中，保留被测电容，并采用一个阻抗为220Q的标准电阻(ZSTD2)作为标准电阻，打开信号源使其输出电压信号V一至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpx2;第八步采用如下方法处理第二、三、四步中的测量结果，计算线缆等效阻抗和被测产品的噪声源内阻抗与专利ZL200610098038.2、CN1996030、CN101082655、CN101093235类似，根据双电流探头原理以及附图l，可以得到被测产品EMI内阻抗为<formula>formulaseeoriginaldocumentpage11</formula>式中，Zx为被测产品的EMI内阻抗，ZIN为电路等效阻抗，Vsig为注入探头电压，Vp为检测探头电压，K为比例因子(复常数)。由式(21)不难可得到<formula>formulaseeoriginaldocumentpage11</formula>不失一般性，设<formula>formulaseeoriginaldocumentpage11</formula>将式(23)带入(22)中，可得<formula>formulaseeoriginaldocumentpage11</formula>将式(24)两边取平方得<formula>formulaseeoriginaldocumentpage11</formula>式中，KKx二K.V一/Vp。两个标准电阻(已知)分别为<formula>formulaseeoriginaldocumentpage11</formula>其幅频特性、相频特性曲线分别如附图2、3所示。根据式(25)、(26)，结合第二、三、四步中测量结果，可得<formula>formulaseeoriginaldocumentpage11</formula>式中<formula>formulaseeoriginaldocumentpage11</formula>由式(27)可得<formula>formulaseeoriginaldocumentpage11</formula>由式(28)可得[116]<formula>formulaseeoriginaldocumentpage12</formula><formula>formulaseeoriginaldocumentpage12</formula>r。、25]其中P真<formula>formulaseeoriginaldocumentpage12</formula>，。根据式(29)、(31)即可得到被测产品的噪声源内阻抗。本发明与双电流探头相比，采用了两个标准电阻，在原理中，考虑到实际中。在lOOKHZ一30Mt{Z频段内，标准电阻很难实现，因此采用标准阻抗进行分析，提高了测量精度，同时获取了EM工内阻抗的幅值和相位。第九步采用麦夸尔特法计算全频段(100KHZ一30M}IZ)的噪声源内阻抗，其步骤为在实际中，噪声源内阻多为纯电阻(R)、纯电容(C)、电感(L)中的两者或三者串联构成，如b。。St电路中的内阻为纯电阻与纯电感的串联。因此，其阻抗方程为Z真R＼j0L+盎(32)]f(真一真真x1；R，L，C)咋．’++<formula>formulaseeoriginaldocumentpage12</formula>，式(33)可表示为<formula>formulaseeoriginaldocumentpage12</formula>(34)式中，X’一[X一1，X、]，Z一[R，L，C]。对其进行N个频点测量，即得到N组数据(X山，xln，yn)。先给定A的一组初始值，记为Z”’，即<formula>formulaseeoriginaldocumentpage12</formula>(35)将f(xT，Z)在Z(0)处按泰勒级数展开，并略去二次及二次以上的项，可得<formula>formulaseeoriginaldocumentpage13</formula>式中，Z=[R，L，C]。这是关于R，L，C的线性函数，式中除R，L，C之外均为已知数，对此用最小二乘法原则，可得<formula>formulaseeoriginaldocumentpage13</formula>其中，d>0称为阻尼因子，当取d=0时，就是高斯_牛顿法，高斯_牛顿法是，麦夸尔特法的特殊形式。<<formula>formulaseeoriginaldocumentpage13</formula>据此，可以得到被测产品的全频段噪声源内阻抗，包括幅值和相位，如附图4所权利要求基于双电阻校准和麦夸尔特法的EMI内阻抗测量方法，其步骤是第一步将被测产品、标准电阻用线缆串联成闭合线路，连接频谱分析仪的检测探头和连接信号源的注入探头套在线缆上；第二步根据第一步，在闭合线路中，移去被测产品和标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vp0；第三步根据第一步，在闭合线路中，移去被测产品，并采用标准电阻a(ZSTD1)作为标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vp1；第四步根据第一步，在闭合线路中，移去被测产品，并采用标准电阻b(ZSTD2)作为标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vp2；第五步根据第一步，在闭合线路中，保留被测产品，并移去标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpx0；第六步根据第一步，在闭合线路中，保留被测产品，并采用标准电阻a(ZSTD1)作为标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpx1；第七步根据第一步，在闭合线路中，保留被测产品，并采用标准电阻b(ZSTD2)作为标准电阻，打开信号源使其输出电压信号Vsig至注入探头，连接检测探头的频谱分析仪得到相应的测量结果，记为Vpx2；第八步计算线缆等效阻抗ZIN和被测产品的噪声源内阻抗被测产品EMI内阻抗为<mrow><msub><mover><mi>Z</mi><mo>·</mo></mover><mi>X</mi></msub><mo>=</mo><mfrac><mrow><mi>K</mi><msub><mover><mi>V</mi><mo>·</mo></mover><mi>sig</mi></msub></mrow><msub><mover><mi>V</mi><mo>·</mo></mover><mi>p</mi></msub></mfrac><mo>-</mo><msub><mover><mi>Z</mi><mo>·</mo></mover><mi>IN</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>式中，ZX为被测产品的EMI内阻抗，Vsig为注入探头电压，Vp为检测探头电压，K为比例因子(复常数)；由式(1)可得到<mrow><mo>|</mo><msub><mover><mi>Z</mi><mo>·</mo></mover><mi>X</mi></msub><mo>+</mo><msub><mover><mi>Z</mi><mo>·</mo></mover><mi>IN</mi></msub><mo>|</mo><mo>=</mo><mo>|</mo><mfrac><mrow><mi>K</mi><msub><mover><mi>V</mi><mo>·</mo></mover><mi>sig</mi></msub></mrow><msub><mover><mi>V</mi><mo>·</mo></mover><mi>p</mi></msub></mfrac><mo>|</mo><mo>=</mo><mfrac><mrow><mo>|</mo><mi>K</mi><mo>|</mo><mo>|</mo><msub><mover><mi>V</mi><mo>·</mo></mover><mi>sig</mi></msub><mo>|</mo></mrow><mrow><mo>|</mo><msub><mover><mi>V</mi><mo>·</mo></mover><mi>p</mi></msub><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>设<mrow><msub><mover><mi>Z</mi><mo>·</mo></mover><mi>X</mi></msub><mo>=</mo><msub><mi>R</mi><mi>X</mi></msub><mo>+</mo><msub><mi>jIm</mi><mi>X</mi></msub></mrow>(3)<mrow><msub><mover><mi>Z</mi><mo>·</mo></mover><mi>IN</mi></msub><mo>=</mo><msub><mi>R</mi><mi>IN</mi></msub><mo>+</mo><msub><mi>jIm</mi><mi>IN</mi></msub></mrow>将式(3)带入(2)中，可得<mrow><mo>|</mo><mrow><mo>(</mo><msub><mi>R</mi><mi>X</mi></msub><mo>+</mo><msub><mi>R</mi><mi>IN</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>j</mi><mrow><mo>(</mo><msub><mi>Im</mi><mi>X</mi></msub><mo>+</mo><msub><mi>Im</mi><mi>IN</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>=</mo><mfrac><mrow><mo>|</mo><mi>K</mi><mo>|</mo><mo>|</mo><msub><mover><mi>V</mi><mo>·</mo></mover><mi>sig</mi></msub><mo>|</mo></mrow><mrow><mo>|</mo><msub><mover><mi>V</mi><mo>·</mo></mover><mi>p</mi></msub><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>将式(4)两边取平方得(RX+RIN)2+(ImX+ImIN)2＝|KKX|2(5)式中，KKX＝K·Vsig/Vp；已知标准电阻a、b分别为ZSTD1＝RSTD1+jImSTD1(6)ZSTD2＝RSTD2+jImSTD2根据式(5)、(6)，可得<mrow><msubsup><mi>R</mi><mi>IN</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>Im</mi><mi>IN</mi><mn>2</mn></msubsup><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup></mrow>(RSTD1+RIN)2+(ImSTD1+ImIN)2＝|KK1|2(7)(RSTD2+RIN)2+(ImSTD2+ImIN)2＝|KK2|2式中，KK0＝K·Vsig/Vp0，KK1＝K·Vsig/Vp1，KK2＝K·Vsig/Vp2；由式(7)可得<mrow><mn>2</mn><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub><msub><mi>R</mi><mi>IN</mi></msub><mo>+</mo><mn>2</mn><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub><msub><mi>Im</mi><mi>IN</mi></msub><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>1</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow><mn>2</mn></msubsup></mrow>(8)<mrow><msub><mrow><mn>2</mn><mi>R</mi></mrow><mrow><mi>STD</mi><mn>2</mn></mrow></msub><msub><mi>R</mi><mi>IN</mi></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>Im</mi></mrow><mrow><mi>STD</mi><mn>2</mn></mrow></msub><msub><mi>Im</mi><mi>IN</mi></msub><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>2</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow><mn>2</mn></msubsup></mrow>由式(8)可得线缆等效阻抗<mrow><msub><mi>R</mi><mi>IN</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>P</mi><mn>1</mn></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>P</mi><mn>2</mn></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced><mfencedopen='|'close='|'separators=''><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced></mfrac></mrow><mrow><msub><mi>Im</mi><mi>IN</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>其中<mrow><msub><mi>P</mi><mn>1</mn></msub><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>1</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow><mn>2</mn></msubsup><mo>,</mo></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>2</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>;</mo></mrow>根据式(9)即可得到线缆等效阻抗；根据式(5)、(6)，可得(RX+RIN)2+(ImX+ImIN)2＝|KKX0|2(RSTD1+RX+RIN)2+(ImSTD1+ImX+ImIN)2＝|KKX1|2(10)(RSTD2+RX+RIN)2+(ImSTD2+ImX+ImIN)2＝|KKX2|2式中，KKX0＝K·Vsig/Vpx0，KKX1＝K·Vsig/Vpx1，KKX2＝K·Vsig/Vpx2；由式(10)可得<mrow><msub><mi>R</mi><mi>X</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>P</mi><mn>1</mn></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>P</mi><mn>2</mn></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced></mfrac><mo>-</mo><msub><mi>R</mi><mi>IN</mi></msub></mrow><mrow><msub><mi>Im</mi><mi>X</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mfrac><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>P</mi><mn>2</mn></msub></mtd></mtr></mtable></mfenced><mfencedopen='|'close='|'><mtable><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow></msub></mtd></mtr></mtable></mfenced></mfrac><mo>-</mo><msub><mi>Im</mi><mi>IN</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>其中<mrow><msub><mi>P</mi><mn>1</mn></msub><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>1</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>R</mi><mrow><mi>STD</mi><mn>1</mn></mrow><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>Im</mi><mrow><mi>STD</mi><mn>1</mn></mrow><mn>2</mn></msubsup><mo>,</mo></mrow><mrow><msub><mi>P</mi><mn>2</mn></msub><mo>=</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>2</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mrow><mo>|</mo><msub><mi>KK</mi><mn>0</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msubsup><mi>R</mi><mrow><mi>STD</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>Im</mi><mrow><mi>STD</mi><mn>2</mn></mrow><mn>2</mn></msubsup><mo>;</mo></mrow>根据式(9)、(11)即可得到被测产品的噪声源内阻抗ZX；第九步采用麦夸尔特法计算被测产品的全频段噪声源内阻抗，其方法为阻抗方程为<mrow><mi>Z</mi><mo>=</mo><mi>R</mi><mo>+</mo><mi>jωL</mi><mo>+</mo><mfrac><mn>1</mn><mi>jωC</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>由(12)可得<mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>;</mo><mi>R</mi><mo>,</mo><mi>L</mi><mo>,</mo><mi>C</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>|</mo><mi>Z</mi><mo>|</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>ωL</mi><mo>+</mo><mfrac><mn>1</mn><mi>ωC</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><mi>R</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><mi>LC</mi><mo>+</mo><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup><msup><mi>L</mi><mn>2</mn></msup><msup><mi>f</mi><mn>2</mn></msup><mo>+</mo><mfrac><mn>1</mn><mrow><mn>4</mn><msup><mi>π</mi><mn>2</mn></msup><msup><mi>C</mi><mn>2</mn></msup></mrow></mfrac><msup><mi>f</mi><mrow><mo>-</mo><mn>2</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>令f2＝x1，f-2＝x-1，式(13)可表示为y＝f(xT，Z)＝f(x-1，x1；R，L，C)(14)式中，xT＝[x-1，x1]，Z＝[R，L，C]；对其进行N个频点测量，即得到N组数据(x-1n，x1n，yn)；先给定Z的一组初始值，记为Z(0)，即Z(0)＝[R(0)，L(0)，C(0)](15)将f(xT，Z)在Z(0)处按泰勒级数展开，并略去二次及二次以上的项，可得<mrow><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><mi>Z</mi><mo>)</mo></mrow><mo>≈</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>)</mo></mrow><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mfrac><mrow><mo>∂</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><mi>Z</mi><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>Z</mi><mi>m</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>Z</mi><mo>=</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup></mrow></msub><mrow><mo>(</mo><msub><mi>Z</mi><mi>m</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>m</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>式中，Z＝[R，L，C]；这是关于R，L，C的线性函数，式中除R，L，C之外均为已知数，对此用最小二乘法原则，可得<mrow><mi>Q</mi><mo>=</mo><msup><mrow><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><mrow><mo>(</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>)</mo></mrow><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mfrac><mrow><mo>∂</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><mi>Z</mi><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>Z</mi><mi>m</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>Z</mi><mo>=</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup></mrow></msub><mrow><mo>(</mo><msub><mi>Z</mi><mi>m</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>m</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mo>+</mo><mi>d</mi><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><msup><mrow><mo>(</mo><msub><mi>Z</mi><mi>m</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>m</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>其中，d≥0称为阻尼因子；欲使Q值达到最小，令Q分别对R，L，C的一阶偏导数等于零，可得<mrow><mn>0</mn><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>Q</mi></mrow><mrow><mo>∂</mo><msub><mi>Z</mi><mi>k</mi></msub></mrow></mfrac><mo>=</mo><mn>2</mn><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><mrow><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>)</mo></mrow><mo>+</mo><munderover><mi>Σ</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mn>3</mn></munderover><mfrac><mrow><mo>∂</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><mi>Z</mi><mo>)</mo></mrow></mrow><msub><mrow><mo>∂</mo><mi>Z</mi></mrow><mi>m</mi></msub></mfrac><msub><mo>|</mo><mrow><mi>Z</mi><mo>=</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup></mrow></msub><mrow><mo>(</mo><msub><mi>Z</mi><mi>m</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>m</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo></mo></mrow><mo>)</mo></mrow><mfrac><mrow><mo>∂</mo><mi>f</mi><mrow><mo>(</mo><msup><mi>x</mi><mi>T</mi></msup><mo>,</mo><mi>Z</mi><mo>)</mo></mrow></mrow><mrow><mo>∂</mo><msub><mi>Z</mi><mi>k</mi></msub></mrow></mfrac><msub><mo>|</mo><mrow><mi>Z</mi><mo>=</mo><msup><mi>Z</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup></mrow></msub><mo>+</mo><mn>2</mn><mi>d</mi><mrow><mo>(</mo><msub><mi>Z</mi><mi>m</mi></msub><mo>-</mo><msubsup><mi>Z</mi><mi>m</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>式中，k＝1，2，3；可化为以下形式(a11+d)(R-R(0))+a12(L-L(0))+a13(C-C(0))＝a1ya21(R-R(0))+(a22+d)(L-L(0))+a23(C-C(0))＝a2y(19)a31(R-R(0))+a32(L-(0))+(a33+d)(C-C(0))＝a3y其中，i＝1，2，3；j＝1，2，3，可得<mrow><mfencedopen='('close=')'><mtable><mtr><mtd><mi>R</mi></mtd></mtr><mtr><mtd><mi>L</mi></mtd></mtr><mtr><mtd><mi>C</mi></mtd></mtr></mtable></mfenced><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msub><mi>Z</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mn>3</mn></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfencedopen='('close=')'><mtable><mtr><mtd><msubsup><mi>Z</mi><mn>1</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mn>2</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Z</mi><mn>3</mn><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msubsup></mtd></mtr></mtable></mfenced><msup><mo>+</mo></msup><msup><mfencedopen='('close=')'><mtable><mtr><mtd><msub><mi>a</mi><mn>11</mn></msub><mo>+</mo><mi>d</mi></mtd><mtd><msub><mi>a</mi><mn>12</mn></msub></mtd><mtd><msub><mi>a</mi><mn>13</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>21</mn></msub></mtd><mtd><msub><mi>a</mi><mn>22</mn></msub><mo>+</mo><mi>d</mi></mtd><mtd><msub><mi>a</mi><mn>23</mn></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>31</mn></msub></mtd><mtd><msub><mi>a</mi><mn>32</mn></msub></mtd><mtd><msub><mi>a</mi><mrow><mn>33</mn><mo>+</mo><mi>d</mi></mrow></msub></mtd></mtr></mtable></mfenced><mrow><mo>-</mo><mn>1</mn></mrow></msup><mfencedopen='('close=')'><mtable><mtr><mtd><msub><mi>a</mi><mrow><mn>1</mn><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mn>2</mn><mi>y</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>a</mi><mrow><mn>3</mn><mi>y</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>20</mn><mo>)</mo></mrow></mrow>据此，可以得到被测产品的全频段噪声源内阻抗，包括幅值和相位。F2009102629750C00053.tif全文摘要本发明公开了一种基于双电阻校准和麦夸尔特法的EMI内阻抗测量方法。其中，一个电流探头作为注入探头，另一个探头作为检测探头分别测量短路、标准电阻1、标准电阻2、被测阻抗、标准电阻1与被测阻抗串联、标准电阻2与被测阻抗串联时的电路状态量，从而获取被测产品的EMI内抗，并采用麦夸尔特法对测得的EMI内阻抗幅频特性进行处理，最终得到被测产品的EMI内阻抗(复阻抗)，包括幅值与相位。本发明简单、高效，测量精度较高且解决了插入损耗法、单电流探头、双电流探头法的相位缺损问题。通过采用本发明方法能够精确的获取被测产品的噪声源内阻抗，从而为EMI滤波器的设计提供理论依据。文档编号G01R27/02GK101710157SQ20091026297公开日2010年5月19日申请日期2009年12月14日优先权日2009年12月14日发明者姜宁秋,戎融,王恩荣,董颖华,赵阳,陆婋泉,颜伟申请人:南京师范大学