資源簡介
優秀論文及配套源碼。電壓暫降是動態電能質量問題。優秀論文及配套源碼
代碼片段和文件信息
%%%周期延拓時將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
評論
共有 條評論