麻烦好心人帮我求下最大特征值和权向量，谢谢啦！！

如题所述

第1个回答 2012-03-05

VC++6.0下运行的。

#include<iostream>
#include<math.h>
using namespace std;
const double m=1e-5;
const int Max=11; //定义最大阶
const double RI[10]={};
class A
{
private:
unsigned x,y,n;
double B[Max][Max];
double v[Max],w[Max],u[Max],f[Max],h[Max];
public:
void creat();
void strum();
void sum();
void open();
};
void A::creat()
{
cout<<"请输入矩阵的阶数: ";
cin>>n;
x=n;
y=n;
for(unsigned i=0;i<x;i++)
{
cout<<"输入第"<<i+1<<"行 ";
for(unsigned j=0;j<y;j++)
cin>>B[i][j];
}
}
void A::strum()
{
double s,e=0,t=0,b,r;
unsigned i,j,c;
for(c=0;c<n-1;c++)
u[c]=0;
u[n-1]=1;
do{
b=0;
r=0;
for(i=0;i<=n-1;i++)
{
v[i]=0;
w[i]=u[i];
for(j=0;j<=n-1;j++)
v[i]+=B[i][j]*u[j];
if(fabs(v[i])>b)
b=fabs(v[i]);
}
for(i=0;i<=n-1;i++)
{
u[i]=v[i]/b;
if(fabs(u[i]-w[i])>r)
r=fabs(u[i]-w[i]);
}

}
while(r>=m);
cout<<"主特征根 T="<<b<<"\n";

for(i=0;i<=x-1;i++)
{ cout<<"特征向量"<<i+1<<"为: ";
cout<<u[i]<<" ";
cout<<"\n";}
s=(b-n)/(n-1);
cout<<"\n"<<"CI="<<s<<"\n";
cout<<"CR="<<s/RI[n-1]<<"\n";
for(i=0;i<=n-1;i++)
{for(j=0;j<=n-1;j++)
t+=B[i][j] ;}
for(i=0;i<=n-1;i++)
{ h[i]=0;
f[i]=0;
for(j=0;j<=n-1;j++)
h[i]+=B[i][j];
f[i]=h[i]/t;
e+=f[i];}
for(i=0;i<=n-1;i++)
cout<<"\n"<<"权重向量"<<i+1<<"为: "<<f[i]/e;
cout<<"\n";
}
void A::sum()
{ unsigned i,j;
cout<<"\n请输入所要乘的一维数组: ";
for(i=0;i<y;i++)
cin>>v[i];
for(i=0;i<x;i++)
{
u[i]=0;
for(j=0;j<y;j++)
u[i]+=B[i][j]*v[j];
cout<<"\n总排序向量"<<i+1<<"为: "<<u[i];
}
cout<<"\n";

}
void A::open()
{
cout<<"\n输入矩阵的行数: ";
cin>>x;
cout<<"输入矩阵的列数: ";
cin>>y;
for(unsigned i=0;i<x;i++)
{
cout<<"输入第"<<i+1<<"行:";
for(unsigned j=0;j<y;j++)
cin>>B[i][j];
}
}
int main()
{
int k;
A a1;
cout<<"求主特征根等,请按1回车 ";
cout<<"\n求总排序向量,请按2回车 ";
cin>>k;
if(k==1)
{a1.creat();
a1.strum();}
else if(k==2)
{a1.open();
a1.sum();}
else
cout<<"\n\n\n\n抱歉,您的输入有误!\n";
return 0;
}追问

辛苦你了朋友，看来这个矩阵不对，那我的改改！！谢谢啦

第2个回答 2012-03-05

相似回答

麻烦大家,最大特征值和权向量,谢谢 1 1 1 2 3 1 1 1 2 3 1 1 1 2...答：最大特征值是5.1508；对应的权向量是【0.5451； 0.5451 ； 0.5451 ； 0.2259 ； 0.2402】（这个是列向量，写成竖着的形式就OK了）；

{求权向量及最大特征值?}这是我们数学建模老师让我们做的作业,要交的...答：aij=1/aji来表示,由于aij的特点,A称为正互反矩阵,SATTY等人已经证明可以将A的最大特征值对应的特征向量当作因素的权向量,并且给出了一致性检验标准,当CR<0.1时便合服一致性检验。仿照SAATY等人的方法我们也可以来确定满载率、步行距离、公交驾车之比、等车时间之间的可以设C1=满载率,C2=步行距离,C3=公交驾车之...

...最大特征根λmax,CI,RI,CR? 解答请详细一些,谢谢!答：if CR<0.1 %检验矩阵Q的一致性 disp('成对比矩阵Q的一致性可以接受!其特征向量可以作为权向量')；disp('CI=')；disp(CI)；disp('CR=')；disp(CR)；disp('由成对比矩阵Q计算出的权向量w为：')；disp(W)；else disp('成对比矩阵Q不是一致阵，请重新评分')return end ...

怎么用matlab6.5来求矩阵的最大特征值和特征向量以及如何进行一致性检验...答：V是对应特征值的特征向量 归一化处理得到的权向量 >> A=[1 3/2 6 2; 3/2 1 4 4/3;1 1 1 1/3;1/3 3/4 3 1];>> for j=1:4 sum(j)=0;for i=1:4 sum(j)=A(i,j)+sum(j);%将第j列所有元素相加 end end for j=1:4 for i=1:4 A(i,j)=...

大家正在搜

特征向量和特征值特征值与特征向量例题最大特征值如何求矩阵的特征值求特征值的方法怎么感谢好心人什么是特征值特征值性质特征值为0代表什么