It takes too long to run the program and in the end I get an error, could you tell me what the error is?

Question

Freddy 2023-12-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2054719-it-takes-too-long-to-run-the-program-and-in-the-end-i-get-an-error-could-you-tell-me-what-the-error

评论： Torsten 2023-12-1

function f=funcionrr(x)

xini=[0 4*pi];

t=[0 2*pi];

options=odeset('RelTol',1e-13,'AbsTol',1e-13);

[ts,xs]=ode45(@campo_pendulo,t,xini,options);

f=xs(end,1);

end

function px = campo_pendulo(t,x)

px(1)= x(2);

px(2)=sin(x(1))+0.5*sin(t);

px=px';

% Definir los límites de la búsqueda

x_min = 0;

x_max = 2*pi;

x = (x_min + x_max)/2;

while abs(px) > 1e-6

if px > 0

x_max = x;

else

x_min = x;

end

x = (x_min + x_max)/2;

end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Torsten 2023-12-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2054719-it-takes-too-long-to-run-the-program-and-in-the-end-i-get-an-error-could-you-tell-me-what-the-error#answer_1363889

编辑：Torsten 2023-12-1

在 MATLAB Online 中打开

You can't change the input for x in function "campo_pendulo".

So you can delete this part of your code because it doesn't influence the results:

% Definir los límites de la búsqueda
x_min = 0;
x_max = 2*pi;
x = (x_min + x_max)/2;
while abs(px) > 1e-6
    if px > 0
        x_max = x;
    else
        x_min = x;
    end
    x = (x_min + x_max)/2;
end

Further it doesn't make sense to use a while loop here because px is not changed within the loop. So once abs(px) > 1e-6, you will enter the loop, but never exit.

f = funcionrr()

f = 82.5899

function f=funcionrr()