기구 문제의 수치해석 풀이

소개글

응용동역학 수업때 작성한 4절기구의 위치해석에 관한 리포트 입니다.
해석은 매트랩을 이용하여 직접 m-file로 작성한 뉴튼랩슨법을 썼습니다.
주석을 정말 자세하게 썼기 때문에 현재 매트랩수업이나 수치해석 수업을 듣고 계신분에게도 큰 도움이 될것입니다.
작성한 m-file을 사용하는 방법도 예시로 정말 자세하게 썼고 그 결과도 넣었습니다.
비선형은 물론 선형도 그냥 풀립니다.

목차

1. M-file 작성
2. 결과
3. 결론

본문내용

본문에 있는 m-file의 첫부분입니다.

function root=NRV(eq, x, xi, es, maxit)
%뉴튼랩슨벡터법
%eq는 연립방정식벡터와 해벡터로써 sym(`[;;]`)으로 표현(syms는 에러), inline은 pretty, jacobian같은 대수연산이 안됨.
%x는 변수로써 x1,x2 문자를 써야함.!!<<-- 해결 과제... -_-
%xi는 초기해값으로써 [;;]으로 표현, 라디안 값.
%es는 희망하는 해의 수렴 정도 (ex 0.01)
%maxit는 희망하는 반복 횟수 (ex 3)
if nargin<5, maxit=50;end %함수내의 입력값 수가 5개 미만이면 maxit를 자동으로 50으로 설정
if nargin<4, es=0.001;end %함수내의 입력값 수가 4개 미만이면 es를 자동으로 0.1%로 설정(그래도 maxit=50 을 넘어설수는 없음)
iter=0; %현재 반복횟수를 0으로 초기화
ea=0; %아래 if 구문에서 xi=0 일때 ea=0 이 적용되게끔함

.
.
.
.
이것을 실행시키는 방법은 다음과 같고
>>NRV(sym(`[5*cos(pi/3)+9*cos(x1)+7*cos(x2)-10;5*sin(pi/3)+9*sin(x1)+
7*sin(x2)]`), sym(`[x1;x2]`),[20/180*pi;250/180*pi])

결과는
[ - 15/2 + 9 cos(x1) + 7 cos(x2) ]
[ ]
[ 1/2 ]
[5/2 3 + 9 sin(x1) + 7 sin(x2)]
iter
4

ea
4.0748e-006

root
0.2907
4.5513

degree
16.6538
260.7725

처럼 나옵니다.


이런식으로 한줄 한줄 주석을 넣었습니다.
수치해석, 매트랩, 응용 동역학, 기구학 수업에 큰 도움이 될것 입니다.

참고 자료

없음