| 查看: 1695 | 回复: 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  |
回复此楼» 猜你喜欢
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求助comsol中各物理场的调整参数的意义。
已经有5人回复- 关于菲克第二定律中的数学求解问题！
已经有4人回复
- 焊接中的扩散问题
已经有4人回复
- 对流-扩散方程中加入导热方程的问题
已经有3人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】求助拉格朗日方法追踪颗粒扩散
已经有8人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】不仅控制方程里含有未知参数，边界条件里也有
已经有5人回复
- 【求助】高维PDE的求解
已经有11人回复
- 【求助】求问一个寻找最优解问题，是泛函问题么
已经有22人回复
ya634984094
铜虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 74.7
- 散金: 4
- 红花: 3
- 帖子: 259
- 在线: 38.5小时
- 虫号: 1695458
- 注册: 2012-03-16
- 专业: 流体力学
2楼2012-04-03 10:54:59
ya634984094
铜虫 (小有名气)
- 应助: 39 (小学生)
- 金币: 74.7
- 散金: 4
- 红花: 3
- 帖子: 259
- 在线: 38.5小时
- 虫号: 1695458
- 注册: 2012-03-16
- 专业: 流体力学
4楼2012-04-03 12:34:11