微分方程稳定性matlab程序
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
例子:
构建函数
function dx=stability(t,x)
dx=zeros(2,1);
dx(1)=-2+x(2)+4*(x(1)^2+x(2)^2)*x(1); dx(2)=x(1)+2*(x(1)^2+x(2)^2)*(2-x(2));执行以下代码:
[t,x]=ode45('stability',[0 6],[4 8]); plot(t,x(:,1),'r-');hold on;
plot(t,x(:,2),'bo-');
从图像可知(0,2)是不稳定的平衡点。