e^ax sinbx 함수의 부정적분 구하는 방법과 그 확장(선형대수학을 이용한 적분방법)
보통 이 적분은 부분적분법을 써서 구하는 방법이 일반적입니다.
그런데 다음과 같이 이 함수를 직접 적분하지 않고 부정적분을 구하는 방법이 있습니다.
e^ax sinbx를 먼저 미분해보죠.
e^ax(asinbx+bcosbx) 가 되겠죠? 이번엔 e^ax cosbx를 미분해봅시다.
e^ax(acosbx-bsinbx) 네요?
아하! 둘다 sinbx, cosbx를 가지고 있네요! 그렇다면 이것들을 적당한 수를 곱해서 더해주면 e^ax sinbx를 얻을 수 있을겁니다.
그러니까,
e^ax(Asinbx+Bcosbx)를 미분한 식이 e^ax sinbx가 되는 A, B를 찾을 수 있다는 겁니다.
음 위 식의 도함수는 e^ax { (aA -bB)sinbx + (aB+bA)cosbx } 이니까 이것이 e^ax sinbx와 같기 위해서는,
aA-bB=1, aB+bA=0이어야 합니다. 이건 일차방정식이죠.
(a -b)(A)=(1)
(b a)(B) (0)
(괄호가 위아래로 이어져있습니다.)
그럼 위 행렬의 역행렬을 구해서 양변에 곱해주면 A, B를 구할 수 있어요
그 값은 A=a/(a^2+b^2), B=-b/(a^2+b^2) 겠죠? 한편 e^axcosbx는 위 일차방정식의 우변 벡터가 0, 1인 경우니까 이 경우에는 A=b/(a^2+b^2), B=a/(a^2+b^2)가 됩니다.
따라서 e^ax sinbx의 부정적분은 e^ax/(a^2+b^2) * (asinbx-bcosbx)입니다. 그리고 e^ax cosbx의 부정적분은 e^ax/(a^2+b^2) * (bsinbx+acosbx)입니다.
적분을 미분을 통해서 구한셈이죠.
e^ax sinbx coscx 라는 복잡한 함수의 적분을 구할 때에도 위와 같은 방법을 이용할 수 있습니다.
이때는 물론 4*4행렬의 역행렬을 구해야 겠죠.
원리가 궁금하신분들은 아래 내용을 참조하세요.
이 글의 요지는 sin, cos, exponential의 x의 계수가 서로 다른 일차곱으로 이루어진 함수는 선형대수학을 이용하면 얼마든지 부정적분을 구할 수 있다는 이야깁니다.
미분은 선형사상입니다. 한편 선형사상f의 image에 있는 원소 b의 역상은 f(a)=b인 a에 대해 a+Ker(f) 입니다. 우리가 부정적분을 하는 과정은 바로 a를 구하는 과정이었죠. (미분 operator의 kernel은 상수함수이므로 a+Ker(f)는 부정적분 + 적분상수가 되는 것입니다)
위 적분법의 아이디어는 T-invariant subspace(사상 T에 대해 불변인 부분공간)의 개념을 응용한 것입니다.
exponential, cos, sin과 같은 함수의 곱들을 기저로 갖는 벡터공간을 생각하자는 겁니다. 얘네들은 미분을 해도 결국엔 같은 집합에 속하기 마련이죠. 이 부분공간안에서 미분사상은 다름아닌 bijection이 됩니다. (1이 기저로 있지 않으므로 상수함수가 이 공간안에 존재할 수 없고, 따라서 미분함수는 injection. surjective는, dimension이 같으므로...) 종합하자면 미분사상이 기저를 어디로 보내는지를 통해서 "가역"행렬을 직접 구할 수 있으므로 원하는 함수의 부정적분을 역행렬을 곱해서 구할 수 있다는 이야기입니다.
여기서 질문(ㅋㅋㅋ...)
역행렬을 어떻게하면 쉽게 구할 수 있을까요?
0 XDK (+0)
유익한 글을 읽었다면 작성자에게 XDK를 선물하세요.
-
왤케동안이지
-
막 의치한 아닌 일반학과애들이 신기하게보나요 의치로 안간 수능 엄청 못본사람들도...
-
슬슬 먼가 전조가 보인다
-
이제 그만 봐야겠따
-
똥테 달성! 4
영광굴비는 식상하고 이 영광의 역광을 그대들에게!
-
지인선님.. 2
제가 이렇게 뻘글만 쓰는 놈이어도 팔로우 계속 해주시다니..
-
누가 위라고 생각
-
???
-
19명 뽑는 경외시 문과라인 최종 94명 지원했습니다 점공 현재 기준 37명 중...
-
김지영이나 매가에 또선생인가 그분도 괜찮다고 하는 거 같은데 그래서 누구 들을까
-
음
-
나 카리나랑 맞팔이야! 13
오르비에서...응...
-
이륙했군 0
굿
-
사실 내가 메디컬에 절박하지 않은 것도 있지만 솔직히 지역인재 없애면 지방의대에...
-
아으 그낭 좀 해라 ㅡㅡ
-
나랑 잘 맞는 시험지로 제발...
-
영광굴비! 11
다들 고마워요! 이제 오르비를 떠날 때가...
-
[속보] 합참 “北 미사일 1100여㎞ 비행 후 탄착” 1
속보 [헤럴드경제=강문규 기자] 합참 “北 미사일 1100여㎞ 비행 후 탄착”
-
사촌이나 만날까 2
나보다 전역 한 달 빨리해서 비슷할텐데
-
물2 공부한 지 2일 차인데, 10번대 문제들이 풀림. 확실히 다른 화1이나 생명1 보단 할 만함.
-
서울대병원 산부인과·응급의학과 '0명'…빅5 전임의 지원 저조 1
국내 '빅5' 병원 전임의 1차 채용 결과 지원자가 모집 인원의 절반에 못 미친...
-
미적에 이골이 난 모태 이과로 밝혀져...
-
ㄹㅇ임
-
중앙대 약학대학 신입생 카페 가입 안내 안녕하세요, 중앙대학교 약학대학 제 41대...
-
본좌는 심심하다 6
빨리 재밌는 뻘글을 싸보거라
-
8명 뽑는데 65명 지원함. 현재 23명 들어왔고 12등임. 최근 3년동안 27명...
-
10억이 잊혀지질 안네 10
10억.. 10억.....
-
점공 필요성을 못느끼는 쌉실수인거임? 아니면 스나허수인거임?
-
2칸 스나 0
24명 뽑고 103명 지원 점공 23/40 연고 공대인데 가능성 있을까요? .....
-
지금 확통 쎈b 푸는데 기초실력 기르기에 좋은것같아서 수1수2도 쎈b 사서 풀까요?...
-
아 똥 존나쌌다 11
시원
-
김승리 커리 타는 중입니다 화작 내신이 1학기에 있어서 화작을 겨울방학에 할...
-
화장실 갔다오는데 둘 중 한 명이 siuuuuu 외치면서 스카에서 나오다가 날...
-
그래도 21수능 물2처럼 블랭크는 안떴잖아 한잔해~
-
옯 이미지 상 선택과목 29
서로 서로 써줘봅시다! 은근 재미있다구!
-
ㄹㅇ
-
육개장 1
-
좀 그래
-
어지간해선 인강 보지 말고 독학+현강하셈 인강수강생들 생각보다 수학 꽤 못해서...
-
당장 버리길 열심히 공부했으면 7-80퍼센트는 늘 먹던맛, 당해년도 6.9변형,...
-
ㄴ몰 너무 어려워서 몰티져스땡김
-
일은 시발같이 시키면서 월급 하루 이틀 매 달 쳐밀리네 개새끼들이
-
2시간전에 헌혈해씀
-
라는 미친 생각을 해보았서요.
-
. 2
배고픔
-
옆자리 대화 개어지럽네 대단한 인생을 사시는듯
-
마! 상남자 아잉교!
-
다닐수도 있을 것 같아서요..!!
슈뢰딩거님답게
수리 게시판에는 금지된 수학문제 질문이 아닌 글을 쓰시면서
또 수학적으로 좋은 글을 써주셨군요!
제가 방금전에 이 문제가지고 다시 생각해봤는데
cos,sin 함수들이 n번 곱해진 함수공간에서의 미분행렬을 점화식으로 찾은것 같아요!
먼저 위의 삼각함수들이 n번 곱해진 함수공간(이걸 n-TFP space(n-삼각함수곱공간)라고 부를게요 편의상...)의 기저는 2^n개입니다. 이 기저의 순서를 잘 매겨서 {v_i,n}이라 합시다. 그리고 x의 계수들을 r_i로 매길게요.
그러면 이제 n+1-TFP space의 기저를 잡는 방법이 문제가 되는데,
v_i,n+1=v_i,n*sin(r_n+1x) (i<=2^n 일때)
v_i,n+1=v_i-2^n,n*cos(r_n+1x) (i>2^n 일때)로 정의하면,
T_n(n-TFP space에서의 미분사상)에 대해서
T_n+1(v_i,n+1)=(T_nv_i,n)sin(r_n+1x)+r_n+1v_i+2^n,n+1 (i<=2^n 일때)
=(T_nv_i-2^n,n)sin(r_n+1x)+r_n+1v_i-2^n,n+1 (i>2^n)가 됩니다.
따라서 n에 관한 미분행렬 T의 점화식을 다음과 같이 얻습니다.
T_n+1= T_n -r_n+1*I_n(n by n 항등행렬)
r_n+1*I_n T_n
역행렬 점화식 까지 찾아냈다면 더 기분좋을 것 같은데 생각대로 안되네요 그건...
차라리 삼각함수 덧셈정리를 n번 써서 전개 하는게 나으려나 -_-...
T_2 구해서 sinaxsinbx 시험삼아 적분해봤더니 -a/(a^2-b^2)cosaxsinbx +b/(a^2-b^2)sinaxcosbx 가 나오네요