数值计算课程设计(陕西科技大学).docVIP

1、 经典四阶龙格库塔法解一阶微分方程 1.1、 算法说明 龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。由于此算法精度高，采取措施对误差进行抑制，所以其实现原理也较复杂。该算法是构建在数学支持的基础之上的。 开始，利用下面的计算方法生成近似序列 (1-1) 1.2、 经典四阶龙格库塔法解一阶微分方程算法流程图 1.3、 经典四阶龙格库塔法解一阶微分方程程序调试 图1-2 经典四阶龙格库塔法解一阶微分方程程序调试 1.4、 经典四阶龙格库塔法解一阶微分方程代码 #include iostream #include iomanip using namespace std; //f为函数的入口地址，x0、y0为初值，xn为所求点，step为计算次数 double Runge_Kuta( double (*f)(double x, double y), double x0, double y0, double xn, int step ) { double k1,k2,k3,k4,result; double h=(xn-x0)/step; if(step=0) return(y0); if(step==1) { k1=f(x0,y0); k2=f(x0+h/2, y0+h*k1/2); k3=f(x0+h/2, y0+h*k2/2); k4=f(x0+h, y0+h*k3); result=y0+h*(k1+2*k2+2*k3+k4)/6;} else { double x1,y1; x1=xn-h; y1=Runge_Kuta(f, x0, y0, xn-h,step-1); k1=f(x1,y1); k2=f(x1+h/2, y1+h*k1/2); k3=f(x1+h/2, y1+h*k2/2); k4=f(x1+h, y1+h*k3); result=y1+h*(k1+2*k2+2*k3+k4)/6;} return(result);} int main() {double f(double x, double y); double x0,y0; double a,b; // int step; cout请输入初值x0,y0:; cinx0y0; cout请输入区间：; cinab; //double x0=0,y0=1; double x,y,step; int i; cout请输入步长：; cinstep; //step=0.1; cout.precision(10); for(i=0;i=(b-a)/step;i++) {x=x0+i*step; coutsetw(8)xsetw(18)Runge_Kuta(f,x0,y0,x,i)endl;} return(0); } double f(double x, double y) { double r; r=(x-y)/2; return(r); 2、 高斯列主元法解线性方程组 2.1、 算法说明 首先将线性方程组做成增光矩阵，对增广矩阵进行行变换。 对第元素，在第i列中，第i行及以下的元素选取绝对值最大的元素，将该元素最大的行与第i行交换，然后采用高斯消元法将新得到的消去第i行以下的元素。一次进行直到。从而得到上三角矩阵。 再对得到的上三角矩阵进行回代操作，即可以得到方程组的解。 2.2、 高斯列主元算法流程图  2.3、 高斯列主元程序调试 对所编写的高斯列主元程序进行编译和链接，然后执行得如下所示的窗口，我们按命令输入增广矩阵的行数为3，输入3行5列的增广矩阵,运行界面为： 图2-2 高斯列主元程序调试 2.4、 高斯列主元算法代码 #includeiostream #includecmath using namespace std; void load(); const N=20; float a[N][N]; int m; int main() { int i,j; int c,k,n,p,r; float x[N],l[N][N],s,d; cout下面请输入未知数的个数m=; cinm; coutendl; cout请按顺序输入增广矩阵a:endl; load(); for(i=0;im;i++) { for(j=i;jm;j++) c=(fabs(a[j][i])fabs(a[i][i]))?j:i; /*找列最大元素*/ for(n=0;nm+1;n++) {s=a[i][n]; a[i][n]=a[c][n]; a[c][n]=s;}