24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangluyoa

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 40
帖子: 41
在线: 14.2小时
虫号: 867634
注册: 2009-10-10
性别: GG
专业: 食品加工技术

[求助] 求解扩散方程问题

有一个方程如图
源码是
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

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助comsol中各物理场的调整参数的意义。已经有5人回复
关于菲克第二定律中的数学求解问题！已经有4人回复
焊接中的扩散问题已经有4人回复
对流-扩散方程中加入导热方程的问题已经有3人回复
【500金币本版活动】“经验力场分子模拟+第一性原理计算问题讨论互助及资源共享” 已经有82人回复
【求助】求助拉格朗日方法追踪颗粒扩散已经有8人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】不仅控制方程里含有未知参数，边界条件里也有已经有5人回复
【求助】高维PDE的求解已经有11人回复
【求助】求问一个寻找最优解问题，是泛函问题么已经有22人回复

1楼2012-04-02 23:43:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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

CODE:

D0=3.55e-6;

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:5000;

s=sum(8./(pi*(2*k+1)).^2.*exp(-D*(pi*(2*k+1)).^2[color=red]*t[/color]/(d^2)));

F=sum((m-s').^2);

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