《数值分析》上机题.docx

资源描述

《数值分析》上机题.docx

《《数值分析》上机题.docx》由会员分享，可在线阅读，更多相关《《数值分析》上机题.docx（74页珍藏版）》请在冰豆网上搜索。

《数值分析》上机题.docx

《数值分析》上机题

《数值分析》上机实习题

一、秦九韶算法

1．源程序

#include

//输出文件指针

FILE*stream=NULL;

typedeffloat*pFloat;

//和用户交互输入多项式各系数

boolInputParam（int&n,pFloat&a,float&x）

{

cout<<"秦九韶算法\n";

cout<<"请输入最高次数n:

cin>>n;

if（n<=0）returnfalse;

a=newfloat[n+1];

for（inti=0;i<=n;i++）

{

cout<<"请输入a["<

cin>>a[i];

}

cout<<"请输入x的值:

cin>>x;

returntrue;

}

//按照格式打印多项式到文件中

voidOutputTheMultinomial（intn,pFloata,floatx）

{

fprintf（stream,"此多项式为:

\n"）;

fprintf（stream,"P（x）="）;

for（inti=n;i>=2;i--）

{

fprintf（stream,"（%f）*x**%d+",a[i],i）;

}

fprintf（stream,"（%f）*x+（%f）\n",a[1],a[0]）;

fprintf（stream,"x=%f",x）;

}

//根据多项式各系数及给定的x值计算

floatQinJiuShao（intn,pFloata,floatx）

{

floatr=a[n];

for（inti=n-1;i>=0;i--）

{

r=x*r+a[i];

}

returnr;

}

//计算第一题

voidCalculate1（）

{

fprintf（stream,"秦九韶计算结果1\n"）;

intn=0;

float*a=NULL;//指向存放多项式系数的指针

floatx=0;

charnext='Y';

//和用户交互进行输入并计算

while（next=='Y'||next=='y'）

{

if（InputParam（n,a,x））

{

//输出多项式

OutputTheMultinomial（n,a,x）;

//输出结果

fprintf（stream,"\n计算结果为:

P（%f）=%f\n",x,QinJiuShao（n,a,x））;

if（a）

{

delete[]a;

a=NULL;

}

cout<<"是否进行下一次运算?

（Y/N）:

cin>>next;

}

//计算第二题

voidCalculate2（）

{

fprintf（stream,"秦九韶计算结果2\n"）;

//多项式系数

floata[]={2.0f,-3.36f,2.81f,-2.74f,-3.27f,0.22f};

intn=5;

floatx=0;

//按每次递增0.05计算

for（inti=0;i<=20;i++）

{

x=0.8+0.05*i;

//输出多项式

OutputTheMultinomial（n,a,x）;

//输出结果

fprintf（stream,"\n计算结果为:

P（%f）=%f\n",x,QinJiuShao（n,a,x））;

}

voidmain（）

{

//打开输出结果文件

if（（stream=fopen（"result.txt","w"））==NULL）

{

printf（"打不开输出文件\n"）;

return;

}

//通过用户输入计算第一题

Calculate1（）;

//通过循环生成x值计算第二题

Calculate2（）;

//关闭文件

if（stream）

fclose（stream）;

}

2．输出结果

秦九韶计算结果1

此多项式为:

P（x）=（7.540000）*x**4+（11.080000）*x**3+（3.820000）*x**2+（0.440000）*x+（-0.480000）

x=3.250000

计算结果为:

P（3.250000）=1262.865112

此多项式为:

P（x）=（2.790000）*x**4+（9.850000）*x**3+（14.150000）*x**2+（5.380000）*x+（7.240000）

x=3.250000

计算结果为:

P（3.250000）=823.586670

此多项式为:

P（x）=（9.360000）*x**4+（12.690000）*x**3+（14.390000）*x**2+（0.750000）*x+（-0.940000）

x=3.250000

计算结果为:

P（3.250000）=1633.377319

此多项式为:

P（x）=（3.450000）*x**4+（-2.910000）*x**3+（3.710000）*x**2+（-6.750000）*x+（-2.380000）

x=3.250000

计算结果为:

P（3.250000）=299.878632

此多项式为:

P（x）=（12.780000）*x**4+（14.350000）*x**3+（17.190001）*x**2+（1.340000）*x+（-1.720000）

x=3.250000

计算结果为:

P（3.250000）=2102.631592

此多项式为:

P（x）=（4.790000）*x**4+（5.380000）*x**3+（-2.860000）*x**2+（7.310000）*x+（4.550000）

x=3.250000

计算结果为:

P（3.250000）=717.187134

此多项式为:

P（x）=（15.650000）*x**4+（17.580000）*x**3+（21.700000）*x**2+（2.780000）*x+（1.340000）

x=3.250000

计算结果为:

P（3.250000）=2589.083984

此多项式为:

P（x）=（8.340000）*x**4+（-7.750000）*x**3+（4.530000）*x**2+（-9.290000）*x+（5.750000）

x=3.250000

计算结果为:

P（3.250000）=687.826416

秦九韶计算结果2

此多项式为: