untitled.asv

电能质量扰动信号

clc;clear;
%震荡暂态oscillatory transients信号
t=0:300;
w0=100*pi;
t1=64;t2=128;%暂态时间2个T
y1=stepfun(t,t1);y2=stepfun(t,t2);
y3=y1-y2;
A=randint(1,1,[100 800])/1000;
c=randint(1,1,[500 1000])/100;
a=randint(1,1,[1000 1500])/100;
data=sin(w0*0.000625*t)+A*exp(-c*(t*0.000625-t1*0.000625)).*sin(a*w0*0.000625*t).*y3;%产生震荡暂态信号
subplot(2,1,1);plot(t,data);
t=0:0.000625:0.2;
t1=0.04;t2=0.0025;
y1=stepfun(t,t1);y2=stepfun(t,t2);
y3=y1-y2;
data1=sin(w0*t)+A*exp(-c*(t-t1)).*sin(a*w0*t).*y3;
subplot(2,1,2);plot(t,data1);