典型数值算法C语言程序设计.docVIP

典型数值算法C++语言程序设计 典型数值算法的C++语言程序设计 1 经典四阶龙格库塔法解一阶微分方程组 1.1 经典四阶龙格库塔法解一阶微分方程组算法分析 经典4阶龙格?库塔公式为 ， （1-1） 其中, ， ， ， ， ， ， （1-2） ，1.2 经典四阶龙格库塔法解一阶微分方程组算法流程图 图1-1 经典四阶龙格库塔法解一阶微分方程组流程图 1.3 经典四阶龙格库塔法解一阶微分方程程序调试运行图　　　例如，应用本程序求解微分方程组 ，其中初值为，求解区间为[0,0.2]。 程序调试运行图如下： 图1-2 经典四阶龙格库塔法解一阶微分方程程序调试运行图 1.4 经典四阶龙格库塔法解一阶微分方程源程序代码 //经典四阶龙格库塔法解一阶微分方程 #include #include //控制输出格式 using namespace std; void Runge_Kutta double *fdouble t,double x, double y,double *gdouble t,double x, double y ,double initial[3], double resu[3],double h double f1,f2,f3,f4,g1,g2,g3,g4,t0,x0,y0,x1,y1;t0initial[0]; //定义初值x0initial[1]; 　　y0initial[2];f1ft0,x0,y0; //四阶龙格库塔法算法　 g1gt0,x0,y0;f2ft0+h/2, x0+h*f1/2,y0+h*g1/2;　 g2gt0+h/2, x0+h*f1/2,y0+h*g1/2; 　 f3ft0+h/2, x0+h*f2/2,y0+h*g2/2;　 g3gt0+h/2, x0+h*f2/2,y0+h*g2/2;f4ft0+h, x0+h*f3,y0+h*g3;　 g4gt0+h, x0+h*f3,y0+h*g3;x1x0+h*f1+2*f2+2*f3+f4/6; 　　y1y0+h*g1+2*g2+2*g3+g4/6;resu[0]t0+h;resu[1]x1; //迭代下一点 　　resu[2]y1; int main double fdouble t,double x, double y; double gdouble t,double x, double y; double INITIAL[3],RESU[3]; double a,b,H; double t,step; int i; cout cout 经典四阶龙格库塔法 endl; cout coutPlease input the initial t0,x0,y0:; cinINITIAL[0]INITIAL[1]INITIAL[2]; coutPlease input the interval:; cinab; coutPlease input step:; cinstep; coutsetiosflagsios::rightsetiosflagsios::fixedsetprecision10;Hb-a/step; cout INITIAL[0]setw18INITIAL[1]setw18INITIAL[2]endl; fori0;istep;i++ Runge_Kutta f,g ,INITIAL, RESU,H; cout RESU[0]setw18RESU[1]setw18RESU[2]endl; INITIAL[0]RESU[0];INITIAL[1]RESU[1];INITIAL[2]RESU[2]; return0; double fdouble t,double x, double y double dx; dxx+2*y; returndx; double gdouble t,double x, double y double dy; dy3*x+2*y; returndy; 2 高斯列主元法解线性方程组 2.1 高斯列主元法解线性方程组算法分析 　　　高斯列主元法解线性方程组过程是：先将线性方程组做成增广矩阵，再对增广矩阵进行初等行变换。对元素，在第i列中，第i行及以下的元素选取绝对值最大的元素，将该最大元素的行与第i行交换，然后采用高斯消元法将新得到的消去第i行以下的元素。这样一次进行直到。从而得到上三角矩阵。再对得到的上三角矩阵进行回代操作，即可以得到方程组的解。 2.2 高斯列主元法解线性方程组流程图　　　　　　　　　　　　　　　　　　　　① ① 图2-1 高斯列主元法解线性方程组流程图 2.3 高斯列主元法解线性方程组程序调试图 图2-2 高斯列主元法解线性方