MATLAB and Numerical Analysis

Size: px

Start display at page:

Download "MATLAB and Numerical Analysis"

재희 후
5 years ago
Views:

1 School of Mechanical Engineering Pusan National University

2 Review 무명함수 >> fun x^2 + y^2; % ff xx, yy = xx 2 + yy 2 >> fun(3,4) >> ans = 25 시작 x=x+1 If문 >> if a == b >> c=d; >> end 만약 a=b 라면, c=d로계산하시오 for문 >> for i = 1:5 >> x(i+1)=x(i)+1; >> end 순서도와동일 i=5 출력

Review nlinfit >> fun = @(p,x) I(1)*exp(-p*x); ff tt = II 0 ee μμtt >> t = [0:2:10]; >> I = [0.9726 0.8001 0.6582 0.

3 Review nlinfit >> fun I(1)*exp(-p*x); ff tt = II 0 ee μμtt >> t = [0:2:10]; >> I = [ ]; >> b0 = 1; >> mu = nlinfit(t,i,fun,b0); >> figure,plot(t,i, ob ); hold on; >> plot(t,fun(mu,t), r );

4 Bisection 추정된근의값 x r = x l + x 2 u

5 Bisection 다음의방법들을사용하여 ff xx = ee xx xx 의근을구하고, 참백분율상대오차를하나의그래프에함께그려비교하시오. 반복횟수에따른오차의특성을살피고, 각방법들간의차이점에대해논하시오. 참값으로 를사용. iteration 함수및초기치설정 xx tttttttt = ff xx = ee xx xx xx ll = 0, xx uu = 1 yy ll = ff(xx ll ), yy uu = ff(xx uu ) xxxx ii = xxxx ii + xxxx ii 2 yyyy ii = ff(xxxx ii ) 만약 yyyy ii yyyy ii < 0 이면, xxxx ii+1 = xxxx ii 만약 yyyy ii yyyy ii > 0 이면, xxxx ii+1 = xxxx ii 최종결과 xxxx eeeeee, eeee eeeeee 획득 xxxx ii+1 = xxxx ii xxxx ii+1 = xxxx ii eeee ii = 100 xx tttttttt xxxx ii xx tttttttt

6 Bisection coding

8 False position 다음의방법들을사용하여 ff xx = ee xx xx 의근을구하고, 참백분율상대오차를하나의그래프에함께그려비교하시오. 반복횟수에따른오차의특성을살피고, 각방법들간의차이점에대해논하시오. 참값으로 를사용. iteration 함수및초기치설정 xx tttttttt = ff xx = ee xx xx xx ll = 0, xx uu = 1 yy ll = ff(xx ll ), yy uu = ff(xx uu ) xxxx ii = xxxx ii yyyy ii xxxx ii xxxx ii yyyy ii yyyy ii yyyy ii = ff(xxxx ii ) 만약 yyyy ii yyyy ii < 0 이면, xxxx ii+1 = xxxx ii 만약 yyyy ii yyyy ii > 0 이면, xxxx ii+1 = xxxx ii 최종결과 xxxx eeeeee, eeee eeeeee 획득 xxxx ii+1 = xxxx ii xxxx ii+1 = xxxx ii eeee ii = 100 xx tttttttt xxxx ii xx tttttttt

9 Newton-Raphson f '( x i ) = f ( xi ) 0 x x i i+ 1 x i+ 1 추정된근의값 = x i f ( xi ) f '( x ) i

10 Newton-Raphson 다음의방법들을사용하여 ff xx = ee xx xx 의근을구하고, 참백분율상대오차를하나의그래프에함께그려비교하시오. 반복횟수에따른오차의특성을살피고, 각방법들간의차이점에대해논하시오. 참값으로 를사용. iteration 함수및초기치설정 xx tttttttt = ff xx = ee xx xx ff (xx) = ee xx 1 xx = 0 yy = ff(xx) yyy = fff(xx) eeee 1 = 100 xx tttttttt xx 1 xx tttttttt xx ii+1 = xx ii yy ii yyy 1 yy ii+1 = ff(xx ii+1 ) eeee ii = 100 xx tttttttt xx ii+1 xx tttttttt yyy ii+1 = fff(xx ii+1 ) 최종결과 xx eeeeee, eeee eeeeee 획득

11 Newton-Raphson coding

13 Fixed-point iteration 다음의방법들을사용하여 ff xx = ee xx xx 의근을구하고, 참백분율상대오차를하나의그래프에함께그려비교하시오. 반복횟수에따른오차의특성을살피고, 각방법들간의차이점에대해논하시오. 참값으로 를사용. iteration 함수및초기치설정 xx tttttttt = ff xx = ee xx xx gg(xx) = ee xx xx = 0 yy = gg(xx) eeee 1 = 100 xx tttttttt xx 1 xx tttttttt xx ii+1 = yy ii yy ii+1 = gg(xx ii+1 ) eeee ii = 100 xx tttttttt xx ii+1 xx tttttttt 최종결과 xx eeeeee, eeee eeeeee 획득

14 Secant 다음의방법들을사용하여 ff xx = ee xx xx 의근을구하고, 참백분율상대오차를하나의그래프에함께그려비교하시오. 반복횟수에따른오차의특성을살피고, 각방법들간의차이점에대해논하시오. 참값으로 를사용. iteration 함수및초기치설정 xx tttttttt = δδ = ff xx = ee xx xx xx 1 = 1 yyy = ff(xx 1 ) yyy = ff xx 1 + δδxx 1 eeee 1 = 100 xx tttttttt xx 1 xx tttttttt xx ii+1 = xx ii δδ xx ii yyy ii yyy ii yyy ii yyy ii+1 = ff(xx ii+1 ) yyy ii+1 = ff xx ii+1 + δδ xx ii+1 eeee ii = 100 xx tttttttt xx ii+1 xx tttttttt 최종결과 xx eeeeee, eeee eeeeee 획득

15 Structure & Save Structure Mat 파일로저장

16 Graph plot 1) Load 2) Plotting

17 Graph plot 3) Axis control

18 Graph plot 4) Labeling

19 Graph plot 5) Legend 6) save

20 Graph plot Line style: -, --, -., : Color: k(black), b(blue), r(red), m(magenta), g(green) Mat 파일불러오기 Y 축 log-scale 로 plot 그래프선굵기 한그래프에여러 line plot 격자 Line 의형태및색 x & y 축범위설정 Legend ( 범례 ) x 축범위설정및해당값설정 x 축설명 타이틀 y 축설명 저장

<32303032BEC7BFECC1F62E687770>

<32303032BEC7BFECC1F62E687770> 2001-2002 전남대학교 산악회 2002년 제12집 먼 훗날 당신이 찾으시면 그때에 그 말이 잊었노라. 당신이 속으로 나무라면 무척 그리다 잊었노라. 오늘도 어제도 아니 잊고 먼 훗날 그때도 잊었노라. 지도교수 인사말 산은 인생의 도장이라는 말이 있습니다. 산은 無言 속에서도 自然 의 理致 를 가르쳐주고, 또한 많은 智慧 를 줍니다. 그러나 生活