| 查看: 2402 | 回复: 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 ] |
回复此楼» 猜你喜欢
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有7人回复
-
申请2026年博士
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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