数学建模社区-数学中国

标题: 求教maple画图 [打印本页]

作者: succin    时间: 2010-2-28 22:56
标题: 求教maple画图
求教关于动态分岔的,方程:diff(z(t), t) = (tau(t)+I)*(z(t)-tau(t)), diff(tau(t), t) = epsilon' p4 X7 d, M" ?8 p- u
/ i5 n$ B2 @/ Q# C: Z
其中z(t), tau(t)都是复数,I是虚数i
% `: P8 B  T1 \# `" ?$ f3 L; s5 d, B& d2 ]3 u! I# E
想画纵轴是z(t)的实部, 横轴是tau(t)的图。
8 C" R1 k$ {$ f* w, I9 M0 p2 M( Y- ]" w* Z' V$ G
各位大侠,谁能帮忙, 先在此谢过,有劳了!
作者: snowy1991    时间: 2010-9-3 21:28
有教程。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
作者: haviet    时间: 2011-9-13 20:50
禁用词语   i like the chapter MATH
作者: cwxiong20098    时间: 2012-11-4 23:31
相当于是参数方程了,plot([z(t),tau(t)],t=-1..1);好像是这样的,具体记不清了,maple帮助里有的,查查就知道了
作者: zsy312    时间: 2012-11-9 19:47

  1. ; Y6 ?$ z: |( S* p7 }
  2. sol := dsolve([diff(z(t), t) = (tau(t)+I)*(z(t)-tau(t)), diff(tau(t), t) = epsilon, z(0) = a, tau(0) = b], [z(t), tau(t)]);
    ( z, s" ]! t0 U  F1 c) r9 o# \" |
  3. soli := eval(sol, [a = .1, b = .2, epsilon = 0.1e-1]);
    ! L7 S) U. c) v' e! w
  4. plot([rhs(soli[1]), Re(rhs(soli[2])), t = 0 .. 10]);
    $ S/ Y3 L# ?; ?" b+ x5 A5 _. |
  5. + p& V! O  m& w; Z6 }

  6. 1 f5 W. n) a: r/ ?' _/ K" [
复制代码





欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5