【答案】应助回帖

» 猜你喜欢

孩子确诊有中度注意力缺陷已经有12人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复
2026博士申请-功能高分子，水凝胶方向已经有6人回复
论文投稿，期刊推荐已经有4人回复
硕士和导师闹得不愉快已经有13人回复
请问2026国家基金面上项目会启动申2停1吗已经有5人回复
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？已经有3人回复

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

1楼2012-08-02 15:29:04

hqulhh

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 220.3
红花: 1
帖子: 26
在线: 40.9小时
虫号: 1263377
注册: 2011-04-12
专业: 零件成形制造

引用回帖:

5楼: Originally posted by csgt0 at 2012-08-02 17:21:37
你没发现你的式子里没有P么，大概这个情况，你再改改吧。需要注意不能出现让K=0

L0=1;
=ode45(@fun,0:100,);
L=y(:,1);
G=y(:,2);
delta=(L-L0)/L0
end

function dy=fun(t,y)
L0=1;
alpha=0.01;
delt ...

用plot画出G值不变，假设G初值为10，怎么回事？
假设各初值L0=100;alpha=14.6;deltaT=10;C=2;r=16;A=8;B=2;
v=0.28;D=1e-12;Q=315;R=8.31;
帮我看下。谢啦

10楼2012-08-03 21:12:30

查看全部 16 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
hqulhh: 金币+2, ★★★很有帮助 2012-08-02 16:54:13
jjdg: 金币+1, 感谢参与 2012-08-02 18:25:43

看这个方程就差不多，其中y(1)对应你的L，y(2)对应G

CODE:

[t,l]=ode45(@fun,[0 1],ones(2,1))

function dy=fun(t,y)

dy=zeros(2,1);

dy(1)=y(1)*t*y(2)^2;

dy(2)=5/y(2);

end

showmethemoney

2楼2012-08-02 16:29:29

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

★
jjdg: 金币+1, 感谢参与 2012-08-02 18:25:33

其他的参数都有，所以方程中所有L都用y(1), dL/dt用dy(1), G都用y(2),dG/dt用dy(2)代替就可以。

showmethemoney

3楼2012-08-02 16:32:27

hqulhh

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 220.3
红花: 1
帖子: 26
在线: 40.9小时
虫号: 1263377
注册: 2011-04-12
专业: 零件成形制造

引用回帖:

3楼: Originally posted by csgt0 at 2012-08-02 16:32:27
其他的参数都有，所以方程中所有L都用y(1), dL/dt用dy(1), G都用y(2),dG/dt用dy(2)代替就可以。

还有个Ρ是时间t的函数 Ρ=kt

4楼2012-08-02 16:53:28