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;
}