matlab怎么建立微分方程,二元二阶非线性微分方程组用matlab求解
你用一下ode45试试
d2θ/ds2=cosθ-sinθ两边同时乘dθdθ/dsd(dθ/ds)=cosθdθ-sinθdθ(dθ/ds)^2=2sinθ+2cosθ+cds=dθ/√(2sinθ+2cosθ+c)这个积分不能用初等函数表示:在https://baike.baidu/item/超越积分/5292075?fr=aladdin第8个,所以方程的解也不能用初等函数表示.
首先可以肯定的告诉你,matlab本身有个pde工具箱2113,可以求解一般的典型二阶线性偏5261微分方程.但是,对于非线性的偏微分方程,我估计很4102难直接用matlab的工具箱或者函数命令得到结果.因为这涉及到一1653个数学处理的问题,你首先要处理这个非线版性问题,或者用数值的方法处理,然后再用权MATLAB编程计算.
matlab怎么建立微分方程
能画一个当然也能画多个了. 画第二个之前用plot on命令,这样后面再绘图的时候就. 另外很多微分方程可以直接解出来的,用dsolve函数. 比如你给的这个例子就可以,.
用matlab求解微分方程组并画图,其实现过程:1、首先应根据已知微分方程组,编写其函数 dy(1)=-0.10982*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3)); dy(2)=1441.5074*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3));2、用ode45函数求解其微分方程组的数值解 [t,y]=ode45(@ode_fun,tspan,y0)3、然后用plot绘图函数,绘出t-Cg和t-θ关系曲线图4、编程后运行结果
先建立一个函数文件function,将所要求解的微分方程写入此函数文件,另建一个m文件或在命令窗口调用ode45函数,得到微分方程的数值解,后用plot函数做图就ok了
用matlab解微分方程组
使用ode45解微分方程组,可以用下列格式来求解.[t,x]=ode45(odefun,tspan,x0) 式中:odefun——给复定微分方程组的自定义函数;tspan——自变量t的范围,制即【t0,tf】;x0——x的初值向量,即x1(0),x1'(0),x2(0),x2'(0),x3(0),x3'(0) 由于没有给出自zd变量t的范围和x的初值,所以不能得到其x1,x2,x3的值.
那y1,y2呢? 那个在y(3)是y的三次导的意思. 上面所有的y1,y2我认为也是求导的意思. 你可以做数值解, 下面做数值解: 在Matlab下输入:editzhidao_rk4_4.m,然后将下.
function hh X0=[17,195,1500,4,62,200]; TSPAN=[0 4]; [t,X] = ode45(@myfun,TSPAN,X0) plot(t,X) legend('x','y','z','u','v','w','Location','NorthWest') xlabel('t'),ylabel('X') function .
matlab求解常微分方程
dsolve('D3y+D2y/x-Dy/(x^2)=x','y(500)=0,Dy(500)=0','x'),用这个能够得到y的一个含参数的解,然后根据x=0时,y为有限值来确定参数的值,之后就能解出相应的y'' 和y' ,后面那个貌似用MATLAB不太好实现,最好是自己来分析.
>> y=dsolve('D2y+Dy=cos(x)','y(0)=2','Dy(0)=1','x')y =sin(x)/2 - cos(x)/2 - exp(-x)/2 + 3
有符号法和数值法.一般使用dsolve和ode命令,比如: 符号法 dsolve('df = f + sin(t)') ans = c4*exp(t) - sin(t)/2 - cos(t)/2 数值法 sol = ode45(@vdp1,[0 20],[2 0]);x = linspace(0,20,100);y = deval(sol,x,1); plot(x,y);
matlab解多元微分方程组
如果是常微分方程,可以用dsolve函数.该函数可以解单变量常微分方程或者多变量常微分方程组,所以5个变量也不在话下.调用格式如下:[y1,,yN]=dsolve(eqns).
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.
初始条件y(t)=0是不是有问题啊,如果是这样,直接解出y=0,是不是y(0)=0啊?clc;cleary=dsolve('d2y=g-k*s*v/m','dy(0)=0,y(0)=0')y1=diff(y,'t')v=solve(y1-'v','v')y=subs(y)结果:y =1/2*(g-k*s*v/m)*t^2y1 =(g-k*s*v/m)*tv =g*t/(k*s*t+m)*my =1/2*(g-k*s*g*t/(k*s*t+m))*t^2
matlab解高阶微分方程
解常微分方程用dsolve实现 在命令窗口输入 [x,y,z]=dsolve('dx=y+x*z,dy=z,dz=x') 方程可以替换为你的啊
1、你的ode函数写的有问题,应该是function f=tuoyuan(t,x)f=[x(3); x(4); (-5*x(1)/sqrt((x. 你确定微分方程没错吗?怎样见得理论上是一个椭圆呢?我知道下面的微分方程 x＂=.
这个方程的解为y=4/3cos(2t)-2sint+1/6(t)3-1.9(t)2+at+b注:(t)2是的t平方;同样(t)3是的立方;a b 是任意常数