Milne格式MATLAB,matlab考试题

)

)(,(12

h k x y h x f k m m ++=)(2

)(211k k h x y y m

m ++=+?????))

(,(1m m x y x f k =5. 设计M 文件计算:

x=0:0.1:10

当sum>1000时停止运算，并显示求和结果及计算次数。

clear

x=0:0.1:10; s=1; sumx=0;

while (sumx<=1000)

sumx=sumx+x(s)*x(s)-2*x(s); s=s+1; end s sumx

二、实验内容：

M 文件设计：仿真时间0～1秒，计算步长h=0.1，初始值y(0)=1，对连续系统

y x dx

dy

+= 采用以下方法建立仿真模型： ① 前向Euler: n n n h f y y +=+1

② 后向Euler: ③ 梯形法： ④ 改进Euler:

⑤ 经典RK4：

)

2(100

2

i i i x x sum -=∑=1

1+++=n n n hf y y ()

112

++++=n n n n f f h

y y ()321151623----+-+

=n n n n n f f f h

y y ))

(,(1m m x y x f k =)2

,2(12k h

y h x f k m m ++

=()43211226

k k k k h

y y m m ++++

=+????

?????

)

,(34h k y h x f k m m ++=)2

,2(23k h y h x f k m m ++

=