24小时热门版块排行榜    

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

信彼南山

木虫 (著名写手)

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的回帖

信彼南山

木虫 (著名写手)

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

射天狼yjx

新虫 (初入文坛)

请问楼上如何解决这个函数无法调用的问题?
5楼2011-11-08 17:08:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 射天狼yjx 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:03 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 3n8v2C8RimXI 2026-09-13 7/350 2026-09-16 02:35 by DgNGHc3h5tPl
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:26 by zNcEhTcH7ihq
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 40ms4Wlo4umh 2026-09-14 3/150 2026-09-15 19:23 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 15:02 by CgyNCDVNhGVg
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 14:25 by VMDqgnlKimtZ
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:22 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
信息提示
请填处理意见