| 查看: 1182 | 回复: 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为任意值时都能得到具体的结果？？？ 多谢了 |
回复此楼» 猜你喜欢
时间戳又变了
已经有16人回复
-
科研孤儿太难了
已经有9人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
今天系统多次维护，明天很可能放榜！
已经有10人回复
-
放榜前的不淡定
已经有12人回复
-
言之凿凿，明天官网会通知放榜？
已经有16人回复
-
重要消息，中午系统在维护
已经有12人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有20人回复
-
今天基金会出结果吗？20260819
已经有15人回复
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
