jcst.net
当前位置:首页 >> 在mAtlAB中用隐式梯形积分法对二阶微分方程(例如... >>

在mAtlAB中用隐式梯形积分法对二阶微分方程(例如...

>> 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))/(...

这个我也不是很懂,仍然是在学习中。 如果是用前时间状态推导后一个时间步长的状态,可以找一下显示欧拉法和隐式欧拉法,我目前学的就是这些,只学懂一些皮毛,仍然在头疼中。

1.新建一个m文件,编写隐式Euler法的程序: function [x,y]=Implicit_Euler(odefun,xspan,y0,h,varargin) % 隐式Euler公式求解常微分方程 % 输入参数: % ---odefun:微分方程的函数描述 % ---xspan:求解区间[x0,xn] % ---y0:初始条件 % ---h...

dsolve('D2y=x','x') ans = x^3/6 + C2*x + C3 D2y的意思就是y的二阶微分项 不明白你的问题什么意思,要输入的话直接定义符号变量输入... syms D2x x D2x=x; 如果是矩阵,那就用矩阵表示...

一、理论解,用dsolve函数(在command window 中输入doc dsolve可以查看帮助) 示例: 代码: [x,y]=dsolve('D2x+2*Dx=x+2*y-exp(-t)','Dy=4*x+3*y+4*exp(-t)') 二、数值解,用ode45,或ode23, ode15s其他函数。(在command window 中输入doc od...

1、对于解析值,你可以用dsolve()函数求解。如求微分方程x*y''+x﹡(y')^2-y'=0的解析解 >> syms y(x) >>Dy = diff(y);D2y = diff(y, 2); >>dsolve(x*D2y+x*(Dy)^2-Dy==0,'x') 运行结果 2、对于数值解,你可以查看二阶常微分方程这个例子,http...

参考代码: m=10600; k=31986497; c=8.9E4;dy = @(t,y)[y(2); (k*y(1)-c*y(2)/abs(y(2))*y(2)^2)/m];ode45(dy,[0 5],[0 3.66]) 其中,蓝线为x,绿线为x'。

这个我也不是很懂,仍然是在学习中。 如果是用前时间状态推导后一个时间步长的状态,可以找一下显示欧拉法和隐式欧拉法,我目前学的就是这些,只学懂一些皮毛,仍然在头疼中。

例 y’= - y+x+1,y(0) = 1 首先建立M-文件 (weif.m) function f = weif(x,y) f=-y+x+1; 求解:[x,y]=ode45('weif',[0,1],1) 再如: 建立文件: function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -...

用matlab的dsolve函数计算该二阶微分方程的解析值是有一定的难度。只能用数值分析的方法来求解其数值解。当然求其数值解,必须还得已知a、θ的具体数值。该问题可以用matlab的ode45函数求得。求解方法如下: 1、建立微分方程的自定义函数,odefun...

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