很好的小波分析谐波的程序，奇异点捕获

很好的小波分析谐波的程序，可用于奇异点捕获，调整量提取等

clear;
%此程序主要用于标准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点处相位突变


% 周期性微小变化幅值处理中采用log
for i=100:100:1024
   %z（i）=z（i）+0.01;
    
end
%个别突变点变化
%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 

 %时域波形
%subplot（411）;
%plot（z）;
%axis（[01024-1.51.5]）;
%title（‘（a）‘）;
%Ylabel（‘幅值‘）;Xlabel（‘点数（时间）‘）;
%subplot（31

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        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\