본문 바로가기

2017/09/22

【MATLAB】 분수 스칼라 연산 문제풀이 스칼라 연산에서 분수 계산문제 풀이입니다. [문제] 아래의 두 식 1)과 2)를 MATLAB으로 계산하시오. 1) 2) [문제풀이] 1) 수학적인 개념을 잘 이해하고 matlab에서 작성해서 계산해야함을 유의합니다. 우선 1)번식은 아래와 같이 작성하고 계산할 수 있습니다. 2) 위 문제보다 조금더 복잡하지만 아래와 같이 작성해서 계산하면 됩니다. 특별한 설명보다 자세히 보실 수 있도록 띄어쓰기로 구분해뒀으니 참고하시길 바랍니다. 이상!
【MATLAB】 화학반응 Arrhenius 방정식 문제풀이 - Chapra의 응용수치해석 3판 문제 2.20 [문제] 화학반응의 온도에 대한 의존성은 다음과 같은 Arrhenius 방정식으로 계산할 수 있다. 여기서 k는 반응속도(s^-1), A는 진동(또는 주파수) 인자, E는 활성에너지(J/mol), R은 기체상수[8.314 J/(mole*K)], 그리고 Ta는 절대온도(K)이다. 혼합물은 E = 1 x 10^5 J/mol 이며 A = 7 x 10^16 이다. MATLAB을 사용하여 온도 253에서 325K까지 반응속도의 값을 구하라. subplot을 사용하여 (a) Ta에 대한 k(녹색 선)와, (b) 1/Ta에 대한 log10k(붉은색 선)의 그래프를 나란히 그려라. semilogy 함수를 적용하여 (b)를 구하라. 두 subplot의 축에 라벨과 제목을..
【MATLAB】 문제 2.18 풀이 - Chapra의 응용수치해석 3판 문제 2.18 [문제] 다음의 MATLAB 명령어를 입력하면 무엇이 화면에 출력되는가? (a) (b) [문제풀이] (a) 순서대로 풀이하면 이렇습니다. 먼저, x에 5를 대입합니다. 그리고 세제곱을 해주어서 결과를 냅니다. 그리고 y = 8-x에 x를 넣어주어 결과를 8 - 5 이므로 y는 3이 됩니다. 여기서 중간에 쓰여진 x^3은 y=8-x 식에서 사용되지 않습니다. 그 이유는 변수가 이용되지 않았다는 것과 두번째줄의 결과는 세미콜론을 붙이지 않으면 알 수 있듯이 ans = 125로 나올 뿐입니다. 위 상태가 세미콜론(;)을 붙이지 않고 계산해본 결과입니다. 위에서 알 수 있듯이 x^3을 하면 ans변수에 125라는 값이 저장되는 것입니다. 만약에 y=8-ans라고..
【MATLAB】 매개변수 방정식, 3차원 원추형 나선 그래프 문제풀이 - Chapra의 응용수치해석 3판 문제 2.17 [문제] 다음의 매개변수 방정식은 원추형 나선(conical helix)을 나타낸다. 범위 t = 0에서 6π까지 증분 Δt = π/64로 x, y 및 z의 값을 계산하라. subplot을 사용하여 위쪽 구획에 (x,y)의 2차원 그래프(붉은색 실선)를, 아래쪽 구획에 (x,y,z)의 3차원 그래프(청록색 실선)를 그려라. 두 그래프 모두 축에 라벨을 표시하라. [문제풀이] 우선, t가 0부터 6π까지 증분을 Δt = π/64으로 아래와 같이 2가지 표현중 하나를 사용하면 됩니다. 1) 콜론 표기법 2) linspace() 저는 이번문제에서 좀더 표시하기 쉬운 1번 콜론 표기법으로 계속 풀어보겠습니다. 다만 2번 방법으로 하려면 마지막에 사이 갯수를 입력..
【MATLAB】 Maclaurim 급수전개로 cosine 문제풀이 - Chapra의 응용수치해석 3판 문제 2.15 [문제] Maclaurim 급수전개로 cosine을 나타내면 다음과 같다. MATLAB을 이용하여 의 항이 포함될 때까지의 급수전개에 대한 그래프(검정색 점선)와 함께 cosine의 그래프(실선)를 그려라. 급수전개를 계산하는데 내장함수 factorial을 사용하라. x축의 범위를 0에서 3π/2까지로 놓아라. [문제풀이] 우선, 쉬운 cosine 함수부터 그려보겠습니다. 먼저, 가장 위에 clear 명령으로 혹시모를 변수중첩을 방지했습니다. 그리고 x값은 linspace(x1,x2)를 이용해서 0부터 3π/2까지 범위로 100개의 인자를 사용했습니다. linspace의 인자 갯수의 defualt는 100개입니다. 그리고 5번째줄에서 y1 = cos(x)..