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

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

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

1、compan函数：生成多项式系数矩阵的友矩阵
2、用法说明
例如，为了求多项式的x3-7x+6的伴随矩阵，可使用命令： p=[1,0,-7,6];
compan§；

compan的c++源码实现

多项式系数矩阵的友矩阵
输入为多项式的系数向量，输出为该多项式的友矩阵。

Matrix Matlab2c::compan(Matrix& a)
{
	if(a.row!=1)
		throw_logic_error("多项式系数矩阵行数不为1");
	Matrix p(a.column-1,a.column-1);
	if(Matlab2c::FloatNotEqual(a(0),0))
		for (int i=0;i<p.row;i++)
			for (int j=0;j<p.column;j++)
			{
				if(i==0)
					p(j)=(a(j+1)/a(0));
				else if (i-j==1)
					p(i,j) =1;
				else
					p(i,j) =0;
			}
	else
		throw_logic_error("向量第一个元素为0，错误");
	return p;
}

compan函数的使用测试

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

int main()
{
	double a[]={1，0，-7，6};
	Matrix aa=Matrix(1,4,a);
	Matrix bb=Matlab2c::compan(aa);
	cout<<bb.toString()<<endl;

	system("pause");
	return 0;

}