본문 바로가기

행렬2

[Julia] 기본1 - 6. 행렬식 해 구하기 링크의 내용을 공부하며 제멋대로 번역한 내용입니다. 기본편은 이 글이 마지막 글이며, 링크에 보간법, 수치적분법, FFT 관련 내용도 있으니 참고바랍니다. 수치해석을 하려면 이제 행렬을 푸는 방법을 알아야한다. 얼마나 빠르고 정확하게 행렬식의 해를 구할 수 있는가는 수치해석의 처음과 끝이다. FEM을 풀던 FDM, FVM, 등 뭐를 풀던 똑같다. EXAMPLE PROB. LAPLACIAN EQU. 예제 라플라시안 방정식 자 격자를 M+1개로 나눈다. 즉 포인트는 M+2개가 되겠다. (0,1,2..., M, M+1) 그리고 Uxx를 2차 중앙 차분을 하면 아래와 같다. 각 경계조건은 알파와 베타이니까, u0 = alpha, u_(m+1) = beta로 놓고 행렬식을 정리하면 다음과 같다. 자 여기서 h는 .. 2021. 10. 11.
[Julia] 기본1 - 5. 행렬 가지고 놀기! 링크의 내용을 공부하며 제멋대로 번역한 내용입니다. 일단 수치해석을 하려면 사용하는 언어의 행렬사용법을 찐하게 알고 있어야한다. 효율이 생명인데 행렬(어레이)자체를 비효율적으로 계산하고 있으면 도루묵이기 때문이다. INITIALIZING ARRAYS 행렬정의! 직접 정의는 여느 프로그래밍언어와 비슷하다. A = [1 2 3; 1 2 4; 2 2 2] A = [1 2 3; 1 2 4; 2 2 2] 요런 식이다. 근데 백터랑 행렬이랑 좀 객체 유형이 달라진단다. 아래 예시를 보자 V1 = [4.0, 5, 6] # 3-element VECTOR (Float64) V2 = [4.0; 5; 6]# 3-element VECTOR (Float64) M1 = [4.0 5 6]# 1x3 Matrix(Float64) 콤마.. 2021. 10. 11.