24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1893  |  回复: 17

cooooldog

铁杆木虫 (著名写手)

ส็็็

[交流] 一个看起来简单的不太容易的定积分问题已有9人参与

求如下定积分的解析形式,希望解析形式中不出现变量,即结果只是常量的函数形式。
比如,椭圆积分表达不符合要求。



此外:可积性是显然的:
一个看起来简单的不太容易的定积分问题



[ Last edited by cooooldog on 2017-7-15 at 18:34 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

ส็็็็็็็็็็็็็็็็็็็็
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Edstrayer

版主 (著名写手)

方寸斗室小天地正气迷漫大世界


小木虫: 金币+0.5, 给个红包,谢谢回帖
函数在区间[0,1]上连续,当然可积,但是精确的积分值却难于计算,只能求近似的数值解。

青葱岁月圣诞夜,浪漫歌舞迎新年。
4楼2015-08-25 04:36:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主


小木虫: 金币+0.5, 给个红包,谢谢回帖
1)如果能搞出x^2/sqrt(1+x^5)的傅里叶变换解析式,或拉普拉斯变换解析式,次积分可以有解析表达。
2)有限项级数近似。
3楼2015-08-24 21:48:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

铸剑人

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
在零点展开成无穷级数,泰勒或者傅立叶级数都可以,然后积分就可以了。也是一个求极限的问题。

[ 发自小木虫客户端 ]
绝尘。
6楼2015-08-25 05:28:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
3楼: Originally posted by feixiaolin at 2015-08-24 21:48:40
1)如果能搞出x^2/sqrt(1+x^5)的傅里叶变换解析式,或拉普拉斯变换解析式,次积分可以有解析表达。
2)有限项级数近似。


我尝试符号计算傅立叶变换,然后得到一个很长的很复杂的式子。

到底有多长呢?
我尝试贴图但是传不上来。

把公式排成一列放长安街上,站在公式的尾部,一眼望过去看不到公司的头部
ส็็็็็็็็็็็็็็็็็็็็
7楼2015-08-25 07:23:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫


小木虫: 金币+0.5, 给个红包,谢谢回帖
int('x*sin(x)*(x^5+1)^(-1/5)',0,1)

ans =

pi^(1/2)*gamma(4/5)*5^(1/2)*2^(1/5)*(-1)^(1/5)*(1/2/pi^(1/2)*csc(2/5*pi)*sin(1/10*pi)*gamma(4/5)^2/gamma(7/10)^2*hypergeom([],[1/5, 2/5, 1/2, 3/5, 7/10, 7/10, 9/10, 11/10, 6/5],-1/10000000000)-1/24/pi*hypergeom([1],[2/5, 3/5, 7/10, 4/5, 9/10, 9/10, 11/10, 6/5, 13/10, 7/5],-1/10000000000)*tan(3/10*pi)*sin(1/10*pi)*gamma(3/5)/gamma(9/10)+1/720*2^(1/5)/pi*hypergeom([],[3/5, 4/5, 9/10, 11/10, 11/10, 13/10, 7/5, 3/2, 8/5],-1/10000000000)*sin(1/10*pi)*csc(1/5*pi)*csc(2/5*pi)/csc(3/10*pi)*gamma(7/10)-1/60480*2^(1/5)/pi*hypergeom([],[4/5, 11/10, 6/5, 13/10, 13/10, 3/2, 8/5, 17/10, 9/5],-1/10000000000)*csc(2/5*pi)^2*cos(1/5*pi)/csc(3/10*pi)*gamma(7/10)+1/30240000*hypergeom([],[6/5, 13/10, 7/5, 3/2, 3/2, 17/10, 9/5, 19/10, 2],-1/10000000000)*sec(1/10*pi)*cot(1/5*pi)*gamma(3/5)/gamma(9/10)+2*i*pi*(1/4*2^(1/5)/pi^(3/2)*csc(2/5*pi)*sin(1/10*pi)*gamma(9/10)/gamma(7/10)*(-1)^(1/5)*hypergeom([],[1/5, 2/5, 1/2, 3/5, 7/10, 7/10, 9/10, 11/10, 6/5],-1/10000000000)-1/48*2^(4/5)*(-1)^(2/5)/pi^(3/2)*hypergeom([1],[2/5, 3/5, 7/10, 4/5, 9/10, 9/10, 11/10, 6/5, 13/10, 7/5],-1/10000000000)*tan(3/10*pi)*sin(1/10*pi)/gamma(4/5)+1/720*2^(1/5)/pi^2*csc(1/5*pi)*sin(1/10*pi)*sin(3/10*pi)*gamma(7/10)*(-1)^(3/5)*hypergeom([],[3/5, 4/5, 9/10, 11/10, 11/10, 13/10, 7/5, 3/2, 8/5],-1/10000000000)-1/120960*2^(1/5)/pi^2*csc(2/5*pi)*cos(1/5*pi)*gamma(7/10)*(-1)^(4/5)*hypergeom([],[4/5, 11/10, 6/5, 13/10, 13/10, 3/2, 8/5, 17/10, 9/5],-1/10000000000)))

以上为matlab code...有耐心自己去搞定
板凳要做十年冷文章不发一个字
10楼2015-08-25 08:23:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wurongjun

专家顾问 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
椭圆积分!
只能数值解!

>> quad('x.*sin(x).*(1+x.^5).^(1/2)',0,1)
ans =
    0.3494
善恶到头终有报,人间正道是沧桑.
2楼2015-08-24 16:49:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

俄罗斯棕熊

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主你能别开玩笑吗?

[ 发自手机版 http://muchong.com/3g ]
5楼2015-08-25 05:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
5楼: Originally posted by 俄罗斯棕熊 at 2015-08-25 05:20:30
楼主你能别开玩笑吗?

你好像很懂。 解释下。

就算用椭圆积分也行,表示出来
ส็็็็็็็็็็็็็็็็็็็็
8楼2015-08-25 07:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主

引用回帖:
7楼: Originally posted by cooooldog at 2015-08-25 07:23:48

我尝试符号计算傅立叶变换,然后得到一个很长的很复杂的式子。

到底有多长呢?
我尝试贴图但是传不上来。

把公式排成一列放长安街上,站在公式的尾部,一眼望过去看不到公司的头部...

1)不行,就选2)。

[ 发自小木虫客户端 ]
9楼2015-08-25 08:03:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cooooldog 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见