事件检测代码，比较简单，仅供参考.rar_nilm_nilm事件检测_事件_事件检测_事件检测；NILM

function S = eventdetection(X) %X表示待检测的负荷采样信号序列(假定为有功功率) %S表示返回的事件标记序列，低频情况下的标记点为事件点，高频情况下的标记点为事件发生时段的中心点 %事件标记序列中，1表示开启事件，-1表示关闭事件 %此事件检测方法只是对事件的初步检测，无法区分用电器事件种类 a=size(X,2);%返回序列X的采样点个数 a S=zeros(1,a); SS=zeros(1,a); b=20;%第一次事件检测阈值，可视事件发生情况自行调整 bb=100;%第二次事件检测阈值，可视事件发生情况自行调整 %============================================================================= %1、低频采样情况下的事件检测（忽略了事件的发生/持续时间） for i=2:a%事件检测(第一次，利用相邻采样点差值进行检测，较为粗糙，针对低频(30秒~1分钟/次的采样)序列) if X(i)-X(i-1)>b%相邻采样点差值超过阈值 S(i-1)=1;%负荷开启事件 end if X(i)-X(i-1)<-b S(i-1)=-1;%负荷关闭事件 end end for j=2:a%事件检测(第二次，针对第一次检测的事件发生可疑点进行二次检测，较为精细，针对低频) if (S(j)==1)&&(j>3)&&(j<a-3)%负荷开启事件 c=(X(j)+X(j-1)+X(j-2))/3;%选前三个点求平均值 cc=(X(j+1)+X(j+2)+X(j+3))/3;%选后三个点求平均值 if abs(c-cc)<bb S(j)=0;%如未超过阈值，判定未发生事件，抹去原标记点 end end if (S(j)==-1)&&(j>3)&&(j<a-3)%负荷关闭事件 c=(X(j)+X(j-1)+X(j-2))/3; cc=(X(j+1)+X(j+2)+X(j+3))/3; if abs(c-cc)<bb S(j)=0; end end end %============================================================================ % %============================================================================ % %2、高频采样情况下的事件检测（未 忽略事件的发生/持续时间） % for i=1:a-10%事件检测(第一次，较为粗糙，针对高频(1~10秒/次的采样)序列) % Y=X(i:i+10);%每十一个点放入一个数组中检测一下（计算量较大） % SS(i)=var(Y,1);%计算十一个点的方差 % if i>1 % if SS(i)-SS(i-1)>b%方差变化超过阈值（可根据实际情况来设） % S(i+5)=1;%将十一个点的中间点对应标号进行标记（这里暂不区分开启和关闭事件） % end % end % end % % for j=2:a%事件检测(第二次，针对第一次检测的事件发生可疑点进行二次检测，较为精细，针对高频，默认事件发生/持续时间为7个采样点) % c=0;cc=0;f=6; % for e=1:10 % c=c+X(j+1-e-f/2);%标记点前3点旁边的10个点求平均值 % cc=cc+X(j+e+f/2);%标记点后3点旁边的10个点求平均值 % end % c=c/10;%标记点前3点旁边的10个点求平均值 % cc=cc/10;%标记点后3点旁边的10个点求平均值 % if abs(c-cc)<bb % S(j)=0;%如未超过阈值，判定未发生事件，抹去原标记点 % end % if abs(c-cc)>bb%如超过阈值 % if cc-c<0 % S(j)=-1;%通过前后点平均值大小来判定是开启事件还是关闭事件 % end % end % end % %============================================================================ %图像展示一下 figure(1); subplot(2,1,1);plot(X); subplot(2,1,2);plot(S); end