매트랩1과제물

소개글

매트랩1과제물에 대한 보고서 자료입니다.

목차

문제 :　f(x)=1/x 에 대하여 x=3.5 에서의 1계 및 2계도함수 (f'(3.5),f''(3.5))를 Numerical Differentiation 으로 구하라.


문제 :　다음 비선형 방정식의 x>0의 구간에 존재하는 2개의 해를 3가지 방법(Bisection method, New-ton-Rapson method, secant method ) 로 구하라 。

본문내용

ts required\')
end
if nargin<4|isempty(es)
es=0.0001;
end
if nargin<5|isempty(maxit)
maxit=50;
end
iter=0;
while(1)
xrold=xr;
xr=xr-func(xr)/dfunc(xr);
iter=iter+1;
if xr~=0, ea=abs((xr-xrold)/xr); end
if ea<=es|iter>=maxit,break, end
root1(end+1)=xr;
ea1(end+1)=ea;
iter1(end+1)=iter;
end
end
★secant. m 파일
function [root1,ea1,iter1]=secant(func,xl,xu,es,maxit,varargin)
root1=[];
ea1=[];
iter1=[];
x=[];
if nargin<3
error(\'at least 3 input arguments required\')
end
if nargin<4|isempty(es)
es=0.0001;
end
if nargin<5|isempty(maxit)
maxit=50;
end
iter=0;
x(end+1)=xl;
x(end+1)=xu;
while(1)
xrold=x(end);
x(end+1)=x(end)-func(x(end))/dfunc(x(end),x(end-1),func);
iter=iter+1;
if x(end)~=0, ea=abs((x(end)-xrold)/x(end)); end
if ea<=es|iter>=maxit,break, end
root1(end+1)=x(end);
ea1(end+1)=ea;
iter1(end+1)=iter;
end
end
function u= dfunc(y,yold,func1)
u=(func1(yold)-func1(y))/(yold-y);
end
function [root1,ea1,iter1]=secant(func,xl,xu,es,maxit,varargin)
root1=[];
ea1=[];
iter1=[];
x=[];
if nargin<3
error(\'at least 3 input arguments required\')
end
if nargin<4|isempty(es)
es=0.0001;
end
if nargin<5|isempty(maxit)
maxit=50;
end
iter=0;
x(end+1)=xl;
x(end+1)=xu;
while(1)
xrold=x(end);
x(end+1)=x(end)-func(x(end))/dfunc(x(end),x(end-1),func);
iter=iter+1;
if x(end)~=0, ea=abs((x(end)-xrold)/x(end)); end
if ea<=es|iter>=maxit,break, end
root1(end+1)=x(end);
ea1(end+1)=ea;
iter1(end+1)=iter;
end
end
function u= dfunc(y,yold,func1)
u=(func1(yold)-func1(y))/(yold-y);
end
★결과값

키워드

매트랩1과제물,   매트랩,   비선형 방정식,   2계도함수