| 查看: 1471 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
关于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这个数值之前已经出现了四次，都计算成功了，在第五次计算时报错。 希望论坛中的各位高人给予指点。 |
回复此楼» 猜你喜欢
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有4人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有7人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
-
让我中一个面上吧！
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
我不理解！
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 不知道怎么回答审稿人的问题，求助大家,(关于分子筛MCM-41的）
已经有10人回复
- 请教各位关于研究计划的一些问题？在线等。。。
已经有10人回复
- 有关几个有机的问题（对应LG WADE第五版）
已经有35人回复
- origin 8.0 对数函数拟合
已经有13人回复
- matlab的函数调用问题
已经有3人回复
- 关于支持向量机与核函数的问题
已经有9人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- 问一个关于函数调用的小问题
已经有3人回复
- 请教如何调用Module里面的函数，谢谢
已经有6人回复
- 菜鸟问题：VC++调用函数关闭子窗口的问题
已经有4人回复
- 关于LANL2DZ加极化函数的问题
已经有15人回复
- ★★★关于三甲基氯硅烷保护羟基的问题★★★
已经有4人回复
- 请教关于origin中自定义函数非线性回归拟合问题！
已经有10人回复
3楼2011-11-05 21:57:25
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2011-11-05 10:28:09
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
4楼2011-11-07 21:05:42
5楼2011-11-08 17:08:42
投票:
