matlab四五阶龙格-库塔方法计算示例

3、有一初始温度为1400K 的铁球投入到无限大的水池中,水温为310K ,假设铁

球温度随时间的变化规律为下列表达式,试用四五阶 R-K (龙格-库塔)方法计算铁球在35min ,95min 的温度。

]310)[K/min -(T -0.04e 310)-0.001(T dt

dT

步骤1:重新打开或直接输入clear

matlab四五阶龙格-库塔方法计算示例

步骤2:创建m 文件 左上角File>>New>>M-File

matlab四五阶龙格-库塔方法计算示例

步骤3:

在m文件编辑器里输入function temp =temp1(t,T)

temp=-0.04*exp(0.001*(T-310))*(T-310)

matlab四五阶龙格-库塔方法计算示例

步骤4:保存m文件,注意文件名是函数名.m File>>Save Workspace As

matlab四五阶龙格-库塔方法计算示例

步骤5:保存m文件,注意文件名是<函数名>.m 不然等会调不出来函数。

matlab四五阶龙格-库塔方法计算示例

步骤6:在matlab 右边命令栏输入

[t,T]=ode45('temp1',0:5:95,1400)

plot(t,T)

xlabel('time t0=0 tt=95');

ylabel('T values T(0)=1400');

matlab四五阶龙格-库塔方法计算示例

步骤7:回车运行就会出来结果了

matlab四五阶龙格-库塔方法计算示例

T t

1.0e+003 *

1.4 0

0.989 5

0.787 10

0.663 15

0.5788 20

0.5187 25

0.474 30

0.4401 35

0.4139 40

0.3935 45

0.3773 50

0.3644 55

0.3541 60

0.3458 65

0.3391 70

0.3337 75

0.3293 80

0.3258 85

0.3229 90

0.3205 95

相关文档
  • 龙格库塔法matlab

  • 龙格库塔法matlab程序

  • 四阶龙格库塔法

相关推荐: