回 帖 发 新 帖 刷新版面

主题:如何将一列数输出为一个矩阵的一列

我现在求一个矩阵的特征值,求10次,每次求完矩阵特征值都有所不同,
eigenvalue=eig(h,w=eigenvector)  
这样可以求得特征值的平均值,我想把每次求的特征值放在另一个矩阵里,按列存放,然后对每行求平均作为这个特征值的平均值。因为特征值是一系列数,可以用reshape函数把它转化为一列数不?

do j=1,4*m*n
     err(j,:)=reshape((/eigenvalue(j)/),(j,1))  
   end do
  write(1,*) eigen(j,nx)
 end do     

回复列表 (共3个回复)

沙发

按列存放
do j=1,4*m*n
   err(:,j)=eigenvalue(:)
end do

还是那句, 兄弟先看看书, 跟着谢谢书上面的代码, 这样比较扎实.

板凳

欲速则不达!
给我的感觉是:楼主问的多,思考的少。自我给其修改如何调用 IMSL 的 Eig 后,我看不出楼主在编程方面有什么进展(不好意思,恕我直言)。一些向量是行还是列不断地纠缠,然后就是矩阵的行列转置。这些基本的东西,自己不会先看看基本的教程书吗?不停地、反复地、接二连三的问,对你帮助不大吧(从你发文我看到了这一点)。
放下手头的工作,花一个把星期好好看看书,跟着书本做做习题,比什么都有效。
祝你早日提高编程水平。

3 楼

哈哈,多谢啦

我来回复

您尚未登录,请登录后再回复。点此登录或注册