基于C语言的科学与工程计算实验报告（非线性方程求根、线性方程组、函数插值、最小二乘数拟合、数值积分等）

约定不等于承诺〃 2024-04-01 09:50 68阅读 0赞

摘 要  
实验一：非线性方程求根  
\------二分法+牛顿法+牛顿下山法+弦截法  
实验二：线性方程组的直接解法  
\------改进的Cholesky方法  
实验三：线性方程组的迭代解法  
\------雅可比+SOR+共轭梯度法  
实验四：函数插值  
\------三次样条插值+拉格朗日插值法  
实验五：最小二乘数拟合  
\------最小二乘法+共轭梯度法+三次样条插值+拉格朗日插值法  
实验六：数值积分  
\------龙格公式  
关键词：二分法、牛顿法、牛顿下山法、共轭梯度法、改进的Cholesky方法、雅可比、SOR、共轭梯度法、三次样条插值、拉格朗日插值、龙格公式  
目 录  
实验一、非线性方程求根 5  
一、实验内容： 5  
二、算法基本思想及复杂度分析： 5  
三、源程序及注释： 5  
四、运行输出结果： 13  
五、调试和运行程序过程中产生的问题、采取的措施及获得的相关经验教训： 16  
实验二、线性方程组的直接解法 18  
一、实验内容： 18  
二、算法基本思想及复杂度分析： 18  
三、源程序及注释： 18  
四、运行输出结果： 21  
五、调试和运行程序过程中产生的问题、采取的措施及获得的相关经验教训： 22  
实验三、线性方程组的迭代解法 23  
一、实验内容： 23  
二、算法基本思想及复杂度分析： 23  
三、源程序及注释： 24  
四、运行输出结果： 38  
五、调试和运行程序过程中产生的问题、采取的措施及获得的相关经验教训： 43  
实验四、函数插值 44  
一、实验内容： 44  
二、算法基本思想及复杂度分析： 44  
三、源程序及注释： 44  
四、运行输出结果： 47  
五、调试和运行程序过程中产生的问题、采取的措施及获得的相关经验教训： 48  
实验五、最小二乘拟合 49  
一、实验内容： 49  
二、算法基本思想及复杂度分析： 49  
三、源程序及注释： 50  
四、运行输出结果： 57  
五、调试和运行程序过程中产生的问题、采取的措施及获得的相关经验教训： 68  
实验六、数值积分 69  
一、实验内容： 69  
二、算法基本思想及复杂度分析： 69  
三、源程序及注释： 69  
四、运行输出结果： 72  
五、调试和运行程序过程中产生的问题、采取的措施及获得的相关经验教训： 72  
心得体会 73  
实验一、非线性方程求根验内容：  
二、算法基本思想及复杂度分析：  
（一）二分法  
将有根区间对分，并找出根所在的小区间，然后再对该小区间对分，依次类推，直到有根区间的长度足够小为止。  
（二）牛顿法  
牛顿法是一种线性化方法，其基本思想是将非线性方程f(x)=0逐步归结为某种线性方程来求解。  
（三）牛顿下山法  
要求每一步迭代满足下降条件（保证全局有效）。也就是由于牛顿法的复杂度依赖于x0的选取，所以在牛顿法中加入了下山因子，本文转载自http://www.biyezuopin.vip/onews.asp?id=16948且逐次减半。  
（四）弦截法  
当f(x)比较复杂的时候，f’(x)会比较难算，该方法就是用差商来取代牛顿法中的  
f’(x)。

