如何模擬二元高斯正態分布數據

⑴ matlab如何產生二維正態分布隨機數

實現的方法和詳細的操作步驟如下：

1、第一步，基於實驗獲得一條數據，要確定是否服從嘩豎分布，將其轉換為直方圖，一些數據如圖所示，轉到下面的步驟。

⑵ 如何實現二維高斯擬合 matlab

你可以參考這個例子。用matlab 將一組數據進行正態分布擬合，你可以在基礎上修改，畢備實現過程如下友叢：

x=[。。。];y=[。。。];z=[。。。];

x=[xy];y=z;

func=@（手告毀a,x）二維高斯表達式

a0=[0 0 0 0];

[a,r] = nlinfit(x,y,func,a0) %函數擬合得到其正態分布參數，

式中：μ1=a（1）,σ1=a（2）,μ2=a（3）,σ2=a（4）

r是殘值（即z與z1的差值）

⑶ 怎麼用MATLAB產生2維或者多維的高斯分布數據

x=randn（m，n）就是二維的高斯分布函數，例如：
x=randn（5，6）便產生5行6列的二維的 高斯分布函數：
ans =

-0.4326 1.1909 -0.1867 0.1139 0.2944 0.8580
-1.6656 1.1892 0.7258 1.0668 -1.3362 1.2540
0.1253 -0.0376 -0.5883 0.0593 0.7143 -1.5937
0.2877 0.3273 2.1832 -0.0956 1.6236 -1.4410
-1.1465 0.1746 -0.1364 -0.8323 -0.6918 0.5711
三維的：randn（m,n,p);依此類推。
例如randn（2，畢侍敏3，4）：
ans(:,:,1) =

-0.3999 0.8156 1.2902
0.6900 0.7119 0.6686

ans(:,:,2) =

1.1908 -0.0198 -1.6041
-1.2025 -0.1567 0.2573

ans(:,:,3) =

-1.0565 -0.8051 0.2193
1.4151 0.5287 -0.9219

ans(:,:,4) =

-2.1707 -1.0106 0.5077
-0.0592 0.6145 1.6924

當然也可以自己構造兩個一維的數據，再自己把之相乘，合成一個手枝二維的數據，但是最簡單的方法，還是直接用函數randn（），
給出的是均值為0，方差為1的高談大斯正態函數的分布的數值。

⑷ 用matlab進行蒙特卡洛模擬，模擬服從二維正態分布，求高手幫忙，十分感謝

用randn()可以生成高斯分布的隨機數。
不過只有先生成隨喚梁機數才有方差和均值，反過來可難了。

也只能使均值和方差近似等於0和6，
用下面方法：
x=randn(1,100)*sqrt(6);
while abs(mean(x))>=0.01 | abs(var(x)-6)>=0.01
x=randn(1,100)*sqrt(6);
end

這樣生成的隨機數列x，平均值約等於0，方差越等於6，誤差不超過櫻凱0.01。
當然你也可以把精度調得更高一些，但很可能很久都找不到合適的隨機數。

mean()是求平均值，var()是求方差，
取隨機數之後乘以sqrt(6)是因為randn()取出的隨機數方差在1左右，所以要放大一下。

還有隨機數的個數我這里是100個，你可以脊鏈喚隨意修改。

⑸ 將數據擬合二維高斯分布並輸出擬合後的高斯分布參數

1、檢查一下x的數據輸入是否正確，應按下列方式輸入

x=[1.1,2.1; 3.1,4.1; 。。。。。。];

2、調整a0的值，可以先悉亂冊用rand（）函數作為睜宏a0的初值，即陪蠢
a0=rand(1,5)
3、檢查一下x(:,1)和x(:,2)中，是否用的是冒號

⑹ matlab中二元正態分布函數

二元正態分布又名高斯分布（英語：Gaussian distribution, 採用德國數學家卡爾·弗里德里希·高斯的名字冠名），是一個在數學、物理及工程等領域都非常重要舉棚的概率分布，由於這個分布函數具有很多非常漂亮的性質，使得其在諸多涉及統計科學離散科學等領域的許多方面都有著重大的影響力。

比如圖像處理中最常用的濾波器類型為Gaussian濾波器（也就是茄咐所顫答純謂的正態分布函數）。

⑺ (三)高斯變數的條件模擬

1.模擬的范圍

