using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
namespace ConsoleApp
{
class Program
{
static void Main(strin
如果M是一個密集的m×n矩陣並且v是一個n分量向量，那麼產品u = Mv是由u[i] = sum(M[i,j] * v[j], 1 <= j <= n)給出的m分量向量。一個簡單的實現該乘法是 allocate m-component vector u of zeroes
for i = 1:m
for j = 1:n
u[i] += M[i,j] * v[j]
e
我不明白爲什麼Schur的分解不適用於複雜矩陣。 我的測試程序是： M <- matrix(data=c(2-1i,0+1i,3-1i,0+1i,1+0i,0+1i,1+0i,1+1i,2+0i), nrow=3, ncol=3, byrow=FALSE)
M
S <- Schur(M)
S
(S$Q)%*%(S$T)%*%(solve(S$Q))
結果是： > M
[,1]
我想找到了，如果我能以某種方式找到一種封閉形式的符號矩陣指數矩陣： n=3;
a = symbols(['a'+str(1+k) for k in range(n)], real=True);
Ts = symbols('T_s',real=True,positive=True);
A = Matrix([zeros(1,n),eye(1,n),a])
然而 expm(A)
似乎不工