1.//实验一（一）牛顿下山法  
    2.//x1=x0-f(x0)/f'(x0)  
    3.  
    4.#include <iostream>  
    5.#include <iomanip>  
    6.#include <cmath>  
    7.using namespace std;  
    8.  
    9.//f(x)=x^2-3x+2-e^x  
    10.double f(double x){
          
    11.    return x*x-3*x+2-exp(x);  
    12.}  
    13.  
    14.//f'(x)=2x-3-e^x  
    15.double f_(double x){
          
    16.    return 2*x-3-exp(x);  
    17.}  
    18.  
    19.int main(){
          
    20.    double lambda=1.0;  //下山因子  
    21.  
    22.    double x[100]={
        0};  
    23.    x[0]=0;  
    24.    x[1]=x[0]-f(x[0])/f_(x[0]);  
    25.  
    26.    //循环并判断精度  
    27.    int k=1;  
    28.    cout<<"x0: "<<x[0]<<endl;  
    29.    while(fabs(x[k]-x[k-1])>=1e-8){
         //设置精度  
    30.        x[k+1]=x[k]-f(x[k])/f_(x[k]);  
    31.        while(fabs(f(x[k+1]))>=fabs(f(x[k]))){
          
    32.            lambda/=2.0;  //下山因子逐次减半  
    33.            x[k+1]=x[k]-lambda*f(x[k])/f_(x[k]);  
    34.        }  
    35.        cout<<"x"<<k<<": "<<setprecision(10)<<x[k]<<endl;  
    36.        k++;  
    37.    }  
    38.    cout<<"x"<<k<<": "<<setprecision(10)<<x[k]<<endl;  
    39.    return 0;  
    40.}  
    
    4. 弦截法
    1.//实验一（一）弦截法  
    2.//x[k+1]=x[k]-f(x[k])/(f(x[k])-f(x[k-1]))*(x[k]-x[k-1])  
    3.  
    4.#include <iostream>  
    5.#include <iomanip>  
    6.#include <cmath>  
    7.using namespace std;  
    8.  
    9.//f(x)=x^2-3x+2-e^x  
    10.double f(double x){
          
    11.    return x*x-3*x+2-exp(x);  
    12.}  
    13.  
    14.int main()  
    15.{
          
    16.    double x[100]={
        0};  
    17.    x[0]=0;     //近似根  
    18.    x[1]=1;     //近似根  
    19.    int k=1;    //计数  
    20.  
    21.    cout<<"x0: "<<setprecision(10)<<x[0]<<endl;  
    22.  
    23.    while(fabs(x[k]-x[k-1])>=1e-8){
         //设置精度  
    24.        cout<<"x"<<k<<": "<<setprecision(10)<<x[k]<<endl;  
    25.        x[k+1]=x[k]-f(x[k])/(f(x[k])-f(x[k-1]))*(x[k]-x[k-1]);  
    26.        k++;  
    27.    }  
    28.    cout<<"x"<<k<<": "<<setprecision(10)<<x[k]<<endl;  
    29.    return 0;  
    30.}

![在这里插入图片描述][475e931eea3f46cc8113ea9a8b613b32.png_pic_center]  
![在这里插入图片描述][421e2556cb2542c799297d79fda394cc.png_pic_center]  
![在这里插入图片描述][13d5beda944d4bf4af7624bd86a0bcf8.png_pic_center]  
![在这里插入图片描述][c79030038b6141ceaf85950663d0f5d8.png_pic_center]  
![在这里插入图片描述][06036136a2e544d3bf58f4cbac2fb995.png_pic_center]  
![在这里插入图片描述][7ea28854b6d24cc8827aaeed6d6c30bb.png_pic_center]  
![在这里插入图片描述][a6169602534b4935b6822eac5244641d.png_pic_center]  
![在这里插入图片描述][6588c50a917e40b9af27f7adc1bb2513.png_pic_center]  
![在这里插入图片描述][658bf10bc5264d15975723498aea3768.png_pic_center]  
![在这里插入图片描述][474b501f67ac4fdc9b9e4552a5270932.png_pic_center]  
![在这里插入图片描述][7eacf3b1fe0545ceb1896e2850f5b5fa.png_pic_center]  
![在这里插入图片描述][530e5064613047f080c9648563628fb5.png_pic_center]  
![在这里插入图片描述][24a6df082d4f4b55975b83a227eb0c33.png_pic_center]  
![在这里插入图片描述][7b188a52b0234e16a7bf73153cc0750e.png_pic_center]

[475e931eea3f46cc8113ea9a8b613b32.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/420498349ca747768bc1093d074a57ff.png
[421e2556cb2542c799297d79fda394cc.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/3cdafbf1fca6465890f7296a5c5feae6.png
[13d5beda944d4bf4af7624bd86a0bcf8.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/37c4dc562bb14c9ebf9e4652ed1d46d8.png
[c79030038b6141ceaf85950663d0f5d8.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/570290ce109e4adfbf8c66753cc83fce.png
[06036136a2e544d3bf58f4cbac2fb995.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/e77af76ef9ef47d6975163a8b71758c9.png
[7ea28854b6d24cc8827aaeed6d6c30bb.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/f7d6534f936c4b7f842e72a4125d7e76.png
[a6169602534b4935b6822eac5244641d.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/ca37dceb2a394f0fbda2ad21b99fa4fd.png
[6588c50a917e40b9af27f7adc1bb2513.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/44eba2c7c0224940bf8e8f3e50903543.png
[658bf10bc5264d15975723498aea3768.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/0f2a1c1d3822479ebb6b80a6891979c7.png
[474b501f67ac4fdc9b9e4552a5270932.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/20539b5ccb7544d89cd02c099c69de90.png
[7eacf3b1fe0545ceb1896e2850f5b5fa.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/eec0a33eb40c48a78e795100cae0080f.png
[530e5064613047f080c9648563628fb5.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/05415911130c4c5e833658ba3bcf2f93.png
[24a6df082d4f4b55975b83a227eb0c33.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/bde814f5ddc94d63ad85cd6354b45a8d.png
[7b188a52b0234e16a7bf73153cc0750e.png_pic_center]: https://image.dandelioncloud.cn/pgy_files/images/2024/04/01/f376cf38309144deae740ee1e267b9ed.png