t=a:h:b;

2 次查看（过去 30 天）
Mos_bad
Mos_bad 2017-9-8
编辑： James Tursa 2017-9-9
I need to code an iteration for the below algorithm:
% a,h,b,s are constants.
t=a:h:b;y(0)=s;
f(k)(x,t)=t(k)+y(k);
y(k+1)=y(k)+h*f(k);
please help me how to use for loop for this iteration.
  1 个评论
James Tursa
James Tursa 2017-9-8
编辑：James Tursa 2017-9-8
I assume that the statement y(0)=s is meant to imply that the initial value of y at the initial value of t (which is "a" in this case) is "s"?
This is just numerical integration using the Euler method.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

James Tursa
James Tursa 2017-9-8
编辑：James Tursa 2017-9-9
An outline of the loop:
s = whatever
a = whatever
b = whatever
h = whatever
y = zeros(size(t));
y(1) = s;
for k=1:numel(t)-1
f = ________; % <-- you fill the blanks in based on your formula from above
y(k+1) = _______; % <-- you fill the blanks in based on your formula above
end
If you want to save each of the derivative values, then you can use f(k) in the loop instead of just f. (And also initialize f = zeros(size(t)) prior to the loop)

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by