品質管制
場景：
新品 70% 會通過品檢，未通過者回生產線修正。
修正品 90% 會通過品檢 ，未通過者回生產線再修正 。
再修品 60% 會通過品檢 ，未通過者報廢 。
當前庫存如下：
新品 = 400
修正品 = 60
再修品 = 40
求算可出貨數量
馬可夫轉移狀態圖
馬可夫轉移陣列
myEq <- array(0, c(3, 3)) myEq[1, 1] = 1 myEq[1, 2] = 0 myEq[1, 3] = 0 myEq[2, 1] = -0.1 myEq[2, 2] = 1 myEq[2, 3] = 0 myEq[3, 1] = 0 myEq[3, 2] = -0.3 myEq[3, 3] = 1 N = solve(myEq) myEqVal <- array(0, c(3, 1)) myEqVal[1, 1] = 1 myEqVal[2, 1] = 1 myEqVal[3, 1] = 1 T = N %*% myEqVal myEqVal <- array(0, c(3, 2)) myEqVal[1, 1] = 0.6 myEqVal[1, 2] = 0.4 myEqVal[2, 1] = 0.9 myEqVal[2, 2] = 0 myEqVal[3, 1] = 0.7 myEqVal[3, 2] = 0 B = N %*% myEqVal print(N) print(T) print(B)