matlab解二阶微分方程

陈信帝3699matlab 求解二阶微分方程并画出y与t函数曲线 -
官池学19434201408 ______ 用matlab 求解二阶微分方程并画出y与t函数曲线,可以用dsolve函数和plot函数来实现.给出的二阶微分方程存在着几个疑点:1、表达式中的z是否是笔误,还是y.如是z,其表达式是什么?2、求解二阶微分方程的数值解,必须有两个初始条件...

陈信帝3699matlab 二阶微分方程 求解 -
官池学19434201408 ______ >> clear>> syms a b c d e;>> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e') y =(2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))...

陈信帝3699如何用matlab求解二阶非线性偏微分方程? -
官池学19434201408 ______ 首先可以肯定的告诉你,matlab本身有个pde工具箱2113,可以求解一般的典型二阶线性偏5261微分方程.但是,对于非线性的偏微分方程,我估计很4102难直接用matlab的工具箱或者函数命令得到结果.因为这涉及到一1653个数学处理的问题,你首先要处理这个非线版性问题,或者用数值的方法处理,然后再用权MATLAB编程计算.

陈信帝3699matlab 中二阶常微分方程的数值解法
官池学19434201408 ______ odefun=@(t,x)[x(2);3*x(2)-2*x(1)+1];[t,y]=ode45(odefun,[0:0.01:2],[1 0]);plot(t,y)[t y] 结果 y(0.5000)=0.7896 y= dsolve('D2y-3*Dy+2*y=1','Dy(0)=0','y(0)=1');>> yy =exp(t) - exp(2*t)/2 + 1/2>> feval(@(t)exp(t) - exp(2*t)/2 + 1/2,0.5)ans = 0.7896

陈信帝3699matlab二阶微分方程已知函Y=(1.082X)^1/2 用四次牛顿前插公式计算f( - 207) 求matlab源程序 最好能列出差分表 -
官池学19434201408 ______[答案] >> clear >> syms a b c d e; >> y=dsolve('a*D2y+b*Dy+c*y=0','y(0)=d','Dy(0)=e') y = (2*a*e + b*d + d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b - (b^2 - 4*a*c)^(1/2)))/(2*a))*(b^2 - 4*a*c)^(1/2)) - (2*a*e + b*d - d*(b^2 - 4*a*c)^(1/2))/(2*exp((t*(b + (b^2 - 4*a*c)^(1/2)))/(2*...

陈信帝3699请问这个二阶微分方程如何用matlab求解 -
官池学19434201408 ______ 使用DSolve[]

陈信帝3699matlab 用ode45解二阶微分方程:t''=sin(t) 麻烦给程序,谢谢 -
官池学19434201408 ______ 我来解答你的关键问题:将之前主程序计算出来的系数用global cof定义,再在ode45的子函数里同样定义global cof,这样就能把值计算出来并且传递进去了.

陈信帝3699用Matlab编程求解 二阶微分方程:4*d^2y(t)/dt^2+y(t)=dx(t)/d(t) - 0.5x(t) -
官池学19434201408 ______ MATLAB提供了dsolve命令可以用于对符号常微分方程进行求解.语法: dsolve('eq','con','v') %求解微分方程 dsolve(...

陈信帝3699matlab中求二阶微分方程 的解析解和数值解,并在同一图形中画出解析解和数值解的图形进行比较 -
官池学19434201408 ______ y=dsolve('(1+x^2)*D2y=2*x*Dy','y(0)=1','Dy(0)=3',x) y = x*(x^2 + 3) + 1 数值解法: 首先编写函数myfun.m function dy=myfun(x,y) dy=zeros(2,1) dy(1)=y(2); dy(2)=2*x*y(2)/(1+x^2) 调用ode45: [x45,y45]=ode45('myfun',[0 10],[1 3]) plot(x45,y45,'+') hold on x1=0:0.5:10; y1=subs(y,x1); plot(x1,y1,'o')