Group Abstract Group Abstract

Message Boards Message Boards

Solve a System of equations using matrix inverse

Posted 5 years ago

Hello, I am attempting to solve a system of equations by using the matrix inverse. I was able to create a system and separate the system into two matrices, but I am having issues generating the matrix inverse. I have attached the nb file containing my code.

Attachments:
6 Replies

display a linear equation as a matrix other than the way I used it below?

What exactly do you mean?

POSTED BY: Hans Dolhaine

You wrote one equation twice in system1. And C is a Mma Symbol. This runs without error

Clear[x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, \
x15, x16, x17, x18, x19, x20, x21, x22, x23]
equation1 = 
  634.14 x1 + 97.54 x10 - 216.48 x11 - 486.49 x12 + 12.46 x13 + 
    649.97 x14 + 877.29 x15 - 755.88 x16 + 365.27 x17 - 442.49 x18 + 
    736.68 x19 + 268.69 x2 - 387.29 x20 + 742.82 x21 - 348.74 x22 + 
    98.25 x23 - 939.96 x3 - 313.98 x4 + 672.09 x5 + 348.98 x6 + 
    754.7 x7 - 535.21 x8 - 72.71 x9 == 102.79;
equation2 = 
  57.49 x1 + 336.13 x10 - 306.59 x11 - 151.18 x12 + 393.28 x13 - 
    229.18 x14 + 21.16 x15 - 878.84 x16 - 733.05 x17 + 567.96 x18 - 
    596.09 x19 - 343.58 x2 + 553.37 x20 - 741.77 x21 + 550.17 x22 + 
    269.73 x23 + 854.11 x3 + 84.37 x4 + 995.49 x5 - 521.77 x6 + 
    477.37 x7 - 865.91 x8 + 908.69 x9 == -762.13;
equation3 = 
  768.64 x1 - 866.74 x10 - 322.2 x11 - 315.52 x12 + 627.03 x13 + 
    392.89 x14 - 477.87 x15 + 936.09 x16 - 379.37 x17 - 681.48 x18 + 
    997.81 x19 - 55.4 x2 - 160.5 x20 + 737.81 x21 - 886.85 x22 - 
    291.88 x23 - 987.62 x3 + 324.21 x4 - 770.56 x5 - 878.62 x6 + 
    880.97 x7 + 832.53 x8 + 916.11 x9 == 767.06;
equation4 = 
  852.33 x1 + 93.77 x10 - 308.41 x11 - 939.14 x12 - 618.9 x13 + 
    573.08 x14 + 61. x15 - 462.59 x16 - 241.62 x17 - 675.85 x18 - 
    38.36 x19 - 990.88 x2 - 783.38 x20 - 296.37 x21 + 921.63 x22 - 
    649.89 x23 - 435.02 x3 + 288.33 x4 + 576.8 x5 + 511.03 x6 - 
    77.95 x7 - 564.5 x8 - 298.31 x9 == 112.21;
equation5 = 
  646.65 x1 - 840.26 x10 - 106.99 x11 - 445.78 x12 + 505.43 x13 + 
    191.41 x14 - 912.5 x15 - 520.01 x16 - 23.7 x17 + 903.8 x18 - 
    631.25 x19 - 525.66 x2 + 530.8 x20 + 431.07 x21 + 66.8 x22 + 
    303.66 x23 - 336.69 x3 - 69.38 x4 - 15.55 x5 - 641.89 x6 - 
    57.87 x7 + 696.96 x8 + 409.45 x9 == 214.6;
equation6 = -552.61 x1 + 912.29 x10 - 975.39 x11 - 870.38 x12 + 
    921.44 x13 - 698.54 x14 + 862.28 x15 + 329.94 x16 + 988.37 x17 + 
    374.54 x18 + 637.58 x19 - 441.51 x2 - 42.11 x20 - 138.73 x21 + 
    690.18 x22 + 917.89 x23 - 142.99 x3 + 343.93 x4 - 58.83 x5 + 
    965.62 x6 + 945.54 x7 + 295.82 x8 - 428.23 x9 == 788.16;
equation7 = 
  38.29 x1 + 403.04 x10 - 690.25 x11 - 235.05 x12 + 636.74 x13 - 
    146.01 x14 - 945.1 x15 + 178.11 x16 + 838.91 x17 + 921.68 x18 - 
    676.65 x19 + 183.46 x2 - 820.45 x20 + 414.96 x21 - 406.99 x22 - 
    222.9 x23 + 910.97 x3 - 374.53 x4 - 899.46 x5 - 774.27 x6 + 
    995.92 x7 + 551.71 x8 + 56.73 x9 == -514.84;
