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

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

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

1、hankel函数：生成Hankel矩阵。Hankel矩阵的构成方式是：先指定第一列和最后一行，矩阵中其他所有元素都与其左下角相邻位置的元素相等，即H(i,j)=H(i+1,j-1)。

2、用法说明

 1．H=hankel(c,r);   生成非对称的Hankel矩阵。第一列元素为c，最后一行元素为r，若c的最后一个元素与r的第一个元素不相等，交叉位置的元素取c的最后一个元素。
 
 2．H=hankel(c);    生成对称的Hankel矩阵。第一列元素为c，次对角线下方所有元素均为零。

hankel的c++源码实现

生成Hankel方阵
输入为两个行向量，输出为其Hankel矩阵。

Matrix Matlab2c::hankel(Matrix& a,Matrix& b)
{
	Matrix p(a.column,b.column);
	int i,j;
	if (a(a.column-1)!=b(0))
	{
		throw_logic_error("第一个向量最后一个元素和第二个向量第一个元素不相等");
	}
	for (i=0;i<p.row;i++)
		for (j=0;j<p.column;j++)
			if ((i+j)<(a.column))
				p(i,j)=a(i+j);
			else
				p(i,j)=b(i+j+1-a.column);
	return p;
}
Matrix Matlab2c::hankel(Matrix& a)
{
	Matrix p(a.column,a.column);
	int i,j;
	for (i=0;i<p.row;i++)
		for (j=0;j<p.column;j++)
			if ((i+j)<(a.column))
				p(i,j)=a(i+j);
			else
				p(i,j)=0;
	return p;
}

hankel函数的使用测试

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

int main()
{
	double a[]={1,2,3};
	double b[]={3,4,5,6,7};
	Matrix aa=Matrix(1,3,a);
	Matrix bb=Matrix(1,5,b);
	Matrix cc1=Matlab2c::hankel(aa,bb);
	cout<<cc1.toString()<<endl;
	Matrix cc2=Matlab2c::hankel(aa);
	cout<<cc2.toString()<<endl;

	system("pause");
	return 0;

}