经过上面的处理后，利用matlab的ode45来求解常微分方程
建立导函数文件fun.m，内容为：
function dzdt=fun(t,z)
a=1;b=-2;c=3;d=-4;m=5;%这些参数我随意设的，根据需要自己修改
dzdt=[a*z(1)*(1-z(1)/m)+b*z(1)*z(2);
c*z(2)+d*z(1)*z(2)];
设定积分区间
TSpan=[0,10]';%自己确定
设定初值
Z0 = [1,-1]';%自定
用ode45求解
[T,Z]=ode45(@fun,TSpan,Z0);
作图
plot(T,Z(:,1),T,Z(:,2))
蓝色的是x(t)，绿色的是y(t)
跪谢大神，太感谢了！
大神你能再帮我看看 知道 y,x,t数据，怎么用matlab优化出参数不？？跪谢跪谢啊！急用啊！！
什么叫优化出参数？不是太懂
...dx\/dt=ax(1-x\/m)+bxy dy\/dt=cy+dxy 要画出x y 分别与t关系啊 求 ...
function dzdt=fun(t,z)a=1;b=-2;c=3;d=-4;m=5;%这些参数我随意设的，根据需要自己修改 dzdt=[a*z(1)*(1-z(1)\/m)+b*z(1)*z(2);c*z(2)+d*z(1)*z(2)];设定积分区间 TSpan=[0,10]';%自己确定 设定初值 Z0 = [1,-1]';%自定 用ode45求解 [T,Z]=ode45(@fun...
用matlab确定dx\/dt=ax-bxy,dy\/dt=-cy+dxy
额 ：[y(t)=C1*exp(-c*t),x(t)=0,[diff(x(t),$(t,2))=-(b*x(t)^3*m*t-e*x(t)^2-x(t)^2*d*diff(x(t),t)+x(t)^3*d*a+x(t)^3*d*e*t+x(t)*c*diff(x(t),t)-x(t)^2*c*a-x(t)^2*c*e*t-diff(x(t),t)^2)\/x(t),y(t)=-(diff(x(t...
微分方程组 dx\/dt=-axy dy\/dt=-bxy 和 dx\/dt=-ay dy\/dt=-bx 如何求解...
如下
微分方程组 (1) dx\/dt=-axy dy\/dt=-bxy 和 (2) dx\/dt=-ay dy\/dt=-bx...
(1) dx\/dt=-axy dy\/dt=-bxy 那么dy\/dx=b\/a 解得y=(b\/a)x+C1 (下面是用t表示）dx\/dt=-ax((b\/a)x+C1)=-bx^2-aC1x dx\/x(bx+aC1)=-dt,解得(1\/aC1)lnx-(b\/aC1)ln(ax+aC1)=-t+lnC2 y(t)=(b\/a)x(t)+C1 其中x(t)由上面的隐函数决定）(2) dx\/dt=-...
全微分方程是什么?
3. 食饵-掠食者的生态奥秘: Lotka-Volterra方程 $\\frac{dx}{dt}=ax-bxy,\\frac{dy}{dt}=-cy+dxy$，虽然是非线性，但它描述了生态系统的动态平衡，尽管解析解复杂，但数值方法却能揭示其微妙变化。4. 热传导的温度分布: $\\frac{\\partial u}{\\partial t}=k\\frac{\\partial^2u}{\\partial x^...
已知一阶二元微分方程组dx\/dt=ax-bxy;dy\/dt=cx+dxy;并且已知多组t,x...
经过上面的处理后，利用matlab的ode45来求解常微分方程 建立导函数文件fun.m，内容为：function dzdt=fun(t,z)a=1;b=-2;c=3;d=-4;m=5;%这些参数我随意设的，根据需要自己修改 dzdt=[a*z(1)*(1-z(1)\/m)+b*z(1)*z(2);c*z(2)+d*z(1)*z(2)];设定积分区间 TSpan=[0,10...
