 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2545 | 回復(fù): 5 | |||
| 【懸賞金幣】回答本帖問題，作者dsmj1995將贈送您 10 個金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
請教如何提高MATLAB函數(shù)norminv(p,0,1)的使用精度
|
|||
|
請教各位大佬一個MATLAB問題：如何提高函數(shù)norminv(p,0,1)的使用精度？ norminv為一維正態(tài)分布的累積分布函數(shù)的逆函數(shù)，假設(shè) p = 1-10^(i) , i = 1,...,.N， 因為matlab的運算精度為雙精度16位有效數(shù)字 ，因此當(dāng)N>16的時候，norminv(p,0,1)返回的都是正無窮大�，F(xiàn)在編程涉及到的p比較小，所以想請教一下如何使得p<1-10^16時，也能使得該函數(shù)得出較為準(zhǔn)確的結(jié)果？ |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有96人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有5人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有2人回復(fù)
|
編制的MATLBA函數(shù) LeftTruncatedNormrnd及其驗證代碼如下，經(jīng)過驗證感覺這種方法還不錯。有路過大佬或蟲友如果覺得有什么不對，或者更高效的方法，還請多指教。 主函數(shù)： function [X,accept_rate] = LeftTruncatedNormrnd(ul,Ns) alpha = (ul+sqrt(ul^2+4))/2; Ns_accept = 0; while Ns_accept==0 % 1. Generate translated exponential distribution Z = exprnd(1/alpha,Ns,1)+ul; % 2. Compute coeffcient of accept rate po = exp(-(Z-alpha).^2/2); % 3. Accept or reject U = rand(Ns,1); index = find(U<=po); X = Z(index,1); % 4. other Ns_accept = length(index); accept_rate = Ns_accept/Ns; end return 驗證代碼： xc_max = max(X); xc = ul  xc_max-ul)/1000:xc_max;dx = xc(2)-xc(1); n1 = hist(X,xc); pdf_sim = n1/Ns_accept/dx; pdf_norm = normpdf(xc,0,1)/normcdf(ul,'upper'); plot(xc,pdf_sim) hold on plot(xc,pdf_norm) |
4樓2020-09-06 14:18:34
kingspin
木蟲 (著名寫手)
- 應(yīng)助: 254 (大學(xué)生)
- 金幣: 4472.5
- 紅花: 74
- 帖子: 1635
- 在線: 315.5小時
- 蟲號: 2585875
- 注冊: 2013-08-07
- 性別: GG
- 專業(yè): 聚合物共混與復(fù)合材料
★
獨孤神宇: 金幣+1, 鼓勵交流 2020-09-07 08:20:11
獨孤神宇: 金幣+1, 鼓勵交流 2020-09-07 08:20:11
|
做數(shù)值計算的第一步就應(yīng)該是把問題歸一化，從而在計算質(zhì)量最高的范圍求解 發(fā)自小木蟲IOS客戶端 |
2樓2020-08-29 20:16:43
|
該問題已經(jīng)解決。實際上，我本來的目的是為生成服從左截斷的正態(tài)分布變量，X~N(mu,sigma,u-), 式中mu,sigma為標(biāo)準(zhǔn)正態(tài)變量的均值與方差，u-為左截斷的位置。通常的的方法是生成均勻分布變量 p~U(p1,1)，其中p1 = normcdf(u-,mu,sigma)，然后再根據(jù)正態(tài)函數(shù)的累積分布函數(shù)將p1轉(zhuǎn)化目標(biāo)的截斷正態(tài)分布變量X， 即是X = norminv(p1,mu,sigma)。其中normcdf，norminv都是MATLAB內(nèi)置函數(shù)，分別是正態(tài)分布函數(shù)的累積分布函數(shù)與其逆函數(shù)。這種做法在u-與mu比較接近的時候，效率還是可以的。但是當(dāng)u->>mu，即是p1非常接近于1，如同我在帖子中所描述的問題，由于數(shù)值計算的精度所限制這種做法是行不通的。 關(guān)于當(dāng)p1非常接近于1時，對應(yīng)的逆函數(shù)值的求解問題，我在鏈接上 https://www.mathworks.com/matlab ... ty-is-close-to-zero 看到兩種解決方法：一種是在MATLAB采用符號變量sym，一種是通過迭代的方法。這里不再贅述。下面說我的問題怎么解決，相應(yīng)的MATLAB代碼如下： Beta = u-；即是上面問題描述的截斷的位置 % % method#1：經(jīng)典方法 U1 = rand(1); X = norminv(U1+(1-U1).*normcdf(Beta)); % % method#2： 采用matlab的sym變量 U1 = rand(1); p = (U1+(1-U1).*(1-sym(normcdf(Beta,'upper')))); X = double(norminv(p)); % % method#3：文獻(xiàn)鏈接https://arxiv.org/pdf/0907.4010.pdf，提供了一種采用平動指數(shù)分布函數(shù)作為抽樣函數(shù)的Accept-Reject模擬方法，根據(jù)此編制了MATLAB函數(shù) LeftTruncatedNormrnd X = LeftTruncatedNormrnd(Beta,Ns);， |
3樓2020-09-06 14:14:19
5樓2020-09-06 14:20:39
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料277分求調(diào)劑 +7 | 飯飯星球 2026-03-04 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +8 | 土土小蟲 2026-03-03 | 10/500 |
|
|
|
[考研] 320材料與化工，求調(diào)劑 +6 | 鶴遨予卿 2026-03-04 | 8/400 |
|
|
[考研] 267化工調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +5 | 青科11 2026-03-02 | 5/250 |
|
|
|
[考研] 成績276，專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
|
[考研] 0703化學(xué) 一志愿東北大學(xué) 279分求調(diào)劑 +4 | 聶樂平 2026-03-03 | 4/200 |
|
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
|
[考研] 0856材料工程，初試313調(diào)劑 +7 | 賣個關(guān)子吧 2026-03-03 | 7/350 |
|
|
|
[考研] 268求調(diào)劑 +10 | 簡單點0 2026-03-02 | 14/700 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
|
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|