back start next


[start] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] [167] [168] [169] [170] [171] [172] [173] [174] [175] [176] [177] [178] [179] [180] [181] [182] [183] [184] [185] [186] [187] [188] [ 189 ] [190] [191] [192] [193] [194] [195] [196] [197] [198] [199] [200] [201] [202] [203] [204] [205]


189

Matrix Solution to Linear Equations and Maikov Chains

DIRECT SOLUTION AND CONVERGENCE METHOD

Before computer programs offered readj solutions, problems such as Markov chains were solved in a direct manner, by algebraically manipulating the equations. This direct solution requires an understanding of simple mafrix arithmetic, and very careful attention to calculating the numbers correctly. The convergence method is now easier although it requires many more calculations. Without the computer, we would never even consider using this approach-with a computer, it is the best choice.

GENERAL MATRIX FORM

A mafrix is a rectangular arrangement of elements into rows and columns A matrix A is said to be m x n (pronounced "m by n") if there are m rows and n columns in A

Certain properties of a maUix make it a valuable tool for solving simultaneous linear equations. These elementarj mafrix operations, called transformations, allow you to alter the rows (which will represent equations) without changing the solution. There are three basic row operations:

1. Multiplication or division of all elements of the row by any number.

2. Interdianging of any two rows (and consequently of all rows)

3. The addition or si4)traction of the elements of one row with the corresponding elements of another.

To relate the maUix to simultaneous linear equations, consider a three-equation example, where only the coefficients, a, remain as unknowns:

The ihm. tkiucnijry n ii. ur



.JUL..



tDVš«hcr Pniblbilitie F.rxiiM « ,

. -ill tqujj (I . l)di 1 1 . tbt cgmiiaiu

(Hii!)

Computer Program Direct Solution

The following FORTRAN program accepts a sjstem of 10 equations in 10 variables (unknowns) and applies the maUix method of solution. This is done using the method of Gaussian elimination, discussed in the previous section.



[start] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] [167] [168] [169] [170] [171] [172] [173] [174] [175] [176] [177] [178] [179] [180] [181] [182] [183] [184] [185] [186] [187] [188] [ 189 ] [190] [191] [192] [193] [194] [195] [196] [197] [198] [199] [200] [201] [202] [203] [204] [205]