need plot a graph using Euler's method for both numerical and analytical

Question

Mpho Mokgotlo 2022-10-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1838848-need-plot-a-graph-using-euler-s-method-for-both-numerical-and-analytical

编辑： Jan 2022-10-30

% Euler's Method
% Initial conditions and setup
h = (10);  % step size
x = (0):h:(60);  % the range of x
y = zeros(size(x));  % allocate the result y
y(1) = (2);  % the initial y value
n = numel(y);  % the number of y values
for i=1:n-1
    pi = 3.1416
    a = 2
    r = 1
    f = a + (2 - 0.4*a) * x / r^2 * pi
    y(i+1) = y(i) + h * f;
end
plot (x,y,'b')
grid on

4 个评论
显示 2更早的评论隐藏 2更早的评论

John D'Errico 2022-10-29

编辑：John D'Errico 2022-10-29

在 MATLAB Online 中打开

Why do you think you need to define a value for pi? And worse, the value you chose is not even very accurate? Of course you are using Euler's method, with an extremely coarse step size. So 5 digits fo pi is probably adequate. But learn tol leave pi alone. It already exists. So there is no need to define it as a variable.

format long g
pi
ans = 
          3.14159265358979

And that is the correct value out to that many digits.

But you seem to be plotting the numrical solution. So what is your problem? What is the differential equation? Can you integrate that analytically? Surely, since f is LINEAR in x, you know the integral, since the integral of x is x^2/2? Or you could use dsolve in the symbolic toolbox?

Do you need to know how to add a second plot to the existing figure?

help hold
 HOLD   Hold current graph
    HOLD ON holds the current plot and all axis properties, including 
    the current color and linestyle, so that subsequent graphing commands
    add to the existing graph without resetting the color and linestyle.
    HOLD OFF returns to the default mode whereby PLOT commands erase 
    the previous plots and reset all axis properties before drawing 
    new plots.
    HOLD, by itself, toggles the hold state.
    HOLD does not affect axis autoranging properties.
 
    HOLD ALL is the same as HOLD ON. This syntax will be removed in 
    a future release. Use HOLD ON instead.
 
    HOLD(AX,...) applies the command to the Axes object AX.
 
    Algorithm note:
    HOLD ON sets the NextPlot property of the current figure and
    axes to "add".
    HOLD OFF sets the NextPlot property of the current axes to
    "replace".
 
    See also ISHOLD, NEWPLOT, FIGURE, AXES.

    Documentation for hold
       doc hold

Mpho Mokgotlo 2022-10-30

Thank you for the help. tht ting is when I try to plot the graph it gives out error message like : unable to perform assignment because the left and right sides have a different number of element, so I was asking what should I do or change and how

Jan 2022-10-30

@Mpho Mokgotlo: If you mention an error, post a copy of the complete error message. Seeing some parts of the message only is less useful, e.g. the readers do not know, in which line the error occurs.

请先登录，再进行评论。

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

Answer 1

Alan Stevens 2022-10-30

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1838848-need-plot-a-graph-using-euler-s-method-for-both-numerical-and-analytical#answer_1087408

编辑：Jan 2022-10-30

在 MATLAB Online 中打开

f = a + (2 - 0.4*a) * x(i) / r^2 * pi;  % x(i) not just x.
%  Also take note of the other comments above.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jan 2022-10-30

I've moved this to an answer, because it solves a problem.

请先登录，再进行评论。

Answer 2

Jan 2022-10-30

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1838848-need-plot-a-graph-using-euler-s-method-for-both-numerical-and-analytical#answer_1087403

在 MATLAB Online 中打开

h = 10;  % step size
x = 0:h:60;  % the range of x
y = zeros(size(x));  % allocate the result y
y(1) = 2;  % the initial y value
n = numel(y);  % the number of y values
a = 2;
r = 1;
for i = 1:n-1
    f      = a + (2 - 0.4 * a) * x(i) / r^2 * pi;
    %                             ^^^
    y(i+1) = y(i) + h * f;  % I guess this line fails
end
plot (x,y,'b')
grid on

A cleaned version of your code. Alan Stevens' hint is included (please accept his answer).

You can solve such problems using the debugger. Type in the command window:

dbstop if error

and run the code again. Matlab stops at the problem and you can check the values of the used variables.

size(y(i+1))         % A scalar, of course
size(y(i) + h * f)   % Not a scalar
size(f)              % Same size as x!

You do not want x in the line to get the value of f, but only the element x(i).

Do you really mean: ... / r^2 * pi or .../ (r^2 * pi) ?

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

请先登录，再进行评论。

need plot a graph using Euler's method for both numerical and analytical

4 个评论
显示 2更早的评论隐藏 2更早的评论

回答（2 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

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

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

need plot a graph using Euler's method for both numerical and analytical

4 个评论 显示 2更早的评论隐藏 2更早的评论

回答（2 个）

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

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

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

WeChat

4 个评论
显示 2更早的评论隐藏 2更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

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