為使模擬區域的3個主方向與變差函數的3個主方向一致,有必要按前述方向將礦區的坐標系統按順時針方向旋轉10°,鑒於計算機硬碟空間和計算機時的限制,我們模擬的區域只限制在1538m水平標高至1588m水平標高之間,在沿走向方向上的長度為425m,垂直於礦體走向方向上的寬度為120m的長方體空間范圍內,在此范圍之內包含了3層探礦坑道,即PD1538,PD1554及PD1584,條件模擬的網格設置為沿走向方向節點間隔為2.5m,垂直於礦體走向方向上的節點間隔為1m,在鉛垂方向上為2.5m,這樣總共模擬的結點數為408000個,這主要是考慮到一方面采樣的密度,另一方面前人的工作經驗,uLsfer於1985年提出如果結點間隔大於變差函數變程的1/4,那麼變差函數的變程就很難驗證要模擬的數據。在進行條件模擬過程中,估計鄰域的確定必須要考慮到變差函數的變程大小,同時還考慮了在每一個模擬點,至少有兩個信息樣在估計鄰域內,所以,在3個方向上估計鄰域的半徑分別設置為在主方向為40m,走向方向為20m,垂直礦體走向方向為8m。

2.條件模擬的結果

對於任何一個模擬來說,隨機數發生器的根必須取一個較大的值,對於不同的根就會產生不同的現實,另外,轉向帶的數目也是模擬的關鍵參數之一。本次研究中它們分別設置為125793(軟體設置的缺少值)和1000(通常應設置為大於40的正整數)。

眾所周知,一個模擬只是區域化變數變異性的一個現實,如果我們進行了大量的模擬,則它們的結果雖各不相同,但總有一些代表著礦體的真實變化性,這樣地質學家或采礦工程師就可以結合自己對礦床的認識從中選取一個能真正反映礦體內變化特徵的模擬結果。鑒於時間因素,本次研究中,只進行了3個條件模擬,而且所使用的參數都相同,其結果之間存在著一定差異,下表就給出了它們的統計結果。

坑道樣品條件模擬的結果統計表

相應的實驗變差函數曲線見下圖,從圖中可以看出這3組條件模擬的實驗變差函數在各方向上表現出既有相同之處又有不同之處,首先在水平面內,模擬1與模擬2的結果十分相近,而模擬1的變程略有增大,在沿走向方向上,模擬2和模擬3幾乎完全一致,而模擬1的基數略有下降,而在沿傾向方向上,模擬1和模擬2的結果十分吻合,而模擬3卻表現出嚴格遞增的趨勢。

通過與信息樣的實驗變差函數曲線的比較可以發現,模擬1和模擬2的實驗曲線都在一定程度上偏離了原始樣品的曲線,只有模擬3的實驗曲線與原始樣品的曲線較吻合,雖然在沿傾向方向上的實驗曲線的變程較短,但主要是由於模擬區在該方向的范圍較小,無法達到變程的要求,所以,在實驗曲線上表現為嚴格遞增區。

所以單純從變差函數的角度考慮,我們認為模擬較好地反映了礦脈在1538m標高及1584m標高之間的變異性。

(3)條件模擬結果與坑道樣品揭示的礦體幾何形態的比較

在地質勘探工作中,沒有任何一種方法比坑道更能准確地揭示礦體的邊界、確定礦體的連續性。因而,與坑道樣品所揭示的結果進行比較,將會更有力地說明條件模擬結果的優劣。下面分別繪出的1538m標高,1554m標高及1584m標高處的條件模擬平面圖及坑道采樣平面圖。從圖中可以看出:

1)條件模擬所揭示的礦體密度遠遠大於坑道樣品揭示的礦脈寬度,但是在脈體內部出現大量的「空白點」(模擬值很低),說明並不是有黑影的位置上部是礦體,這需要地質人員根據野外勘查結果加以區分,並且使礦體連續化。

2)在沿礦體的走向方向上模擬所揭示的礦脈連續性較差,尤其是1538m標高上,在礦脈的北段表現極為零亂,而野外工作人員卻能根據自己的認識連續地將礦體往南延續,所以說條件模擬並不能加入人為因素而只是一種統計意義上的結果。

3)在1554m標高上模擬所揭示的礦脈長度將遠大於坑道所揭示的礦脈長度,在模擬區域的北半部的坑道內沒有採集任何樣品,而在其下部的1538m坑道內及上部1584m標高的坑道內均有樣品揭示有礦脈存在,根據統計意義,在其中部也應有礦脈存在,所以這種差異可能由於礦脈變化太復雜,在一個很短的距離范圍內礦脈突然尖滅,也可能由於坑道遺漏了確定存在的礦脈所以有待於進一步驗證。

