Question

MATLAB迭代算法问题

dy=2.112*exp(19.837-13636/y)-0.005*y+3.201，
dy为y的增量，每次算一次都要叠加上去。
y的初值是640，终值为670。
求出迭代次数及依次的解。

Answer 1

y =

640.0000
640.4870
640.9794
641.4776
641.9816
642.4917
643.0081
643.5309
644.0605
644.5970
645.1408
645.6920
646.2509
646.8178
647.3930
647.9769
648.5697
649.1719
649.7837
650.4055
651.0379
651.6811
652.3356
653.0021
653.6808
654.3724
655.0775
655.7967
656.5306
657.2799
658.0454
658.8279
659.6281
660.4471
661.2858
662.1452
663.0265
663.9308
664.8596
665.8143
666.7963
667.8074
668.8494
669.9244
671.0345
迭代次数 i=45 （初始为640算作第一次）
代码：
y=zeros(100,1);
y(1,1)=640;
i=1
while i<100 & y(i,1)<=670
y(i+1,1)=y(i,1)+2.112*exp(19.837-13636/y(i,1))-0.005*y(i,1)+3.201
i=i+1
end

Answer 2

clear all;
clc;
y=640;
i=1;
while y<=670;
if ceil(y)==670 % 判断y是否等于670，是跳出while循环，否就需执行
return;
else
dy=2.112*exp(19.837-13636/y)-0.005*y+3.201;
y=y+dy;
mat(i)=dy; % 存储每一次的dy值
i=i+1; % i的终值为迭代次数
end
end

追问

要附带输出各阶段的y值，这点怎么办？

追答

在y=y+dy；下一步加一个 ymat(i)=y；