应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 基础数学 » 【求助】帮我看看这个积分怎么积？
北京石油化工学院2026年研究生招生接收调剂公告
161/2返回列表12下一页
查看: 1429  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

rising007

木虫 (正式写手)

[交流] 【求助】帮我看看这个积分怎么积？

帮我看看这个积分，其中a，b，r0，Cs，Cs2，Ys，Ds均为常数

谢谢

[ Last edited by laizuliang on 2009-3-19 at 12:21 ]
回复此楼

» 猜你喜欢
1楼 2009-03-18 16:18:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hongyan3278
2楼2009-03-18 16:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黄凯

禁虫 (小有名气)
本帖内容被屏蔽
3楼2009-03-18 19:06:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)
★ ★ ★
点点风(金币+3,VIP+0):感谢参与讨论，欢迎常来理工学版~ 5-28 17:51
>> help erf
ERF Error function.
    Y = ERF(X) is the error function for each element of X.  X must be
    real. The error function is defined as:

      erf(x) = 2/sqrt(pi) * integral from 0 to x of exp(-t^2) dt.

下面是我用计算机算的，你看看我写的式子和你写的一样吗，要是不一样你就写成这种形式，用计算机算就方便了，看上去很复杂啊，如果是实际应用，建议采用数值计算，这么复杂的式子你也不可能代进数去算吧，呵呵
>> a=sym('a')

a =

a


>> b=sym('b')

b =

b


>> fun=sym('4*pi*(r0-x)^2*(cs+(cs2-cs)*erf(x/2/sqrt(Ds*l))/erf(Ys/2/sqrt(Ds)))*r0/(r0-x)')

fun =

4*pi*(r0-x)^2*(cs+(cs2-cs)*erf(x/2/sqrt(Ds*l))/erf(Ys/2/sqrt(Ds)))*r0/(r0-x)


>> y=int(fun,a,b)

y =

2*pi^(1/2)*r0*(2*exp(1/4*a^2/Ds/l)*cs2*erf(1/2*b/(Ds*l)^(1/2))*r0*b*pi^(1/2)*exp(1/4*b^2/Ds/l)-exp(1/4*a^2/Ds/l)*cs2*erf(1/2*b/(Ds*l)^(1/2))*b^2*pi^(1/2)*exp(1/4*b^2/Ds/l)+4*exp(1/4*a^2/Ds/l)*cs2*(Ds*l)^(1/2)*r0-2*exp(1/4*a^2/Ds/l)*cs2*(Ds*l)^(1/2)*b+2*exp(1/4*a^2/Ds/l)*cs2*Ds*l*pi^(1/2)*erf(1/2*b/(Ds*l)^(1/2))*exp(1/4*b^2/Ds/l)-2*exp(1/4*a^2/Ds/l)*cs*erf(1/2*b/(Ds*l)^(1/2))*r0*b*pi^(1/2)*exp(1/4*b^2/Ds/l)+exp(1/4*a^2/Ds/l)*cs*erf(1/2*b/(Ds*l)^(1/2))*b^2*pi^(1/2)*exp(1/4*b^2/Ds/l)-4*exp(1/4*a^2/Ds/l)*cs*(Ds*l)^(1/2)*r0+2*exp(1/4*a^2/Ds/l)*cs*(Ds*l)^(1/2)*b-2*exp(1/4*a^2/Ds/l)*cs*Ds*l*pi^(1/2)*erf(1/2*b/(Ds*l)^(1/2))*exp(1/4*b^2/Ds/l)+2*exp(1/4*a^2/Ds/l)*cs*b*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)*exp(1/4*b^2/Ds/l)*r0-exp(1/4*a^2/Ds/l)*cs*b^2*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)*exp(1/4*b^2/Ds/l)-2*exp(1/4*b^2/Ds/l)*cs2*erf(1/2*a/(Ds*l)^(1/2))*r0*a*pi^(1/2)*exp(1/4*a^2/Ds/l)+exp(1/4*b^2/Ds/l)*cs2*erf(1/2*a/(Ds*l)^(1/2))*a^2*pi^(1/2)*exp(1/4*a^2/Ds/l)-4*exp(1/4*b^2/Ds/l)*cs2*(Ds*l)^(1/2)*r0+2*exp(1/4*b^2/Ds/l)*cs2*(Ds*l)^(1/2)*a-2*exp(1/4*b^2/Ds/l)*cs2*Ds*l*pi^(1/2)*erf(1/2*a/(Ds*l)^(1/2))*exp(1/4*a^2/Ds/l)+2*exp(1/4*b^2/Ds/l)*cs*erf(1/2*a/(Ds*l)^(1/2))*r0*a*pi^(1/2)*exp(1/4*a^2/Ds/l)-exp(1/4*b^2/Ds/l)*cs*erf(1/2*a/(Ds*l)^(1/2))*a^2*pi^(1/2)*exp(1/4*a^2/Ds/l)+4*exp(1/4*b^2/Ds/l)*cs*(Ds*l)^(1/2)*r0-2*exp(1/4*b^2/Ds/l)*cs*(Ds*l)^(1/2)*a+2*exp(1/4*b^2/Ds/l)*cs*Ds*l*pi^(1/2)*erf(1/2*a/(Ds*l)^(1/2))*exp(1/4*a^2/Ds/l)-2*exp(1/4*b^2/Ds/l)*cs*a*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)*exp(1/4*a^2/Ds/l)*r0+exp(1/4*b^2/Ds/l)*cs*a^2*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)*exp(1/4*a^2/Ds/l))/erf(1/2*Ys/Ds^(1/2))/exp(1/4*b^2/Ds/l)/exp(1/4*a^2/Ds/l)


>> simplify(y)

ans =

-2*pi^(1/2)*r0*(-2*exp(1/4*(a^2+b^2)/Ds/l)*cs2*erf(1/2*b/(Ds*l)^(1/2))*r0*b*pi^(1/2)+exp(1/4*(a^2+b^2)/Ds/l)*cs2*erf(1/2*b/(Ds*l)^(1/2))*b^2*pi^(1/2)-4*exp(1/4*a^2/Ds/l)*cs2*(Ds*l)^(1/2)*r0+2*exp(1/4*a^2/Ds/l)*cs2*(Ds*l)^(1/2)*b-2*exp(1/4*(a^2+b^2)/Ds/l)*cs2*Ds*l*pi^(1/2)*erf(1/2*b/(Ds*l)^(1/2))+2*exp(1/4*(a^2+b^2)/Ds/l)*cs*erf(1/2*b/(Ds*l)^(1/2))*r0*b*pi^(1/2)-exp(1/4*(a^2+b^2)/Ds/l)*cs*erf(1/2*b/(Ds*l)^(1/2))*b^2*pi^(1/2)+4*exp(1/4*a^2/Ds/l)*cs*(Ds*l)^(1/2)*r0-2*exp(1/4*a^2/Ds/l)*cs*(Ds*l)^(1/2)*b+2*exp(1/4*(a^2+b^2)/Ds/l)*cs*Ds*l*pi^(1/2)*erf(1/2*b/(Ds*l)^(1/2))-2*exp(1/4*(a^2+b^2)/Ds/l)*cs*b*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)*r0+exp(1/4*(a^2+b^2)/Ds/l)*cs*b^2*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)+2*exp(1/4*(a^2+b^2)/Ds/l)*cs2*erf(1/2*a/(Ds*l)^(1/2))*r0*a*pi^(1/2)-exp(1/4*(a^2+b^2)/Ds/l)*cs2*erf(1/2*a/(Ds*l)^(1/2))*a^2*pi^(1/2)+4*exp(1/4*b^2/Ds/l)*cs2*(Ds*l)^(1/2)*r0-2*exp(1/4*b^2/Ds/l)*cs2*(Ds*l)^(1/2)*a+2*exp(1/4*(a^2+b^2)/Ds/l)*cs2*Ds*l*pi^(1/2)*erf(1/2*a/(Ds*l)^(1/2))-2*exp(1/4*(a^2+b^2)/Ds/l)*cs*erf(1/2*a/(Ds*l)^(1/2))*r0*a*pi^(1/2)+exp(1/4*(a^2+b^2)/Ds/l)*cs*erf(1/2*a/(Ds*l)^(1/2))*a^2*pi^(1/2)-4*exp(1/4*b^2/Ds/l)*cs*(Ds*l)^(1/2)*r0+2*exp(1/4*b^2/Ds/l)*cs*(Ds*l)^(1/2)*a-2*exp(1/4*(a^2+b^2)/Ds/l)*cs*Ds*l*pi^(1/2)*erf(1/2*a/(Ds*l)^(1/2))+2*exp(1/4*(a^2+b^2)/Ds/l)*cs*a*erf(1/2*Ys/Ds^(1/2))*pi^(1/2)*r0-exp(1/4*(a^2+b^2)/Ds/l)*cs*a^2*erf(1/2*Ys/Ds^(1/2))*pi^(1/2))*exp(-1/4*(a^2+b^2)/Ds/l)/erf(1/2*Ys/Ds^(1/2))


>> pretty(y)

             /            2                               2
      1/2    |           a                  1/2          b
  2 pi    r0 |2 exp(1/4 ----) cs2 %3 r0 b pi    exp(1/4 ----)
             \          Ds l                            Ds l

                     2                             2
                    a            2   1/2          b
         - exp(1/4 ----) cs2 %3 b  pi    exp(1/4 ----)
                   Ds l                          Ds l

                       2                                  2
                      a              1/2                 a              1/2
         + 4 exp(1/4 ----) cs2 (Ds l)    r0 - 2 exp(1/4 ----) cs2 (Ds l)    b
                     Ds l                               Ds l

                       2                               2
                      a               1/2             b
         + 2 exp(1/4 ----) cs2 Ds l pi    %3 exp(1/4 ----)
                     Ds l                            Ds l

                       2                              2
                      a                 1/2          b
         - 2 exp(1/4 ----) cs %3 r0 b pi    exp(1/4 ----)
                     Ds l                           Ds l

                     2                            2
                    a           2   1/2          b
         + exp(1/4 ----) cs %3 b  pi    exp(1/4 ----)
                   Ds l                         Ds l

                       2                                 2
                      a             1/2                 a             1/2
         - 4 exp(1/4 ----) cs (Ds l)    r0 + 2 exp(1/4 ----) cs (Ds l)    b
                     Ds l                              Ds l

                       2                              2
                      a              1/2             b
         - 2 exp(1/4 ----) cs Ds l pi    %3 exp(1/4 ----)
                     Ds l                           Ds l

                       2                           2
                      a              1/2          b
         + 2 exp(1/4 ----) cs b %1 pi    exp(1/4 ----) r0
                     Ds l                        Ds l

                     2                            2
                    a        2      1/2          b
         - exp(1/4 ----) cs b  %1 pi    exp(1/4 ----)
                   Ds l                         Ds l

                       2                               2
                      b                  1/2          a
         - 2 exp(1/4 ----) cs2 %2 r0 a pi    exp(1/4 ----)
                     Ds l                            Ds l

                     2                             2
                    b            2   1/2          a
         + exp(1/4 ----) cs2 %2 a  pi    exp(1/4 ----)
                   Ds l                          Ds l

                       2                                  2
                      b              1/2                 b              1/2
         - 4 exp(1/4 ----) cs2 (Ds l)    r0 + 2 exp(1/4 ----) cs2 (Ds l)    a
                     Ds l                               Ds l

                       2                               2
                      b               1/2             a
         - 2 exp(1/4 ----) cs2 Ds l pi    %2 exp(1/4 ----)
                     Ds l                            Ds l

                       2                              2
                      b                 1/2          a
         + 2 exp(1/4 ----) cs %2 r0 a pi    exp(1/4 ----)
                     Ds l                           Ds l

                     2                            2
                    b           2   1/2          a
         - exp(1/4 ----) cs %2 a  pi    exp(1/4 ----)
                   Ds l                         Ds l

                       2                                 2
                      b             1/2                 b             1/2
         + 4 exp(1/4 ----) cs (Ds l)    r0 - 2 exp(1/4 ----) cs (Ds l)    a
                     Ds l                              Ds l

                       2                              2
                      b              1/2             a
         + 2 exp(1/4 ----) cs Ds l pi    %2 exp(1/4 ----)
                     Ds l                           Ds l

                       2                           2
                      b              1/2          a
         - 2 exp(1/4 ----) cs a %1 pi    exp(1/4 ----) r0
                     Ds l                        Ds l

                     2                            2  \     /             2
                    b        2      1/2          a   |   / |            b
         + exp(1/4 ----) cs a  %1 pi    exp(1/4 ----)|  /  |%1 exp(1/4 ----)
                   Ds l                         Ds l / /   \           Ds l

                  2  \
                 a   |
        exp(1/4 ----)|
                Ds l /

                 Ys
  %1 := erf(1/2 -----)
                  1/2
                Ds

                    a
  %2 := erf(1/2 ---------)
                      1/2
                (Ds l)

                    b
  %3 := erf(1/2 ---------)
                      1/2
                (Ds l)
