[매틀랩]수치해석 매틀랩 소스
- 최초 등록일
- 2006.06.21
- 최종 저작일
- 2006.05
- 10페이지/ 한컴오피스
- 가격 1,500원
소개글
토막 선형 보간함수 ,3차 spline 함수 ,사다리꼴 법칙과 Simpson rule을 사용하여 매틀랩을 풀어본다
목차
1. 다음 데이터를 고려해보자.
(a) 토막 선형 보간함수 l(x)를 구하라.
(b) 위 자료에 보간하고 비 매듭 경계조건(not-a-knot-boundary) 식 (5.78)의 s(z1)=f(z1), s(z2)=f(z2)을 만족시키는 3차 spline 함수 s(x)를 구하라. 이 경우 n=3, x1=1, x2=3,x3=5, z1=2, z2=4임을 주의하라. 1≤x≤5일 때, s(x)와 l(x)의 그래프를 그려라.
2. 사다리꼴 법칙과 Simpson rule을 사용하여 아래의 적분에 대해 n=2,4,8...512일 때, Tn(f)값의 표를 작성하여라. 또, 오차가 감소하는 비율을 구하라.
본문내용
(a) 토막 선형 보간함수 l(x)를 구하라.
l1x=((x1-x)*y0 +(x-x0)*y1)/(x1-x0)
l2x=((x2-x)*y1 +(x-x1)*y2)/(x2-x1)
l3x=((x3-x)*y2 +(x-x2)*y3)/(x3-x2)
l4x=((x4-x)*y3 +(x-x3)*y4)/(x4-x3)
>> clear all; %% initializing
>>
>> x0=1; y0=3; %% 위에 나타난 값을 대응시켜준다
>> x1=2; y1=1; %% 위에 나타난 값을 대응시켜준다
>> x2=3; y2=2; %% 위에 나타난 값을 대응시켜준다
>> x3=4; y3=3; %% 위에 나타난 값을 대응시켜준다
>> x4=5; y4=2; %% 위에 나타난 값을 대응시켜준다
>>
>> x=sym(`x`);
>> l1x=((x1-x)*y0+(x-x0)*y1)/(x1-x0); %% 책에 나오는 공식에 대입
>> l2x=((x2-x)*y1+(x-x1)*y2)/(x2-x1); %% 책에 나오는 공식에 대입
>> l3x=((x3-x)*y2+(x-x2)*y3)/(x3-x2); %% 책에 나오는 공식에 대입
>> l4x=((x4-x)*y3+(x-x3)*y4)/(x4-x3); %% 책에 나오는 공식에 대입
>>
>> xx=1:5; %% x값의 범위를 설정
>> y=[3 1 2 3 2] %% x값에 대응하는 y값 지정
참고 자료
없음