| 查看: 2754 | 回复: 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人回复
-
为什么网上很多人说本周 12号出结果
已经有6人回复
-
听说今天filecode变了
已经有49人回复
-
帮忙看看fileCode
已经有12人回复
-
FileCode能看出啥？
已经有17人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有6人回复
-
filecode
已经有9人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有23人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有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
100