Question

这是Euler在MATLAB上的显式方法的代码

function [t,x] = meuler(f, intervalo, x0, N)
h = (intervalo(2)-intervalo(1))/N;
t = intervalo(1):h:intervalo(2);
x(:,1) = x0(:);
for i = 1:N
    x(i+1,:) = x(i,:) + h*(f(t(i),x(i,:)));
end

它需要在另一个文件中定义的函数f。当我尝试使用以下函数运行代码时

function f = funccorazon(t,x)
f1 = x(2);
f2 = 16*x(1)+4*sin(2*t);
f=[f1;f2];

我收到此错误

>> meuler(@funccorazon, [0 2*pi], [0 2], 1000)
Attempted to access y(2); index out of bounds because numel(y)=1.

Error in funccorazon (line 2)
f1 = y(2);

我不知道为什么。显然，当我使用ode45求解微分方程时，似乎没有任何问题。任何帮助都将不胜感激。谢谢！

Answer 1

错误的原因是您需要使用行向量而不是列向量。

首先，在函数调用中，将输入x0定义为行向量：

meuler(@funccorazon, [0 2*pi], [0,2], 1000)

第二，在funccorazon函数中，类似地定义输出f：

f=[f1,f2];

最后，在meuler函数中进行以下更改：

x = zeros(N,2);
x(1,:) = x0;

Euler在MATLAB上的显式方法通过示例引发错误

1 个答案: