//=======================================================================
//函數說明
//函數名稱：Correlation
//函數功能：計算最小二乘法擬合的多項式的相關系數
//使用方法：int?M------?擬合多項式的階數（已知條件）
//??????????double?*b---?擬合曲線的系數排列順序為由高階到低階（已知條件）
//			double?*x---?結點x軸數據（已知條件）
//			double?*y---?結點y軸數據（已知條件）
//			double?*Yg--?結點估計值，個數為m（過程變量）
//			int?m------?結點個數（已知條件）
//注意事項：多項式階數最高為10多項式的形式為?y?=?a0?+?a1x?+a2x2
//=======================================================================
#include?
double?Correlation（int?M?double?*b?double?*x?double?*y?double?*Yg?int?m）
{
	int?ij;
	double?power[10]tempYavrQrQtR;?//最大10階

	Yavr?=?0.0;	?????????//求y的平均值
	for（j=0;j	{
		Yavr?+=?y[j];
		Yg[j]?=?0.0;?????//初始化Yg
	}
	Yavr?=?Yavr/m;

	Qr?=?0.0;
	Qt?=?0.0;
	for（j=0;j	{
		temp?=?1.0;
		for（i=0;i		{
			temp?*=?x[j];
			power[i]?=?temp;??
			Yg[j]=?Yg[j]?+?b[M-i-1]*power[i];
		}
		Yg[j]?=?Yg[j]?+?b[i];

		Qr?=?Qr?+?（Yg[j]?-Yavr）*（Yg[j]?-Yavr）;??//求Qr
		Qt?=?Qt?+?（y[j]?-?Yavr）*（y[j]?-?Yavr）;??//求Qt
	}

	R?=?sqrt（Qr/Qt）;

	return（R）;
}

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????目錄??????????0??2005-07-07?12:17??算法子程序

?????文件???????1244??2005-07-07?12:14??算法子程序\Correlation（多項式相關系數的計算方法（多項式形式1））.c

?????文件???????1341??2005-07-07?12:16??算法子程序\Correlation（多項式相關系數的計算方法（多項式形式2））.c

?????文件???????2570??2005-07-07?12:12??算法子程序\PolyFit（最小二乘法曲線擬合）.C

?????文件???????2305??2005-07-07?12:15??算法子程序\Spline（三次樣條插值（自然邊界條件））.c

-----------?---------??----------?-----??----

?????????????????7460????????????????????5