 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (5937)
- >導(dǎo)師招生 (1194)
- >考研 (1143)
- >文獻(xiàn)求助 (431)
- >蟲(chóng)友互識(shí) (423)
- >休閑灌水 (190)
- >考博 (151)
- >碩博家園 (132)
- >招聘信息布告欄 (131)
- >基金申請(qǐng) (131)
- >論文投稿 (107)
- >博后之家 (88)
- >綠色求助(高懸賞) (43)
- >SciFinder/Reaxys (36)
- >論文道賀祈福 (32)
- >找工作 (32)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2601 | 回復(fù): 5 | |||
| 【懸賞金幣】回答本帖問(wèn)題，作者dsmj1995將贈(zèng)送您 10 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
請(qǐng)教如何提高M(jìn)ATLAB函數(shù)norminv(p,0,1)的使用精度
|
|||
|
請(qǐng)教各位大佬一個(gè)MATLAB問(wèn)題：如何提高函數(shù)norminv(p,0,1)的使用精度？ norminv為一維正態(tài)分布的累積分布函數(shù)的逆函數(shù)，假設(shè) p = 1-10^(i) , i = 1,...,.N， 因?yàn)閙atlab的運(yùn)算精度為雙精度16位有效數(shù)字 ，因此當(dāng)N>16的時(shí)候，norminv(p,0,1)返回的都是正無(wú)窮大�，F(xiàn)在編程涉及到的p比較小，所以想請(qǐng)教一下如何使得p<1-10^16時(shí)，也能使得該函數(shù)得出較為準(zhǔn)確的結(jié)果？ |
回復(fù)此樓» 猜你喜歡
2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有71人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有5人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有25人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有27人回復(fù)
|
驗(yàn)證代碼更正，需先保存主函數(shù)LeftTruncatedNormrnd： ul = 2 Ns = 1e5; % method#1 % U1 = rand(Ns,1); % X = norminv(U1+(1-U1).*normcdf(ul)); % % % method#2 % U1 = rand(Ns,1); % p = (U1+(1-U1).*(1-sym(normcdf(ul,'upper')))); % X = double(norminv(p)); % % method#3 [X,accept_rate] = LeftTruncatedNormrnd(ul,Ns); close all 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_rate)/dx; pdf_norm = normpdf(xc,0,1)/normcdf(ul,'upper'); plot(xc,pdf_sim) hold on plot(xc,pdf_norm) legend('sim','target') |
6樓2020-09-06 14:48:54
kingspin
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 254 (大學(xué)生)
- 金幣: 4472.5
- 紅花: 74
- 帖子: 1635
- 在線: 315.5小時(shí)
- 蟲(chóng)號(hào): 2585875
- 注冊(cè): 2013-08-07
- 性別: GG
- 專業(yè): 聚合物共混與復(fù)合材料
★
獨(dú)孤神宇: 金幣+1, 鼓勵(lì)交流 2020-09-07 08:20:11
獨(dú)孤神宇: 金幣+1, 鼓勵(lì)交流 2020-09-07 08:20:11
|
做數(shù)值計(jì)算的第一步就應(yīng)該是把問(wèn)題歸一化，從而在計(jì)算質(zhì)量最高的范圍求解 發(fā)自小木蟲(chóng)IOS客戶端 |
2樓2020-08-29 20:16:43
|
該問(wèn)題已經(jīng)解決。實(shí)際上，我本來(lái)的目的是為生成服從左截?cái)嗟恼龖B(tài)分布變量，X~N(mu,sigma,u-), 式中mu,sigma為標(biāo)準(zhǔn)正態(tài)變量的均值與方差，u-為左截?cái)嗟奈恢�。通常的的方法是生成均勻分布變�?p~U(p1,1)，其中p1 = normcdf(u-,mu,sigma)，然后再根據(jù)正態(tài)函數(shù)的累積分布函數(shù)將p1轉(zhuǎn)化目標(biāo)的截?cái)嗾龖B(tài)分布變量X， 即是X = norminv(p1,mu,sigma)。其中normcdf，norminv都是MATLAB內(nèi)置函數(shù)，分別是正態(tài)分布函數(shù)的累積分布函數(shù)與其逆函數(shù)。這種做法在u-與mu比較接近的時(shí)候，效率還是可以的。但是當(dāng)u->>mu，即是p1非常接近于1，如同我在帖子中所描述的問(wèn)題，由于數(shù)值計(jì)算的精度所限制這種做法是行不通的。 關(guān)于當(dāng)p1非常接近于1時(shí)，對(duì)應(yīng)的逆函數(shù)值的求解問(wèn)題，我在鏈接上 https://www.mathworks.com/matlab ... ty-is-close-to-zero 看到兩種解決方法：一種是在MATLAB采用符號(hào)變量sym，一種是通過(guò)迭代的方法。這里不再贅述。下面說(shuō)我的問(wèn)題怎么解決，相應(yīng)的MATLAB代碼如下： Beta = u-；即是上面問(wèn)題描述的截?cái)嗟奈恢?br /> % % 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，提供了一種采用平動(dòng)指數(shù)分布函數(shù)作為抽樣函數(shù)的Accept-Reject模擬方法，根據(jù)此編制了MATLAB函數(shù) LeftTruncatedNormrnd X = LeftTruncatedNormrnd(Beta,Ns);， |
3樓2020-09-06 14:14:19
|
編制的MATLBA函數(shù) LeftTruncatedNormrnd及其驗(yàn)證代碼如下，經(jīng)過(guò)驗(yàn)證感覺(jué)這種方法還不錯(cuò)。有路過(guò)大佬或蟲(chóng)友如果覺(jué)得有什么不對(duì)，或者更高效的方法，還請(qǐ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 驗(yàn)證代碼： 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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +3 | share_joy 2026-03-16 | 3/150 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[文學(xué)芳草園] 伙伴們，祝我生日快樂(lè)吧 +17 | myrtle 2026-03-10 | 26/1300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 妮妮ninicgb 2026-03-15 | 9/450 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂(lè)的香蕉 2026-03-11 | 4/200 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[論文投稿]
投稿問(wèn)題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | 柯淮然 2026-03-10 | 8/400 |
|