主题:如何将一列数输出为一个矩阵的一列
我现在求一个矩阵的特征值,求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
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