matlab2c使用c++实现matlab函数系列教程-var函数

全栈工程师开发手册 （作者：栾鹏）

matlab2c动态链接库下载
matlab库函数大全
matlab2c基础教程
matlab2c开发全解教程

matlab2c调用方法：

1、下载动态链接库
2、将Matlab2c.dll拷贝到exe同目录下
3、将Matlab2c.h、Matlab2c.lib放到项目头文件目录下
4、在cpp文件中引入下面的代码

#include "Matlab2c.h"
#pragma comment(lib,"Matlab2c.lib")  
using namespace Matlab2c;

---

matlab中var函数简介

1、var函数：

求样本方差

2、用法说明

D=var(X) %var(X)=,若X为向量，则返回向量的样本方差。
D=var(A) %A为矩阵，则D为A的列向量的样本方差构成的行向量。
D=var(X, 1) %返回向量（矩阵）X的简单方差（即置前因子为的方差）
D=var(X, w) %返回向量（矩阵）X的以w为权重的方差

var的c++源码实现

返回按列方差向量，输入为行向量时，返回行向量的方差
输入为一个矩阵，输出为每列的方差向量。
实现的为matlab中D=var(X, 1) 函数

Matrix Matlab2c::var(Matrix& a)
{
	Matrix p,m;
	int i,j;
	double sum=0;
	m = mean(a);
	if(a.row==1)  //若数据矩阵为行向量
	{
		p=Matrix(1,1);
		for (j=0;j<a.column;j++)
			sum +=(a.data[j]-m.data[0])*(a.data[j]-m.data[0]);
		p.data[0] = sum/a.column;
		return p;
	}
	p=Matrix(1,a.column);
	for (j=0;j<a.column;j++)
	{
		sum=0;
		for (i=0;i<a.row;i++)
		{
			sum +=(a(i,j)-m.data[j])*(a(i,j)-m.data[j]);
		}
		p.data[j]=sum/a.row;
	}
	return p;
}

var函数的使用测试

#include "Matlab2c.h"
#pragma comment(lib,"Matlab2c.lib")  
using namespace Matlab2c;

int main()
{


	double a[]={1,2,3,4,5,6,1,3,4};
	Matrix aa=Matrix(3,3,a);
	Matrix cc1=Matlab2c::var(aa);
	cout<<cc1.toString()<<endl;


	system("pause");
	return 0;
}