24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1819  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见