%%%周期延拓時將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