matlab怎么画出dp/dt=(207p^0.17-P)/2.7的相图？

Question

区间p从0-20

Answer 1

根据问题的已知条件，对于此微分方程，通过变形可以得到dt/dp形式的方程，然后用dsolve函数解答其解析值，最后用plot函数绘出其相位图。求解过程如下：

>> syms t(p)

>> Dt=diff(t);

>> t=dsolve(Dt==2.7/(207*p^0.17-p),t(0)==0)；%给出的问题缺初始条件

>> p=0:20;

>> t=eval(t);

>> plot(t,p)

>> xlabel('t'), ylabel('p')

运行结果