| 查看: 2765 | 回复: 6 | |||
lqglz木虫 (正式写手)
|
[求助]
Matlab(或者Mathematica) 如何进行复杂函数的二重积分？
|
|
变量两个：theta和t 函数如下：f(theta, t)=(2.391e-7)./(0.00105+t).^4.83.*exp(-4.*((theta-10.94).^2+3.362e-3)./(0.00105+t))./(3.8825e-4+t).^5.*1.021e-4.*exp(-0.5.*((theta-10.91).^2+2.381e-3)./(3.8825e-4+t)), 公式也可见附件1图片。 需要对式中t和theta分别在（0, inf）和（-inf, inf）进行二重积分.  1.jpg [ Last edited by lqglz on 2012-10-23 at 22:11 ] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有49人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有21人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有5人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有8人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有7人回复
-
filecode
已经有10人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab高手帮忙做个自定义复杂函数拟合，含有复数，
已经有25人回复
- MATLAB函数速查手册
已经有211人回复
- matlab的函数调用问题
已经有3人回复
- 已知一随机信号的功率谱密度函数，如何利用Matlab产生其单次时域样本
已经有6人回复
- 问个关于matlab里find函数的问题！
已经有4人回复
- matlab做函数图象的问题
已经有6人回复
- 如何用MATLAB拟合指数函数
已经有3人回复
- 【求助】修改MATLAB支持向量机核函数程序遇到的疑问？求助
已经有5人回复
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
2楼2012-10-27 01:39:26
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
3楼2012-10-27 01:40:49
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
4楼2012-10-27 01:41:58
walk1997
金虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4676.2
- 红花: 22
- 帖子: 1066
- 在线: 798.1小时
- 虫号: 416039
- 注册: 2007-06-29
- 性别: GG
- 专业: 粒子物理学和场论
|
这个积分有点意思..... 我觉得估计你最好需要表达式里面有些数值更高精度的，我假设你里面的数值都是解析的 无穷精度 （不过检查了下，用小数 算出来结果一样） 最后2个输出结果应该是对了.... 1.74992*10^9 没有求倒数 --------------------------------------------------------------------------- f[theta_, t_] := (2391*10^-10)/(105/100000 + t)^(483/100)* Exp[-4*((theta - 1094/100)^2 + 3362*10^-6)/(105/100000 + t)]/(38825*10^-8 + t)^5*1021*10^-7* Exp[-1/2*((theta - 1091/100)^2 + 2381*10^-6)/(38825*10^-8 + t)]; f[1, 10.91] NIntegrate[f[theta, t], {theta, -Infinity, Infinity}, {t, 0, Infinity}] NIntegrate[f[theta, t], {theta, -100, 100}, {t, 0, 10}] NIntegrate[ f[theta, t], {theta, -Infinity, -10, -5, -1, 0, 1, 5, 1091/100, 20, Infinity}, {t, 0, 1, 5, 10, 20, 100, Infinity}, PrecisionGoal -> 12, MaxRecursion -> 100] NIntegrate[ f[theta, t], {theta, -Infinity, 1091/100, 1094/100, Infinity}, {t, 0, Infinity}] －------------------------------- 3.15879*10^-39 2.22875*10^-19 2.33261*10^-20 1.74992*10^9 1.74992*10^9 |
5楼2012-10-27 02:23:06
lqglz
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 5451.8
- 散金: 50
- 红花: 1
- 帖子: 946
- 在线: 197.3小时
- 虫号: 754370
- 注册: 2009-04-21
- 性别: GG
- 专业: 岩土力学与岩土工程
6楼2013-04-12 22:27:27
lqglz
木虫 (正式写手)
- 应助: 8 (幼儿园)
- 金币: 5451.8
- 散金: 50
- 红花: 1
- 帖子: 946
- 在线: 197.3小时
- 虫号: 754370
- 注册: 2009-04-21
- 性别: GG
- 专业: 岩土力学与岩土工程
7楼2013-04-12 22:27:52
