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

微机继电保护实验报告 项目名称： 微机电流保护算法 姓 名： 学 号： 班 级： 实验时间： 2016.5.30 实验地点： 实验楼六楼 实验成绩： 一、实验目的 1．熟悉MATLAB桌面和命令窗口； 2．通过编写滤波程序、电流幅值计算程序以及电流保护动作判据程序，了解微机保护工作原理。 3．定性分析各种算法的优缺点。 实验步骤 1.打开MATLAB软件，建立新的编辑器； 2.依据材料所给的算法案例，设计编写滤波算法程序，要求滤除直流分量、二次谐波分量，并画出幅频特性曲线； 3.用“load”函数导入I_DATA.txt短路电流数据，对其进行滤波处理。该电流数据第一列为采样时间，第二列为采样电流，在t=0.5s时发生短路故障，短路前正常运行。并比较处理前后的波形； 4．编写微机保护算法程序，包括电流幅值算法和动作判据算法。并利用该程序对步骤2处理后的数据进行计算，观察保护的动作情况。电流保护的整定值为:425A。 三、实验模型及程序 1、设计编写滤波算法程序，要求滤除直流分量、二次谐波分量。并画出幅频特性曲线； clc;clear;%设置减法滤波器的传递函数系数 a1=1;b1=[1 0 0 0 0 0 0 0 0 0 0 0 -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); %滤波效果仿真 %模拟输入参数 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); for jj=13:m(2); Y(jj)=(Va1(jj)-Va1(jj-12))/2; end %输入波形 figure; plot(t1,Va,-r,t1,Va1,-b,t1,Y,-g); xlabel(t/s);ylabel(v/V); grid on 2、用“load”函数导入I_DATA.txt短路电流数据，对其进行滤波处理。该电流数据第一列为采样时间，第二列为采样电流，在t=0.5s时发生短路故障，短路前正常运行。并比较处理前后的波形； load(C:\Documents and Settings\Administrator\桌面\I_DATA.txt.txt); clc; a=load(C:\Documents and Settings\Administrator\桌面\I_DATA.txt.txt); %设置减法滤波器的传递函数系数 % a1=1;b1=[1 0 0 0 0 0 0 0 0 0 0 0 -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); IA=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)); for jj=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); subplot(212) plot(t1,Y,-g); xlabel(t/s);ylabel(v/V); 3、编写微机保护算法程序，包括电流幅值算法和动作判据算法。并