equation8 = 
 190.75 x1 + 238.13 x10 - 984.05 x11 - 531.96 x12 - 919.46 x13 + 
   412.37 x14 - 379.33 x15 + 24.95 x16 + 823.48 x17 - 811.69 x18 - 
   228.43 x19 + 152.52 x2 + 468.03 x20 - 472.1 x21 + 466.73 x22 - 
   722.04 x23 - 556.3 x3 + 232.74 x4 + 231.01 x5 + 422.89 x6 + 
   158.29 x7 + 534. x8 + 542.32 x9 == -233.44; equation9 = 
 870.91 x1 + 556. x10 - 366.9 x11 - 144.72 x12 - 472.31 x13 + 
   30.68 x14 - 969.56 x15 - 403.76 x16 + 476.46 x17 - 869.8 x18 - 
   458.33 x19 - 705.33 x2 - 157.96 x20 + 212.32 x21 - 612.48 x22 - 
   739.55 x23 - 892.3 x3 - 138.97 x4 + 201.55 x5 - 293.1 x6 - 
   189.66 x7 - 240.01 x8 - 
   346.22 x9 == -963.15; equation10 = -778.6 x1 + 385.63 x10 - 
   900.46 x11 + 28.64 x12 + 270.83 x13 + 434.3 x14 - 314.46 x15 - 
   669.93 x16 + 122.05 x17 + 7.24 x18 + 615.19 x19 - 946.84 x2 + 
   433.12 x20 + 779.5 x21 - 849.43 x22 - 334.84 x23 + 580. x3 - 
   493.3 x4 - 526.83 x5 - 665.79 x6 - 352.95 x7 + 219.66 x8 - 
   560.62 x9 == -95.48
equation11 = 
 236.42 x1 - 15.5 x10 + 525.42 x11 - 177.3 x12 + 899.17 x13 + 
   622.3 x14 - 238.89 x15 + 292.48 x16 + 561.03 x17 + 29.5 x18 - 
   406.12 x19 + 343.69 x2 - 703.95 x20 - 443. x21 - 142.38 x22 - 
   208.98 x23 + 54.93 x3 + 206.79 x4 - 787.09 x5 - 151.16 x6 - 
   986.95 x7 + 273.85 x8 + 522.39 x9 == 155.73; equation12 = 
 680.18 x1 + 708.45 x10 - 664.91 x11 - 946.99 x12 - 392.22 x13 + 
   92.91 x14 - 63.53 x15 + 703.83 x16 - 320.08 x17 - 268.2 x18 + 
   422.39 x19 - 419.69 x2 - 420.1 x20 + 172.18 x21 + 328.86 x22 - 
   870.12 x23 - 982.71 x3 - 480.25 x4 - 825.52 x5 - 248.07 x6 + 
   652.84 x7 - 942.37 x8 - 244.09 x9 == 802.85;
equation13 = 
  50.63 x1 - 486.19 x10 - 346.43 x11 + 194.1 x12 + 791.08 x13 - 
    356.58 x14 + 566.14 x15 + 952.08 x16 - 602.24 x17 - 695.1 x18 - 
    928.78 x19 + 257.52 x2 - 542.2 x20 + 787.83 x21 + 519.67 x22 - 
    239.3 x23 + 253.41 x3 + 858.31 x4 + 775.11 x5 - 33.93 x6 + 
    949.52 x7 - 501.49 x8 - 932.67 x9 == 756.75;
equation14 =
 -282.05 x1 + 422.61 x10 + 768.42 x11 + 221.72 x12 - 283. x13 + 
   951.33 x14 + 5.37 x15 + 786.37 x16 - 20.16 x17 + 591.71 x18 + 
   821.73 x19 + 988.3 x2 + 720.57 x20 - 636.77 x21 + 583.91 x22 + 
   190.09 x23 - 70. x3 - 582.82 x4 - 584.18 x5 - 927.73 x6 - 
   66.66 x7 + 334.4 x8 - 314.91 x9 == -777.53; equation15 = 
 347.84 x1 + 342.35 x10 + 672.1 x11 + 403.15 x12 + 961.69 x13 + 
   428.25 x14 + 183.05 x15 + 855.27 x16 - 835.89 x17 - 169.98 x18 - 
   492.42 x19 + 108.05 x2 + 663.31 x20 + 26.26 x21 + 916.49 x22 - 
   176.56 x23 - 772.21 x3 + 188.46 x4 - 883.49 x5 + 244.16 x6 - 
   880.37 x7 + 0.84 x8 + 233.13 x9 == -275.74;
equation16 = -689. x1 + 44.41 x10 - 94.25 x11 - 712.75 x12 - 
   365.06 x13 - 80.49 x14 - 161.44 x15 + 242.22 x16 + 90.75 x17 - 
   986.12 x18 - 819.92 x19 - 296.6 x2 - 58.88 x20 - 372.61 x21 + 
   272.23 x22 + 66.93 x23 + 360.46 x3 - 976.77 x4 - 776.85 x5 + 
   897.3 x6 - 47.97 x7 - 234.66 x8 - 711.09 x9 == 786.52; equation17 =
  983.17 x1 + 863.82 x10 + 371.05 x11 - 719.86 x12 - 415.52 x13 - 
   589.89 x14 + 12.68 x15 - 88.53 x16 - 980.51 x17 - 178.96 x18 - 
   937.55 x19 - 257.57 x2 + 347.72 x20 + 212.37 x21 + 33.56 x22 + 
   851.85 x23 - 9.23 x3 - 634.54 x4 + 573.32 x5 + 281.41 x6 + 
   915.45 x7 - 140.7 x8 - 
   540.64 x9 == -68.35; equation18 = -518.67 x1 - 934.33 x10 + 
   399.3 x11 - 185.26 x12 + 740.53 x13 + 213.06 x14 + 158.16 x15 + 
   659.67 x16 + 216.62 x17 - 993.05 x18 + 89.73 x19 - 619.61 x2 + 
   603.47 x20 + 812.34 x21 + 488.95 x22 - 461.56 x23 - 587.24 x3 + 
   399.91 x4 - 837.61 x5 + 964.37 x6 - 255. x7 + 628.13 x8 - 
   231.17 x9 == 758.41;
