全棧工程師開發手冊 （作者：欒鵬）

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中prod函式簡介

1、prod函式：

矩陣按列相乘求積

2、用法說明

B =prod(A)返回矩陣按列相乘求積的結果。
如果A是一個向量，將返回一個包含A各元素相乘求積的結果，一個double值。
如果A是一個矩陣，將返回一個行向量，包含A每一列向量相乘求積的結果。

prod的c++原始碼實現

返回按列求積後的行向量，輸入為行向量，返回行積（1x1矩陣）
輸入為一個矩陣，輸出為每列的乘積向量。返回一個行向量

Matrix Matlab2c::prod(Matrix& a)
{
    Matrix p;
    int i,j;
    double prod=1;
    if(a.row==1)  //若資料矩陣為行向量
    {
        p=Matrix(1,1);
        for (j=0;j<a.column;j++)
            prod=prod*a.data[j];
        p.
 
data[0] = prod;
        return p;
    }
    p=Matrix(1,a.column);
    for (j=0;j<a.column;j++)
    {
        prod=1;
        for (i=0;i<a.row;i++)
        {
            prod =prod*(a(i,j));
        }
        p.data[j]=prod;
    }
    return p;
}

CMatrix Matlab2c::prod(CMatrix& a)
{
    Complex prod
 
=1;
    if(a.row==1)  //若資料矩陣為行向量
    {
        CMatrix p = CMatrix(1,1);
        for (int j=0;j<a.column;j++)
            prod=prod*a.data[j];
        p.data[0] = prod;
        return p;
    }
    CMatrix p = CMatrix(1,a.column);
    for (int j=0;j<a.column;j++)
    {
        prod=1;
        for (int i=0;i<a.row;i++)
        {
            prod =prod*(a(i,j));
        }
        p.data[j]=prod;
    }
    return p;
}

prod函式的使用測試

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


    system("pause");
    return 0;
}