| 查看: 2248 | 回复: 14 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
dengwei_19木虫 (著名写手)
铁齿金不换诚实可靠小郎君
|
[求助]
现在做到高斯函数的误差函数一步，又是一个方程不知道用matlab怎么写
|
|||
求w，求详解！谢谢！！！ |
回复此楼» 本帖@通知
» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有15人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有14人回复
-
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够
已经有17人回复
-
基金系统什么内容也没有
已经有10人回复
-
学科评审组评审是指会评吗？
已经有5人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
dengwei_19: 回帖置顶 2012-04-13 15:38:54
|
fun=@(x)-5.4*(pi^(1/2)*erf(10*2^0.5./x))/2-2.8*(pi^(1/2)*erf(40*2^0.5./x))/2+2.6; 这个是将你的积分方程式转为误差函数，-(pi^(1/2)*erf(10*2^0.5./x))/2，(pi^(1/2)*erf(40*2^0.5./x))/2，是你的两个积分 w（即x）的解使得fun=0 x=0.1:0.1:100; y=feval(fun,x); plot(x,y) 这个是w（即x）在0到100之间变化时fun的取值，用来在做图猜测w大概在什么位置，fun=0 [x,feval,flag]=fsolve(fun,80)是用来求精确数值解 x解即w feval 取x的值 fun的函数值越接近0越好 flag 判断解的情况 1 FSOLVE converged to a root. 2 Change in X too small. 3 Change in residual norm too small. 4 Computed search direction too small. 0 Too many function evaluations or iterations. -1 Stopped by output/plot function. -2 Converged to a point that is not a root. -3 Trust region radius too small (Trust-region-dogleg) or Regularization parameter too large (Levenberg-Marquardt). -4 Line search failed. |
dengwei_19
11楼2012-04-13 15:23:53
onesupeng
金虫 (职业作家)
- 数学EPI: 17
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2116.7
- 散金: 9437
- 红花: 92
- 帖子: 4585
- 在线: 1305小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
2楼2012-04-13 13:47:46
dengwei_19
木虫 (著名写手)
铁齿金不换诚实可靠小郎君
- 应助: 103 (高中生)
- 贵宾: 0.209
- 金币: 1797
- 散金: 906
- 红花: 17
- 帖子: 2059
- 在线: 404.6小时
- 虫号: 242330
- 注册: 2006-04-11
- 性别: GG
- 专业: 玻璃材料
不会吧。。。。。没办法求吗？ 3楼2012-04-13 13:53:05
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
dengwei_19: 金币+6, ★★★很有帮助 2012-04-13 14:36:22
dengwei_19: 回帖置顶 2012-04-13 15:40:01
感谢参与，应助指数 +1
dengwei_19: 金币+6, ★★★很有帮助 2012-04-13 14:36:22
dengwei_19: 回帖置顶 2012-04-13 15:40:01
|
>> syms x >> syms w;int(exp(-x^2),0,w) ans = (pi^(1/2)*erf(w))/2 >> syms w;int(exp(-x^2),0,-w) ans = -(pi^(1/2)*erf(w))/2 =========== >> fun=@(x)-5.4*(pi^(1/2)*erf(10*2^0.5./x))/2-2.8*(pi^(1/2)*erf(40*2^0.5./x))/2+2.6; >> x=0.1:0.1:100; >> y=feval(fun,x); >> plot(x,y) >> [x,feval,flag]=fsolve(fun,80) Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. x = 81.4787 feval = -9.0775e-010 flag = 1 |
» 本帖已获得的红花（最新10朵）
4楼2012-04-13 14:19:25
30