| 查看: 2437 | 回复: 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 ] |
回复此楼» 猜你喜欢
上海工程技术大学张培磊教授团队招收博士生
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有5人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 专业: 粒子物理学和场论
4楼2012-10-27 01:41:58
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
- 专业: 粒子物理学和场论
|
这个积分有点意思..... 我觉得估计你最好需要表达式里面有些数值更高精度的，我假设你里面的数值都是解析的 无穷精度 （不过检查了下，用小数 算出来结果一样） 最后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