24小时热门版块排行榜    

查看: 2260  |  回复: 9

beefly

专家顾问 (职业作家)

地沟油冶炼专家

[交流] 物理学中尚未解决的积分,数学家们有啥办法? 已有6人参与

在固体物理中会涉及以下形式的三重积分,其中被积函数是一个一般的二次多项式,积分范围也涉及二次多项式。想得到积分的严格表达式,或者高精度的数值积分结果。

目前只有二维情况(去掉维度z)的解析积分是知道的,但是尚未见到三重积分的解析表达式,在实际计算中遇到这类问题,只能做近似,例如简化g(x,y,z),把积分区域均匀分割成8块,每一块用一个线性多项式g'(x,y,z)=a1'+a2'*x+a3'*y+a4'*z去拟合g(x,y,z),但是在计算金属、半导体的时候,这么做会遇到问题。

本人也做过一些不成功的尝试。例如,把g(x,y,z)替换为阶梯函数,乘到f(x,y,z)上去,这样积分范围简单了,但是被积函数太复杂,无法进行严格的积分。于是改做数值积分,用了从5阶到8阶的高斯积分,结果均不理想,还不如上面g(x,y,z)线性近似的方法准。

物理学中尚未解决的积分,数学家们有啥办法?
001.png
回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

至尊木虫 (职业作家)

科苑小木虫


小木虫: 金币+0.5, 给个红包,谢谢回帖
主要难度在于,g是四维的,画不出来,因此积分线不好定
板凳要做十年冷文章不发一个字
2楼2016-04-28 16:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enlearner

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这个问题其实不难,我说一下框架吧。(其实这个问题在代数学中已经解决了关键的部分)
解决这个问题的关键是积分区域g(x,y,z)的规则表达,也就是代数学中的“二次型的标准化,具体做法如下:
第一步:通过坐标旋转把交叉项(也就是xy,xz,yz项)去掉,一次项不用管,自然累加就可以了。三维的旋转角有两个,计算有一定有难度,但一定能算出来,而且后面还要用---具体的内容百度或者查(二次型或二次曲面方面的数学书)
第二步:通过配方法将上一步所得方程化为标准型,结果不外乎三种:椭球(包含球),双曲型,抛物型。(如果是虚球,刚此积分无意义,不与考虑)
第三步:将前面0<x<1,0<y<1,0<z<1在旋转的情况下算出新的取值范围(只与旋转角有关)
第四步:将被积函数f(x,y,z)在旋转下用新的未知量替换,积分变量同步替换。
第五步:最标准的三重积分解法,本科二年级难度。

补充:算法二:
第一步:同上。
第二步:用球坐标变换或柱坐标变换将被积函数和积分区域同步变换。以下内容相信本科的任何一本高数书上都有。
看起来困难,其实最关键的就是第一步---二次型的标准化。两种算法都要用到这个。但这个内容已经是古典数学中较为成熟的东西,早已经解决。如果第一步解决,其余步骤只是顺其自然。
希望对你有所帮助。
期待交流
3楼2016-04-29 10:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujm2000

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
g(x,y,z)的定义有问题

发自小木虫Android客户端
暴走族
4楼2016-04-29 16:05:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

花痕逸韵

铁虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
题目有不严谨的地方,首先你不能确定它围成的空间是不是可以积分,这样就像要通解是不是有点过早。其次,如果参数的不确定,也可能导致积分不出来,比如混沌出现。

发自小木虫Android客户端
5楼2016-04-29 16:23:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
3楼: Originally posted by enlearner at 2016-04-29 10:41:05
这个问题其实不难,我说一下框架吧。(其实这个问题在代数学中已经解决了关键的部分)
解决这个问题的关键是积分区域g(x,y,z)的规则表达,也就是代数学中的“二次型的标准化,具体做法如下:
第一步:通过坐标旋转 ...

第一步有人做过,如J. Phys.: Condens. Matter. 3, 6721-6742 (1991),其中最后三种一次曲面已经解决。不过原文用的是四面体+二次曲面g(x,y,z),再往下就做不下去了。我想改用立方体+二次曲面g(x,y,z),原以为能简化,但是一样很复杂。
物理学中尚未解决的积分,数学家们有啥办法?-1
001.png

beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
6楼2016-04-29 20:51:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

引用回帖:
5楼: Originally posted by 花痕逸韵 at 2016-04-29 16:23:05
题目有不严谨的地方,首先你不能确定它围成的空间是不是可以积分,这样就像要通解是不是有点过早。其次,如果参数的不确定,也可能导致积分不出来,比如混沌出现。

...

如果立方体与g(x,y,z)无重合,积分结果就是0。参数除了知道是有限的实数之外,无法事先预测。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
7楼2016-04-29 20:56:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghlsx

铁杆木虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by enlearner at 2016-04-29 10:41:05
这个问题其实不难,我说一下框架吧。(其实这个问题在代数学中已经解决了关键的部分)
解决这个问题的关键是积分区域g(x,y,z)的规则表达,也就是代数学中的“二次型的标准化,具体做法如下:
第一步:通过坐标旋转 ...

牛啊
不愧是数学专业的
数理真是基础科学啊
8楼2016-04-29 22:33:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

凯旋在子夜

金虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
拉格朗日(夹逼)定理

发自小木虫Android客户端
9楼2016-04-29 22:58:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

花痕逸韵

铁虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
10楼2016-04-29 23:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 beefly 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:23 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 8/400 2026-08-01 18:22 by Kddn1pB6pAzX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:17 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:11 by Kddn1pB6pAzX
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:03 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:51 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:42 by Kddn1pB6pAzX
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 19/950 2026-08-01 15:18 by 40862026
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见