オリジナル：https://blog.csdn.net/weixin_40583722/article/details/79892289
まず、原理を描きます：
CWT三つの機能は、ウェーブレットツールボックスを使用する必要がある（）、centfrq（）、scal2frq（）
COEFS = CWT（S、スケール、 「wname」）
連続ウェーブレットのwname、ウェーブレット名により、スケールのためにSは入力信号である関数、スケールを変換します。
FREQ = centfrq（「wname」）
という名前wnameウェーブレット母の中心周波数を見つけるための機能。
F = scal2frq（A、「wname 」、DELTA）
この関数は、Aは、wname、名前、DELTAによってウェーブレットサンプリング周期スケールされる実際の周波数スケールに変換することができます。
スケールと周波数との関係
FS / *スケールを聞かせ、fsはサンプリング周波数であり、Fcは中心周波数をウェーブレット、対応する実際の周波数FaがFaを=たFcでは、明らかに、サンプリング定理によれば、周波数範囲にスケイログラムは、（0でありますFS / 2）、スケールの範囲は、INFが無限を表す（2 * FC、INF）であるべきです。実際には、取るためにちょうど十分な大きさの規模。
スケールの配列を決定するために、
上記の式から分かるように、配列への変換器の周波数は、以下の形式をとる必要があり、算術、スケール配列の配列である：C / totalscal、C /（totalscal-1）、...、C / 2、C請求、配列の長さスケールtotalscal（通常プログラムPRE）信号をウェーブレット変換するときに使用される、Cは定数です。対応totalscal実際の周波数スケールcは/ FS / 2は、次いで、C = 2 *のFc * totalscal配列こうして必要な寸法を得て得ることができるべきです。
時間周波数スケッチ
ウェーブレット基盤と規模を決定した後、我々はCWTウェーブレット係数coefs用途を見つけることができる（係数が複素数を法とする）、及び実際の周波数にscal2frq規模配列は関数imagesc（T、で、F、最後の時系列Tを配列パイロットは、ウェーブレット図を描くことができるようになり、F、ABS（coefs））。
第二に、適用例は、
ウェーブレット周波数を図にプロットされたときに具体例を例示するために以下に与えられます。シミュレートされた信号は、それぞれ、合成信号の50Hzのと100Hzの2つの正弦波成分の周波数で撮影しました。
周波数解析％ウェーブレット
CLC
すべてクリア
閉じるすべての
％原信号
FS = 1000;
F1 = 50;
F2 = 100;
T = 0 :. 1 / FS :. 1;
S = SIN（2 * PI * F1 * T）+ SIN（ * F2 * PI * 2 T）;
図
プロット（T、S）
％連続ウェーブレット変換
; wavename 'はcmor3-3' =
totalscal = 256;
のFc = centfrq（wavename）;％ウェーブレットの中心周波数
C = 2 *のFc * ; totalscal
scals = C ./（1：totalscal）;
F = scal2frq（scals、wavename ,. 1 / FS）;％コンバート周波数スケールに
coefs = CWT（S、scals、 wavename）;％は、 連続ウェーブレット係数シーク
図を
関数imagesc （T、F、ABS（coefs））;
SET（GCA、 'YDIR'、 '通常'）
カラーバー;
は、xlabel（ '時間T / S'）;
（ '周波数F /ヘルツ'）ylabelの;
タイトル（「ウェーブレットを使用して図頻度「）。
説明：
この例では、モレット複素ウェーブレットの最良の選択は、他の効果を分析しないウェーブレットとモレットウェーブレットの中心周波数と帯域幅パラメータを大きく、図周波数濃度の頻度が良く反映します。
他の時間-周波数解析との比較、短い時間として時間周波数プロットは、フーリエ変換
要約：
ウェーブレット変換の分解能は、高周波、高解像度で自動的に調整することができるのでとき、高周波は、良い効果をウェーブレット。
周波数解析％
CLC
すべてクリア
すべて閉じる
％元の信号
FS = 1000;
F1 = 50;
F2 = 100;
T = 0 :. 1 / FS :. 1;
S = SIN（2 * PI * F1 * T）+ SIN（2 PI * F2 * T *）;関数randn％+）（1、長さ（T）;
％S =チャープ（T、20,0.3,300）;
S =チャープ（T、20,1,500、 'Q'）;
図
プロット（T、S）
％の頻度図連続ウェーブレット変換
wavename 'をcmor3-3を' =;
totalscal = 256;
のFc = centfrq（wavename）;ウェーブレット％の中心周波数
C = 2 *のFc totalscal;
scals = C /。 （1：totalscal）;
F = scal2frq（scals、wavename ,. 1 / FS）;％スケール周波数変換
coefs = CWT（S、scals、 wavename）を;％は、 連続ウェーブレット係数シーク
図
関数imagesc（T、F、ABSを（ coefs））;
SET（GCA、 'YDIR'、 '通常'）
カラーバー;
は、xlabel（ '時間T / S'）。
ylabelの（「周波数F /ヘルツ」）;
タイトル（「周波数図をウェーブレット」）;
周波数ビュー％短時間フーリエ変換
図
スペクトログラム（S、256250256、FS）を、
周波数解析ツールボックス短い％フーリエ変換
F = 0：FS / 2;
TFR = tfrstft（S「）;
TFR TFR =（1：床（長さ（S）/ 2）、:);
図
関数imagesc（T、F、ABS（TFR） ）;
SET（GCA、 'YDIR'、 '通常'）
カラーバー;
は、xlabel（ '時間T / S'）;
ylabelの（ '周波数F /ヘルツ'）;
タイトル（ '周波数ビューSTFT'） ;
