ZKX's LAB

欧拉 龙格-库塔与精确解绘图 matlab 里面 powergui选项中simulation type中三个选项什么区别 phashor continuous discrete

2020-09-27知识10

从欧拉方法、改进欧拉方法、2阶龙格-库塔方法、4阶龙格-库塔方法中选择一种方法,每一步从精确解出发计算出下一 欧拉方法 ;nbsp;yn+1=yn+h·f(xn,yn),xn=x0+n·h. ;nbsp;改进欧拉方法 ;nbsp;nbsp;nbsp;nbsp;nbsp;2阶龙格-库塔方法 ;nbsp;yn+1=yn+hk2,k1=f(xn,yn)。

matlab 里面 powergui选项中simulation type中三个选项什么区别 phashor continuous discrete 尤其是continuous,计算机就是计算机,根本不可能是连续的,怎么会有。

自然常数 e 的两个表达式为什么是等价的,有什么联系? 谢邀. 要讲明白这个事,简单也不简单.有必要先回顾一下历史. 所谓自然常数 通常称为\"自然对数的底\",事实上,因为在历史上,对数函数先于指数函数。?www.zhihu.com 注: 推导

急求用MATLAB用龙格库塔和外推法解一阶微分方程 f=inline('-y+x+1','x','y');微分方程的右边项dx=0.05;x方向步长xleft=0;区域的左边界xright=10;区域的右边界xx=xleft:dx:xright;一系列离散的点n=length(xx);点的个数y0=1;(1)欧拉法Euler=y0;for i=2:nEuler(i)=Euler(i-1)+dx*f(xx(i-1),Euler(i-1));end(2)龙格库塔法RK=y0;for i=2:nk1=f(xx(i-1),RK(i-1));k2=f(xx(i-1)+dx/2,RK(i-1)+k1*dx/2);k3=f(xx(i-1)+dx/2,RK(i-1)+k2*dx/2);k4=f(xx(i-1)+dx,RK(i-1)+k3*dx);RK(i)=RK(i-1)+dx*(k1+2*k2+2*k3+k4)/6;endEuler和Rk法结果比较plot(xx,Euler,xx,RK)hold on精确解用作图syms xrightsolve=dsolve('Dy=-y+x+1','y(0)=1','x');求出解析解rightdata=subs(rightsolve,xx);将xx代入解析解,得到解析解对应的数值plot(xx,rightdata,'r*')legend('Euler','Runge-Kutta','analytic')

#微分方程#欧拉#龙格库塔法#近似算法#欧拉常数

随机阅读

qrcode
访问手机版