%%%%%%%%%%%%%%%%%%%%%%%%用改進的morlet_wavelet做信號（用文件打開）的小波變換%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%得到小波變換譜圖并提取小波脊%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%求加權平均頻率%%%%%%%%%%%%%%%%%%%%%%%%%
clear?all;
close?all;
%%%%%%%%%%%%%%選擇初始值%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
nn=1024;
dt=0.001;
df=1/（nn*dt）;
sigma=3;????
beta=0;
tao=0.5;
mmax=400;?????????????%輸入頻率范圍
%%%%%%%%%%%%%%%%%%讀原始信號%%%%%%%%%%%%%%%%%%
fid=fopen（‘signal‘‘r‘）;
[f]=fread（fid[nn1]‘float‘）;
fclose（fid）;
%%%%%%%%%%%%%%%%%%進行小波變換%%%%%%%%%%%%%%%%%%%%
fw=fft（f）;
p=（2*tao/pi）^0.25*（4*（exp（-sigma^2/（2*tao））-exp（-3*sigma^2/（8*tao）））*（cos（beta*sigma））^2+...
????1-exp（-sigma^2/（2*tao）））^（-0.5）;%physical?wavelet中的參數p
q=（2*tao/pi）^0.25*（4*（exp（-sigma^2/（2*tao））-exp（-3*sigma^2/（8*tao）））*（sin（beta*sigma））^2+...
????1-exp（-sigma^2