在 Python 中实现矩形矩阵的伪谱

问题描述

我想在 Python 中实现以下算法，该算法计算矩形（非方形）矩阵的伪谱。 T 矩阵是舒尔分解的上三角。Band_Qr 是第二个 qr 分解。网格参数必须是复数序列。

有什么帮助吗？到目前为止，我的努力是这样的：

import numpy as np
import matplotlib.pyplot as plt
from scipy.linalg import svdvals,schur


A = np.array([[0,2,4],[0,1,[1,4]])
A
def RECTSPA(A,grid):
    m,n = A.shape
    I = np.eye(n)
    if m >= 2*n:
        a = A[n+1:m,:]
        Q,R = np.linalg.qr(a)
        T,U = schur(a,output='real')
        R[0:n]
        T = I 
    else:
        S2,T2,Q,Z = linalg.qz(A[m-n+1,:],I[m-n+1,:])
    for z in range(1,grid):
        q,r = np.linalg.qr(z*T-S2)
        psevals = np.min(svdvals(r))
    return(psevals)

但它不起作用。有什么帮助吗？

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

decomposition matrix matrix python