clear;
%此程序主要用于標(biāo)準(zhǔn)HARMONIC分解

%單頻及諧波
for?i=0:1023;
????z（i+1）=1*sin（2*pi*50*i/1024）;%+25*sin（2*pi*100*i/1024）;
??
?%????if?（z（i+1）>=60）?
?%??????z（i+1）=55;
?%??end????
??
end

figure（1）
%?plot（i（0:1023）z（1:1024）‘:‘）
plot（0:1023z）
%z（512）=1*sin（2*pi*50*i/1024+5/180*pi）;%512點(diǎn)處相位突變


%?周期性微小變化幅值處理中采用log
for?i=100:100:1024
???%z（i）=z（i）+0.01;
????
end
%個(gè)別突變點(diǎn)變化
%z（200）=z（200）+1;
%幅值突變???幅值處理中no?log

for?i=512:1023;
%z（i）=1.4*sin（2*pi*50*i/1024）;??%幅值突變???幅值處理中no?log
??z（i）=sin（2*pi*50*i/1024）+1*sin（2*pi*100*i/1024）;???%頻率突變???幅值處理中no?log
end



AA=hmapdn（z）;%求諧波小波均方值
for?i=1:8
???for?j=1:1024
??%??AA（ij）=1;
??end
end
%?AA?

?%時(shí)域波形
%subplot（411）;
%plot（z）;
%axis（[01024-1.51.5]）;
%title（‘（a）‘）;
%Ylabel（‘幅值‘）;Xlabel（‘點(diǎn)數(shù)（時(shí)間）‘）;
%subplot（31

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件????????1376??2007-09-13?11:06??HarmoniousWavelet\harmonicA.m
?????文件?????????444??2013-03-01?10:57??HarmoniousWavelet\hmapdn.m
?????文件?????????341??2007-07-08?21:30??HarmoniousWavelet\hwtdn.m
?????文件?????????329??2007-07-08?21:29??HarmoniousWavelet\ihwtdn.m
?????文件???????34641??2013-02-28?16:36??HarmoniousWavelet\剖面圖2.fig
?????目錄???????????0??2013-02-28?16:36??HarmoniousWavelet\