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

版权声明:本文为博主原创文章,转载请注明来源。开发合作联系luanpenguestc@sina.com https://blog.csdn.net/luanpeng825485697/article/details/77803403

全栈工程师开发手册 (作者:栾鹏)

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;
}
阅读更多

没有更多推荐了,返回首页