| 5 | 1/1 | 返回列表 |
| 查看: 1249 | 回復: 4 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
關于erfinv函數(shù)調(diào)用錯誤的問題
|
||
|
在 for 語句中調(diào)用erfinv函數(shù),程序如下: N=0; X_sum=0.0; Y_sum=0.0; X2_sum=0.0; XY_sum=0.0; for d=0.001:dp_step:dp_max1+0.001; N=N+1; dx=log10(d); X_sum=X_sum+dx; X2_sum=X2_sum+dx*dx; y_erf=2*Derf_2(N)-1; if abs(abs(y_erf)-1)<=0.00001; y_erf=y_erf/abs(y_erf)*0.99999; end dy=erfinv(y_erf); Y_sum=Y_sum+dy; XY_sum=XY_sum+dx*dy; end 運行過程中出錯提示如下 ??? Error using ==> erfinv Y must be real. Error in ==> LC3 at 118 dy=erfinv(y_erf); 需要說明的兩個情況: 1. 我所編寫的程序是根據(jù)一個參考程序改寫的。整體程序結構與參考程序相同,只是更換了部分函數(shù)的表達式。參考程序中也包含上述調(diào)用erfinv函數(shù)的程序,賦相同初值的情況下,參考程序運行無誤,改寫程序出現(xiàn)上述錯誤。 2. 我驗算過y_erf的取值區(qū)間,當y_erf=-0.99999時,程序報錯,但-0.99999這個數(shù)值之前已經(jīng)出現(xiàn)了四次,都計算成功了,在第五次計算時報錯。 希望論壇中的各位高人給予指點。 |
木蟲 (著名寫手)
木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 8/400 |
|
|
[考研] 085602 307分 求調(diào)劑 +5 | 不知道叫什么! 2026-03-26 | 5/250 |
|
|
[考研] 考研化學308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
[考研] 0703化學/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 342求調(diào)劑 +3 | 加油a李zs 2026-03-26 | 3/150 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 材料277求調(diào)劑 +5 | min3 2026-03-24 | 5/250 |
|
|
[考研] 化學調(diào)劑一志愿上海交通大學336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
[考研] 生物技術與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
[考研] 0703化學調(diào)劑,求導師收 +7 | 天天好運來上岸?/a> 2026-03-24 | 7/350 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 291求調(diào)劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|