1. 關於matlab如何求協方差的問題~幫助啊...
使用cov函數分別計算協方差。獨立運行時候與同時運行時候得到的協方差可能會不同,取決於A與B是否獨立。不可以計算三組數之間的協方差,只能兩兩之間計算。
2. matlab股票協方差
pp概述
原則上,所有圖象處理都是圖像的變換,而本章所謂的圖象變換特指數字圖象經過某種數學工具的處理,把原先二維空間域中的數據,變換到另外一個"變換域"形式描述的過程。例如,傅立葉變換將時域或空域信號變換成頻域的能量分布描述。
任何圖象信號處理都不同程度改變圖象信號的頻率成分的分布,因此,對信號的頻域(變換域)分析和處理是重要的技術手段,而且,有一些在空間域不容易實現的操作,可以在頻域(變換域)中簡單、方便地完成。
Pp
如上所述,圖象變換是將 維空間圖象數據變換成另外一組基向量空間(通常是正交向量空間)的坐標參數,我們希望這些離散圖象信號坐標參數更集中地代表了圖象中的有效信息,或者是更便於達到某種處理目的。下圖描述了數字圖象處理中空域處理與變換域處理的關系。
pp
圖象變換的實質就是將圖象從一個空間變換到另一個空間,各種變換的不同之處關鍵在於變換的基向量不同。以下給出幾種不同變換基向量的變換示例。
例如,由直角坐標系變化到極坐標系,見下圖
pp
同樣,一幅彩色圖象可以按照某種准則,分解成若干個基本色彩分量圖象的和。
傅立葉變換可以將一維信號從時間域變換到頻率域,例如下圖,一個正弦信號經過傅立葉變換後,得到它的頻率分布零頻(直流分量)和基頻。
一維傅立葉變換的定義:
一維傅立葉反變換定義:
F(u)包含了正弦和餘弦項的無限項的和,u稱為頻率變數,它的每一個值確定了所對應的正弦-餘弦對的頻率。
根據尤拉公式
傅立葉變換系數可以寫成如下式的復數和極坐標形式:
其中:
傅立葉譜(幅值函數)為
相角為
能量譜為
pp
連續二維函數的傅立葉變換對定義
二維函數的傅立葉正變換
二維函數的傅立葉逆變換
二維函數的傅立葉譜
二維函數的傅立葉變換的相角
二維函數的傅立葉變換的能量譜
pp
2離散傅立葉變換
由於實際問題的時間或空間函數的區間是有限的,或者是頻譜有截止頻率。至少在橫坐標超過一定范圍時,函數值已趨於 而可以略去不計。將 和 的有效寬度同樣等分為 個小間隔,對連續傅立葉變換進行近似的數值計算,得到離散的傅立葉變換定義。
其中,一維離散傅立葉正變換
一維離散傅立葉逆變換
pp
二維離散傅立葉變換:對於 圖象
對於 圖象
pp
1.3離散傅立葉變換的性質
性質1:可分離性
二維傅立葉變換可分解成了兩個方向的一維變換順序執行。
pp
性質2:平移性
空間域平移:
頻率域平移:
pp
當 時有:
可以簡單的用 乘以 將 的傅立葉變換的原點移動到相應 頻率方陣的中心。
(圖)
pp
性質3:周期性及共軛對稱性
離散的傅立葉變換和它的反變換具有周期為 的周期性:
傅立葉變換也存在共軛對稱性:
pp
性質4:旋轉性質
平面直角坐標改寫成極坐標形式:
做代換有:
如果 被旋轉 則 被旋轉同一角度。即有傅立葉變換對:
pp
(圖)
性質5:線性性質
如果:
則有:
pp
性質6: 與圖象均值的關系
二維圖象灰度均值定義:
而傅立葉變換變換域原點的頻譜分量:
所以有:
即 數值 倍於圖象灰度均值。
Pp
性質7:圖象拉普拉斯運算元處理後的傅立葉變換
圖象拉普拉斯運算元處理的定義:
則圖象拉普拉斯運算元處理後的傅立葉變換對為:
pp
性質8:卷積與相關定理
卷積定理 一維序列的卷積運算定義為:
當
則有
注意在用傅立葉變換計算卷積時, 由於函數被周期化,為了保證卷積結果正確,計算過程中兩個序列長度N1,N2都要補零加長為N1+ N2-1。二維圖象序列卷積定理的定義和計算過程與一維情況相同。*為卷積符號。
pp
相關定理:
一維、二維兩個離散序列的相關可以寫作
則有相關定理
pp
4快速傅立葉變換
由一維傅立葉變換入手,換一種表示方法
pp
定義:
則:
因為:
pp
傅立葉變換的快速計算示意圖:
(圖)
pp
一維傅立葉變換:
其逆變換為: R
則有:
對於二維情況:
pp
§2離散餘弦變換(DCT)
從第一節內容我們可以看到,傅立葉變換是用無窮區間上的復正弦基函數和信號的內積描述信號中總體頻率分布,或者是將信號向不同頻率變數基函數矢量投影。實際上,基函數可以有其它不同類型,相當於用不同類型基函數去分解信號(圖象)。餘弦變換是其中常用的一種。
pp
設離散序列 ,為一離散序列,根據下式延拓成偶對稱序列 :
其中 。 是關於 為中心的偶對稱序列如下圖所示。
(圖)
pp
以 代入在 范圍內作 點的傅立葉變換:
pp
餘弦變換的變換核為:
表示成矩陣形式為:(其中各列模為1)
pp
定義偶餘弦變換(EDCT)和逆變換為:
pp
二維餘弦變換:
二維餘弦變換具有可分離性:
表示成矩陣形式:
pp
餘弦變換可以利用傅立葉變換實現:
將 延拓為:
則有:
藉助傅立葉變換計算餘弦變換的步驟:
1)把 延拓成 ,長度為 ;
2)求 的 點的FFT;
3)對 各項乘上對應的因子 ;
4)取實部,並乘上因子 ;
5)取 的前 項,即為 的餘弦變換。
Pp
餘弦反變換:
首先延拓 ,
反變換,
pp
§3 正弦變換
一維正弦變換核
一維正弦變換
二維正弦變換核
二維正弦變換
pp
§4 沃爾什-哈達瑪變換
沃爾什-哈德瑪(Walsh-Hadamard)變換的變換核是一類非正弦的正交函數(Walsh函數),例如方波或矩形波。與正弦波頻率相對應,這種非正弦波形可用"列率"(單位時間內波形通過零點數平均值的一半)描述。Walsh函數可以由Rademacher函數構成,Rademacher函數集是一個不完備的正交函數集,Rademacher函數有兩個自變數 和 ,用 表示。
pp
Rademacher函數波形圖和矩陣表示
(圖)
pp
用Rademacher函數構造沃爾什函數:
其中: 表示 所選用的二進制位數
是Rademacher函數
是 的自然二進制的位序反寫後的第 位數字,
例: 用三位二進制碼, ,求
pp
(圖)
pp
Walsh函數的矩陣形式
其變換核矩陣有遞推關系:(直積)
沃爾什-哈達瑪變換定義:
一維沃爾什-哈達瑪變換可表示成矩陣形式:
pp
例:
二維沃爾什-哈達瑪變換:
,
其中 , 階數相同。
pp
例:
另外有
pp
§5 奇異值分解
基於矩陣奇異值分解的二維酉變換:
任何 矩陣 可以分解成:
和 是正交矩陣, 是矩陣 的特徵值。
則:
稱為矩陣 的奇異值。
Pp
對矩陣 作奇異值分解:
令:
則;
是 或 的特徵值所構成的對角陣。
矩陣的特徵矢量
矩陣的特徵矢量
pp
矩陣 可用奇異值分解的級數展開式表示成:
是 矩陣的第 行第 列陣元,即 分別是 矩陣的第 列矢量。
(圖)
pp
§6 K_L變換
K_L變換又稱為Hotelling變換和主成分分析。
當變數之間存在一定的相關關系時,可以通過原始變數的線性組合,構成為數較少的不相關的新變數代替原始變數,而每個新變數都含有盡量多的原始變數的信息。這種處理問題的方法,叫做主成分分析,新變數叫做原始變數的主成分。例如人臉圖象可表示成:
(圖)
pp
主成分分析與線性回歸的比較:
設有 個觀測點 ,散布如圖所示,線性回歸的問題是要找一條對 個點 的擬合直線 ,使偏差平方和最小。
主成分的基本思想是,先對 個點 求出第一條"最佳"擬合直線,使得這 個點到該直線的垂直距離的平方和最小,並稱此直線為第一主成分。然後再求與第一主成分相互獨立(或者說垂直)的,且與 個點 的垂直距離平方和最小的第二主成分。
Pp
(圖)
有 幅圖象 ,大小為 。每幅圖象表示成向量:
向量的協方差矩陣定義為:
其中:
令 和 是 的特徵向量和對應的特徵值。
特徵值按減序排列,
變換矩陣的行為 的特徵值,則變換矩陣為:
對應第 個特徵向量的第 個分量。
K_L變換定義為:
變換後,有:
pp
K_L變換的計算步驟:
1. 求協方差矩陣 ;
2. 求協方差矩陣的特徵值 ;
3. 求相應的特徵向量
4. 用特徵向量 構成變換矩陣 ,求 。
Pp
K-L變換的一種快速演算法:
輸入圖像樣本集合為: ,每一個樣本圖的大小為 , 可以用 維的向量 來表示(即把原圖像按行連到一起構成 維向量)。它也可看作 維空間的一個點,稱此空間為原始圖像空間S。實際上樣本圖像具有較大的相似性的, 因此,全部樣本圖象不會
充滿整個 維的空間,只是會聚集
在圖像空間的一個相對狹小子空間內。
Pp
(圖)
樣本圖象在原始圖像空間中的分布
如果以樣本圖集的總體協方差矩陣為主成分分析的產生矩陣,則所有樣本圖像的總體協方差矩陣為:
式中, , 是全體樣本圖像均值。 滿足下面的方程
是矩陣 的特徵向量, 是對應的特徵值。根據主成分分析理論,得到一個從原始圖像空間到新特徵空間的線性變換 。 是由 的特徵向量構成的變換矩陣。
Pp
但是,直接求矩陣 的特徵值和特徵向量很困難。如果樣本圖象個數 不太多,可以先計算出 維矩陣 的特徵值 和特徵向量 。因為
左乘矩陣 ,得
那麼 就是矩陣 的特徵向量。
Pp
-主成分空間的基。根據主成分分析,可以選擇 個較大特徵值對應的特徵向量(主成分),構造新的 維主成分空間 。每一幅圖象在此空間的投影對應一個 維向量 ,它們就是低維新特徵向量(主成分)。
Pp
小結
傅立葉變換(FFT) 具有快速演算法,數字圖象處理中最常用。需要復數運算。可把整幅圖象的信息很好地用若干個系數來表達。
餘弦變換(DCT) 有快速演算法,只要求實數運算。在相關性圖象的處理中,最接近最佳的K_L變換,在實現編碼和維納濾波時有用。同DFT一樣,可實現很好的信息壓縮。
正弦變換(DST) 比快速DCT快一倍。只需實數運算,可導出快速的K_L變換演算法。在實現編碼和濾波時有用。具有很好的信息壓縮效果好。
沃爾什-哈達瑪變換(WHT) 在數字圖象處理的硬體實現時有用。容易模擬但很難分析。在圖象數據壓縮、濾波、編碼中有應用。信息壓縮效果好。
K_L變換(KLT) 在許多意義下是最佳的。無快速演算法。在進行性能評估和尋找最佳性能時有用。對小規模的向量有用,如彩色多譜或其他特徵向量。對一組圖象集而言,具有均方差意義下最佳的信息壓縮效果。
奇異值分解(SVD) 對任何一幅給定的圖象而言,具有最佳的信息壓縮效果。無快速演算法。設計有限沖激響應(FIR)濾波器時,尋找線性方程的最小范數解時有用。潛在的應用是圖象恢復,能量估計和數據壓縮。
>> A=[0,0,0;2,0,2;]
A =
0 0 0
2 0 2
>> v = diag(cov(A))'
v =
2 0 2 說實在的我也看不懂啦,。
3. Matlab中自協方差函數
xcorr2函數是對矩陣做相關的
自相關就是xcorr2(x,x)這種的,自己跟自己來一下就行了
比如:
a=ones(5,3);
b=xcorr2(a,a);
4. 請教用MATLAB如何做協方差模型分析
在方差分析中協變數必須是連續性變數,否則結果會出現錯誤。不過在你的實驗中,性別應該作為混雜因素來處理,在實驗設計階段可以採用限制、匹配、隨機化的方法以避免其產生混雜作用;如果其混雜作用已經產生,即實驗數據已得出,則只能通過分層分析或多因素分析中的Logistic 回歸分析來解決了。
5. 如何用matlab求協方差矩陣特徵值
先寫出協方差矩陣s,再調用eig(s)這個庫函數,
調用方法:[ev,ed]=eig(s). ed為特徵值矩陣,ev特徵向量矩陣,
排列順序:從低階到高階.
s=[2291.333134019342523.3331245.3332482;
1340956.666715961401.333883.33331480;
193415964281.6671436.66716631945.667;
2523.3331401.3331436.6672984.66712362800.667;
1245.333883.333166312368431343;
248214801945.6672800.66713432729.667]
[ev,ed]=eig(s)
6. matlab 計算協方差報錯
您好!
cov計算UINT8 會超限。
希望對您有幫助!
7. 怎麼求協方差矩陣啊,R,matlab,或者excel都可以。
顯然你這個問題不在於matlab,而是你要先搞清楚這個模型。顯然國債、市場、規模和成長是幾個隨機變數,這些變數最終作用在某一個指標上,你得搞清楚這四個變數的協方差矩陣具體是怎麼計算的,然後再用matlab做計算。
就打個比方,你讓matlab去計算某一個函數f(x),而你自己都不知道這個函數的具體內容,matlab怎麼給出算出來?
8. MATLAB 怎麼計算協方差
>> x=rand(1,5);
>> y=2*rand(1,5);
>> cov(x,y) %計算協方差
ans =
0.1079 -0.0225
-0.0225 0.6148
協方差分析是建立在方差分析和回歸分析基礎之上的一種統計分析方法。 方差分析是從質量因子的角度探討因素不同水平對實驗指標影響的差異。一般說來,質量因子是可以人為控制的。 回歸分析是從數量因子的角度出發,通過建立回歸方程來研究實驗指標與一個(或幾個)因子之間的數量關系。但大多數情況下,數量因子是不可以人為加以控制的。
9. MATLAB 怎麼計算協方差
>> x=rand(1,5);
>> y=2*rand(1,5);
>> cov(x,y) %計算協方差
ans =
0.1079 -0.0225
-0.0225 0.6148