clzm.net
当前位置:首页 >> mAtlAB解微分方程 >>

mAtlAB解微分方程

楼主是要解析解吧 dsolve('Dy=A/V-B/V*sqrt(2*g*y)','y(0)=h0','t') ans = (2^(1/2)*A + 2^(1/2)*A*lambertw(0, (2^(1/2)*B*g^(1/2)*exp(-(A + B^2*g*(t/V - (A*(log(-(exp(-(A + 2^(1/2)*B*(g*h0)^(1/2))/A)*(2*B*(g*h0)^(1/2) + 2^(1/2)*A))/(...

dsolve()求解微分方程后,使用ezplot()作图,参考代码: >> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间

问题1、function函数中未用到x,为何要将它设置为“输入变量” 因为微分方程的解 y=f(x),在求解时是用到变量x的。 问题2、输出结果中Y有三列,分别表示什么意思? 输出结果中Y有三列数值,分别表示y,y',y"的值,即y——y(1),y'——y(2),y"——y(3)

好像没什么特别好的办法: >> syms t q y; q=4*t^3;>> syms Dy>> eq=Dy-(q-2*t);>> y=dsolve(char(subs(eq)),'y(0)=5','t')y = t^4-t^2+5 也就是要先移项把微分方程写成等号右端为零的形式(需额外定义Dy),然后subs代入,再求解。

%欧拉法解一阶常微分方程 %例子dy/h=-y+x+1 %f=inline('-y+x+1','x','y'); %微分方程的右边项 f = inline('x-2*y','x','y'); y0 = 2; %初始条件 h = 0.025; %步长 xleft = 0; %区域的左边界 xright = 1; %区域的右边界 x = xleft:h:xright; n = ...

分段函数的一个简单方法就是分了几段就写几个函数: 定义一个主调函数、两个子函数 最后一个主调程序调用 function dx=fun(t,x) dx=zeros(2,1); f=; dc(t)=diff('f','t');%f就是你的c(t) dx(1)=x(2)*dc(t); dx(2)=x(1)+6*t^2;%0

在Matlab下输入:edit,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function dxdt=zhidao_rk4_5(t,x) global a b c d r x1=x(1); x2=x(2); x3=x(3); x4=x(4); dxdt=[ a*(x2-x1)+x4; d*x1-x...

说明 求解常微分方程一般使用ode系列函数(其中最常用的是ode45),基本做法就是把微分方程用函数表示出来,然后调用ode函数求解即可。对于本题而言,使用匿名函数描述微分方程更为方便一些。 参考代码 k0=8; k1=1; K=1; k2=5; n=2; % 常数定义d...

改微分方程的符号特解为: ys=dsolve('D2y+6*Dy+9*y-10*sin(x)','y(0)=0','D(y)=0','x') 不知道你这道题是不是还要求微分方程的数值解,自变量的范围也没给定,需要画图也没说明。

在Matlab下输入:edit zhidao_feiying.m,然后将下面两行百分号之间的内容,复制进去,保存 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function y=zhidao_feiying(t,x) a=1; b=2; c=3; d=4; %%比方说 %f=a*x+y-z; %g=b*sin(x*y)-c*cos...

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