Scribd
Upload
4K views

Shooting Method For Boundary Value Problem Solving

This document discusses using the shooting method to solve boundary value problems (BVPs) in MATLAB. It presents the shooting method as a continuous method for solving BVPs numerically by converting them into initial value problems. An example is provided of using the Newton-Raphson method to solve the system of nonlinear equations that arise in the shooting method formulation. Plots of the solution are generated to verify that the BVP is accurately solved.

Uploaded by

minnntheino
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
4K views

Shooting Method For Boundary Value Problem Solving

This document discusses using the shooting method to solve boundary value problems (BVPs) in MATLAB. It presents the shooting method as a continuous method for solving BVPs numerically by converting them into initial value problems. An example is provided of using the Newton-Raphson method to solve the system of nonlinear equations that arise in the shooting method formulation. Plots of the solution are generated to verify that the BVP is accurately solved.

Uploaded by

minnntheino
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

Shooting Method for Boundary Value Problem Solving

Boundary Value Problem Solving method Shooting


method Continuous Method Compare
matlab bvp4c, bvp5c

Propblem solve problem

ODE BVP Solve

..

Example

to=0, tf=1, y1o=y2=1, y1f=2, y2f=3;

the system of nonlinear equations

function F = dEqs(t,y) % Differential equations.


F = zeros(4,1);
F(1) = y(3)^2;
F(2) = sqrt(y(4))-y(1);
F(3) = y(2)+y(1);
f(4) = y(1)-y(3);
Initial Con

function y = inCond(u) % Initial conditions; u(1)


y = [1 1 u(1) u(2)]; % and u(2) are unknowns.
Initial condition ODE Integrate the system of nonlinear
equations

function r = residual(u) % Bounday residuals.


global XSTART XSTOP
r = zeros(length(u),1);
x = XSTART;
[xSol,ySol] = ode45(@dEqs,[x XSTOP],inCond(u));
[m,n] = size(ySol);
r(1) = ySol(m,1)-2;
r(2) = ySol(m,2)-3;

function u(1) and u(2) equations က r(1) and r(2)

two unknowns and two equations solving the system of


nonlinear equations solving system of nonlinear equations F(x)=0
Matlab fsolve function fsolve
built-in trust region method ပါ။ built-in algorithm
options Gauss-newton or Levenberg-Marquardt change
root finding
newton-raphson method newton-raphson method
Advanced Engineering mathematics
Newton-rapson method matlab

function root = newtonRaphson2(func,x,tol)


% Newton-Raphson method of finding a root of simultaneous
% equations fi(x1,x2,...,xn) = 0, i = 1,2,...,n.
% USAGE: root = newtonRaphson2(func,x,tol)
% INPUT:
% func = handle of function that returns[f1,f2,...,fn].
% x = starting solution vector [x1,x2,...,xn].
% tol = error tolerance (default is 1.0e4*eps).
% OUTPUT:
% root = solution vector.
if nargin == 2; tol = 1.0e4*eps; end
if size(x,1) == 1; x = x'; end % x must be column vector
for i = 1:150
[jac,f0] = jacobian(func,x);
if sqrt(dot(f0,f0)/length(x)) < tol
root = x; return
end
dx = jac\(-f0);
%dx=-1e-3;
x = x + dx;
if sqrt(dot(dx,dx)/length(x)) < tol*max(abs(x),1.0)
root = x; return
end
end
error('Too many iterations')

function [jac,f0] = jacobian(func,x)


% Returns the Jacobian matrix and f(x).
h = 1.0e-4;
n = length(x);
jac = zeros(n);
f0 = feval(func,x)
for i =1:n
temp = x(i);
x(i) = temp + h;
f1 = feval(func,x);
x(i) = temp;
jac(:,i) = (f1 - f0)/h;
end

newton-raphson jacobian function finite difference method


.. Newton-raphson function two
unknowns two equations function

function shoot4nl
% Shooting method for nonlinear 4th-order boundary
% value problem in Example.
global XSTART XSTOP % Make these params. global.
XSTART = 0; XSTOP = 1; % Range of integration.
u = [-1 1]; % Trial values of u(1)
% and u(2).
x = XSTART;
u = newtonRaphson2(@residual,u)

function run y3 and y4 initial value

u=
-1.3777
11.4519

y3 and y4
Boundary Value Problem

[xSol,ySol] = ode45(@dEqs,[x XSTOP],inCond(u));


y1=ySol(:,1);
y2=ySol(:,2);
plot(xSol,y1)
grid on
figure
plot(xSol,y2)

Fig(1) change of y1

ႈ
ႈ
fig(2).change of y2
fig(1) and fig(2)
Newton-raphson method matlab က fsolve function
u
u=
-1.2250 11.0197
y3 and y4 initial ODE integrate

Fig(3).change of y1

Fig(4).change of y4

fig(3) and fig(4)


solving system of nonlinear equations
system of nonlinear equations

(1) steepest desent method

(2) bisection method

(3) newton method

(4) newton-raphson method

(5) gauss-newton method

(6) quasi-newton method

(7) Levenberg-Marquardt method

(8) dogleg method


problem method
matlab fsolve
strong
convergence Myanmar Astronomy
presentation for Shanghai BVP convergence

initial sensitive Problem post


က continuous method convergence
core2duo prosessor, RAM 4Gb Run
solve အ run run
method powell’s hybrid method

bluephoenix
‘
’
‘
’

‘ ’ post
။

Mtssnrty

You might also like