Faddeev_Leverrier法确定矩阵特征方程系数
来源:互联网 发布:网络层的数据帧结构 编辑:程序博客网 时间:2024/05/16 12:28
//Faddeev_Leverrier法确定矩阵特征方程系数
#include <iostream>
using namespace std;
class faddeev
{
private:
int i, j, k, l, n;
double sum1, sum2;
double *alpha;
double **a, **b, **c;
public:
void solution();
double trace();
~faddeev()
{
delete[] alpha;
for (i = 0; i < n; i++)
{
delete[] a[i];
}
delete[] a;
for (i = 0; i < n; i++)
{
delete[] b[i];
}
delete[] b;
for (i = 0; i < n; i++)
{
delete[] c[i];
}
delete[] c;
}
};
void main()
{
faddeev eigen;
eigen.solution();
}
//Faddeev_Leverrier法
void faddeev::solution()
{
cout << "输入矩阵阶数:";
cin >> n;
a = new double*[n];
for (i = 0; i < n; i++)
{
a[i] = new double[n];
}
b = new double*[n];
for (i = 0; i < n; i++)
{
b[i] = new double[n];
}
c = new double*[n];
for (i = 0; i < n; i++)
{
c[i] = new double[n];
}
alpha = new double[n+1];
for (i = 0; i < n; i++)
for (j = 0; j < n; j++)
{
cout << "\n输入a[" << i << "][" << j << "] = ";
cin >> a[i][j];
}
alpha[0] = 1.0;
for (i = 0; i < n; i++)
for (j = 0; j < n; j++)
{
b[i][j] = a[i][j];
}
alpha[1] = trace();
for (k = 2; k <= n; k++)
{
for (i = 0; i < n; i++)
for (j = 0; j < n; j++)
{
if (i == j)
{
c[i][j] = b[i][j] - alpha[k-1];
}
else
{
c[i][j] = b[i][j];
}
}
for (i = 0; i < n; i++)
for (j = 0; j < n; j++)
{
sum2 = 0.0;
for (l = 0; l < n; l++)
{
sum2 += a[i][l]*c[l][j];
}
b[i][j] = sum2;
}
alpha[k] = trace() / k;
}
cout << "\n特征多项式的系数是:" << endl;
for (i = 1; i <= n; i++)
{
cout << "\nalpha[" << i << "] = " << alpha[i] << endl;
}
}
//计算矩阵迹的函数
double faddeev::trace()
{
sum1 = 0.0;
for (i = 0; i < n; i++)
for (j = 0; j < n; j++)
{
if (i == j)
{
sum1 += b[i][j];
}
}
return(sum1);
}
- Faddeev_Leverrier法确定矩阵特征方程系数
- hdu3802 矩阵快速幂&特征方程&降幂
- 特征方程
- 不共线三点确定特征平面并求解平面系数
- 根据国防科大论文确定的特征点坐标拟合平面方程-拟合优度分析
- matlab之自定义方程系数解方程
- 读书笔记: 变系数波方程
- 最小二乘法求解直线方程系数
- 两点确定直线方程
- 递归方程的特征方程
- 高斯消元法求矩阵系数
- 数列的特征方程
- 特征根方程
- 特征根方程
- 数列特征方程
- 特征方程及其应用
- l1范数最小化求解系数方程_正交匹配追踪法(orthogonal matching pursuit)
- 使用牛顿法确定逻辑斯谛回归(Logistic Regression)最佳回归系数
- Struts2.x+JFreeChart搭建框架出错(一)
- Linux/UNIX脚本编程学习笔记(0728)
- __stdcall, __cdecl, __fastcall, __declspec几种调用约定及区别
- ttl值的含义
- oracle初始参数文件问题,求解,谢谢
- Faddeev_Leverrier法确定矩阵特征方程系数
- Householder变换将对称矩阵化为三对角形式
- poj 2762
- Jacobi法确定对称矩阵特征值与特征向量
- 乘幂法计算矩阵绝对值最大特征值
- Google Hacking
- POJ 3485 Highway(贪心算法)
- VC使用API连接mysql方法
- poj2299 - Ultra-QuickSort (求逆序数)