24小时热门版块排行榜    

查看: 1470  |  回复: 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的回帖

射天狼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的回帖
查看全部 5 个回答

信彼南山

木虫 (著名写手)

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

信彼南山

木虫 (著名写手)

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

射天狼yjx

新虫 (初入文坛)

请问楼上如何解决这个函数无法调用的问题?
5楼2011-11-08 17:08:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 03:31 by xxniao123
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 怎么查啊 +5 huang1991js 2026-08-26 5/250 2026-08-26 16:09 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见