條件模擬變數實驗變差函數曲線圖

4)通過對3個模擬結果與坑道樣品比較發現,無論是礦體連續性上還是礦體的空間位置與形態上模擬3均優於模擬1和模擬2。

⑻ 高斯正態分布

正態分布，一般都只會講公式，怎麼證明的就不提了。我遇到這傢伙有兩個地方：一個是高中數學課上；一個是本科《誤差理論和測量平差》課上。找點資料，想自己推到一下如何得到高斯正態分布的公式。

正態分布（Normal distribution），也稱「常態分布」，又名高斯分布（Gaussian Distribution）。我習慣合起來叫 高斯正態分布 。（剛好和linux Distribution： linux發行版 的英文單詞一樣）

期望（平均數）：μ
標准差 ，
方差 為。
當 和 時候稱為： 標准正態分布 。

matlab繪制正態分布概率密度函數圖像的命令為normpdf，normpdf函數的調用格式為normpdf(x,mu,sigma),其中mu為0，sigma為1時，為標准正態分布悄蠢。

在高斯分布中有三個數學符號，先來解釋這個三個數學符號的含義，然後再說明這個公式的推導思路和推導方法。
三個符號 在數學上分別叫做平均值（又稱數學期望），標准差，自然數。即：
平均值（又稱數學期望）：
標准差:
自然數:

對於數據:

平均數：

語言解釋： 平均數 就是所有數加起來的和除以數據個數n。
數學的含義是：數據中間位置的具體數值。

詳細說明方差方差和標准差之前，先復習一下關於 勾股定理 （在西方又稱 畢達哥拉斯定理 ）和 平面兩點間距離公式 。
在直角三角形中，對於邊長a，b，c有如下關系：

即

在平面坐標系x-o-y下對任意兩點 間的距離D有:

通過勾股定理和平面兩點間距離公式可以看出，型如
表示的含義為兩個之間的距離。數值越小，證明兩個之間越近。

一組數據，平均數是這個數據的中心，那麼就可以用其他數據到平均數的距離來衡量數據和平均數的遠近關系。即這組數據是聚攏一些呢，還是分散一些呢。
方差

因為距離D是需要開方的，所以 方差的含義是距離的平方 。對開方後的方差稱為標准差 。

假設有兩組數據：

說明兩組數據的中間值數值一樣，且都為零。平均值可以諒解為含悉此數組中的中心位置。

即 說明：
A組數據之間的距離較小，數據較聚攏；
B組數據之間的距離較大，數據較分散；

從 歐拉公式看出，把字母e定義成自然數，和歐拉是有直接關系的。倒不太相信網路里說的 歐拉選這個字母的原因，不太可能是因為這是他自己名字Euler的首字母，因為他是個很謙虛的人，總是恰當地肯定他人的工作。

其實從這個公式還是不太能看出來e=2.71828，一開始誰會想到這個式子就極限就是自然數e呢。

但我們可以從對數和指數的關系來聯系e是怎麼來的。

利用對數運算性質中的化平方為相乘的特性，我們知道自然數在對數運算中是最常用的底數。

對於對數運算：
對於指數求導
那麼如果 就好了，a等於多少，才會使得 呢？
恰巧a等於自然數e的時候，lne=1.
於是，可以將a=e帶入指數求導公式：

對函數求導後依舊是其本身，這是一個很好的性質。

е主啟老陪要出現在涉及增長的地方，比如說經濟增長、人口增長、放射性衰變等，可以說е代表了自然率之美。
比如某個市人口為120萬人，每年的人口增長率為20%:
一年後人口：100萬+100萬x20%=100萬（1+20%）=120萬
兩年後人口：120萬+120萬x20%=120萬（1+20%）=100萬（1+20%）（1+20%）=

三年後人口：=
四年後人口：=
X年後人口：=

當人口增長率不可能一直保持20%，因為生存空間有限，增長率應該是隨著時間而降低的。假設增長率和時間X成反比，即增長率為

那麼上述人口增長的數學模型可以抽象為：

當我們想知道很多年後的人口增長，即時間X趨向無窮 的人口時候即可得極限：

為什麼正態分布如此常見
為什麼數據科學家都鍾情於最常見的正態分布？