matlab解常微分方程范例
管桑秋3362用MATLAB解二阶常微分方程里,两个边界条件,其中一个边界在无穷大用MATLAB解二阶常微分方程里,两个边界条件,一个在一个已知x1处为0,另一个... -
欧强策19478053644 ______[答案] 没法设置无穷大处的边界条件的,貌似 不过可以设置一个很大的数作为边界 比如y(10000)=1 然后就用matlab自带的ode45等解就可以了,可以参考一下以下的帖子
管桑秋3362matlab求解常微分y=dsolve('D3y - 3*D2y - Dy*y=0','y(0)=1 Dy(0)= - 1 D2y(0)=0','t')错误如下? Error using ==> dsolveat offset 84, `=` unexpectedError in ==> fun1 at ... -
欧强策19478053644 ______[答案] 第一,你有语法错误 y=dsolve('D3y-3*D2y-Dy*y=0','y(0)=1','Dy(0)=-1','D2y(0)=0','t') 第二,没有显式解 只能把初始条件去掉可得隐式解
管桑秋3362高手进,用matlab求解复杂常微分方程的数值解 -
欧强策19478053644 ______ 是怎样复杂的常微分方程?一般常微分方程的数值解可以用ode()函数.求解方法:y0=[a,0]; %初值条件 [t,y]=ode45(@(t,y)odefun(t,y,a),[0 20],y0) %数值解 plot(t,y(:,1),'*',t,y(:,2),'.-'),grid on,hold on %绘图 自定义函数的调用 function dy=odefun(t,y,a); dy=【y(2);(1/a)*sqrt(1+y(1)^2)】;
管桑秋3362常微分方程求解,急用!y''+ay+by^3=0,a,b为常数.这样一个常微分方程怎么求解?毕业设计急用!matlab能求出数值解吗? -
欧强策19478053644 ______[答案] 你学过数理方程没有?其实将其做个拉普拉斯变换就可以了,令y=t,且K(P)=∫te^-pt dt=1/p^2 则将两边都进行变换,得到p^2*1/p^2-1+a/p^2+b*6/p^4=0 得到 a*p^2+6b=0 则p=√-6b/a 再进行逆变换得到y的值 y=1/2∏i∫√-6b/a*e^pt dp=√-6b/a*e^t
管桑秋3362求助过路的matlab大神,老师留作业:用两种欧拉方法解常微分方程方程是 20y＂+y'+0.5y=5sin(3x) 其中 h=0.1,y'(0)=1,y＂(0)= - 1 -
欧强策19478053644 ______[答案] 欧拉方法的matlab 先定义函数euler function [x,y]=euler(fun,x0,xfinal,y0,n); if nargin
管桑秋3362matlab 解常微分方程
欧强策19478053644 ______ dsolve('M*Dv=-N*(f+b*sin(C*atan(B*(1-w*R/v)-E*(B*(1-w*R/v)-atan(B*(1-w*R/v))))))','t') dsolve('I*Dw=N*(f+b*sin(C*atan(B*(1-w*R/v)-E*(B*(1-w*R/v)-atan(B*(1-w*R/v))))))-T','t')
管桑秋3362在MatLab里面用隐式欧拉法(backward euler)解决常微分方程.初学matlab 好多都不会,知道的帮下忙 -
欧强策19478053644 ______ 1.新建一个m文件,编写隐式Euler法的程序: function [x,y]=Implicit_Euler(odefun,xspan,y0,h,varargin) % 隐式Euler公式求解常微分方程 % 输入参数: % ---odefun:微分方程的函数描述 % ---xspan:求解区间[x0,xn] % ---y0:初始条件 % ---h:迭...
管桑秋3362matlab能解什么样的微分方程 -
欧强策19478053644 ______ MATLAB 求解常微分方程的函数是 dsolve().应用此函数可以求得常微分方程(组)的通解,以及给定边界条件(或初始条件)后的特解.常微分方程(组)可以是高阶的.不能解偏微分方程但可以解非常系数的微分方程(组). 注:MATLAB 之所以有如此强大的符号运算功能,完全要归功于 MAPLE.Maths 公司拥有MAPLE 的内核后,MATLAB的功能得到了巨大的增强.
管桑秋3362【SOS】在matlab中四阶Runge - Kutta法求解常微分方程
欧强策19478053644 ______ 没试过matlab,算这玩意太慢了,有fortran版的要不,有兴趣的话可以参考一下. 代码: SUBROUTINE runge_kutta() !关于Runge-Kutta方法,该方法是用来解形如y'=f(t,y)的常微分方程的 !经典的4阶R-K方法的公式如下: ! Yn+1 = Yn + ...