>>

[ Last edited by fspdlh on 2009-3-18 at 20:41 ]
赞一下(10人) 回复此楼
4楼2009-03-18 20:27:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ying7666

银虫 (正式写手)
★
点点风(金币+1,VIP+0):感谢参与讨论，欢迎常来理工学版~ 5-28 17:51
我晕，原来不是初等函数。求近似吧。明白楼上的做法了。
赞一下(10人) 回复此楼
5楼2009-05-28 16:41:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuishao_小

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
见识了,我还没用过,符号积分. 这位仁兄, 能不能对结果再简化一下, 弄一个大家,看起来容易接受的形式啊.
赞一下(1人) 回复此楼
6楼2009-06-04 17:20:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

栾钰琨

木虫 (正式写手)
真牛呀
回复此楼
7楼2009-06-05 20:14:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duzhuqi

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
太复杂了哈，只要知道erf(x)是误差函数，知道它的表达式，结合广义积分还是很好求解得哈
赞一下(1人) 回复此楼
8楼2009-06-06 22:11:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuewu_liu

银虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
用高斯积分进行近似计算，很方便的
赞一下(1人) 回复此楼
yuewu
9楼2009-06-08 17:34:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujin41658

木虫 (著名写手)

执法队伍新手
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
学习了，不过不太明白哦~~~~~~
赞一下(1人) 回复此楼
欢迎前辈指导！
10楼2009-06-09 14:09:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 rising007 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +6 双马尾痞老板2 2026-04-01 6/300 2026-04-01 23:33 by chaolymer
[考研] 266求调剂 +7 阳阳哇塞 2026-04-01 7/350 2026-04-01 22:27 by barlinike
[考研] 270求调剂 +7 小杰pp 2026-03-31 8/400 2026-04-01 21:36 by 七度不信任
[考研] 339求调剂，想调回江苏 +7 烤麦芽 2026-03-27 10/500 2026-04-01 21:35 by 495374996
[考研] 085601材料工程找调剂 +20 oatmealR 2026-03-29 21/1050 2026-04-01 21:00 by lijunpoly
[考研] 070300一志愿211，312分求调剂院校 +14 小黄鸭宝 2026-03-30 14/700 2026-04-01 20:19 by 赖春艳
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
[考研] 一志愿厦门大学化学工程（专硕）-数二英二406分-求调剂 +4 厦大化工 2026-04-01 4/200 2026-04-01 17:54 by JYD2011
[考研] 350求调剂 +7 阿佳～ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
[考研] 318求调剂 +3 笃行致远. 2026-03-31 3/150 2026-03-31 20:27 by 求调剂zz
[考研] 生医工0831调剂求推荐 +4 小熊睿睿_s 2026-03-27 6/300 2026-03-31 17:12 by 记事本2026
[考研] 340求调剂 +4 希望如此i 2026-03-31 4/200 2026-03-31 16:40 by 690616278
[考研] 求化学调剂 +12 wulanna 2026-03-28 12/600 2026-03-31 16:38 by 690616278
[考研] 353求调剂 +3 江上枫_26 2026-03-28 3/150 2026-03-31 15:53 by jp9609
[考研] 085601 329分调剂 +6 yzsa12 2026-03-31 6/300 2026-03-31 15:23 by yanflower7133
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +13 辰沐5211314 2026-03-26 13/650 2026-03-31 14:37 by 记事本2026
[考研] 303求调剂 +7 DLkz1314. 2026-03-30 7/350 2026-03-30 16:05 by shuang5186
[考研] 一志愿上海理工能源动力（085800）310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
信息提示
关闭
请填处理意见
关闭