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

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

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

1、wilkinson函数：生成Wilkinson特征值测试阵
2、用法说明

 A = wilkinson(4) 函数产生4维Wilkinson特征值测试阵

wilkinson的c++源码实现

生成Wilkinson特征值测试阵
输入为矩阵的维数，输出相应大小的Wilkinson矩阵。

Matrix Matlab2c::wilkinson(int n)
{
	int i,j;
	Matrix p(n,n);
	if (n==1)
		p.data[0]=0;
	else if (n>1)
	{
		if (n%2==0)
		{
			for (i=0;i<p.row;i++)
				for (j=0;j<p.column;j++)
					if (i-j==1 || j-i==1)
						p(i,j) = 1;
					else if(i==j)
						p(i,j) =std::fabs((n-1)/2-i+0.5);
					else
						p(i,j) =0;
		}
		else
		{
			for (i=0;i<p.row;i++)
				for (j=0;j<p.column;j++)
					if (i-j==1 || j-i==1)
						p(i,j) = 1;
					else if(i==j)
						p(i,j) =std::abs(n/2-i);
					else
						p(i,j) =0;
		}
	}
	return p;
}

wilkinson函数的使用测试

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

int main()
{
	Matrix cc=Matlab2c::wilkinson(4);
	cout<<cc.toString()<<endl;

	system("pause");
	return 0;
}