equation19 = 
  333.16 x1 + 779.73 x10 - 559.82 x11 - 781.88 x12 - 465.08 x13 - 
    367.24 x14 - 742.58 x15 + 432.99 x16 + 215.89 x17 + 218.75 x18 + 
    230.71 x19 - 324.23 x2 - 545.6 x20 - 283.21 x21 + 676.47 x22 + 
    18.24 x23 - 475.34 x3 + 886.15 x4 - 887.09 x5 - 398.95 x6 - 
    660.04 x7 - 89.91 x8 + 264.35 x9 == 891.45;
equation20 = 
 913.14 x1 + 578.23 x10 + 544.77 x11 - 281.47 x12 + 45.8 x13 + 
   818.8 x14 - 27.56 x15 - 853.88 x16 - 734.34 x17 - 887.11 x18 - 
   426.76 x19 - 105.89 x2 - 286.53 x20 - 485.57 x21 - 762.6 x22 + 
   137.72 x23 - 270.55 x3 + 179.56 x4 + 109.98 x5 + 300.58 x6 - 
   578.92 x7 + 135.87 x8 + 270.43 x9 == 863.38; equation21 = 
 244.79 x1 + 891.92 x10 - 208.48 x11 - 199.29 x12 + 25.62 x13 + 
   510.89 x14 + 300.53 x15 + 953.35 x16 + 960. x17 + 341.88 x18 - 
   780.57 x19 - 772.53 x2 - 347.5 x20 + 328.68 x21 - 946.35 x22 + 
   877.18 x23 + 475.96 x3 - 991.6 x4 - 636.24 x5 + 37.6 x6 - 
   233.62 x7 + 851.72 x8 - 65.36 x9 == -932.19;
equation22 = -728.22 x1 + 571.07 x10 + 523.25 x11 - 47.01 x12 - 
    214.58 x13 + 35.19 x14 + 642.52 x15 - 125.15 x16 + 374.71 x17 - 
    116.44 x18 + 705.87 x19 - 355.76 x2 + 308.85 x20 + 275.03 x21 + 
    401.73 x22 + 424.5 x23 + 382.19 x3 - 804.71 x4 + 448.04 x5 - 
    606.38 x6 + 675.33 x7 - 360.01 x8 - 136.95 x9 == -562.9;
equation23 = 
  516.12 x1 - 189.52 x10 - 543.08 x11 - 610.05 x12 - 906.02 x13 - 
    388.31 x14 - 294.28 x15 + 377.96 x16 - 433.34 x17 + 96.14 x18 - 
    118.19 x19 - 577.27 x2 - 842.95 x20 - 124.59 x21 - 400.7 x22 + 
    695.05 x23 - 207.36 x3 + 919.4 x4 + 705.66 x5 + 921.18 x6 - 
    16.24 x7 + 783.9 x8 + 599.5 x9 == -805.68;
system1 = {equation1, equation2, equation3, equation4, equation5, 
   equation6, equation7, equation8, equation9, equation10, 
   equation11, equation12, equation13, equation14, equation15, 
   equation16, equation17, equation18, equation19, equation20, 
   equation21, equation22, equation23};
mSystem = 
  Normal@CoefficientArrays[
    system1, {x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, 
     x14, x15, x16, x17, x18, x19, x20, x21, x22, x23}];
B = -mSystem[[1]];
A = mSystem[[2]];
(*times1Style@Column[{TraditionalForm@DisplayForm[GridBox[{{"\
\[Piecewise]",Column[system1]}}]],
Row[{"A = ",MatrixForm[A],"; b = ",MatrixForm[B]}]
}]*)
Ainv = Inverse[A];
X = Ainv.B
POSTED BY: Hans Dolhaine

Thank you! I have another question as I am ew to using Wolfram Mathematica. Is there another way for me to display a linear equation as a matrix other than the way I used it below?

(*times1Style@Column[{TraditionalForm@DisplayForm[GridBox[{{"\
\[Piecewise]",Column[system1]}}]],
Row[{"A = ",MatrixForm[A],"; b = ",MatrixForm[B]}]
}]*)

Your matrix A is not a square matrix, so it has no inverse

Dimensions[A]
POSTED BY: Hans Dolhaine
Posted 5 years ago
POSTED BY: Rohit Namjoshi

Thank you so much for the assistance. I'm brand new to Mathematica, so I didn't even notice I had equation 13 in there twice, in my system1 statement.

Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard