优秀论文及配套源码。电压暂降是动态电能质量问题

优秀论文及配套源码。电压暂降是动态电能质量问题。优秀论文及配套源码

%%%周期延拓时将y信号进行周期平移
close all
tic
clear;
%%原始信号 
f=50;      % 基波信号频率
fs=1600;   % 采样频率
N=512;     % 采样点 
S=zeros（NN）;%生成零阶矩阵
n=0:N-1;
%x=load（‘3.txt‘）;
%x=sin（2*pi*f*n/fs）;% 标准正弦
x=sin（2*pi*f*n/fs）-（（n>=150）&（n<=300））.*（0.5*sin（2*pi*f*n/fs））;%短路
%x=sin（2*pi*f*n/fs）+0.3*sin（6*pi*f*n/fs）+0.2*sin（10*pi*f*n/fs）;
%x=sin（2*pi*f*n/fs）-（（n>=64）&（n<=512））.*（（0.2-n/（5*（N-1）））.*sin（2*pi*f*n/fs））;%感应电动机
% x=sin（2*pi*f*n/fs）-（（n>=100）&（n<=（512）））.*（0.15*sin（2*pi*f*n/fs））+（0.1）*sin（4*pi*f*n/fs）.*（（n>=100）&（n<=（512）））+（0.1）*sin（8*pi*f*n/fs）.*（（n>=100）&（n<=（512）））+（（n>=150）&（n<=512））.*（（n-150）/3000）.*sin（2*pi*f*n/fs）;
%变压器
%Diejia_D = awgn（x45）;%加入高斯噪音



%%%原始信号的离散傅立叶变换函数的产生
X=[fft（xN）fft（xN）]; 
%%%窗函数的产生
W=zeros（N-1N）;
S=zeros（NN）;
m=[0:1:N-1];%列向量0，1，2......
for n=1:N-1;
   W=exp（（-2*pi^2*m.^2）/（n*n））;%高斯窗函数
   S（n+1:） = ifft（X（（n+1）:（N+n）） .* W）; %S变换，快速傅立叶变换逆变换
end
S0=（1/N）.*sum（fft（xN）2）;
S（1:）=S0;
for n=1:N/2
    S1（n:）=S（n:）;
end
ST=abs（S1）;%绝对值

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      51637  2010-06-17 17:28  ceshi.mdl

     文件       3664  2010-06-17 11:29  tiaoshi.m

     文件    1062912  2010-11-18 16:21  基于S变换的电压暂降检测研究.doc

----------- ---------  ---------- -----  ----

              1118213                    3