기초전산유체역학 과제1
- 최초 등록일
- 2021.05.17
- 최종 저작일
- 2018.03
- 4페이지/ 압축파일
- 가격 1,000원
소개글
"기초전산유체역학 과제1"에 대한 내용입니다.
목차
1. Finished the attached tridiagonal.m
2. Use the tridiagonal.m. Modify the assignment01.m if necessary.
1) k = 0.2 f(0) = 1,f(L) = 6
2) k = 0.3 df/dx(0) = -1, f(L) = 4
3) k = 0.3 f(0) = 2, df/dx(L) = 0
4) k = 0.4 df/dx(0) = -1, df/dx(L) = 1 (If this causes error, why?)
본문내용
1. Finished the attached tridiagonal.m
% Tridiagonal Solver
function u = tridiagonal(a, b, c, d)
n = length(a);
a(1) = 0; c(n) = 0;
P = zeros(size(a));
Q = zeros(size(a));
u = zeros(size(a));
% Finish the two lines below
P(1) = -c(1)/b(1); Q(1) = d(1)/b(1);
% Finish the for loop
for i = 2:n
P(i) = -c(i)/(b(i)+a(i)*P(i-1));
Q(i) = (d(i)-a(i)*Q(i-1))/(b(i)+a(i)*P(i-1));
end
% Finish the line below
u(n) = Q(n);
% Finish the for loop
j=0;
for j = (n - 1):(-1):1
u(j)=P(j)*u(j+1)+Q(j);
end
참고 자료
없음
압축파일 내 파일목록
assignment 1.hwp
assignment 1.txt
tridiagonal.txt