tic
clearclose?all?
%%?無隨機取，取1到9個連續整周波擾動，在十個整周波上，（含噪），200個測試集全部模型?
?
w=100*pi;??
t=0:1/6400:0.2;%取1028?2^10=1024?
t1=0.071;?
t2=0.142;?
%pure?sin??正常正弦

figure（1）;
for?c1=1:1:200;?
subplot（411）;?
?u0（c1:）=sin（w*t）;?
?plot（tu0）;?
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘正常電壓波形‘）;??
end

%%?swell?電壓驟升
for?c1=1:1:200;?
subplot（412）;?
u（c1:）=sin（w*t）.*（1+0.5*（heaviside（t-t1）-heaviside（t-t2）））;?
plot（tu）;
ylabel（‘幅值‘）;
xlabel（‘t/s‘）;
title（‘電壓驟升波形‘）;?
end?

%%?sag?電壓驟降

for?c1=1:1:200;?
v（c1:）=sin（w*t）.*（1-0.8*（heaviside（t-t1）-heaviside（t-t2）））;?
subplot（413）;
plot（tv）;
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘電壓驟降波形‘）;?
end?

%%?interupt??電壓中斷
for?c1=1:1:200;?
vi（c1:）=sin（w*t）.*（1-1*（heaviside（t-t1）-heaviside（t-t2）））;?
subplot（414）;
plot（tvi）;?
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘電壓中斷波形‘）;?
end?
%%?homonics??諧波（間）
figure（2）;
for?c1=1:1:200;?
?vv（c1:）=sin（w*t）+[0.1*sin（3*w*t）+0.1*sin（5*w*t）+0.1*sin（7*w*t）].*（heaviside（t-t1）-heaviside（t-t2））;?
?subplot（411）;?
plot（tvv）;?
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘諧波波形‘）;?
end?

%%?pluse?暫態脈沖

for?c1=1:1:200;?
z（c1:）?=sin（w*t）.*（1-2*（heaviside（t-0.071）-heaviside（t-0.07115625））-2*（heaviside（t-0.142）-heaviside（t-0.14215625）））;??
%?[z11t3]?=?gensig（‘pulse‘t10.20.00015625）;%1.2??0.0011673152?
%?z1（c1:）=sin（w*t）+[0.4*z11]‘;%幅度0.1-0.8?
subplot（412）;
plot（tz）;
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘暫態脈沖波形‘）;?
end?
%%%%%[u??t]=gensig（type??tau?tf??ts）
%%%%%Type—輸入信號類型（sin正弦信號；?square?方波信號?；pulse周期脈沖信號）
%%%%%Tau—輸入信號的周期??tf—信號的持續時間??ts—信號時間的步長
%%%%%??[u??t]—輸出信號歷史記錄

%%?osc?暫態振蕩
for?c1=1:1:200;?
zz（c1:）=sin（w*t）+0.4*exp（10*（t1-t2））*sin（20*w*t）.*（heaviside（t-t1）-heaviside（t-t2））;%取20比較象?osc.m?
subplot（413）;
plot（tzz）;
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘暫態振蕩波形‘）;?
end?

%%?flicker?電壓閃變
for?c1=1:1:200;
??g（c1:）=sin（w*t）.*[1+（0.1*（sin（5*w*t）+0.5*sin（10*w*t）+（1/3）*sin（15*w*t）））.*（1*（heaviside（t-t1）-heaviside（t-t2）））];
?subplot（414）;
plot（tg）;
?ylabel（‘幅值‘）;
?xlabel（‘t/s‘）;
?title（‘電壓閃變波形‘）;?
end

%%?小波分解與能量系數提取
??for?c1=1:1:200;?
??
[c0l0]=wavedec（u0（c1:）10‘db4‘）;
[Ea0Ed0（c1:）]?=?wenergy（c0l0）;

[culu]=wavedec（u（c1:）10‘db4‘）;?
[EauEdu（c1:）]?=?wenergy（culu）;
?EDu（c1:）?=?Edu（c1:）-Ed0（c1:）;

[cvlv]=wavedec（v（c1:）10‘db4‘）;?
[EavEdv（c1:）]?=?wenergy（cvlv）;?
?EDv（c1:）?=?Edv（c1:）-Ed0（c1:）;

[cvilvi]=wavedec（vi（c1:）10‘db4‘）;?
[EaviEdvi（c1:）]?=?wenergy（cvilvi）;?
EDvi（c1:）?=?Edvi（c1:）-Ed0（c1:）;

[cvvlvv]=wavedec（vv（c1:）10‘db4‘）;??
[EavvEdvv（c1:）]?=?wenergy（cvvlvv）;??
EDvv（c1:）?=?Edvv（c1:）-Ed0（c1:）;

[czlz]=wavedec（z（c1:）10‘db4‘）;?
[EazEdz（c1:）]?=?wenergy（czlz）;?
EDz（c1:）?=?Edz（c1:）-Ed0（c1:）;

[czzlzz]=wavedec（zz（c1:）10‘db4‘）;?
[EazzEdzz（c1:）]?=?wenergy（czzlzz）;?
?EDzz（c1:）?=?Edzz（c1:）-Ed0（c1:）;

[cglg]=wavedec（g（c1:）10‘db4‘）;?
[EagEdg（c1:）]?=?wenergy（cglg）;?
EDg（c1:）?=?Edg（c1:）-Ed0（c1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????8723??2012-03-22?15:14??pq_pcasvm\buhanzaoxiaobonengliangcha.m

?????文件???????9270??2012-03-24?23:07??pq_pcasvm\hanzaoxiaobonengliangcha002.m

?????文件???????9698??2012-06-22?21:24??pq_pcasvm\hanzaoxiaobonengliangcha01.m

?????文件???????4640??2012-06-29?14:43??pq_pcasvm\PCASVM.m

?????文件???????2637??2012-05-21?14:43??pq_pcasvm\pcv1.m

?????文件???????6483??2012-06-22?21:26??pq_pcasvm\xiaobonengliangcha.m

?????目錄??????????0??2012-07-23?09:39??pq_pcasvm

-----------?---------??----------?-----??----

????????????????41451????????????????????7