24小时热门版块排行榜

返回列表

dengmin155

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2.3
帖子: 13
在线: 5.1小时
虫号: 1113369
注册: 2010-10-03
专业: 光学

[求助] 这个程序怎么改啊指针

这个程序是用来进行矩阵对角化的，但是，编译好后，一旦输入有零的矩阵元后，就会发生指令使得内存不可读的错误，如果矩阵中没有零元素，就很好用，这是采用雅可比方法编写的~求大侠相助！！

#include
#include
#include
#include
class jacobi
{private:
int flag,i,iteration,j,k,n,b,h;
double csa,eps1,eps2,eps3,p,q,s1,s2,sna,spq,sum,temp;
double  *eigen,**a,**t;
public:
jacobi()
{
iteration=0;flag=1;
eps1=1e-10;eps2=1e-10;eps3=1e-10;

}
void solution ();
~jacobi()
{
   delete [] eigen;
   for (i=0;i; }
   delete [] a;
   for (i=0;i; }
delete [] t;
}
};

void main()
{ jacobi eigenpair;
eigenpair.solution();

}

void jacobi::solution()
{ cout<<"输入矩阵阶数：";
cin>>n;
eigen=new double [n];
for (i=0;i=0.0; }
a=new double*[n];
for(i=0;i=new double [n];}
t=new double*[n];
for (i=0;i=new double [n];}
for (i=0;i for (j=0;j {t[j]=0.0;a[j]=0.0; }
for (i=0;i for (j=i;j {
cout<<"\n输入a["< cin>>a[j];
a[j]=a[j];
}
cout<<"输入矩阵结束"< s1=0.0;sum=0.0;
for (i=0;i {
t=1.0;
s1+=a*a;
sum+=a;

}
do
{
iteration++;
for(i=0;i for(j=(i+1);j {
if(fabs(a[j])[j]=0.0; a[j]=0;}
else
{
q=fabs(a-a[j][j]);
if (q>eps2)
{ b=i;h=j;
p=2*a[j]*q/(a-a[j][j]);
spq=sqrt(p*p+q*q);
csa=sqrt((1+q/spq)/2);
sna=p/(2*csa*spq);
}
else
{
csa=0.707107;sna=csa;
}
for(k=0;k { temp=t[k];
t[k]=temp*csa+t[k][j]*sna;
t[k][j]=temp*sna-t[k][j]*csa;

}

         temp=a;
   a=temp*csa*csa+a[h][h]*sna*sna+2*a[h]*csa*sna;
         a[h][h]=temp*sna*sna+a[h][h]*csa*csa-2*a[h]*csa*sna;
         a[h]=0;a[h]=0;
   for(k=0;k    if((k!=b)&&(k!=h))
   {
   temp=a[k];
      a[k]=temp*csa+a[h][k]*sna;
      a[h][k]=temp*sna-1*a[h][k]*csa;//这里有点问题

   }

for (k=0;k if((k!=b)&&(k!=h))
{
temp=a[k];
a[k]=temp*csa+a[k][h]*sna;//有问题，这里要注意a[h]的索引标的问题
a[k][h]=temp*sna-1*a[k][h]*csa;
}


}

}
s2=0.0;
for (i=0;i {
eigen=a;
s2+=eigen*eigen;
}
if (fabs(1-s1/s2) else { s1=s2;}

}while(flag==1);

for (i=0;i {
cout<<"\n特征值["<< }
cout<<"\n收敛于"< cout<<"\n输入矩阵的迹="< sum=0.0;
for(i=0;i; }
cout<<"\n特征值之和="< cout<<"\n特征向量是:"< cout.precision (4);
for(j=0;j {
   cout<<"\n特征向量["<    for (i=0;i[j]<
   cout< }
system("pause";
}