1樓:
樓下是正解,cwt後用real和imag這兩個函式即可。
用matlab如何做morlet連續復小波分析?
2樓:
標準格式為coefs=cwt(s,32,'cmor','plot');即對訊號s進行最大尺度為32的連續小波變換並作圖,由於使用復morlet小波(cmor)所以就是連續復小波變換。
matlab用連續小波變換後的小波係數的資料怎麼提取
3樓:李辰正傳
close all ; %關閉當前及所有工作區中所有開啟的資料庫,表,和索引,並選擇 1 號工作區.
clc;
clear all ; % 從記憶體中釋放所有的記憶體變數和陣列以及所有使用者自定義的選單欄,選單,和視窗的定義.
%訊號濾波
l1=filter(b,1,c3); %c3 8-12hz濾波
l2=filter(b,1,c4); %c4 8-12hz濾波
figure;
subplot(2,1,1);plot(1:1152,c3(:,2),'r',1:1152,l1(:,2),'k');
xlabel('時間');ylabel('幅度');legend('c3 濾波前','c3濾波後');
subplot(2,1,2);plot(1:1152,c4(:,2),'r',1:1152,l2(:,2),'k');
xlabel('時間');ylabel('幅度');legend('c4 濾波前','c4濾波後');
title('8-12hz濾波後的結果');
c3_shang=-(l1.^2).*log(l1.^2);
c4_shang=-(l2.^2).*log(l1.^2);
figure;
plot(1:1152,c3_shang(:,1),'r',1:1152,c3_shang(:,2),'k');
legend('左手c3','右手c3');
title('左手c3的熵和右手c3的熵的圖');
figure;
plot(1:1152,c3_shang(:,1),'r',1:1152,c4_shang(:,1),'k');
legend('左手c3','左手c4');
title('左手c3的熵和左手c4的熵的圖');
%由上圖可以發現左手c3和右手c3的熵在600-1000點之間差別較大,因此在下面的分類中,選用600-1000之間的點進行分類
y_train=(y_train-1.5)*2;
model= svmtrain(y_train(1:90),c3(600:1000,1:90)','-s 0 -t 2 -c 1 -g 0.1');
[predictlabel,accuracy] = svmpredict(y_train(91:140),c3(600:1000,91:140)',model)
在matlab中cwt函式可以用mexihat小波,morlet小波嗎
求matlab的morlet小波分析**,我有近50年的風速資料,求大神啊
哪位好人懂matlab,幫忙看看這幅圖,為什麼5年以下尺度上的週期在小波實部圖上體現不出來呢?**如下:
4樓:
**沒啥明顯問題,出現問題的可能1.原始訊號中可能本身就沒有尺度5以下的資訊,這和你原始訊號的特徵還有采樣多少有關。2.
使用的contourf函式,在matlab中所有與二維網格化有關的問題都解決的不好,contourf函式對實際資料繪製的等值線圖有時就是慘不忍睹,對於數值較小或較大的不連續區域可能繪不出等值線,而尺度5以下小波係數分佈都比較細碎或數值比較其他區域相差較大,這樣就繪不出等值線嘍!
另外,還想對這種週期分析談一些個人看法。對於cwt中的尺度是與所用小波函式的數學性質有關的,與所用的小波函式的中心頻率有關,根本不是上面的這種理解。連續小波變換對尺度的定義根本不是時間尺度,它是沒有量綱的(還扯什麼年),上面的程式沒錯,但對小波分析的理解是錯誤的。
cwt的尺度到底與時間如何對應要通過所用的小波函式的中心頻率、取樣點數和經歷時間計算(就是取樣頻率)。一直都見到有人這麼搞,其實一直都是錯的,沒辦法審稿人也只是混飯吃的棒槌,搞得錯誤文章滿天飛,禍害後來人吶!
怎麼用matlab實現小波變換???急!!!
MATLAB中用Simulink測電壓,怎麼不是正弦波形,而
可能是解算器選擇有問題 也就是simulation中的configuration parameters中solver選第一個 離散 再看看 如果還不行,那雙擊powergui 彈出一個選單 頂上不選中continuous 而選中第二個 然後close 這下就正常了 置於為什麼,你每 一次 系統應該顯...
在matlab中,用兩種方法求函式sinx x cos2x
注意matlab裡表示式的格式,利用syms定義符號,然後利用diff函式求符號導數 y diff x,n 求函式x的n階導數matlab裡執行如下 syms x y dyy sin x x 2 4 x 3 dy diff y,4 執行結果有點長,這是因為你的函式求導表示式比較複雜。做個測試,希望有...
關於小matlab程式設計的小波變換
呵呵,這個程式是告訴你用morlet小波進行連續小波變換時的程式,相當於是內matlab中的底層程式,容function wt,freqbins,scales cwt morlet sig,winlen,nlevel 這個是定義了這樣的一個函式,其中sig,winlen,nlevel這幾個變數為輸入...