sympy로 미분방정식 푸는 방법

Question

정확한 문제는 말씀 못드리고, 대신에 함수나 식을 조금씩 바꿔서 설명드리자면fx&#43;fxy&#43;fyy&#61;0과 같은 partial 미분 방정식이 있고 f가 a*sin(b*x&#43;c*y)와 같은 꼴로 주어졌을 때 계수 a, b, c를 찾는 문제입니다 a,b,c&#61;symbols(&#34;a b c&#34;) x,y&#61;symbols(&#34;x y&#34;) f&#61;a*sin(b*x&#43;c*y) eq&#61;Eq(f.diff(x)&#43;f.diff(x).diff(y)&#43;f.diff(y,2),0) solve(eq,[a,b,c],dict&#61;True)을 했지만 계속 실행만 되고 결과가 나오지 않습니다조언해주시면 감사하겠습니다!

anydog · Answer

실행만 되고 결과가 나오지 않는다면 컴퓨터 연산속도가 느린 겁니다. 정의하신 f와 eq대로라면 solver 결과는 프로그램 돌려볼 것도 없이 a&#61;0입니다.이런걸 원하시는 건 아닐테니 a를 지우고 다시 연산해보면 에러가 발생합니다. a, b, c &#61; sp.symbols(&#34;a b c&#34;) x, y &#61; sp.symbols(&#34;x y&#34;) f &#61; sp.sin(b*x &#43; c*y) eq &#61; sp.Eq(f.diff(x) &#43; f.diff(x).diff(y) &#43; f.diff(y,2), 0) result &#61;sp.solve(eq, [b, c]) > No algorithms are implemented to solve equation -b*c*sin(b*x &#43; c*y) &#43; b*cos(b*x &#43; c*y) - c**2*sin(b*x &#43; c*y)참고로 위 에러는 방정식에 해가 없거나 해를 기호로 표현이 불가능할 때 발생합니다. (eg. x-cos(x)&#61;0)

초보보 · Answer

감사합니다 저 경우는 제가 막 지어낸 식이라 해가 없나봐요

초보보 · Answer

한가지 더 질문드려도 될까요? 원래 문제는 sech가 들어간 더 복잡한 식이었고 solve를 실행했을 때 계속 돌아가기만 해서 나름대로 식을 정리한 후 다시 solve를 해보니까 잘 되더라구요.그런데 식 정리과정에서 (tanh2, A2 등등은 A^2같은 제곱 표현입니다)이 나왔는데, 분명 1-tanh^2를 인수로 가짐에도 1-tanh^2로 나눠지지도 않고(/) collect 되지도 않아서 그냥 제가 직접 나눈 식을 새로 입력했습니다. 왜 이러는지 혹시 아시나요?

anydog · Answer

무슨 의미로 말씀하신건지 잘 모르겠습니다.이 경우에도 원 질문처럼 solver가 공통인수인 1-tanh2(...) &#61; 0을 풀었을 것 같네요.코드를 질문주실 땐 가급적 코드를 올려주세요.복붙이 안되니 재현할 마음이 안생기네요.

초보보 · Answer

네 답변 감사합니다!

sympy로 미분방정식 푸는 방법

5개의 답변

#equation 목록