sagst.asv

电能质量扰动信号

%电压突降（voltage sag）的ST变换
clc;clear;
%电压突降sag
t=0:600  %(600/32=18.75个周期,每个周期可采样32个点，因为基频50赫兹，采样频率1.6KHZ);
w0=100*pi;%w0=2*pi/T=100*pi,T=0.02;
t1=64;t2=224;%起止时间设定为第3个周波到第7个周波；
y1=stepfun(t,t1);y2=stepfun(t,t2);
y3=y1-y2;
a=randint(1,1,[100 900])/1000;
data=(1-a*y3).*sin(w0*0.000625*t);%%%%%产生SAG信号
noise=awgn(data,30,'measured');%产生30dB的高斯白噪声
data_noise=data+noise;%将产生的30dB高斯白噪声迭加到SAG信号中
%figure;
%no noise
%以下代码输出相应的图片
subplot(2,4,1);c
[st_matrix,st_time,st_freq]=st(data);%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));

%add noise 30 dB
%subplot(2,4,5);plot(t,data_noise);title('电压突降 add noise')
%[st_matrix_noise,st_time_noise,st_freq_noise]=st(data_noise);
%subplot(2,4,6);mesh(st_time,st_freq,abs(st_matrix_noise));
%subplot(2,4,7);contourf(st_time,st_freq,abs(st_matrix_noise));
%subplot(2,4,8);plot(st_time,abs(st_matrix_noise));
%figure;contourf(abs(st_matrix_noise));title('电压突降等高线')
%figure;contourf(st_time,st_freq,abs(st_matrix));
ss=abs(st_matrix);
[m,n]=size(ss);
ss_t=zeros(1,n);
for i=1:n
    sum=0;
    for j=1:m
        sum=sum+ss(j,i).*ss(j,i);
    end
    mean=sum/m;
    ss_t(1,i)=mean;
end
figure;
subplot(3,1,1);
plot(ss_t);title('S矩阵各列的平方和均值');