[em1][em1][em1][em2][em2][em2]
在做矩阵相乘，遇到几个问题：一直不知道该怎么解决
eigenvector：矩阵特征向量矩阵，大小为（4mn,4mn），m=4=n,由外部量给出
traneigv：eigenvector的逆矩阵
temp:控制前4行（列）
j：控制后面的行（列）
hx,hy:大小为（4mn,4mn）的矩阵
       do temp=1,4
    do j=5,4*m*n

       ea(temp,:)=traneigv(temp,:).x.hx  ！ea是4*4mn矩阵
        ec(temp,:)=traneigv(temp,:).x.hy  ！ec是4*4mn矩阵 
        eb(j,:)=traneigv(j,:).x.hx       ！eb：(4mn-4)*4mn
        ed(j,:)=traneigv(j,:).x.hy       !ed:(4mn-4,4mn)
        ee(:,temp)=eigenvector(:,temp)   !ee:前4个特征向量
        eg(:,j)=eigenvector(:,j)         ！后面特征向量
            
    conduc(temp)=conduc(temp)+(ea(temp,:).x.eg(:,j)).x.((ed(j,:).x.ee(:,temp)))/(eigenvalue(temp)-eigenvalue(j))**2    
          
     end do  
       end do
      con=sum(conduc)

我的目的乘法循环过程：
比如：
     ea(temp,:).x.eg(:,j)：    当temp=1时，ea(1,:)是逆矩阵第一行与hx矩阵相乘得到的结果，是1行4mn列，然后让它跟特征向量矩阵eg(:,j)相乘，得到的结果是1行j列的一维数组，j每改变一个数，就在这行上增加一个数，直到j=4mn，这样一行就有(1，4mn-4)个数，所以最后得到的ea是(4,4mn)的矩阵。 
    ed(j,:).x.ee(:,temp)：还是temp=1,由于ed是特征向量逆矩阵第4行之后的矩阵与hy相乘得到的结果，大小为(4mn-4,4mn)，让它与第1 (temp)个特征向量相乘，应该得到一个(4mn-4,1)的向量，一列向量。
    然后让上面得到的一行矩阵和一列矩阵相乘，得到一个数，我把它存在conduc(temp)里，用temp表示是第几个temp得到的值，以上面temp=1,4为例，conduc里应该有4个数，每个temp对应自己的conduc.
   然后，将conduc求和,所有元素求和，放到con里。

这样描述明白吗？

按照这种思路，我做了下，发现出了很多问题，后来我修改了下：
     do temp=1,i
       do j=i+1,4*m*n         
    ea(temp,:)=traneigv(temp,:).x.hx   
    ec(temp,:)=traneigv(temp,:).x.hy   
    eb(j,:)=traneigv(j,:).x.hx         
    ed(j,:)=traneigv(j,:).x.hy         
    ee(:,temp)=eigenvector(:,temp)     
    eg(:,j)=eigenvector(:,j)           
    
    conduc=conduc+((ea.x.eg).x.((ed.x.ee))-conjg((ea.x.eigenvector).x.(ed.x.eigenvector)))/(eigenvalue(temp)-eigenvalue(j))**2
        end do      
       end do  
    con=sum(conduc)
虽然程序正常运行了，但是数据明显不合理啊。

程序肯定有问题，高手们帮忙支招吧。


[em8][em8][em8][em8]