Gram
Purpose
Generalized rank anihilation method.
Synopsis
- [ord1,ord2,ssq,aeigs,beigs] = gram(a,b,tol,scl1,scl2,out)
Description
GRAM determines the joint invariant subspaces common to the two input matrices a and b, the ratio of their magnitudes ssq, and the response in both modes/orders ord1 and ord2. GRAM assumes that the input matrices a and b are bilinear, i.e. are the summation over outer products.
Inputs are the two response matrices a and b, and the number of factors to calculate or tolerance on the ratio of smallest to largest singular value tol. Optional inputs scl1 and scl2 are scales to plot against when producing plots of the reponse in each mode/order. Optional input out suppresses plotting and printing of results to the command window when set to 0 {default out = 1}.
Outputs are the pure component responses in each mode ord1 and ord2, the table of eigenvalues and their ratios ssq, and the eigenvalues for each matrix aeigs and beigs.