| 查看: 2084 | 回复: 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  |
回复此楼» 猜你喜欢
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天放榜吗？
已经有13人回复
-
filecode，4个jtjc了
已经有9人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有115人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助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
- 专业: 污染控制化学
- 管辖: 计算模拟
7楼2012-04-03 17:58:17
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
