ZKX's LAB

薛定谔方程的解法? 欧拉法与龙格库塔法的区别

2021-03-11知识4

欧拉法,改进欧拉法的一阶微分方程组迭代格式.. 随便找本数值分析或者计算方法的书上都有

急求用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')

薛定谔方程的解法? 欧拉法与龙格库塔法的区别

薛定谔方程的解法? 定态方程所谓势场,就是粒子在其中会有势能的场,比如电场就是一个带电粒子的势场;所谓定态,就是假设波函数不随时间变化。其中,E是粒子本身的能量;v(x,y,z)是描述势。

数值分析的内容简介 《数值分析(高校教材)》系统地阐述了数值分析的基本知识,介绍了各种数值计算方法,全书共分十三章。第一章介绍数值计算的基本概念和误差分析的知识;。

什么是欧拉方法(Euler's method)? 孤光一点萤:数值常微分方程-欧拉法与龙格-库塔法 ? zhuanlan.zhihu.com 比较了一下欧拉法与二阶龙格-库塔法,希望能有所帮助。如有错误疏漏,还请指正。这里直接复制粘贴。

Matlab/simulink中,什么叫Oder45和Ode23bt算法? ode45是基于四点法和五点法的解微分方程数值解的方法,ode23等也一样,都是基于已知点“预测”下一个点的函数值的方法,不同的算法“预测”的方法不一样。比较著名的“预测。

两周速成参加数学建模美赛可行吗? ?www.zhihu.com?www.zhihu.com?www.zhihu.com ?www.zhihu.com 欢迎大家关注 别说2周,3天都行,先挖个坑,过几天来补,这几天家里忙。开始填坑!先说一下3天是指任何。

Matlab二阶非线性常微分方程的求解 d2θ/ds2=cosθ-sinθ两边同时乘dθdθ/dsd(dθ/ds)=cosθdθ-sinθdθ(dθ/ds)^2=2sinθ+2cosθ+Cds=dθ/√(2sinθ+2cosθ+C)这个积分不能用初等函数表示:在https://baike.baidu.com/item/超越积分/5292075?fr=aladdin第8个,所以方程的解也不能用初等函数表示。

#欧拉法与龙格库塔法的区别

随机阅读

qrcode
访问手机版