1樓:李辰正傳
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小波變換的係數和恢復問題
2樓:塗顏帥
cs是個好東西,首先非零個數可以直接用find, length( find(a~=0) ) 就是a中非零元素的個數。
求解1範數有工具包的,l1-magic.
你要得到右圖,第一步需要把小波基寫成矩陣phi,假設要分解的訊號是y, 利用l1magic 求解 y=a*phi*x , a是測量矩陣,如果你只是想用小波分解y,a取1就好了。 得到的x才是稀疏的,否則直接小波分解,得到的係數一般不稀疏
多看看壓縮感知的基礎,l1magic 也可以適當了解他的用法,對你肯定有幫助
matlab 一維連續小波變換原理 50
3樓:
不知是如何個不一樣法? cwt理論積分式你又用滴是哪個,卷積計算有沒有啥處理?每個小波函式滴積分又是咋計算滴?
你問滴這問題可能會有很多涉及到滴其它問題,別想通過這種問答解決,如果你滴基礎很好可以參看matlab滴幫助文件有cwt計算類似差分一樣滴積分做差滴公式,其**你可以解讀一下cwt函式滴**,都是開源滴程式,祝順利!
求matlab 連續小波變換大神!!
4樓:
cwt重構個啥?可以計算得到該尺度cwt變換結果對應的頻率,但不能得到該尺度對應頻率段的訊號,要得到訊號只能用dwt,重構dwt的結果係數可以得到訊號。
請問一維時域訊號的連續小波變換用matlab怎麼做 5
5樓:匿名使用者
1,把第乙個圖的資料在matlab裡的workplace搞出來
用matlab如何與c 連線,用MATLAB如何與C 連線
matlab 是當前應用最為廣泛的數學軟體,具有強大的數值計算 資料分析處理 系統 分析 圖形顯示甚至符號運算等功能 1 利用這一完整的數學平台,使用者可以快速實現十分 複雜的功能,極大地提高工程分析計算的效率 2 3 但與其他高階程式 3 相比,matlab 程式 是一種解釋執行程式,不用編譯等預...
matlab求斜率曲線,求教用matlab求曲線的斜率或切線的方法
斜率即導數。離散資料的話,很簡單,按定義求,後乙個離散資料減去當前乙個 或前乙個,看具體應用 除以取樣間隔,通過這樣來近似當前點的斜率。下面的答覆diff是對表示式進行微分。你有表示式嗎?有就diff,但從你的描述來看,是離散資料,這在數值 時經常遇到。這樣的話不能diff,除非先做曲線擬合。但是擬...
用Matlab畫傾斜的圓平面,matlab裡怎樣畫乙個圓麵,要求是有座標軸的,圓麵裡面填充黑色
主要需要 bai空間圓的引數方程 du啊。syms f1 t f2 t r 6 x linspace sqrt 24 0 linspace 0,sqrt 24 f1 t t sqrt 2 r zhi2 3 t 2 2 f2 t t sqrt 2 r 2 3 t 2 2 y1 eval f1 x z1...