sagunitst.m

电能质量扰动信号

%将“持续不同时间长度”的电压凹陷信号进行归一化
clc;clear;
%产生持续时间不同的电压突降sag
t=1:640;
w0=100*pi;
a=randint(5,1,[100 900])/1000;

t1=randint(5,1,[32,160]);
t2=randint(5,1,[192,288]);
data=zeros(5,640);
y1=zeros(5,640);y2=zeros(5,640);
for i=1:5
    y1(i,:)=stepfun(t,t1(i));
    y2(i,:)=stepfun(t,t2(i));
    y3(i,:)=y1(i,:)-y2(i,:);   
    data(i,:)=(1-a(i)*y3(i,:)).*sin(w0*0.000625*t);

    subplot(2,4,1);plot(t,data(i,:));title('电压突降 no noise')
    [st_matrix,st_time,st_freq]=st(data(i,:));%st of the sag signal

    subplot(2,4,2);mesh(abs(st_matrix));
    subplot(2,4,3);contourf(abs(st_matrix));
    subplot(2,4,4);plot(st_time,abs(st_matrix));
        
     if (t2(i)-t1(i))>32
        x=(t2(i)-t1(i))/32;
        
     end
    data(i,:)=(1-a(i)*y3(i,:)).*sin(w0*0.000625*t/x);
    subplot(2,4,5);plot(t,data(i,:));title('电压突降 no noise')
    [st_matrix,st_time,st_freq]=st(data(i,:));%st of the sag signal

    subplot(2,4,6);mesh(abs(st_matrix));
    subplot(2,4,7);contourf(abs(st_matrix));
    subplot(2,4,8);plot(st_time,abs(st_matrix));   
    pause
    
end
%%%对持续时间不同的凹陷进行归一化