본문 바로가기

컴퓨터

MATLAB을 이용한 알기 쉬운 수치해석) 3장 6번 미분 함수에서 2*sqrt(x)가 아닌 sqrt(2*x)로 해야한다. 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 29.
MATLAB을 이용한 알기 쉬운 수치해석) 3장 5번 %main.m i = newton("current", "diff_current", 1, 0.01); disp(i); %current.m function f = current(t) f = 9*exp(-t)*sin(2*pi*t) -3.5; end %diff_current.m function f = diff_current(t) f = -9*exp(-t)*sin(2*pi*t) + 9 *exp(-t) * cos(2*pi*t) * 2*pi; end 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 29.
MATLAB을 이용한 알기 쉬운 수치해석) 3장 3번 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 28.
MATLAB을 이용한 알기 쉬운 수치해석) 3장 2번 x0 = 0; cnt = 0; while (1) x1 = x0 - feval("func", x0) / feval("diff_func", x0); x0 = x1; cnt = cnt + 1; fprintf("%d회 결과 : %.3f \n", cnt, x0); if cnt == 3 break; end end function f = func(x) f = exp(x) - cos(x) + x^3 - 4*x - 2; end function f = diff_func(x) f = exp(x) + sin(x) + 3*x^2 - 4; end 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 28.
MATLAB을 이용한 알기 쉬운 수치해석) 2장 8번 x = 3; zero = func(1); first = func(1) + first_diff_func(1) * (x-1); second = func(1) + first_diff_func(1) *(x-1) + second_diff_func(1) *(x-1)/2; third = func(1) + first_diff_func(1) + second_diff_func(1)/2 + third_diff_func(1) *(x-1)/6; fourth = func(1) + first_diff_func(1) + second_diff_func(1)/2 + third_diff_func(1) *(x-1)/6 + fourth_diff_func(1) * (x-1)/24; rel = func(x); zero_error = abs((r.. 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 28.
MATLAB을 이용한 알기 쉬운 수치해석) 2장 7번 % x - 1 = 1 이므로, zero = func(1); first = func(1) + first_diff_func(1); second = func(1) + first_diff_func(1) + second_diff_func(1)/2; third = func(1) + first_diff_func(1) + second_diff_func(1)/2 + 25; rel = func(2); first_error = abs((rel-first)/rel*100); second_error = abs((rel-second)/rel*100); third_error = abs((rel-third)/rel*100); fprintf("%f %f %f %f %f \n", zero, first, second, third, rel.. 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 28.
MATLAB을 이용한 알기 쉬운 수치해석) 2장 6번 x = pi/6; first = x; second = x; third = x - x^3/6; first_error = abs((sin(x) - first)/sin(x)) * 100; second_error = abs((sin(x)-second)/sin(x)) * 100; third_error = abs((sin(x)-third)/sin(x)) * 100; fprintf("1차 근사값 : %f \n", first_error); fprintf("2차 근사값 : %f \n", second_error); fprintf("3차 근사값 : %f \n", third_error); 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 28.
MATLAB을 이용한 알기 쉬운 수치해석) 2장 5번 rel_result = func(1.37); err_result = floor(rel_result * 100) / 100; fprintf("오차 : %f %% \n", (rel_result-err_result)/rel_result*100); trans_result = func2(1.37); err_result = floor(trans_result * 100) / 100; fprintf("변형한 오차 : %f %% \n", (trans_result-err_result)/trans_result*100); function f = func(x) f = x^3 - 7 * x^2 + 8*x - 0.35; end function f = func2(x) f = ((x-7)*x+8)*x-0.35; end 차이가 없는데... 컴퓨터/MATLAB을 이용한 알기 쉬운 수치해석 2021. 5. 28.

티스토리툴바