본문 바로가기

수치해석

【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)..
【MATLAB】 Manning 방정식 개수로 행렬문제 풀이 - chapra의 응용수치해석 3판 문제 2.10 [문제] Manning 방정식은 직사각형 단면의 개수로 내에서 물의 속도를 계산하는데 사용된다. 여기서 U는 속도(m/s), S는 수로의 기울기, n은 거칠기 계수, B는 폭(m), 그리고 H는 깊이(m)이다. 다음의 데이터는 5개의 수로에서 얻은 것이다. 거칠기 계수 n 수로의 기울기 S 수로의 폭 B 수로의 깊이 H 0.035 0.0001 10 2 0.020 0.0002 8 1 0.015 0.0010 20 1.5 0.030 0.0007 24 3 0.022 0.0003 15 2.5 이 값을 행렬에 저장하는데, 각 행은 각 수로를 표현하고 각 열은 각 매개변수를 표현하도록 하라. 매개변수 행렬 내의 값을 기초로 속도의 열벡터를 계산하는 MATLAB 문장을..
【MATLAB】 간단한 수학계산 문제풀이 [ MATLAB 간단한 수학계산 문제풀이 ] 1번. 2번. matlab에서 곱하기 기호는 * 입니다. 3번. 나누기 기호는 슬래쉬(/) 입니다. 역슬래쉬는 반대로 나눠줍니다. 4번. 괄호안에 4+3을 먼저 수행하여 7이되고 더하기보다 곱하기가 먼저이므로 2*7을 해서 14가 되고 마지막으로 3+14를 통해 답이 17이 나옵니다. 5번. 곱하기와 나누기는 우선순위가 같기때문에 좌측부터 수행이 됩니다ㅏ. 6번. 소수점계산도 가능합니다. 7번. 8번. 더하기보다 지수승이 먼저 계산됩니다. 9번. sqrt()는 루트(root)를 붙여주는 기능입니다. 괄호안에 숫자가 루트 내에 들어가게 됩니다. help sqrt하게 되면 위와 같이 'square root'라는 설명을 보실 수 있습니다. 10번. matlab에서는..
【MATLAB】 Hooke의 법칙, 포텐셜에너지 문제풀이 (for, max) [문제] 힘 F(N)이 스프링을 압축하기 위해 작용할 때, 변위 x(m)는 Hooke의 법칙에 의해 다음과 같이 표현할 수 있다. 여기서 k는 스프링 상수(N/m)이다. 스프링에 저장된 포텐셜에너지 U(J)는 다음의 식으로 계산할 수 있다. 5개의 스프링으로 실험을 하여 다음과 같은 데이터를 수집하였다. 힘 F(N) 14 18 8 9 13 변위 x(m) 0.013 0.020 0.009 0.010 0.012 MATLAB을 이용하여 F와 x를 벡터로 저장하라. 그리고 스프링 상수와 포텐셜에너지의 벡터를 계산하라. max 함수를 이용하여 최대 포텐셜에너지를 구하라. [문제풀이] 우선, 위 표를 벡터화해보면 다음과 같습니다. 다음으로 스프링 상수를 구해보겠습니다. 5개의 스프링을 사용했다고 했으니 위 F와 x의..
【MATLAB】 표준정규확률 밀도함수 그래프 문제풀이 [문제] 표준정규확률 밀도함수는 다음과 같이 종 모양의 곡선으로 나타낼 수 있다. MATLAB을 이용하여 z가 -5에서 5까지 이 함수의 그래프를 그려라. 수평축은 z, 수직축은 주파수(frequency)로 놓아라. [문제풀이] 우선, 위 함수를 matlab에서 이해할 수 있는 수식으로 적어주는 것이 중요한 듯 합니다. 위에서 2번줄이 문제에서 알려준 공식을 적은 것입니다. 루트는 sqrt() 명령어를 사용하고 지수함수 e는 exp() 명령어를 통해 작성했습니다. 그리고 z가 행렬이기때문에 요소끼리 계산하기위에 바로 뒤에 콤마(.)를 붙여주었습니다. 그리고 1번 줄, 2번 줄 뒤에 각각 세미콜론(;)을 붙여주어 값이 보이지 않게 하였습니다. 위와 같이 작성한 후, 실행(F5)를 눌러주면, 위와 같은 그래..