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

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

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中harmmean函数简介

1、harmmean函数：

矩阵按列求调和平均数

2、用法说明

B =harmmean(A)返回矩阵按列求调和平均数的结果。
如果A是一个向量，将返回一个包含A各元素求调和平均数的结果，一个double值。
如果A是一个矩阵，将返回一个行向量，包含A每一列向量求调和平均数的结果。

harmmean的c++源码实现

计算矩阵各列的调和平均数，输入为行向量时返回行向量的调和平均数
输入：矩阵a
输出：矩阵a各列的调和平均数
提示：做调和运算要求矩阵元素不能为0

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

CMatrix Matlab2c::harmmean(CMatrix& a)
{
	Complex m=0;
	if(a.row==1)  //若数据矩阵为行向量
	{
		CMatrix p=CMatrix(1,1);
		for (int j=0;j<a.column;j++)
			m +=1.0/(a.data[j]);
		p.data[0]=(Complex)a.column/m;
		return p;
	} 
	CMatrix p(1,a.column);
	for (int j=0;j<a.column;j++)
	{
		m=1;
		for (int i=0;i<a.row;i++)
			m +=1.0/(a(i,j));
		p.data[j]=(Complex)a.row/m;
	}
	return p;
}

harmmean函数的使用测试

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


int main()
{
	double a[]={1,2,3,4,5,6};
	Matrix aa=Matrix(1,6,a);
	Matrix bb=Matrix(2,3,a);
	Matrix cc=Matlab2c::harmmean(aa);
	cout<<cc.toString()<<endl;
	Matrix dd=Matlab2c::harmmean(bb);
	cout<<dd.toString()<<endl;


	system("pause");
	return 0;
}