网页资讯视频图片知道文库贴吧地图采购
进入贴吧全吧搜索

 
 
 
日一二三四五六
       
       
       
       
       
       

签到排名:今日本吧第个签到,

本吧因你更精彩,明天继续来努力!

本吧签到人数:0

一键签到
成为超级会员,使用一键签到
一键签到
本月漏签0次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行补签。
连续签到:天  累计签到:天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
05月18日漏签0天
mathematica吧 关注:19,972贴子:74,300
  • 看贴

  • 图片

  • 吧主推荐

  • 游戏

  • 10回复贴,共1页
<<返回mathematica吧
>0< 加载中...

【请指点】数值计算重积分总是报错(两个明确积分限、一个变限)

  • 只看楼主
  • 收藏

  • 回复
  • 微分方程
  • 还未摆脱
    8
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
具体代码放在了楼下


  • 微分方程
  • 还未摆脱
    8
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
数学符号变成了代码:
[CapitalTheta] = m^2 X Cos[ [Sigma]]^2 Cot[ [Theta]]^2;
R[ [Epsilon]_, [Sigma]_, F_] := ( [Epsilon]^2 - 1) r^4 +
2 M r^3 + (a^2 ( [Epsilon]^2 - 1) - X Sin[ [Sigma]]^2 - F) r^2 +
2 M (F + (Sqrt[X] Sin[ [Sigma]] - a [Epsilon])^2) r - a^2 F;


2026-05-18 10:58:10
广告
不感兴趣
开通SVIP免广告
  • 微分方程
  • 还未摆脱
    8
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
para0 = {a -> 0.5, M -> 1, [Theta] -> [Pi]/3};
R1 = R[ [Epsilon], [Sigma], F] /. para0;
DR1 = D[R1, r];
trpl = 2 (M + Sqrt[M^2 - a^2]) /. para0
para = {a -> 0.5, M -> 1, [Beta] -> 1,
r -> trpl, [Theta] -> [Pi]/3}


  • 微分方程
  • 还未摆脱
    8
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
[Sigma]m = 0; [Sigma]M = 2 [Pi];
[Epsilon]m = 1; [Epsilon]M = 100;
Fm = 0.01;
Fc[ [Epsilon]_, [Sigma]_] := (-377.54 + 585.47 [Epsilon]^2 +
2.61 Sin[ [Sigma]]^2 - 2.61 [Epsilon]^2 Sin[ [Sigma]]^2)/(
13.68 + 41.78 Sin[ [Sigma]]^2);


  • 微分方程
  • 还未摆脱
    8
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
n = 100;
core[ [Epsilon]_, [Sigma]_,
F_] := ( [Epsilon] Exp[- [Beta] [Epsilon]])/(
2 Sqrt[R[ [Epsilon], [Sigma], F]] Cos[ [Theta]]) /. para;
core1[ [Epsilon]_, [Sigma]_] :=
NIntegrate[
core[ [Epsilon], [Sigma], F], {F, Fm,
Fc[ [Epsilon], [Sigma]]}] /. para;
core1data =
Flatten[Table[{{ [Epsilon], [Sigma]},
core1[ [Epsilon], [Sigma]]}, { [Epsilon], [Epsilon]m,
[Epsilon]M, ( [Epsilon]M - [Epsilon]m)/
n}, { [Sigma], [Sigma]m, [Sigma]M, ( [Sigma]M - [Sigma]m)/
n}], 1];
corefunc = Interpolation@core1data;
NIntegrate[
corefunc[ [Epsilon], [Sigma]], { [Epsilon], [Epsilon]m,
[Epsilon]M}, { [Sigma], [Sigma]m, [Sigma]M}]


  • 微分方程
  • 还未摆脱
    8
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼


请大佬指点报错部分如何修正处理。
将报错点邻域挖去后会出现新的报错点。


  • xzcyr
  • 吧主
    15
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
……特殊字符前面的 \ 呢?百度贴吧又出bug了还是你自己去掉了?如果是后者,你这样反而把正确的代码弄得没法执行了,那个 \ 是代码的一部分,不信你把 \[Theta] 贴回软件试下。


  • xzcyr
  • 吧主
    15
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
然后,你看到的东西不叫“报错”,叫警告信息,它只是提示出错的可能性。遇到警告时请不要六神无主,请冷静地阅读它到底写了什么。考虑到警告信息里给的误差估计不是很大,你遇到的警告很可能是良性的。(当然这个还是要结合问题的物理背景去判断。)如果NIntegrate计算出的结果精度尚可接受,那么不妨直接无视警告。(嫌烦的话就直接用Off或者Quiet关掉。)


登录百度账号

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!
  • 贴吧页面意见反馈
  • 违规贴吧举报反馈通道
  • 贴吧违规信息处理公示
  • 10回复贴,共1页
<<返回mathematica吧
分享到:
©2026 Baidu贴吧协议|隐私政策|吧主制度|意见反馈|网络谣言警示