實驗一 larange插值

很簡單，隨便寫了

#include <iostream>
#include <stdio.h>
#include <cstdio>
using namespace std;

int main()
{
	float l[10],x[10],xx,y[10],sum;     //這裡要用浮點型儲存。
	int n,i,j;
	cout<<"請輸入已知點的個數"<<endl;
	cin>>n;
	cout<<"請輸入要計算的點的數值(x)"<<endl;
	cin>>xx;
	cout<<"已知的點(x)和(y)"<<endl;
	
	for(i=0;i<n;i++)
	{
		l[i]=1.0;                 //賦初值，不然的話會錯錯的
	}

	for(i=0;i<n;i++)
		cin>>x[i]>>y[i];         //輸入，沒什麼好說的
	for(j=0;j<n;j++)
	{
		for(i=0;i<n;i++)
		{
			if(i==j)
			{
				i++;               //因為x[j]不能等於x[i]，小學生都知道的事情。所以遇到i=j的情況就i++好了。
			}
			l[j]=l[j]*((xx-x[i])/(x[j]-x[i]));         //按照指導書的公式敲成程式碼，用兩重迴圈就可以了。
		}
	}
	sum=0.0;        //迭代前賦初值
	for(i=0;i<n;i++)
		sum=sum+l[i]*y[i];       //公式，沒什麼好說的
	printf("%f\n",sum);               
		return 0;
}

 

測試資料

3
11.5
11 0.190809
12 0.207912
13 0.224951

答案：0.199368

2.Newton  

#include <iostream>
#include <stdio.h>
#include <cstdio>
using namespace std;

int main()
{
	float l[10],xx,a[10][10],sum,ca[10];
	int n,i,j,k;
	cout<<"請輸入已知點的個數"<<endl;
	cin>>n;
	cout<<"請輸入要計算的點的數值(x)"<<endl;
	cin>>xx;
	cout<<"已知的點(x)和(y)"<<endl;
	
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=2;j++)
		{
		cin>>a[i][j];      //我這裡用二維陣列來存目的是為了打出差商表
		}
	}

	ca[0]=1;
	for(i=0;i<n;i++)
	{
		ca[i+1]=ca[i]*(xx-a[i+1][1]);     //這裡就是（x-x0）*(x-x1)....（x-xn）的意思
	}
	
	for(j=3;j<=n+1;j++)
	{
		k=1;                                                             
		for(i=j-1;i<=n;i++)
		{
			a[i][j]=(a[i][j-1]-a[i-1][j-1])/(a[i][1]-a[k][1]);       //打差商表，推導過程很簡單。用上幾組資料找出規律就好了
			k++;
		}
	}

	for(i=1;i<=n;i++)
	{
		l[i]=a[i][i+1];          //把f（x0）,f(x1)....賦予l[i]
	}
	
	sum=l[1];
	for(i=1;i<n;i++)
	{
		sum=sum+ca[i]*l[i+1];    //最後加起來得出答案
	}
	printf("%f\n",sum);
	return 0;
}
 

樓主在做這題的時候參考了老師的PPt   ch01第36頁，擦。這組資料用3來計算答案是50.但是，老師的PPT 有錯，這裡錯了。所以答案是52.

這組資料答案是52.       

4
3
-2 17
0 1
1 2
2 17

實驗測試資料

4
0.596
0.4 0.41075
0.55 0.57815
0.65 0.69675
0.80 0.88811  

答案：0.63194.

實驗三

1.二分法

按上面的變成程式碼就可以了，沒什麼好說的

#include <cstdio>
#include <algorithm>
#include <iostream>
#include <string.h>
#include <cstring>
#include <cmath>
#include <math.h>
using namespace std;

double f(double x)
{
	return x*x*x-x-1; //函式
}

int main()
{
	double a=1.0,b=1.5,x,e=0.00001;
	if(f(a)*f(b)<0)
	{
		while(fabs(a-b)>e)
		{
			x=(a+b)/2;
			if(fabs(f(x))<e)
			{
				break;
			}
			else
			{
				if(f(a)*f(x)<0)
				{
					b=x;
				}
				else if(f(x)*f(b)<0)
				{
					a=x;
				}
			}
		}
		x=(a+b)/2;
		printf("%0.5lf\n",x); 
	}
	return 0;
}
 

答案：1.32472   如果是精確到小數點後兩位那麼就是1.32

可是。你大爺 你的答案居然是就算是精確到後兩位也沒有差上0.002吧！無語了，什麼答案。有圖為證，答案錯了。就是這樣。

2.牛頓迭代法

這題很簡單，隨便寫了，過程很簡單。

沒心情去弄成通用的了。如果你想弄也很簡單的。

#include <iostream>
#include <stdio.h>
#include <math.h> 
using namespace std;


double f(double x)
{
	return x*x*x+x*x-3*x-3;
}

double f1(double x)
{
	return 3*x*x+2*x-3;
}

int main()
{
	int m,N;
	double x,xx,n,e;
	cout<<"請輸入初值"<<endl;
	cin>>x;
	cout<<"請輸入根的容許誤差"<<endl;
	cin>>e;
	cout<<"請輸入|f(x)|的容許誤差"<<endl;
	cin>>n;
	cout<<"請輸入迭代次數的容許值"<<endl;
	cin>>N;
	m=0;
	while(f1(x)!=0 || m<=N)
	{
		xx=x-f(x)/f1(x);
		if(fabs(xx-x)==e || f(xx)<n)
		{
			printf("%0.6lf\n",xx);
			break;
		}
		x=xx;
		m++;
	}
	return 0;
}

答案：1.732051

後續更新

福利