| 查看: 1668 | 回复: 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  |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
孩子确诊有中度注意力缺陷
已经有6人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求助comsol中各物理场的调整参数的意义。
已经有5人回复- 关于菲克第二定律中的数学求解问题！
已经有4人回复
- 焊接中的扩散问题
已经有4人回复
- 对流-扩散方程中加入导热方程的问题
已经有3人回复
- 【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享”
已经有82人回复
- 【求助】求助拉格朗日方法追踪颗粒扩散
已经有8人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】不仅控制方程里含有未知参数，边界条件里也有
已经有5人回复
- 【求助】高维PDE的求解
已经有11人回复
- 【求助】求问一个寻找最优解问题，是泛函问题么
已经有22人回复
10楼2016-11-22 18:58:42
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