两个同阶矩阵的Bialternate Product

mxlzhenzhu

两个同阶矩阵的Bialternate Product，定义及其性质如下：

下面的代码可以用来计算两个矩阵的Bialternate Product， 拿走不谢。

1. function AB=Bialternate_Prod(A,B)
   
2. %% Compute the bialternate_Product of matrix A & B
   
3. %% Format:AB=Bialternate_Prod(A,B)
   
4. %% mxl.2019-1-12
   
5. n=size(A,1);
   
6. if n<=2
   
7.     error('Input Error.');
   
8. end
   
9. 
   
10. if size(A,1)~=size(A,2)||size(A,1)~=size(B,1)||size(A,1)~=size(B,2)
    
11.     error('Input Error.');
    
12. end
    
13. 
    
14. 
    
15. 
    
16. 
    
17. 
    
18. [v1,v2]=find(tril(true(n,n),-1));
    
19. [v1,ki]=sort(v1,'ascend');
    
20. v2=v2(ki);
    
21. 
    
22. m=n*(n-1)/2;
    
23. AB=zeros(m,m);
    
24. for loopi=1:m
    
25.     p=v1(loopi);q=v2(loopi);
    
26.     for loopj=1:m
    
27.         r=v1(loopj);s=v2(loopj);
    
28.         
    
29.         M=[A(p,r) A(p,s)
    
30.             B(q,r) B(q,s)];
    
31.         N=[B(p,r) B(p,s)
    
32.             A(q,r) A(q,s)];
    
33.         
    
34.         deta=det(M);detb=det(N);
    
35.         
    
36.         AB(loopi,loopj)=(deta+detb)/2;
    
37.         
    
38.     end
    
39. end
    
40. 
    
41. 
    
42. end
    

复制代码