• lake 2012/02/07 16:11

    방명록으로 질문하지말라고하셨는데...어디다 물어봐야할지 몰라서 이렇게 남겨요!!!

    2계 연랍상미분 해결

    m*x''=(P1-P2)*Ac-b*x'-k*x -1번식

    Vc/be*dP2/dt=sgn(P1-P2)*Cd*x*Lr*(2/ro*|P1-P2|)^1/2+Ac*x' -2번식

    여기서 P1=f/Ap
    f=([t1>=0.5]-[t1>=1]+[t1>=1.5]-[t1>=2]+[t1>=2.5]-[t1>=3]+[t1>=3.5]-[t1>=4]+[t1>=4.5])*2000
    이렇게 f를 주고 시간을 0 에서 5초로 풀려는데요...

    function dx=flow(t,x)
    global ap ac m b k Vc be cd ro Lr
    t1=0:0.001:t;
    n=length(t1);
    f=([t1>=0.5]-[t1>=1]+[t1>=1.5]-[t1>=2]+[t1>=2.5]-[t1>=3]+[t1>=3.5]-[t1>=4]+[t1>=4.5])*2000;
    p1=f/ap;
    dx=[((p1(n)-x(3))*ac-b*x(1)-k*x(2))/m; (sign(p1(n)-x(3))*cd*x(2)*Lr*sqrt(2/ro*abs(p1(n)-x(3)))+ac*x(1))*be/Vc; x(1)];

    이렇게하고 rk4차 풀이법으로 계산을 했는데 계속 값이 NAN 나오네요 ㅠㅠ 어떻게 잘못된거죠??
    도와주세용~~~

    • PinkWink 2012/02/08 12:52

      먼저 모르는 부분이 생기고 그를 해결하는 과정이 재미있으면서도 또 답답하기도 합니다. 그래서 저도 항상 여기저기 돌아다니게 되는데요.

      그런 마음을 알지만, 또 대답을 드리기가 너무 곤란합니다.

      먼저 주어진 질문에 답변을 하기 위해서
      보여주신 코드만 딱보고 제가 답변을 찾을 수 있으면 좋겠지만, 그렇지 못하면 오랫동안 고민할 수도 있기 때문입니다. 현재 직장인으로 점심시간이나 이른 아침에만 짬짬이 블로그를 하기때문인데요.

      현재 코드상에서 제가 답변을 드릴 여건이 되지 못해 죄송합니다.

    • lake 2012/02/08 17:40

      아니에요~ 답변해주신것만으로도 감사합니다..~~!

      ode45로 풀리다가 t=0.5초부터 에러가나더군요ㅠㅠ

      답변해주셔서 감사합니다 ㅠㅠ~!

    • PinkWink 2012/02/09 08:39

      도움이 되지 못해 안타까울 뿐이랍니다.

  • 푸샵 2012/02/04 07:39

    봄이 오는 걸 시샘하는지 동장군의 심술;;;이 이만저만 아니네요. ㅎㅎ
    한파에 무사하시죠? ^^
    그래도 봄은 오나봅니다. 벌써 입춘이니.......ㅋㅋ
    행복하고, 건강한 주말 보내시고 입춘대길(立春大吉)!! 하세요. (⌒▽⌒) 핑크윙크님!~

    • PinkWink 2012/02/05 20:10

      ㅎㅎ 무사하지 않아요.. 지독한 감기에....^^
      그래도 활기찬 한주를 또 보내야죠^^
      푸샵님도 화이팅입니다.^^

  • 푸샵 2012/02/02 07:57

    완전 춥네요. ㅎㅎ 갠적으론 좋지만. ㅋㅋ
    영하 16도에 체감 23도라니...추위 주의하세요. 핑크윙크님!~
    행복한 2월 시작하셨죠?

    • PinkWink 2012/02/02 08:43

      아이쿠.. 제가 요즘 정신없이 보내는 사이 시간이 흘렀네요
      언능언능 포스팅활동도 하고, 답방도 다니고 해야할텐데요^^
      좋은 2월 보내세용^^

  • Like a Pro 2012/01/31 11:41

    오랜만이에요. 이전에 골프어플 관계로 메일 드렸었는데... 기억하실지 모르겠네요.
    우연히 제품출시 소식을 보다가 다리부상 소식을 봤어요. 사실 저도 '11년 5월에 눈길도 아닌(^^) 타일 바닥에서 넘어져서 왼쪽 무릎뼈(슬개골)이 부러졌거든요 ㅠㅠ 아직도 정상은 아니구요. 올해 5월경에 철심을 빼야합니다.
    핑크윙크님의 다리 수술 사진은 저보다 좀 더 심하더군요. 많이 좋아 지셨는지 궁금합니다.
    작년에 출시한 어플은 이제 제가 직접 개발하고 있어요... (자세한 내용은 메일로 드릴께요)
    제 생각을 제가 직접 구현하는게 제일 좋겠더라구요.
    무료버전은 지금도 다운 받을 수 있구요. 유료버전은 지금 앱스토어 심사중입니다.
    http://itunes.apple.com/kr/app/like-a-pro/id494568731?mt=8#
    사용해 보시고 리뷰의견 주시면 많은 도움 될 것 같아요.
    늦었지만 새해 복 많이 받으시고, 더욱 건강하시기 바랍니다.

    • PinkWink 2012/02/01 07:04

      네 오랜만입니다
      좋은 소식을 또 가지고 오셨군요^^
      한번 사용해보겠습니다
      근데
      제가 골프의 골~자도 몰라서 ㅠㅠ 아무튼 사용해볼께요^^

  • 푸샵 2012/01/30 07:56

    1월도 벌써 다 지나가려고 하네요. ㅎㅎ
    이번주엔 한파 소식이 있군요...한주 내내 춥다는데...감기조심하세요. ^^
    이 추위가 지나면 봄이 조금씩 다가오지 않을까 하고 생각해보게 됩니다.
    행복하고, 활기찬 한주 시작하시고, 한파에 감기 조심하세요. (⌒▽⌒) 핑크윙크님!~

    • PinkWink 2012/01/30 08:28

      네.. 차 폐차하고, 처음 걸어봤는데.. 좀 걸으니까 추위는 금새 없어지더군요.. ㅋㅋㅋ
      (이 저질 체력같으니라고.. ㅎㅎㅎㅎ)

1 2 3 4 5 ... 91