Questions about bvp4c solution

조회 수: 7 (최근 30일)
Elia Paini
Elia Paini 2021년 4월 26일
댓글: Elia Paini 2021년 4월 29일
Hi, I solved correctly a system of ODE second order with BCs, with the instruction bvp4c.
However, I have some questions about it.
My code is as follows:
...
xmesh = linspace(1e-5,1.5*1e-3,20);
solinit = bvpinit(xmesh,1e-3*[0,0,0,0,1,1,1,1]);
sol = bvp4c(@bvpfcn,@bcs,solinit);
plot(sol.x,sol.y)
1) When I plot the solution, I should choose in x interval sol.x or xmesh? What is their difference?
2) I don't understand how Matlab uses the second vector [0,0,0,0,1,1,1,1] in bvpinit.
It should represent the initial guess of solution (I have 8 variables), but it seems that its values don't change final solution. Is it right?
I've already read Matlab guide about bvp4c, but I have not found these answers.
Thank you!

채택된 답변

Shadaab Siddiqie
Shadaab Siddiqie 2021년 4월 29일
From my understanding you have few doutes on the arguments of bvp4c function.
  1. 'bvp4c' function not output 'sol.stats' when 'NMax' is exceeded, This is a due to what the 'stats' object returns. Typically, it has values that are associated with a solved problem, rather than an incomplete one.
  2. solinit are the Initial guess of solution, specified as a structure. Use bvpinit to create solinit. Unlike initial value problems, a boundary value problem can have no solution, a finite number of solutions, or infinitely many solutions. An important part of the process of solving a BVP is providing a guess for the required solution. The quality of this guess can be critical for the solver performance and even for a successful computation. For some guidelines on creating a good initial guess, see Initial Guess of Solution. But it is not nessary for the solution to change if the initial conditions are changed.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Boundary Value Problems에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by