morlet小波分析matlab程序

%%%%%%%%%%%%%%%%%%%%%%%%用改进的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