| 查看: 1189 | 回复: 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为任意值时都能得到具体的结果？？？ 多谢了 |
回复此楼» 猜你喜欢
看来今天不会放榜了？
已经有9人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有18人回复
-
感觉是下周放榜了
已经有12人回复
-
我面上完蛋了
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有25人回复
-
什么时候开奖？
已经有7人回复
-
时间戳又变了
已经有18人回复
-
时间戳又变了8-15
已经有27人回复
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
2楼2012-04-12 09:43:00
4楼2012-04-12 09:56:11
