main.cpp

变分独立因子分析C++代码,H.Attias说比独立分量分析要好,但是这个程序分析效果不好,可能是程序问题,也可能是对理论理解不透

#include "matrix.h"
#include "ifa.h"
#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;


////////////////////////////////////////////////////////////////////
void displaymatrix(double** matrix, int nr, int nc)
{
	for(int i=0;i<nr;i++)
	{
		for(int j=0;j<nc;j++)
		{
			cout<<matrix[i][j]<<"  ";
		}cout<<endl;
	}cout<<endl;
}

void display(const CMatrix & matrix)
{
	displaymatrix(matrix.getmatrix(),matrix.getnr(),matrix.getnc());
}
/////////////////////////////////////////////////////////////////////


void main()
{	
	double** a=new double*[3];a[0]=new double[3];a[1]=new double[3];a[2]=new double[3];
	a[0][0]=1;a[0][1]=2;a[0][2]=3;
	a[1][0]=4;a[1][1]=5;a[1][2]=6;
	a[2][0]=7;a[2][1]=8;a[2][2]=0;

	CMatrix y(a,3,3);

	CMatrix sigma=~y*y/100;display(sigma);
	CMatrix x(sigma.getmatrix(),3,1);
	CMatrix mu(y.getmatrix(),3,1);
	display(x);
	display(mu);
	cout<<Gauss(x,0.1*mu,sigma)<<endl;
	symmdeterminant(2*sigma);


	
}