版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (1324)
>虫友互识 (129)
>导师招生 (44)
>考研 (25)
>考博 (23)
>硕博家园 (22)
>论文投稿 (22)
>休闲灌水 (22)
>教师之家 (21)
>基金申请 (17)
>博后之家 (12)
>论文道贺祈福 (12)
>公派出国 (12)
>找工作 (10)
>文献求助 (10)
>电化学 (8)

返回列表

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";
}

回复此楼

» 猜你喜欢

1楼 2011-10-03 21:38:07

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

高校试用期调动怎么一个手续啊？急已经有13人回复
求助有机牛人，有什么办法可以强有力的促进还原消除啊~~ 已经有13人回复
African journal of Biotechnology怎么8月份刊出的文章还没有进SCI数据库啊已经有10人回复
最全的matlab字符串操作总结--一点点的改变让程序变得更快速已经有48人回复
正交试验的最优结果怎么选取已经有5人回复
论文小修，没改不小心就发过去了，怎么办啊已经有27人回复
100BB急求matlab编程中的问题解决，如何修改这个程序呢？已经有4人回复
文章已经接受，但是Elsevier数据库编辑要我改变表的格式，这个怎么办？已经有11人回复
【求助】液相色谱测定血浆问题，怎么做？已经有3人回复
【求助】module 里面封装的程序的形参能是指针类型的吗？【已完结】已经有5人回复
【求助】求将c语言改写为C++，求教高人啊！已经有6人回复
【求助】做程序升温氧化，出峰时抖动需要怎么改进已经有17人回复
【求助】指针型真空泵与固相萃取装置连用时怎么操作已经有6人回复
【求助】热固性塑料成型添加玻纤或者碳纤改性怎么加啊？已经有12人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 dengmin155 的主题更新