clzm.net
当前位置:首页 >> MAtlAB解非线性微分方程组的问题 >>

MAtlAB解非线性微分方程组的问题

没解析解你肯定得用数值解。以下是步骤。 第一步,在Matlab里写m-file,命名微分方程对应的函数: function f=myfun(t,x) f = [-x(1)+x(2)^2; -2*x(2)+x(1)^2]; 这里的x(1)就是你这儿的x, x(2)就是你这儿的y。写完后自动保存,不要对保存目录做...

fun=@(t,x)[36*(x(2)-x(1));-x(1)*x(3)-16*x(1)+28*x(2)-x(4);... x(1)*x(2)-3*x(3);x(1)+0.2] [t,X] = ode45(fun,[0:0.01:20],[0 0 0 0]); plot(t,X)

x y z 分别用 y(1) y(2) y(3) 表示; x' y' z' 分别用 y(4) y(5) y(6) 表示。 dy = @(t,y) [y(4);y(5);y(6);... y(4).*y(4).*y(1)/2 + 2*y(5);... y(5).*y(5).*y(2)/2-2*(y(4)+y(6));... y(6).*y(6).*y(3)/2 + y(5)] y0 = [0 0 0 0.1 0.2 0.3]; ...

参考http://zhidao.baidu.com/question/479644084.html?oldq=1 回答的一样的问题,还有下面的程序 function dy=odefun(t,y) %y(1) x %y(2) x' %y(3) y %y(4) y' %y(5) z %y(6) z' u=9.8*1; dy=zeros(6,1); dy(1)=y(2); dy(2)=-u*y(1)/sqrt(y(1)^...

第一个的t的步长可不是你后来写的以1为步进哦,所以就是你第一个图画的t和p其实不是对应的 clear;clc; gh=@(t,x)[x(2);(-4.5/200000)*(x(2)^2)] y0=[0,240];%初始条件 [t,xx]=ode45(gh,[0,200],y0);%参数区间 p=xx(:,1);%水平速度 plot(t,p) yli...

用MATLAB ode45命令求解,对于分段函数可以用if语句,初值可以任意设定,只是不能全部为0,例如可以设为(0,1,-1,0)。

使用符号运算工具。 首先定义变量符号 syms s x t; 使用字符串定义公式 eq1 = 'diff(x, t) = 0.11*s*x/(0.06*x+s)'; eq2 = 'diff(s, t) = (-1)*0.029*x-0.234*s/(0.006*x+s)-s/(30*s^2+300*s+0.03)'; 接下来是求解 比如现在要求解s,下面的代码...

由题目可知,自变量是t,因变量是x,y,z,化为matlab时是x(1),x(2),x(3),不会有x(4)--x(6),说明函数式有问题。

请贴出具体的非线性微分方程组及要求,因为不同的方程其所用的方法有所不同。

%我在用matlab解n元一次非线性方程组中遇到了困难,不知道如何编写代码,想麻烦您帮我看看%?entry=qb_ihome_tag%n不是某个定值,n可以由自己任取c=t(x)./x;sum(x)=Q;function f=t(x)f=x.^5+x.^6-1./x;end你怎么编???

网站首页 | 网站地图
All rights reserved Powered by www.clzm.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com