计算方法实验报告PC版 —王世儒、王金金

最新推荐文章于 2024-02-19 16:25:40 发布

而濡木染

最新推荐文章于 2024-02-19 16:25:40 发布

阅读量2.9k

点赞数 4

分类专栏： ------------专业知识储备---------- 计算方法 Nutritious_Planet 文章标签：实验报告

本文链接：https://blog.csdn.net/PNAN222/article/details/50452188

版权

这篇实验报告探讨了计算方法中的多项式求和算法，对比了秦九韶算法与普通算法的效率，并通过C语言实现验证。此外，还介绍了用雅克比迭代法和高斯-赛德尔迭代法解线性方程组的步骤与结果，展示了这两种迭代法在求解过程中的迭代次数和精度。最后，报告探讨了数值积分的梯形公式、复化梯形公式和复化辛甫生公式在求解定积分中的应用，以及它们在求解圆面积问题上的精度表现，得出复化新甫生公式更为精确的结论。

摘要由CSDN通过智能技术生成

计算方法实验

—（14级计算机科学与技术2班）

实验仪器和环境：PC-机、codeblocks 10.05、MATLAB R2013a

编程语言：C、C++、MATLAB

实验一：编写程序求多项式的和

实验目的：

比较两种求和方法发计算复杂度，通过比较两种程序的运行时间来实现。

实验内容：

计算多项式：

方法一：一般算法，程序运行过程中，所算乘法次数为，由于所加项共(n+1)项，则加法次数为n，C语言代码实现如下:

#include <iostream>
#include <cstdio>
#include <ctime>

using namespace std;

double a[100005];
int main()
{

   double second1,second2;
   second1=(double)clock()/CLOCKS_PER_SEC;//clock()函数返回值本为毫秒，除以CLOCK_PER_SEC之后单位即变为秒
   int n,x;
   double mul,sum;
   n=10000;//记得给n更换取值看所用时间的差别
   x=1;//防止数据过大超范围，将x初始化为1
   sum=0;
   for(int i=0;i<=n;i++)
    {
       a[i]=1;//防止数据过大超范围，初始化为1
    }
   for(int i=0;i<=n;i++)
    {
       mul=1;
       for(int j=0;j<i;j++)
       {
           mul*=x;
       }
           mul*=a[i];
           sum+=mul;
    }
   cout<<"Pn(x)="<<sum<<endl;
   second2=(double)clock()/CLOCKS_PER_SEC;
   cout<<"second1="<<second1<<endl;
   cout<<"second2="<<second2<<endl;
   cout<<"用时为： "<<second2-second1<<endl;
   return 0;
}

 
n=10000时，输出结果为：
Pn(x)=10001
second1=0
second2=0.465
用时为： 0.465

 
n=100000时，输出结果为：
Pn(x)=100001
second1=0
second2=29.317  用时为： 29.317

方法二：秦九昭算法，已知递推公式为：，u[n]即为所求结果，则整个过程中共进行了n次加法和n次乘法，C语言代码实现如下：

#include <iostream>

#include <cstdio>

#include <ctime>

#include <cstring>

using namespace std;

double a[100005];

double u[100005];

int main()

{

double second1,second2;

second1=(double)clock()/CLOCKS_PER_SEC; //clock()函数返回值本为毫秒，除以CLOCK_PER_SEC之后单位即变为秒