【答案】应助回帖

» 猜你喜欢

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

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

csgt0

荣誉版主 (著名写手)

彩色挂图

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

哦，时间太久都忘了。
matlab不用预先声明和定义的。但是A还是可以定义下，以防出现rou超过0~1的范围出现A为空值而出错的情况。
delta=y(1)-L0/L0;
对每个y就有一个delta，然后得到一个rou，从而计算出dy，符合你的方程。
你把它放在计算rou之前就行了。

showmethemoney

16楼2012-08-09 10:41:18

查看全部 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