function?wcoefs?=?mymorletcwt（SigScalesfcfb）
%============================================================%
%??Continuous?Wavelet?Transform?using?Morlet?function????????????????
%%%%%%%%%%%%%%%%%%%%%%%%輸入%%%%%%%%%%%%%%%%%%%%%%%%%%
%????Sig:?輸入信號??????????????????????????????????????????
%????Scales:?輸入的尺度序列?
%????fc:?morlet小波中心頻率??（默認為2）?
%????fb:?morlet小波帶寬參數???（默認為2）??????
%%%%%%%%%%%%%%%%%%%%%%%%輸出%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%????wcoefs:??morlet小波變換計算結果
%============================================================%
if?（nargin?<=?1）
?????error（‘At?least?2?parameters?required‘）;
end;
if?（nargin??????fb?=?2;
elseif?（nargin??????fc?=?2;
end;
wavsupport=8;???????????????????????????%?默認morlet小波的支撐區為[-88]
nLevel=length（Scales）;??????????????????%?尺度的數目
SigLen?=?length（Sig）;???????????????????%?信號的長度
wcoefs?=?zeros（nLevel?SigLen）;?????????%?分配計算結果的存儲單元?
for?m?=?1:nLevel????????????????????????%?計算各尺度上的小波系數
???
????a?=?Scales（m）;???????????????????????????????????%?提取尺度參數??????????????????????????????
????t?=?-round（a*wavsupport）:1:round（a*wavsupport）;??????????%?在尺度a的伸縮作用下，此時小波函數的支撐區會變為[-a*wavsupa*wavsup]，采樣頻率為1Hz
????Morl?=?fliplr（（pi*fb）^（-1/2）*exp（i*2*pi*fc*t/a）.*exp（-t.^2/（fb*a^2）））;?????%?計算當前尺度下的小波函數，按小波變換的定義這里需要倒置
????temp?=?conv（SigMorl）?/?sqrt（a）;????????????%?計算信號與當前尺度下小波函數的卷積???
????d=（length（temp）-SigLen）/2;??????????????????%?由于卷積計算所得結果的長度可能遠遠大于原信號，只需提取按原信號的長度提取中間部分的系數
????first?=?1+floor（d）;?????????????????????????%?區間的起點
????wcoefs（m:）=temp（first:first+SigLen-1）;???
????
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1784??2011-06-13?20:58??mymorletcwt.m
?????文件????????1435??2011-06-13?20:59??mymorletcwt_exp.m