|
수정된 할선법
clear all
clc
f='1/3*(sin(x).^3)+1/2*(cos(x).^2)-x'
n=0.1; err=0.000001; x0=0.1; q=0.0001; num=1;
fprintf(' 횟수 x(n) x(n+1) 오차 \n');
while num>=err
n=n+1;
x=x0; f0=eval(f);
x=x0+q*x0; f1=eval(f);
x1=x0-(f1*(q*x0)/(f1-f0));
num=abs(x1-x0);
fprintf('%2.0f %2.8f %2.8f %2.8f \n', n,
|
|
|
<수정된 할선법 m파일 코드>
function root = modsec(f,xrNew,offage,es,maxit)
%modsec()함수는 수정된 할선법(modified secant method)을 위한 함수이다.
%입력으로는
% f=근을 구하고자하는 대상이되는 방정식
% xrNew = 초기 가정값
% offage = 변동량
% es=
|
|
|
할선법을 사용하라.
(세 번 반복계산하고, 초기 가정값은 x_i-1 = 2.5 와 x_i = 3.5를 사용하라)
⒟ 수정된 할선법을 사용하라.
(세 번 반복계산하고, 초기 가정값은 x_i = 3.5 와 δ = 0.01을 사용하라)
결론
수학으로 푼 점과 수치해석적으로 푼 관점이
|
|
|
|
초기값 X0, 변동량 δ, 한계추정오차 ε, 주어진 함수 f(x)를 input으로 갖는 다음의 프로그램 작성
→함수 f(x)의 근을 구하는 수정된 할선법을 이용해서 근을 구하는 프로그램 작성, 단 X0,δ,ε는 프로그램 실행시 값을 변경할수 있어야함.
|
|
메뉴펼치기
