- 세상의 모든 계산기 자유(질문) 게시판 질문 ()
[nspire] desolve에서 삼각함수 들어갔을 때 올바른 해가 나오지 않습니다 ㅠㅠ
-
- 2018.09.26 - 17:50 2018.09.25 - 22:21 4658 10
deSolve(60*x''+1500*x=90*sin(2.5*t) and x(0)=0 and x'(0)=0,t,x)
를 입력시
x=−0.04*sin(5.*t)-5242.88*sin(t)*(((sin(t))^(2)-0.25)*(cos(t))^(2)-0.0625)*((sin(0.5*t))^(2)+2.75)*(cos(0.5*t))^(13)+5242.88*cos(t)*(((sin(t))^(2)+0.25)*(cos(t))^(2)-0.3125)*sin(0.5*t)*((sin(0.5*t))^(2)+2.25)*(cos(0.5*t))^(12)-29491.2*sin(t)*(((sin(t))^(2)-0.25)*(cos(t))^(2)-0.0625)*((sin(0.5*t))^(2)-0.23611111111111)*(cos(0.5*t))^(9)+21626.88*cos(t)*(((sin(t))^(2)+0.25)*(cos(t))^(2)-0.3125)*sin(0.5*t)*((sin(0.5*t))^(2)-0.375)*(cos(0.5*t))^(8)+9216.*sin(t)*(((sin(t))^(2)-0.25)*(cos(t))^(2)-0.0625)*(cos(0.5*t))^(7)-4300.8*cos(t)*(((sin(t))^(2)+0.25)*(cos(t))^(2)-0.3125)*sin(0.5*t)*(cos(0.5*t))^(6)-1935.36*sin(t)*(((sin(t))^(2)-0.25)*(cos(t))^(2)-0.0625)*(cos(0.5*t))^(5)+645.12*cos(t)*(((sin(t))^(2)+0.25)*(cos(t))^(2)-0.3125)*sin(0.5*t)*(cos(0.5*t))^(4)+15.36*sin(t)*(((sin(t))^(2)-0.25)*(cos(t))^(2)-0.0625)*((sin(0.5*t))^(2)+11.916666666667)*(cos(0.5*t))^(3)-15.36*cos(t)*(((sin(t))^(2)+0.25)*(cos(t))^(2)-0.3125)*sin(0.5*t)*((sin(0.5*t))^(2)+2.0833333333333)*(cos(0.5*t))^(2)-9.6*sin(t)*(((sin(t))^(2)-0.25)*(cos(t))^(2)-0.0625)*cos(0.5*t)+1.28*cos(t)*(((sin(t))^(2)+0.25)*(cos(t))^(2)-0.3125)*sin(0.5*t)
와 같은 해가 나오네요
deslove(60x'' + 1500x=90sin2.5t and x(0)=0 and x'(0)=0 , t, x )
라는 식을 풀면 비제차 항에 대한 해는 깔끔하게 나오는데 제차 해는
sint 2.5t 에 관해 정리되는게 아니고 삼각함수 공식을 어떻게든 적용해 sint 항으로 정리하는 것 같더군요..
덕분에 답이 매우매우 더러워집니다 ㅠㅠ 맞는지도 모르겠구요.
Asin5t+Bcos5t+Csin2.5t+Dcos2.5t 이런 식으로 깔끔하게 나오는 방법이 있을까요 ?ㅠㅠ
exact로 표기해도 안되네요
우변에 삼각함수가 들어가면 어떤 때는 제대로 나오고 어떤때는 저런식으로 꼬여서 나옵니다 ㅠㅠ
-
1
댓글10
-
세상의모든계산기
2. 결과값을 살펴보면 삼각함수 형식으로 굉장히 길~~게 그리고 근사값으로 도출되었음을 확인할 수 있습니다.
삼각함수 결과가 깔끔하게 정리되지 못하는 현상은 TI-nspire 에서 비교적 흔한 상황입니다.
(이럴 때는 계산기를 약간 거들어주는 작업이 필요할 수 있습니다. 그래도 해결이 안될 수도 있구요)하지만 계수들이 강제로 근사값(approximate result) 로 계산되는 것은 흔하게 보여지는 것은 아닌데 이유는 잘 모르겠습니다.
(아마도 desolve와 결합하여 나타난 현상일 듯 한데 말이죠)
ㄴ 입력식에 2.5t 가 있어서 그랬을 것 같습니다. -
세상의모든계산기
3. (틀린 답 가지고 시간낭비하는 경우를 막기 위해서는) 답이 맞는지를 우선 확인해 보면 좋겠지요.
울프람알파에서 정답을 확인해 봅니다.
http://www.wolframalpha.com/input/?i=Solve+2*x%27%27(t)%2B50*x(t)%3D3*sin(2.5*t),+x(0)%3D0,+x%27(0)%3D0그래프 모양은 대략 이렇구요.
수식으로 표현하면 이렇다는군요.
x(t) = sin(5t)*(0.06 cos(2.5t) - 0.02 cos(7.5t) - 0.04) + (0.02 sin(7.5t) - 0.06 sin(2.5t))*cos(5t)1의 답(TI-nspire)과 3의 답(울프람 알파)을 각각 그래프로 그려 봅니다.
완벽하게 겹치는 걸로 봐선 TI-nspire 의 결과가 오답은 아님을 확인할 수 있습니다.
그러면 원래 긴~~~ 답에서 어떻게 짧은 답을 유도하느냐가 문제가 되겠네요. -
-
2018.09.26 - 15:24 #24094앗 네 답 맞아요!
공수1 에서는
해는 제차해+특수해로 표현되고
제차해(homogeneous solution)가 Acos5t+Bsin5t로 표현되며
특수해(particular solution)는
미정계수법에 의해
Ccos2.5t+Dsin2.5t 로만 표현될 수 있으니
결과는 이 네개의 항 조합으로만 나오는게 당연한데 계산기가 미방을 다른 알고리즘으로 푸나보네요 ㅠㅠ
rst기호부터는 이해가 안되어서 집에가서 찬찬히 따라하면서 이해해보겠습니다!
쉬는날인데도 시간써서 답변주셔서 너무 감사합니다
사진까지 친절하게 첨부해주시고 ㅠㅠㅠ
kwj**** 님의 최근 댓글