1樓:
我是這樣解決的
用1到(X1+1)(Y1+1)(Z1+1)的乙個數t=H*(Y1+1)(Z1+1)+m*(Z1+1)+n +1代表乙個狀態(H,m,n)
然後求這樣乙個(X1+1)(Y1+1)(Z1+1)^2的矩陣.
Module[,
prob[i_, j_] := Module[,
H1 = Quotient[i - 1, d2*d3];
m1 = Quotient[Mod[i - 1, d2*d3], d2];
n1 = Mod[Mod[i - 1, d2*d3], d2];
H2 = Quotient[j - 1, d2*d3];
m2 = Quotient[Mod[j - 1, d2*d3], d2];
n2 = Mod[Mod[j - 1, d2*d3], d2];
Return[
If[H1 == 0, If[H2 == 0 && n2 == 0 && m2 == 0, 1, 0If[H1 - H2 == 1 && m1 == m2 && n1 == n2If[m1 == 0 || n1 == 0, If[m1 == 0 && n1 == 0, 1, 1/2], 1/3If[m1 - m2 == 1 && H1 == H2 && n1 == n2, If[n1 == 0, 1/2, 1/3If[n1 - n2 == 1 && H1 == H2 && m1 == m2If[m1 == 0, 1/2, 1/3], 0transfermatrix = Array[prob, ];
Mathematica 中如何定義 f x y f x f y ?
Attributes f HoldAll f a Plus f Unevaluated a 寫完了才想起來上面這種複雜的寫法其實毫無必要。最簡單的就足夠了 ClearAll f Attributes f HoldAll f x y f x f y f a b c f a f b f c 因為Plus...
mathematica中畫圖函式的區分?
吉聆遠 具體的區別 withparadox2 已經說得很明晰了。但對於 3,當不好求出顯示解時,如 x 2 Sin y Cos x y,不好寫出以上兩種形式時便可以用ContourPlot ImplicitPlot淘汰掉 來繪製,例如 ContourPlot x 2 Sin y Cos x y,我來...
mathematica中怎樣進行矩陣初等行變換
紐澤西管風琴 有三種初等行變換 將矩陣的第行和第行交換 將矩陣的第行乘以乙個係數 將矩陣的第行加上第行的倍。分別對應三個函式 swap m i j Permute m Cycles times m i k MapAt k m i add m i j k MapAt km j m,i 作為例子,定義函...