24小时热门版块排行榜    

查看: 1911  |  回复: 11

阳光男孩羊

铜虫 (小有名气)

[求助] 怎样求解我这个简单积分,哪位大神帮忙,小弟谢过啦 已有2人参与

帮我求求这个式子吧,万分感谢

怎样求解我这个简单积分,哪位大神帮忙,小弟谢过啦
P50804-100310.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

peterflyer

木虫之王 (文学泰斗)

peterflyer


引用回帖:
5楼: Originally posted by 阳光男孩羊 at 2015-08-04 19:58:44
恩恩,那大神数值积分怎么计算啊?我没写的参数都是定值,不用管...

要计算数值,就要先把不定积分改为定积分:Integral{dφ/sqrt[Sin(φ)] , α0, α}=C'- β*t
其中为直线的初始位置的角度;α为t时刻直线的位置的角度。当α=π/2时直线就转到了X轴的位置。
由 t=0时α=α0,故代入方程中后得:C'=0
故:Integral{dφ/sqrt[Sin(φ)] , α0, α}=- β*t
当=π/2时,t=T
T=-1/β*Integral{dφ/sqrt[Sin(φ)] , α0, π/2}
现在用高斯勒让德数值计分方法计算。
令φ=(α0+π/2)/2+(π/2-α0)/2*ξ ,dφ=(π/2-α0)/2*dξ
原积分=-1/β*Integral{(π/2-α0)/2*dξ/Sin[(α0+π/2)/2+(π/2-α0)/2*ξ] ,-1 , 1}
=-(π/2-α0)/(2*β)*Integral{dξ/Sin[(α0+π/2)/2+(π/2-α0)/2*ξ] ,-1 , 1}
=-(π/2-α0)/(2*β)*Integral{dξ/f(ξ) ,-1 , 1}
=-(π/2-α0)/(2*β)*Sum{δk*1/f(ηk), k=1~n}
其中,n为高斯点数量,根据需要选取。原则是计算精度要求高时可多取,反之可少取。δk为插值点ηk处对应的被积函数的值的系数。
比如,n=1时,η=0,δ=2;若n=2,则η1=-1/sqrt(3),δ1=1,η2=1/sqrt(3),δ2=1 ;n=3时,......。具体可参考下面链接:
http://wenku.baidu.com/link?url= ... RZmLWeZFFALlBgoKTvC
6楼2015-08-04 20:34:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
由于φ(t)是逐渐减小的,因此dφ/dt=-ω(t)=-β*sqrt[Sin(φ)] ,这里的β代表前面的系数,由于楼主的图片不清晰,只好用这个代替了。
dφ/sqrt[Sin(φ)]=-ω(t)=-β*dt
Integral{dφ/sqrt[Sin(φ)]}=C'-β*t
C' 为积分常数。
左边的积分中,令t=tan(φ/2)可得:
Integral{dt/sqrt[t*(1+t^2)]}=C-β/sqrt(2)*t
俄国数学家切比雪夫在一百年前就已经证明,这种积分是不能够用初等函数表示的,也就是积不出来的。只能够借助数值积分方法了。
2楼2015-08-04 11:26:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

越接近x轴,角速度越慢并无限向零靠近,如此下去这个直线可能永远也到不了x轴。
3楼2015-08-04 11:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阳光男孩羊

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by peterflyer at 2015-08-04 11:33:42
越接近x轴,角速度越慢并无限向零靠近,如此下去这个直线可能永远也到不了x轴。

不会吧,肯定能算近似值吧
4楼2015-08-04 19:57:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阳光男孩羊

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by peterflyer at 2015-08-04 11:26:44
由于φ(t)是逐渐减小的,因此dφ/dt=-ω(t)=-β*sqrt ,这里的β代表前面的系数,由于楼主的图片不清晰,只好用这个代替了。
dφ/sqrt=-ω(t)=-β*dt
Integral{dφ/sqrt}=C'-β*t
C' 为积分常数。
左边的积分中, ...

恩恩,那大神数值积分怎么计算啊?我没写的参数都是定值,不用管
5楼2015-08-04 19:58:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阳光男孩羊

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by peterflyer at 2015-08-04 20:34:59
要计算数值,就要先把不定积分改为定积分:Integral{dφ/sqrt , α0, α}=C'- β*t
其中为直线的初始位置的角度;α为t时刻直线的位置的角度。当α=π/2时直线就转到了X轴的位置。
由 t=0时α=α0,故代入方程中 ...

大神,我是菜鸟,还是不懂。能帮人帮到底吗?我很急啊,万分感谢。结果到底是什么啊?
7楼2015-08-05 08:57:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


引用回帖:
7楼: Originally posted by 阳光男孩羊 at 2015-08-05 08:57:02
大神,我是菜鸟,还是不懂。能帮人帮到底吗?我很急啊,万分感谢。结果到底是什么啊?...

比如就取两个点,n=2。为简化计算,假设α0=0,积分式里的系数β=1。那么就得到:
原式=-π/4*{1*1/Sin{π/4+π/4*[-1/sqrt(3)]}+1*1/Sin{π/4+π/4*[1/sqrt(3)]}}
      =具体算出数字
8楼2015-08-05 09:24:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阳光男孩羊

铜虫 (小有名气)

引用回帖:
8楼: Originally posted by peterflyer at 2015-08-05 09:24:32
比如就取两个点,n=2。为简化计算,假设α0=0,积分式里的系数β=1。那么就得到:
原式=-π/4*{1*1/Sin{π/4+π/4*}+1*1/Sin{π/4+π/4*}}
      =具体算出数字...

积分里的β是个初始角,不能先赋值吧?我想得到一个含有β的式子
9楼2015-08-12 09:26:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


引用回帖:
9楼: Originally posted by 阳光男孩羊 at 2015-08-12 09:26:52
积分里的β是个初始角,不能先赋值吧?我想得到一个含有β的式子...

这个只是积分的系数,不是你那个β角。
10楼2015-08-12 11:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阳光男孩羊 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 3/150 2026-11-02 01:00 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 6/300 2026-11-02 00:29 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-02 00:15 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-02 00:08 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 19/950 2026-11-02 00:02 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 23:55 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 23/1150 2026-11-01 23:49 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 19:50 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 14:48 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
请填处理意见