rowstandard.cpp

来自「本程序将独立分量分析技术和数字水印技术有机地结合在一起」· C++ 代码 · 共 27 行

#include <math.h>
/*使得矩阵每一行都为单位方差*/
void rowstandard (float *ans, int n, int p)
{

	double tmp[2];
	double tmp1;
	int i, j;
	for (i = 0; i < n; i++) {
		tmp[0] = 0;
		tmp[1] = 0;
		
		for (j = 0; j < p; j++) {
			tmp[0] += (double) ans[i * p + j];
			tmp[1] += ((double) ans[i * p + j]) * ((double) ans[i * p + j]);
		}
		
		tmp[0] = tmp[0] / p;
		tmp1 = (tmp[1] - p * (tmp[0]) * (tmp[0])) / (p - 1);
		
		tmp[1] = sqrt (tmp1);
		for (j = 0; j < p; j++) {
			ans[i * p + j] =
				(float) (((double) ans[i * p + j])  / tmp[1]);
		}
	}
}