| 查看: 2244 | 回复: 9 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
求解扩散方程问题
|
|||||
|
有一个方程如图 源码是 function AA D0=3.55e-19; options=optimset('tolx',1e-9); [D,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@fei,D0,[],[],options) function F=fei(D) C0=0.687; d=0.02; Ceq=20; t=60; C=4.136; m=1-(C0-C)/(C0-Ceq); s=[]; k=0:1000; s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2*t*60/(d^2))); s F=sum((m-s').^2); 而得到的结果怎么也不对 如下 得到的D值怎么都是跟赋的D0值是一样的额 s = 0.9998 s = 0.2158 s = Inf s = 0.2447 s = Inf s = 0.2161 Optimization terminated: norm of the current step is less than OPTIONS.TolX. D = 3.5500e-019 resnorm = 0.0010 residual = 0.0318 exitflag = 2 output = firstorderopt: 7.1546e+005 iterations: 2 funcCount: 6 cgiterations: 1 algorithm: 'large-scale: trust-region reflective Newton' message: [1x77 char] lambda = lower: 0 upper: 0 jacobian = (1,1) 2.2480e+007 到底怎么弄得，小弟是菜鸟，希望有高手跟好心人帮我5555  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
考博求助
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助comsol中各物理场的调整参数的意义。
已经有5人回复
- 关于菲克第二定律中的数学求解问题！
已经有4人回复
- 焊接中的扩散问题
已经有4人回复
- 对流-扩散方程中加入导热方程的问题
已经有3人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】求助拉格朗日方法追踪颗粒扩散
已经有8人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】不仅控制方程里含有未知参数，边界条件里也有
已经有5人回复
- 【求助】高维PDE的求解
已经有11人回复
- 【求助】求问一个寻找最优解问题，是泛函问题么
已经有22人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
臭水沟: 金币+3, 谢谢交流~~ 2012-04-04 11:15:22
zhangluyoa: 金币+4, ★有帮助 2012-04-11 16:35:41
感谢参与，应助指数 +1
臭水沟: 金币+3, 谢谢交流~~ 2012-04-04 11:15:22
zhangluyoa: 金币+4, ★有帮助 2012-04-11 16:35:41
|
Local minimum possible. lsqnonlin stopped because the final change in the sum of squares relative to its initial value is less than the default value of the function tolerance. D = 5.0774e-008 resnorm = 1.1370e-007 residual = 3.3720e-004 exitflag = 3 output = firstorderopt: 39.0251 iterations: 13 funcCount: 28 cgiterations: 0 algorithm: 'large-scale: trust-region reflective Newton' message: [1x458 char] lambda = lower: 0 upper: 0 jacobian = (1,1) 1.1573e+005 |
6楼2012-04-03 17:56:51
ya634984094
铜虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 74.7
- 散金: 4
- 红花: 3
- 帖子: 259
- 在线: 38.5小时
- 虫号: 1695458
- 注册: 2012-03-16
- 专业: 流体力学
2楼2012-04-03 10:54:59
3楼2012-04-03 11:35:39
ya634984094
铜虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 74.7
- 散金: 4
- 红花: 3
- 帖子: 259
- 在线: 38.5小时
- 虫号: 1695458
- 注册: 2012-03-16
- 专业: 流体力学
4楼2012-04-03 12:34:11
