西电常微分方程的数值解报告.docxVIP

计算方法实验题目：班级：学号：姓名：1 实验目的通过编程加深对欧拉法、改进欧拉法和四阶龙格-库塔法求解常系数微分方程的理解观察上述三种方法的计算稳定性和求解精度并比较各种方法利弊2 实验步骤2.1环境配置：VS2013，C++控制台程序2.2添加头文件#include stdio.h#include stdlib.h#include stdafx.h2.3主要模块程序一共分成三层，最底层是数据结构部分，负责存储数据，第二层是交互部分，即多项式方程部分，负责输入输出获得数据，最上层是核心的算法部分，负责处理已获得的数据。具体功能如下：数据结构部分数据结构部分是整个程序的最底层，负责存储部分。因方程系数作为数据元素插入和删除操作较少，而顺序表空间利用率大且查看方便，故此程序选用顺序表保存系数。数据结构文件中写的是有关顺序表的所有基本操作以供其他文件调用。多项式方程部分多项式方程部分是程序的第二层，内容是有关方程的所有函数、构建方程、输出方程、方程求导等等，同时在此文件中获得方程系数和导数系数并储存，同时此文件还负责显示菜单。算法部分此文件负责核心算法，处于整个程序最上层部分，负责二分法、牛顿迭代法和割线法的具体实现过程。调用方程文件的函数，将获得的数据进行处理运算，将结果返回给方程主函数和输出的第二层。主函数负责获取方程系数并显示，算法和方程作为后台程序，顺序表作为存储手段。和方程求根思路相同，整体而言更简单。3 代码3.1主程序部分// numerical_solutions_of_linear_equations.cpp : 定义控制台应用程序的入口点。//#include stdafx.h#include equation.hint _tmain(int argc, _TCHAR* argv[]){GetEquation();while (Exflag)ShowMenu();return 0;}3.2多项式方程部分方程部分头文件#ifndef _EQUATION_H#define _EQUATION_H#include squencelist.hextern datatype Firsty;extern datatype Dx;extern datatype Xl;extern datatype Xh;extern int Exflag;extern int N;extern sequenlist *B;extern sequenlist *D;extern sequenlist *L;datatype Function(sequenlist *A, datatype x, datatype y);void GetEquation(void);void ShowMenu(void);void printfunctionf(sequenlist *A);void Tip(void);#endif方程部分CPP文件#include stdafx.h#include equation.h#include squencelist.h#include algorithm.h//全局变量datatype Firsty=0;datatype Dx = 0;datatype Xl = 0;datatype Xh = 0;int N = 0;int Exflag = 1;sequenlist *B;sequenlist *D;sequenlist *L;//////////////////////////多项式函数///////////////////////////datatype Function(sequenlist *A, datatype x,datatype y){double f = 0.00;f = A-data[1] * pow(x, A-data[2])*pow(y, A-data[3]);return f;}////////////////////////获得给定数据/////////////////////////void GetEquation(void){int j = 0;datatype x = 0;B = InitList();D = InitList();cout 输入方程导数,输入00结束（如y=(-5)x^2*y^6输入-5 2 6 00）: endl;cin x;while (x != 00){for (j = 1; j = 3; j++){if (!Insert(B, x, j))exit(0);cin x;}}cout 输入x的初值和终值: endl;cin XlXh;while (Xl = Xh){cout 输入x的初值和终值: endl;cin Xl Xh;