微机继电保护实验报告微机电流保护算法

《微机继电保护实验报告微机电流保护算法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、微机继电保护实验报告项目名称：微机电流保护算法姓名：学号：班级：实验时间：2016.5.30实验地点：实验楼六楼实验成绩：一、实验目的1．熟悉MATLAB桌面和命令窗口；2．通过编写滤波程序、电流幅值计算程序以及电流保护动作判据程序，了解微机保护工作原理。3．定性分析各种算法的优缺点。二、实验步骤1.打开MATLAB软件，建立新的编辑器；2.依据材料所给的算法案例，设计编写滤波算法程序，要求滤除直流分量、二次谐波分量，并画出幅频特性曲线；3.用“load”函数导入I_DATA.txt短路电流数据，对其进行滤波处理。该电流数据第一列为采样时间，第二列为采样

2、电流，在t=0.5s时发生短路故障，短路前正常运行。并比较处理前后的波形；4．编写微机保护算法程序，包括电流幅值算法和动作判据算法。并利用该程序对步骤2处理后的数据进行计算，观察保护的动作情况。电流保护的整定值为:425A。三、实验模型及程序1、设计编写滤波算法程序，要求滤除直流分量、二次谐波分量。并画出幅频特性曲线；clc;clear;%设置减法滤波器的传递函数系数a1=1;b1=[100000000000-1];f=0:1:600;h1=abs(freqz(b1,a1,f,1200));%由传递函数系数确定传递函数的幅频特性H1=h1/max(h1)

3、;%绘制幅频特性plot(f,H1);xlabel('f/Hz');ylabel('H1')figureplot(f,h1);%滤波效果仿真%模拟输入参数N=24;t1=(0:0.02/N:0.04)m=size(t1);%基波电压Va=100*sin(2*pi*50*t1);%叠加直流分量和4,8次谐波Va1=75*exp(-0.1*t1)+100*sin(2*pi*50*t1)+30*sin(4*pi*50*t1);%采用减法滤波器滤除谐波Y=zeros(1,12);forjj=13:m(2);Y(jj)=(Va1(jj)-Va1(jj-12))/2

4、;end%输入波形figure;plot(t1,Va,'-r',t1,Va1,'-b',t1,Y,'-g');xlabel('t/s');ylabel('v/V');gridon2、用“load”函数导入I_DATA.txt短路电流数据，对其进行滤波处理。该电流数据第一列为采样时间，第二列为采样电流，在t=0.5s时发生短路故障，短路前正常运行。并比较处理前后的波形；load('C:DocumentsandSettingsAdministrator桌面I_DATA.txt.txt');clc;a=load('C:DocumentsandSett

5、ingsAdministrator桌面I_DATA.txt.txt');%设置减法滤波器的传递函数系数%a1=1;b1=[100000000000-1];%f=0:1:600;%h1=abs(freqz(b1,a1,f,1200));%%由传递函数系数确定传递函数的幅频特性%H1=h1/max(h1);%%绘制幅频特性%plot(f,H1);%xlabel('f/Hz');ylabel('H1')%figure%plot(f,h1);%滤波效果仿真%模拟输入参数Ts=a(2,1)-a(1,1);N=0.02/Ts;k=N/2;t1=a(:,1);I

6、A=a(:,2);%t1=(0:0.02/N:0.04)m=size(t1);%基波电压%Va=100*sin(2*pi*50*t1);%叠加直流分量和4,8次谐波%Va1=75*exp(-0.1*t1)+100*sin(2*pi*50*t1)+30*sin(4*pi*50*t1);%采用减法滤波器滤除谐波Y=zeros(1,m(1));forjj=k+1:m(1);Y(jj)=(IA(jj)-IA(jj-k))/2;end%输入波形figure;subplot(211)plot(t1,IA,'r');xlabel('t/s');ylabel('v/V'

7、);subplot(212)plot(t1,Y,'-g');xlabel('t/s');ylabel('v/V');3、编写微机保护算法程序，包括电流幅值算法和动作判据算法。并利用该程序对步骤2处理后的数据进行计算，观察保护的动作情况。电流保护的整定值为:425A。电流保护clc;clear;%数据导入，数据N=200a=load('C:DocumentsandSettingsAdministrator桌面I_DATA.txt.txt');Ts=a(2,1)-a(1,1);%采样周期N=0.02/Ts;%每周期采样点数t=a(:,1)';I=a(

8、:,2)';m=size(t);%数字滤波1，差分滤波Y=zeros(1,m(2