24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1240  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

射天狼yjx

新虫 (初入文坛)

[求助] 关于erfinv函数调用错误的问题

在 for 语句中调用erfinv函数,程序如下:

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. 我所编写的程序是根据一个参考程序改写的。整体程序结构与参考程序相同,只是更换了部分函数的表达式。参考程序中也包含上述调用erfinv函数的程序,赋相同初值的情况下,参考程序运行无误,改写程序出现上述错误。
2. 我验算过y_erf的取值区间,当y_erf=-0.99999时,程序报错,但-0.99999这个数值之前已经出现了四次,都计算成功了,在第五次计算时报错。

希望论坛中的各位高人给予指点。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

哦,原来函数保密
4楼2011-11-07 21:05:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

信彼南山

木虫 (著名写手)

Derf_2是什么东西?
2楼2011-11-05 10:28:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

射天狼yjx

新虫 (初入文坛)


臭水沟(金币+1): 请使用“引用回复”,否则对方看不到你的回复~ 2011-11-05 21:59:59
Derf_2是之前一段程序定义的变量。也是d循环迭代的结果。
同样地程序我用MATLAB2010b计算了一下,错误提示如下:

??? Error using ==> erfinv
Input must be real and full.

请问这是什么意思?
3楼2011-11-05 21:57:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

射天狼yjx

新虫 (初入文坛)

请问楼上如何解决这个函数无法调用的问题?
5楼2011-11-08 17:08:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿武汉理工,总分321,英一数二,求老师收留。 +5 nnnnnnn5 2026-03-25 5/250 2026-03-27 04:42 by wxiongid
[考研] 286求调剂 +3 lim0922 2026-03-26 3/150 2026-03-26 23:07 by 不吃魚的貓
[考研] 281求调剂 +3 亚克西good 2026-03-26 5/250 2026-03-26 19:48 by 不吃魚的貓
[考研] 一志愿211 初试270分 求调剂 +6 谷雨上岸 2026-03-23 7/350 2026-03-26 18:55 by 不吃魚的貓
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 调剂推荐 +4 清酒714 2026-03-26 5/250 2026-03-26 17:49 by fmesaito
[考研] 275求调剂 +9 Micky11223 2026-03-25 11/550 2026-03-26 15:54 by 不吃魚的貓
[考研] 315分求调剂 +5 26考研上岸版26 2026-03-26 5/250 2026-03-26 12:11 by laoshidan
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 285求调剂 +3 AZMK 2026-03-24 3/150 2026-03-25 12:23 by userper
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 305分求调剂(食品工程) +5 Sxy112 2026-03-21 7/350 2026-03-24 12:27 by 544594351
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 材料/农业专业,07/08开头均可,过线就行 +3 呵唔哦豁 2026-03-23 4/200 2026-03-23 22:30 by 汪!?!
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-20 3/150 2026-03-21 16:46 by barlinike
[考研] A区线材料学调剂 +5 周周无极 2026-03-20 5/250 2026-03-20 21:33 by laoshidan
信息提示
请填处理意见