Question

在matlab中，想要画dx/dt的图，已知dx/dt=0.8x(1-x),0<x<1,0<t<25,请问在matlab中怎么输入？要过程详细！

Answer 1

h=dsolve('Dx=0.8*x*(1-x)','x(0)=2')%假设t=0时x=2,即‘x(0)=2'没有初始条件，x(t)函数有无穷个

ezplot(h,[0 1])
ylabel('x')

追问

你写的好像有问题，输的结果显示error，结果图应该是这样的，虽然我有图，但是我需要matlab运行过程。

追答

clc

clear

close all

for j=0.15:0.05:0.95

x0=j;

[t,x]=ode45(@(t,x)0.8*x*(1-x),[0 25],x0);

plot(t,x);

hold on

end

xlabel('t')

ylabel('x')

title('dx/dt=0.8x(1-x)')

或者：准确解非上述解析解

clc

clear

close all

for j=0.15:0.05:0.95

h=dsolve('Dx=0.8*x*(1-x)',['x(0)=',num2str(j)]);

ezplot(h)

hold on

end

xlabel('t(0<t<25)')

ylabel('x(0<x<1)')

title('dx/dt=0.8x(1-x) (0<t<25)')

axis([0 25 0 1])

Answer 2

x0=0.5;%x(0)的值，自己输
[t,x]=ode45(@(t,x)0.8*x*(1-x),[0 25],x0);
plot(t,x)