罗氏线圈与采集器系统模型建模与仿真.doc

1罗氏线圈微分回路传递函数 根据Rogowski线圈的电气特点，将其等效为图（2）所示的电路： 图1 Rogowski线圈等效电路图 其中为线圈绕组和引线总电阻；L为线圈的电感；为线圈的杂散电容；为负载电阻。 对电源部分，有： （2） 对回路部分，由电压、电流方程知 （3） （4） 联立（2）、（3）和（4）式，有 （5） 根据拉氏变换，其传递函数为 （6） （7） 令A=-M，B= LC0 ， C= ， D= 。 根据电子式互感器公司提供的参数如下： 其中C0 作为回路中的等效耦合杂散电容存在不确定性，是互感器公司反推出来的参数。 A=- 3e-3; B=1.875e-9; C= 11.15e-6; D= 1.0096; 从而得到罗氏线圈的系统传递函数系数如下： 2采集器积分回路传递函数 下图是采集器积分回路的等效电路如下图所示： 图1 积分回路等效电路 我们可以推导出其其频响函数为： 其中Rf =10M R =300K C =0.01uf 令A=- ，B=RfC。 将其带入频响函数得到积分器的传递函数原型如下： 3.系统传递函数 如下图所示设传感头的传递函数为H1，采集器的传递函数H2，则在matlab中将H1与H2做卷积运算得到系统的传递函数H。 得到系统传递函数的参数如图所示： 4.系统的幅频特性和相频特性仿真 整个系统仿真的伯德图反映了系统的幅频特性和相频特性 由图可知，-3dB截止频率：1500Hz。则需要采样频率最低为3000Hz；目前采集器为4000Hz采样频率，则保证在40次谐波内，傅氏算法保证精度。如果被采信号含40次以上谐波，则周波80点算法中基波不能保证精度； 对于保护的24点采样，则其采样频率为1200Hz，则不允许被采信号中含600hz以上分量，即不含12次以上谐波。但是1200~1500Hz之间会有信号进来，理论上是会对基波计算有影响的。 5.MATLAB仿真输入激励 在matlab中构建仿真输入激励，分别选取50hz的基波，10倍谐波，25倍谐波，40倍谐波，输入系统，并进行输出仿真。 构造激励的MATLAB代码如下： M =1; f0 = 50; Ts = 1/4000; t=0: Ts :1; n = 0:4000; y=10*sin(2*pi*M*f0*t); simin=[t,y]; 这里选取采样点为 80点，采样周期250us。 其中M表示谐波的倍数，1代表基波，10代表10次谐波。 各次谐波的输入激励波形在离散域的绘图如下： 基波输入波形离散域采样结果 10次谐波输入波形离散域采样结果 25次谐波输入波形离散域采样结果 40次谐波输入波形离散采样结果 可以看出在25次谐波下，波形失真比较严重，40次谐波下，采用4K的采样频率进行采样已经很难还原出原始波形。 6.积分回路后的暂态响应 若单独对采集器施加激励，采集器积分回路暂态响应图如下，选取基波和10次谐波作为参考。 基波采集器积分回路的暂态响应 10次谐波采集器积分回路的暂态响应 7.系统激励输入输出对比 我们利用MATLAB的simulink工具对输入激励进行仿真，以对比其在不同次谐波下输入和输出之间的特性。 其仿真模型如图所示： Simulink仿真模型 我们将输入波形和输出波形同时绘制在一个图纸上，同时选取基波，10次谐波，25次谐波和40次谐波作为对比，来研究传感头与采集器构成的系统对不同频率信号的还原能力。 基波输入输出波形比较(蓝色：输入；绿色：输出) 10次谐波输入输出波形比较(蓝色：输入；绿色：输出) 25次谐波输入输出波形比较(蓝色：输入；绿色：输出) 40次谐波输入输出波形比较(蓝色：输入；绿色：输出) 8．傅里叶变换法分析其相频特性 为了更好的分析整个系统的相频特性，我们通过分别对输入激励和输出激励进行傅里叶变换，来分析激励在通过传感头和采集单元后，其各次谐波的相频变化。