matlab 定义细胞矩阵,【图片】mathematica中有像matlab中定义细胞，每个元都是矩阵的方法吗？_mathematica吧_百度贴吧...

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

Input["n=", n];

n = %

Input["q=", q];

q = %

B = IdentityMatrix[n + 1]

X = Array[x, n + 1];

X[[1]] = 35/1000

For[j = 2, j <= n + 1, j++, X[[j]] = X[[1]] + j*10^-5]

f[r_, q_] := If[q == 1, r, (1 - q^r)/(1 - q)]

h[r_, q_] := Product[f[i, q], {i, 1, r}]

g[r_, q_] := If[r == 0, 1, h[r, q]]

k[n_, i_, q_] :=

If[q == 1, Binomial[n, i], h[n, q]/(h[i, q] h[n - i, q])];

z[q_, x_, n_, i_] :=

k[n, i, q]*x^i*Product[(1 - q^s*x), {s, 0, n - i - 1}]

For[j = 1, j <= n + 1, j++, B[[1]][[j]] = -z[q, X[[2]], n, j - 1]]

For[i = 2, i <= n + 1, i++,

For[j = 1, j <= n + 1, j++,

B[[i]][[j]] =

z[q, X[[i]], n,

j - 1] - (z[q, X[[2]], n, j - 1] +

z[q, X[[1]], n, j - 1])*(z[q, X[[i]], n, 0]/

z[q, X[[2]], n, 0])]]

NumberForm[MatrixForm[B, 16], 16]

%此处位第一个矩阵

Transpose[

prod[n_, r_, q_] :=

Block[{mat, $MaxPrecision = 16, $MinPrecision = 16},

mat[idx_] :=

SparseArray[{Band[{1, 1}] -> 1,

Band[{1, 2}] ->

Table[r^(idx - i) - q^(n - idx), {i, idx, 1, -1}]}, n + 1];

Array[mat, n - 1, 1, Dot]]]

%此处位第二个矩阵

Input["n=", n];

n = %

M = IdentityMatrix[n + 1]

For[i = 1, i <= n + 1, i++, M[[i]][[i]] = Binomial[n, i - 1]]

NumberForm[MatrixForm[M, 16], 16]

此处为第三个矩阵

Input["n=", n];

n = %

Input["q=", q];

q = %

Dq = IdentityMatrix[n + 1]

f[r_, q_] := If[q == 1, r, (1 - q^r)/(1 - q)]

h[r_, q_] := Product[f[i, q], {i, 1, r}]

g[r_, q_] := If[r == 0, 1, h[r, q]]

k[n_, i_, q_] :=

If[q == 1, Binomial[n, i], h[n, q]/(h[i, q] h[n - i, q])];

For[i = 1, i <= n + 1, i++, Dq[[i]][[i]] = k[n, i, q]] NumberForm[

MatrixForm[Dq, 16], 16]

%此处为第四个矩阵