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

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

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

1、median函数：

返回一个数组不同维中的中位数。

2、用法说明

C = median(A)
如果A是一个向量，median(A)返回A中的中位数。
如果A是一个矩阵，median(A)将A的每一列作为一个向量，返回一行向量包含了每一列的中位数。
如果A是一个矩阵，median(A,dim)，dim=1按列求，dim=2按行求

median的c++源码实现

返回各列中位数组成的行向量，
输入为行向量，返回输入的中位数
输入为一个矩阵，输出为每列的中位数向量。
dim=1（默认）按列求中位数，dim=2按行求中位数

函数声明
Matrix median(Matrix& a,int dim=1);
函数定义
Matrix Matlab2c::median(Matrix& a,int dim)
{
  int lp=0;
  Matrix c = Matlab2c::sort(a,dim);   //对没一列或每一行按升序排序
  Matrix p;
  
  if(a.isVector())  //若数据矩阵为向量
  {
    p=Matrix(1,1);
    if ((a.row*a.column)%2==1)
      p.data[0]=c.data[(a.row*a.column/2)];
    else
      p.data[0]=0.5*(c.data[(a.row*a.column/2)])+0.5*(c.data[(a.row*a.column/2-1)]);
    return p;
  }
  if (dim==1)
  {
    p=Matrix(1,a.column);
    if (a.row%2==1)
      for (int j=0;j<p.column;j++)
        p.data[j]=c(a.row/2,j);
    else
      for (int j=0;j<p.column;j++)
        p.data[j]=0.5*c(a.row/2,j)+0.5*c(a.row/2-1,j);
    return p;
  }else
  {
    p=Matrix(a.row,1);
    if (a.column%2==1)
      for (int i=0;i<p.row;i++)
        p.data[i]=c(i,a.column/2);
    else
      for (int i=0;i<p.row;i++)
        p.data[i]=0.5*c(i,a.column/2)+0.5*c(i,a.column/2-1);
    return p;
  }
}

median函数的使用测试

#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::median(aa);
  cout<<cc.toString()<<endl;
  Matrix dd=Matlab2c::median(bb,1);
  cout<<dd.toString()<<endl;
  Matrix ee=Matlab2c::median(bb,2);
  cout<<ee.toString()<<endl;

  system("pause");
  return 0;
}