| 查看: 2088 | 回复: 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  |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有118人回复
-
今天放榜吗？
已经有15人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
filecode，4个jtjc了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天基金会出结果吗？20260819
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助comsol中各物理场的调整参数的意义。
已经有5人回复
- 关于菲克第二定律中的数学求解问题！
已经有4人回复
- 焊接中的扩散问题
已经有4人回复
- 对流-扩散方程中加入导热方程的问题
已经有3人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】求助拉格朗日方法追踪颗粒扩散
已经有8人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】不仅控制方程里含有未知参数，边界条件里也有
已经有5人回复
- 【求助】高维PDE的求解
已经有11人回复
- 【求助】求问一个寻找最优解问题，是泛函问题么
已经有22人回复
5楼2012-04-03 14:58:46
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
