| 查看: 956 | 回复: 3 | |||
[求助]
求助：改变变量的数值为什么结果就计算不出来了？？？
|
|
开始： >> x1=1; >> x2=0.1; >> syms y1 y2; >> t=1/2.*(1-sqrt(x1.*x2)/pi.*int(int(exp(-(y1-1)^2.*x1-(y2-1)^2.*x2).*log2(1+exp(-4.*y1.*x1-4.*y2.*x2)),y1,-inf,inf),y2,-inf,inf)); Warning: Explicit integral could not be found. Warning: Explicit integral could not be found. >> digits(3) >> vpa(t) ans = 0.376 将x2从0.1变成0.2就计算不出结果来了 >> x1=1; >> x2=0.2; >> syms y1 y2; >> t=1/2.*(1-sqrt(x1.*x2)/pi.*int(int(exp(-(y1-1)^2.*x1-(y2-1)^2.*x2).*log2(1+exp(-4.*y1.*x1-4.*y2.*x2)),y1,-inf,inf),y2,-inf,inf)); Warning: Explicit integral could not be found. Warning: Explicit integral could not be found. >> digits(3) >> vpa(t) ans = 0.5 - 0.0712*numeric::int(numeric::int(log(1/exp(4*y1 + (4*y2)/5) + 1)/(exp((y1 - 1)^2 + (y2 - 1)^2/5)*log(2)), y1 = -Inf..Inf), y2 = -Inf..Inf) 求高人指点这是什么情况？？？ 该怎么改程序才能让x2为任意值时都能得到具体的结果？？？ 多谢了 |
回复此楼» 猜你喜欢
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
-
垃圾破二本职称评审标准
已经有17人回复
-
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
EST投稿状态问题
已经有7人回复
2楼2012-04-12 09:43:00
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
3楼2012-04-12 09:45:19
4楼2012-04-12 09:56:11