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

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

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

1、poly函数：

根据根求多项式系数 或求矩阵特征多项式

2、用法说明

poly(A)
当A是一个NN矩阵式，poly(A)命令求出A的特征多项式det(lambdaeye(size(A))-A)
当V是向量时，命令poly(A)生成以V为根的多项式

poly的c++源码实现

根据根求多项式系数 或求矩阵特征多项式（还未实现）
输入：多项式的根向量a
输出：多项式的系数矩阵

Matrix Matlab2c::poly(Matrix& a)
{
	
	int i;
	Matrix p(1,1),b(1,2);
	p.data[0]=1;
	b.data[0] =1;
	if (a.isVector())
	{
		for (i=0;i<a.row*a.column;i++)
		{
			b.data[1] = -(a.data[i]);
			p=conv(p,b);
		}
	}else if (a.row==a.column)   //求方阵的特征多项式
	{
		//p=det(lambda*eye(a.row))-a;，还未实现
	}else
	{
		throw_logic_error("poly函数参数需为方阵或向量");
	}
	
	return p;
}

poly函数的使用测试

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


int main()
{
	double a[]={3};
	Matrix aa=Matrix(1,1,a);
	Matrix cc=Matlab2c::poly(aa);
	cout<<cc.toString()<<endl;

	system("pause");
	return 0;
}