版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

当前主题已经存档。

rising007

版主

应助: 1 (幼儿园)
金币: 2073.7
散金: 561
帖子: 709
在线: 299.8小时
虫号: 654091
注册: 2008-11-14
性别: GG
专业: 钢铁冶金

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

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

谢谢

[ Last edited by laizuliang on 2009-3-19 at 12:21 ]

回复此楼

» 猜你喜欢

08生物与医药专硕初试346找调剂已经有5人回复
江苏科技大学招材料研究生已经有4人回复
0856，269分求调剂已经有5人回复
303分 0807学硕求调剂已经有3人回复
085600 一志愿9 总分351 求调剂学校已经有9人回复
一志愿郑大材料工程290求调剂已经有18人回复
一志愿北京科技，085601总分305求调剂已经有4人回复
0710生物学求调剂已经有3人回复
285求调剂已经有8人回复
273求调剂已经有18人回复

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

超级版主

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

★ ★ ★
点点风(金币+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

超级版主

应助: 0 (幼儿园)
金币: 427.9
帖子: 688
在线: 11.4小时
虫号: 598967
注册: 2008-09-11
专业: 计算机科学的基础理论

★
点点风(金币+1,VIP+0):感谢参与讨论，欢迎常来理工学版~ 5-28 17:51

我晕，原来不是初等函数。求近似吧。明白楼上的做法了。

赞一下(10人)

回复此楼

5楼2009-05-28 16:41:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

cuishao_小

超级版主

应助: 0 (幼儿园)
金币: 5560.7
散金: 510
红花: 2
沙发: 93
帖子: 1402
在线: 289.2小时
虫号: 751467
注册: 2009-04-18
性别: GG
专业: 概率论与随机分析

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

见识了,我还没用过,符号积分. 这位仁兄, 能不能对结果再简化一下, 弄一个大家,看起来容易接受的形式啊.

赞一下(1人)

回复此楼

6楼2009-06-04 17:20:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

栾钰琨

主管区长

应助: 0 (幼儿园)
金币: 1538
帖子: 400
在线: 48分钟
虫号: 765883
注册: 2009-05-08
专业: 机械电子工程

真牛呀

回复此楼

7楼2009-06-05 20:14:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

duzhuqi

管理员

应助: (幼儿园)
在线:
虫号: 716479

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

太复杂了哈，只要知道erf(x)是误差函数，知道它的表达式，结合广义积分还是很好求解得哈

赞一下(1人)

回复此楼

8楼2009-06-06 22:11:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yuewu_liu

超级版主

应助: 3 (幼儿园)
金币: 801.8
红花: 4
帖子: 442
在线: 348小时
虫号: 464760
注册: 2007-11-21
性别: GG
专业: 生物信息学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

用高斯积分进行近似计算，很方便的

赞一下(1人)

回复此楼

yuewu

9楼2009-06-08 17:34:23

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wujin41658

兑换贵宾

执法队伍新手

应助: 40 (小学生)
贵宾: 0.58
金币: 2058
散金: 16
红花: 9
帖子: 1618
在线: 62小时
虫号: 566863
注册: 2008-05-31
性别: GG
专业: 食品加工技术

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

学习了，不过不太明白哦~~~~~~

赞一下(1人)

回复此楼

欢迎前辈指导！

10楼2009-06-09 14:09:08

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 rising007 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 349求调剂 +6	吃的不少 2026-04-01	6/300	2026-04-01 17:55 by JYD2011
[考研] 285求调剂 +5	FZAC123 2026-03-30	5/250	2026-04-01 15:50 by 韩雨涵
[考研] 085600，材料与化工321分调剂 +5	大馋小子 2026-03-27	7/350	2026-04-01 15:13 by oooqiao
[论文投稿] chinese chemical letters英文版投稿求助 120+4	Yishengeryi 2026-03-30	5/250	2026-04-01 14:11 by 陆小果画大饼
[考研] 291求调剂 +3	迷蒙木木 2026-04-01	4/200	2026-04-01 11:07 by 逆水乘风
[考研] 085701环境工程，267求调剂 +17	minht 2026-03-26	17/850	2026-04-01 09:11 by xiayizhi
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +9	手机用户 2026-03-26	9/450	2026-03-31 18:09 by 544594351
[考研] 求调剂生物学 377分 +6	zzll03 2026-03-31	6/300	2026-03-31 17:33 by 唐沐儿
[考研] 求化学调剂 +12	wulanna 2026-03-28	12/600	2026-03-31 16:38 by 690616278
[考研] 生物考研337分求调剂 +4	cgxin 2026-03-30	6/300	2026-03-31 14:18 by 记事本2026
[考研] 吉大生物学326分求调剂 +3	sunnyupup 2026-03-31	3/150	2026-03-31 09:28 by longlotian
[考研] 085601一志愿西北工业大学初试346 +4	085601初试346 2026-03-30	4/200	2026-03-31 07:47 by jp9609
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu锂
[考研] 一志愿南昌大学324求调剂 +9	hanamiko 2026-03-27	9/450	2026-03-30 20:10 by 无际的草原
[考研] 0703化学求调剂 +6	丹青奶盖 2026-03-26	8/400	2026-03-30 18:33 by 探123
[考研] 284求调剂 +14	junqihahaha 2026-03-26	15/750	2026-03-30 14:12 by 探123
[考研] 295求调剂 +5	wei-5 2026-03-26	5/250	2026-03-30 08:34 by 探123
[考研] 药学105500求调剂 +3	Ssun。。 2026-03-28	3/150	2026-03-28 11:24 by lxf170613
[考研] 305求调剂 +5	哇卢卡库 2026-03-26	5/250	2026-03-27 14:01 by laoshidan
[考研] 考研调剂 +9	小蜡新笔 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚的貓