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

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

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

1、kron函数：求两个矩阵的Kronecker积

2、用法说明

 格式 C=kron (A,B)    %A为m×n矩阵，B为p×q矩阵，则C为mp×nq矩阵。
 kron即为Kronecker积，所谓Kronecker积是一种矩阵运算，其定义可以简单描述成：
 X与Y的Kronecker积的结果是一个矩阵：
 X11*Y   X12*Y … X1n*Y
 X21*Y   X22*Y … X2n*Y
 ……
 Xm1*Y   Xm2*Y … Xmn*Y

kron的c++源码实现

矩阵的张量积
两个矩阵两两元素(一个来自a一个来自b)可能的乘机

Matrix Matlab2c::kron(Matrix& a,Matrix& b)
{
	Matrix p(a.row*b.row,a.column*b.column);
	for (int i=0;i<a.row;i++)
		for (int j=0;j<a.column;j++)
			for (int ii=0;ii<b.row;ii++)
				for (int jj=0;jj<b.column;jj++)
				{
					//行i*b.row+ii，列j*b.column+jj
					p(i*b.row+ii,j*b.column+jj)=a(i,j)*b(ii,jj);
				}
	return p;
}

kron函数的使用测试

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


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

	system("pause");
	return 0;
}