Membuat sebuah kontur dapat dilakukan dengan berbagai cara, bisa menggunakan surfer ataupun matlab. kamu tinggal pilih saja mana suka dan mana yang mudah. tetapi jika kamu menggunakan matlab untuk membuat contour silahkan pelajari syntax berikut;
clc;
clear;
%data
i=[1;1;1;1;1;1;1;1;1;1;1];
T=[20.6744;11.8442;33.2533;38.3768;27.5353;54.4092;57.8916;49.1237;56.2729;58.2464;58.1329];
z=[0;-0.1;-0.2;-0.3;-0.4;-0.5;-0.6;-0.7;-0.8;-0.9;-1];
%orde1
G1=[i z];
m1=inv(G1’*G1)*G1’*T;
a1=m1(1,1);
b1=m1(2,1);
h1=a1+(b1*z);
%orde2
z2=z.^2;
G2=[G1 z2];
m2=inv(G2’*G2)*G2’*T;
a2=m2(1,1);
b2=m2(2,1);
c2=m2(3,1);
h2=a2+(b2*z)+(c2*z2);
%orde3
z3=z.^3;
G3=[G2 z3];
m3=inv(G3’*G3)*G3’*T;
a3=m3(1,1);
b3=m3(2,1);
c3=m3(3,1);
d3=m3(4,1);
h3=a3+(b3*z)+(c3*z2)+(d3*z3);
%orde4
z4=z.^4;
G4=[G3 z4];
m4=inv(G4’*G4)*G4’*T;
a4=m4(1,1);
b4=m4(2,1);
c4=m4(3,1);
d4=m4(4,1);
e4=m4(5,1)
h4=a4+(b4*z)+(c4*z2)+(d4*z3)+(e4*z4);
%orde5
z5=z.^4;
G5=[G4 z5];
m5=inv(G5’*G5)*G5’*T;
a5=m5(1,1);
b5=m5(2,1);
c5=m5(3,1);
d5=m5(4,1);
e5=m5(5,1);
f5=m5(6,1);
h5=a5+(b5*z)+(c5*z2)+(d5*z3)+(e5*z4)+(f5*z5);
%orde6
z6=z.^6;
G6=[G5 z6];
m6=inv(G6’*G6)*G6’*T;
a6=m6(1,1);
b6=m6(2,1);
c6=m6(3,1);
d6=m6(4,1);
e6=m6(5,1);
f6=m6(6,1);
g6=m6(7,1);
h6=a6+(b6*z)+(c6*z2)+(d6*z3)+(e6*z4)+(f6*z5)+(g6*z6);
%orde7
z7=z.^7;
G7=[G6 z7];
m7=inv(G7’*G7)*G7’*T;
a7=m7(1,1);
b7=m7(2,1);
c7=m7(3,1);
d7=m7(4,1);
e7=m7(5,1);
f7=m7(6,1);
g7=m7(7,1);
i7=m7(8,1);
h7=a7+(b7*z)+(c7*z2)+(d7*z3)+(e7*z4)+(f7*z4)+(g7*z5)+(i7*z5);
%orde8
z8=z.^8;
G8=[G7 z8];
m8=inv(G8’*G8)*G8’*T;
a8=m8(1,1);
b8=m8(2,1);
c8=m8(3,1);
d8=m8(4,1);
e8=m8(5,1);
f8=m8(6,1);
g8=m8(7,1);
i8=m8(8,1);
j8=m8(9,1);
h8=a8+(b8*z)+(c8*z2)+(d8*z3)+(e8*z4)+(f8*z5)+(g8*z6)+(i8*z7)+(j8*z8);
%contour
figure(1);
w1=[T z h1];
y1=contour(w1);
text_handle=clabel(y1);
figure(2);
w2=[T z h2];
y2=contour(w2);
text_handle=clabel(y2);
figure(3);
w3=[T z h3];
y3=contour(w3);
text_handle=clabel(y3);
figure(4);
w4=[T z h4];
y4=contour(w4);
text_handle=clabel(y4);
figure(5);
w5=[T z h5];
y5=contour(w5);
text_handle=clabel(y5);
figure(6);
w6=[T z h6];
y6=contour(w6);
text_handle=clabel(y6);
figure(7);
w7=[T z h7];
y7=contour(w7);
text_handle=clabel(y7);
figure(8);
w8=[T z h8];
y8=contour(w8);
text_handle=clabel(y8);
keterangan:
text_handle=clabel(…) : untuk memunculkan nilai konturnya
…=contour(…) : plot contournya
Berdarkan data diatas kamu bisa saja meletakkan plot contournya di setiap orde, tetapi untuk memudahkan disini plot konturnya saya satukan di bagian akhir.
Berikut hasil contour untuk Figure 3:
