关于matlab 四阶龙格库塔算法的 上面用的好像不是FFT,应该用的是DCT变换吧。换上DCT应该相似了。
龙格库塔算法应用方向? 就是一种省去方程的求微分的方法,一般是四阶形式
二阶龙格库塔方法 最低0.27元开通文库会员,查看完整内容>;原发布者:龙吟太虚2012-2013(1)专业课程实践论文二阶Runge-Kutta方法董文峰,0818180123,R数学08-1班一、算法理论由改进的Euler方法得到:凡满足条件式有一簇形如上式的计算格式,这些格式统称为二阶龙格—库塔格式。因此改进的欧拉格式是众多的二阶龙格—库塔法中的一种特殊格式。若取,就是另一种形式的二阶龙格-库塔公式。(1)此计算公式称为变形的二阶龙格—库塔法。二级龙格-库塔方法是显式单步式,每前进一步需要计算两个函数值。由上面的讨论可知,适当选择四个参数y0,a,b,n,可使每步计算两次函数值的二阶龙格-库塔方法达到二阶精度。下面以式子(1)为依据利用VC+6.0编译程序进行问题的求解。二、算法框图三、算法程序#include#include/*n表示几等分,n+1表示他输出的个数*/intRungeKutta(doubley0,doublea,doubleb,intn,double*x,double*y,double(*function)(double,double)){doubleh=(b-a)/n,k1,k2;inti;x[0]=a;y[0]=y0;for(i=0;i;i+){x[i+1]=x[i]+h;k1=function(x[i],y[i]);k2=function(x[i]+h/2,y[i]+h*k1/2);y[i+1]=y[i]+h*k2;}return1;}doublefunction(doublex,doubley){returny-2*x/y;}intmain(){inti;doublex[6],y。
二阶的改进欧拉近似算法和二阶龙格-库塔法 精度一样的情况下 哪个更精确啊 后面那个
用龙格库塔算法解二阶常微分方程,利用c++编程 另z=y'原来的方程就可以化简成以z'和y,x的方程和y'=z 带入两个初始条件,就可以进行迭代了
分别用改进的欧拉法和四阶龙格-库塔公式求解微分方程初值问题 分别用改进的欧拉法和四阶龙格-库塔公式求解微分方程初值问题(1)Y'=Y-2X/Y,Y(0)=1,X=[0,1],H=0.1(2)Y'=X2+。
龙格-库塔-芬尔格算法是什么算法?解决什么问题的?? 是求解抄百ODE的算法度。http://wenku.baidu.com/view/bb3bc1ea551810a6f5248633.html
龙格库塔求数值解有什么限制条件? 你好,龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。由于此算法精度高,采取措施对误差进行抑制,所以其实现原理也较复杂。同前几种算法一样,该算法也是构建在数学支持的基础之上的。对于一阶精度的欧拉公式有:yi+1=yi+h*K1 K1=f(xi,yi)当用点xi处的斜率近似值K1与右端点xi+1处的斜率K2的算术平均值作为平均斜率K*的近似值,那么就会得到二阶精度的改进欧拉公式:yi+1=yi+h*(K1+K2)/2 K1=f(xi,yi)K2=f(xi+h,yi+h*K1)下面的具体程序实现同改进的欧拉算法类似,只需作些必要的改动,下面将该算法的关键部分代码清单列出:…for(floatx=0;x;x+0.1){r=x+expf(-x);K1=x-y[i]+1;file:/求K1K2=(x+(float)(0.1/2))-(y[i]+K1*(float)(0.1/2))+1;file:/求K2K3=(x+(float)(0.1/2))-(y[i]+K2*(float)(0.1/2))+1;file:/求K3K4=(x+0.1)-(y[i]+K3*0.1)+1;file:/求K4y[i+1]=y[i]+(float)(0.1*(K1+2*K2+2*K3+K4)/6);file:/求yi+1r=fabs(r-y[i]);file:/计算误差str.Format(\"y[%d]=fr=f\\r\\n\",i,y[i],r);i+;msg+str;}AfxMessageBox(msg);file:/报告计算结果及误差情况