常微分方程matlab实例

鲜福芸4618高手进,用matlab求解复杂常微分方程的数值解 -
王肥晏15028134758 ______ 是怎样复杂的常微分方程?一般常微分方程的数值解可以用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)】;

鲜福芸4618matlab中符号法求常微分方程(x^2+1)y' - (y^2+1)=0 -
王肥晏15028134758 ______[答案] 把自变量x变为t >> y=dsolve('(t^2+1)*Dy-(y^2+1)=0') y =tan(atan(t)+C1) 即 y =tan(atan(x)+C1) 给10分吧,好吗?看看我的回答 总分:16464 ;回答总数:682;题均:16464/682=24.1408

鲜福芸4618matlab中求常微分方程的数值解 -
王肥晏15028134758 ______ t0 = 0; tf = 5; % y'' = (2-3y-2ty')/(1+t^2) dy = @(t,y) [y(2);(2-3*y(1)-2*t*y(2))/(1+t^2)]; [T,Y] = ode45(dy,[t0 tf],[0 1]); plot(T,Y(:,1))

鲜福芸4618matlab如何画常微分方程的解曲线与轨道 -
王肥晏15028134758 ______ 如果你会matlab下的simulink仿真模块的话,直接搭建两个模块就解决了.这种方法可以解决几乎所有微分方程的图像问题. 否则,你只能通过下述办法了,用dsolve求出目标方程的解析解表达式,然后用fplot函数画出上述解析解表达式的图像.或者求出解析解之后,给自变量一些值,求出响应的因变量的一组值,然后用plot函数画图.这两种方法本质都一样,都严重依赖于解析解.不过你们学生做题的话所有的题目都是有解析解的.

鲜福芸4618如何使用matlab用梯形法解常微分方程dy/dx=y/x - y^2/x,y(1)=0.5,求出其数值解 -
王肥晏15028134758 ______[答案] y=dsolve('Dy=y/x-y^2/x','y(1)=0.5','x') y = 1/(1/x + 1) 梯形公式: y(1)=0.5;x(1)=1;h=0.1; for n=1:99 x(n+1)=x(n)+h; y(n+1)=y(n)+h*( y(n)/x(n)-y(n)^2/x(n)); y(n+1)=y(n)+h/2*( y(n)/x(n)-y(n)^2/x(n) +y(n+1)/x(n+1)-y(n+1)^2/x(n+1) ); end plot(x,y,'y')

鲜福芸4618matlab 一阶常微分方程求参数 -
王肥晏15028134758 ______ 正着做:给一个r和xm,通过ode45求得t,x,然后再与你的数据对比.最终选择一个合适的r和xm.就要一直变参数,我觉得比较难.至少来说r和xm的选择范围太大了.而且ode45得到的t,x与你的实验t,x肯定不是同一个t下面的数据,也不好比较...

鲜福芸4618matlab解高阶常微分方程 -
王肥晏15028134758 ______ dY(2) = 2*m/h^2*[U - En]*Y(1); 这个语句的m h等变量 看看都是几维的 只能是一个数值才行

鲜福芸4618matlab中符号法求常微分方程
王肥晏15028134758 ______ 把自变量x变为t >> y=dsolve('(t^2+1)*Dy-(y^2+1)=0') y =tan(atan(t)+C1) 即 y =tan(atan(x)+C1) 给10分吧,好吗?看看我的回答 http://passport.baidu.com/?center&tpl=ik&aid=7&default_tab=4#4,0 总分:16464 ;回答总数: 682;题均: 16464/682=24.1408

鲜福芸4618一类二阶常微分方程的几种解法 -
王肥晏15028134758 ______[答案] 1、引言常微分方程有着深刻而生动的实际背景,它从生产实践与科学技术中产生,而又称为现代科学技术中分析问题与解决问题的一个强有力的工具.人们对二阶及以上微分方程(包括线性、常系数、隐性)的研究,产生了许多理论成果.如胡爱莲[1]...

鲜福芸4618matlab用dsolve解常微分方程并画图, -
王肥晏15028134758 ______ >> y=dsolve('Dy=1-y^2','y(0)=0') y = tanh(t) >> t=-2*pi:0.1*pi:2*pi; >> y=tanh(t); >> plot(t,y)