非线性二元二阶微分方程组
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/23 19:31:13
% by dynamic of Matlab技术论坛% 2010-02-26 18:57:08%首先回答,Matlab可以能求解二元高阶微分
代码:clearallclcf=@(t,y)([y(2);2*y(2)+.15*y(4)^2-10; y(4);.01*y(4)^2]);[t,Y]=ode45(f,
(1)a/c≠b/d方程组有唯一解a/c==b/d==m/n无穷多组解另:a/c==b/d≠m/n方程组无解(2)非零解即x,y均不为0的解3x+my=0(1)4x-11y=0(2)(1)×4-(2)
用MATLABode45命令求解,对于分段函数可以用if语句,初值可以任意设定,只是不能全部为0,例如可以设为(0,1,-1,0).
这个问题很难,
function [t,x]=solv(tspan)% argument% tspan stands for time spanx
clcYd=inline('[(YY(2)-1)*2*pi*50;(1-(1.47*1/2.82)*sin(YY(1)))/11.28]','t','YY');[t,YY]=ode45(Yd,0.3,
matlab里面常使用龙格库塔方法求解常微分方程组,命令是ode45,还有其他一些函数,但是最常用的是ode45,lz可以help一下,很简单的,另外给你一个文档,讲的还是比较详细,希望可以帮到你ht
假设x(1)=x1,x(2)=x1',x(3)=x2;把高阶导数换算成低阶倒数;得到3个方程联立的方程组1:x(1)'=x(2)2:x(2)'=-2*x(2)-2*x(3)3:x(3)'=-x(2)-
F表达式多了中括号,去掉后得到正解.Equationsolved.fsolvecompletedbecausethevectoroffunctionvaluesisnearzeroasmeasured
人家问的就是解析解,不知道上边几个在说什么你这个方程高版本的matlab可以解,只是算出来之后带积分项目:symsabQSTzZmcAPMeESS=dsolve('Dx=a*x','Dy=(b*Q*x
可以用ode45求数值解可以参考我给别人回答的再问:我是想对这个方程组做图。先解这个方程组,才能做出图来是不是?请问针对这个问题,m-file文件的程序该怎么编?再答:%%%%%%%%%%%%%fun
警告的翻译:无法满足积分条件,因此在时间t=1.768731e-001s处无法积分帮你运行了一下,原因是在1.768731e-001s处,有个值趋于负的无穷大,因此无法积分可能的原因:1.方程出错;2
这个只能求数值解吧.求数值解,先要将微分方程转化为差分方程,即离散化.然后自己给的初值,对t进行不断的迭代.就可以求出数值解.可以看看偏微分方程数值解的书,学习离散化过程.再问:谢谢,推荐几本书能否?
使用符号运算工具.首先定义变量符号symssxt;使用字符串定义公式eq1='diff(x,t)=0.11*s*x/(0.06*x+s)';eq2='diff(s,t)=(-1)*0.029*x-0.
当然可以解.通过定义两个变量dx和dy,把二阶方程组改成四个一阶方程就可以了.另外,你的问题不完整,还需要多两个初始条件x'(0)、y'(0).
没解析解你肯定得用数值解.以下是步骤.第一步,在Matlab里写m-file,命名微分方程对应的函数:functionf=myfun(t,x)f=[-x(1)+x(2)^2;-2*x(2)+x(1)^
参考http://zhidao.baidu.com/question/479644084.html?oldq=1回答的一样的问题,还有下面的程序functiondy=odefun(t,y)%y(1)x
用solve可以解决,helpsolve看看用法,很简单的.
程序给你稍稍改动一下,本身没有问题,调试通过,最为关键的是初值的选取.强烈建议买一本程耿东院士的《工程结构优化设计基础》一书,是国内最具权威的关于优化方法入门的书籍,大连理工大学出版社.我根据自己的理