48 BAB 3. APLIKASI ELIMINASI GAUSS PADA MASALAH INVERSIadalah⎡⎢⎣8 219 8547 | 349, 890 2551, 88 159448, 88 | 3316, 570 159448.88 10656457, 88 | 221101, 6⎤⎥⎦9. Masih dalam proses triangularisai, operasi berikutnya (P 3 − (159448, 88/2551, 88)P 2 ) →P 3 . Hasilnya adalah⎡⎢⎣8 219 8547 | 349, 890 2551, 88 159448, 88 | 3316, 570 0 693609, 48 | 13872, 19⎤⎥⎦ (3.11)Seperti catatan yang lalu, saya ingin menyertakan pula notasi masing-masing elemenpada matrik augment sebelum melakukan proses substitusi mundur.⎡⎢⎣8 219 8547 | 349, 890 2551, 88 159448, 88 | 3316, 570 0 693609, 48 | 13872, 19⎤⎡⎥ ⎢⎦ ⇔ ⎣⎤a 11 a 12 a 13 | a 14⎥a 21 a 22 a 23 | a 24 ⎦a 31 a 32 a 33 | a 3410. Terakhir, tentukan konstanta m 1 , m 2 <strong>dan</strong> m 3 yang merupakan elemen-elemen vektorkolom m, dengan proses substitusi mundur. Pertama tentukan m 3m 3 = a 34a 33=13872, 19= 0, 02693609, 48lalu m 2m 2 = a 24 − a 23 m 3 3316, 57 − (159448, 88)(0, 02)= = 0, 05a 22 2551, 88<strong>dan</strong> m 1m 1 = a 14 − (a 12 m 2 + a 13 m 3 ) 349, 89 − [(219)(0, 05) + (8547)(0, 02)= = 21a 11 83.2.1 Script matlab inversi model parabolaPerbedaan utama script ini dengan script inversi model garis terletak pada inisialisasi elemenelemenmatrik kernel. Elemen-elemen matrik kernel sangat ditentukan oleh model matematikayang digunakan. Seperti pada script ini, matrik kernelnya diturunkan dari persamaan parabola.1 clc2 clear all34 z1=5;5 z2=8;6 z3=14;7 z4=21;
3.2. INVERSI MODEL PARABOLA 498 z5=30;9 z6=36;10 z7=45;11 z8=60;1213 T(1,1)=21.75;14 T(2,1)=22.68;15 T(3,1)=25.62;16 T(4,1)=30.87;17 T(5,1)=40.5;18 T(6,1)=48.72;19 T(7,1)=63.75;20 T(8,1)=96;2122 G(1,1)=1;23 G(1,2)=z1;24 G(1,3)=z1^2;25 G(2,1)=1;26 G(2,2)=z2;27 G(2,3)=z2^2;28 G(3,1)=1;29 G(3,2)=z3;30 G(3,3)=z3^2;31 G(4,1)=1;32 G(4,2)=z4;33 G(4,3)=z4^2;34 G(5,1)=1;35 G(5,2)=z5;36 G(5,3)=z5^2;37 G(6,1)=1;38 G(6,2)=z6;39 G(6,3)=z6^2;40 G(7,1)=1;41 G(7,2)=z7;42 G(7,3)=z7^2;43 G(8,1)=1;44 G(8,2)=z8;45 G(8,3)=z8^2;4647 G48 d=T;49 d5051 N=8; %jumlah data52 M=3; %model parameter53 pause545556 %%%%%===========Proses inversi==============57 disp(’Mencari G transpos’)58 for i=1:N59 for j=1:M60 GT(j,i)=G(i,j);61 end62 end63 GT6465 disp(’Perkalian GT <strong>dan</strong> G’)66 for i=1:M