Moh Khairul Imam Assidiqqi
السَّلاَمُ عَلَيْكُمْ وَرَحْمَةُ اللهِ وَبَرَكَاتُهُ
Profile | |
---|---|
Nama Lengkap | Moh Khairul Imam Assidiqqi |
NPM | 1906301210 |
Jurusan | Teknik Mesin |
Alamat Email | moh.khairul@ui.ac.id |
Contents
Perkenalan Singkat
Assalamu'alaikum Wr. Wb. Perkenalkan nama saya Moh Khairul Imam Assidiqqi. Saya adalah mahasiswa Departemen Teknik Mesin angkatan 2019.
Pertemuan Pertama Metode Numerik (11-11-2020)
Pada hari Rabu, 11 November 2020 adalah pertemuan pertama kelas Metode Numerik bersama Pak Dr. Ir. Ahmad Indra Siswantara. Pak Dai memberikan penjelasan terlebih dahulu bagaimana cara menggunakan Wiki AIR. Kemudian beliau menjelaskan apa itu metode numerik dan untuk apa kita mempelajarinya.
Ada 4 tujuan mempelajari metode numerik:
1. Memahami konsep dan prinsip dasar di dalam metode numerik.
2. Dapat menerapkan pemahaman terhadap konsep di dalam permodelan numerik.
3. Mampu menerapkan metode numerik di dalam persoalan keteknikan.
4. Mendapat nilai tambah sehingga menjadi orang yang beradab.
Setelah itu Pak Dai menjelaskan secara singkat tentang Open Modelica.
Tugas 1 Metode Numerik
Sebelum kelas berakhir, Pak Dai memberikan kami tugas untuk mempelajari Open Modelica. Setelah itu kami ditugaskan untuk membuat video tutorial cara menggunakan Open Modelica.
Berikut adalah hasil tugas saya:
Pertemuan Kedua Metode Numerik (18-11-2020)
Pada pertemuan kedua dengan Pak Dai, kami diberikan pelajaran yang sangat berharga bahwa kita harus menjadi orang yang beruntung. Maksudnya adalah kita harus menjalani hari ini lebih baik dari hari kemarin. Dan jangan menjadi orang yang merugi yang hari malah lebih buruk dari hari kemarin.
Setelah itu Pak Dai juga menyuruh kami untuk menayangkan tugas-tugas kami pada minggu pertama. Kami diperbolehkan bertanya apapun yang masih kami belum ketahui. Kemudian Pak Dai mengajarkan kami tentang Class dan Function di OpenModelica dan Pak Dai juga mengajarkan bahwa suatu fungsi bisa dipanggil ke sebuah kelas.
Tugas 2 Metode Numerik
Di akhir kelas, Pak Dai memberikan kami tugas untuk membuat video dimana kami diminta untuk membuat sebuah fungsi berupa persamaan aljabar simultan dengan variabel array kemudian membuat class untuk memanggil fungsi tersebut.
Berikut adalah hasil tugas saya:
function
function Gauss input Real A[5,5]; input Real B[5]; output Real X[5]; algorithm X:=Modelica.Math.Matrices.solve(A,B); end Gauss;
class
class SPL_5_Variabel parameter Real A[5,5]=[2,-2,0,4,2; -2,4,2,0,5; 2,-2,2,-2,0; 0,4,2,4,4; 4,2,2,4,-2]; parameter Real B[5]={2,4,6,8,-12}; Real X[5]; equation X=Gauss(A,B); end SPL_5_Variabel;
Pertemuan Ketiga Metode Numerik (25-11-2020)
Pada pertemuan ketiga, Pak Dai memaparkan tiga aplikasi metode numerik yang sering digunakan dalam menyelesaikan permasalahan teknik, pertama ada Computation Fluid Dynamics (CFD), lalu Finite Element Analysis, dan Metode Stokastik. Setelah itu Pak Dai meminta beberapa mahasiswa untuk mereview tugas-tugas yang telah diberikan sebelumnya. Salah satu teman kami yaitu Josiah memberikan sharing mengenai aljabar simultan dan variabel array.
Kemudian Pak Dai menjelaskan kembali mengenai metode Gauss. Pak dai mendorong kami untuk membuat code sendiri untuk menyelesaikan Gauss Jordan. Salah satu mahasiswa, yaitu Cristo memberikan sharing mengenai code yang telah ia buat untuk Gauss Jordan, berikut code modelica untuk Gauss Jordan dari Cristo:
// Gauss-Jordan Algorithm
// Transforms input matrix A into reduced row echelon form matrix B
// Christopher S.E. November 2020
input Real [:,:] A; // An augmented matrix of m*n
output Real [:,:] B; // Output matrix in reduced row echelon form
// Local variables
protected
Integer h = 1; // Initialize pivot row
Integer k = 1; // Initialize pivot column
Integer m = size(A,1); // Number of rows in matrix
Integer n = size(A,2); // Number of columns in matrix
Integer c = 0; // Index counter
Integer max_row; // Row index of max number in pivot column
Real [:] pivot_column; // Vector containing pivot column data
Real [:] pivot_row; // Vector containing backwards pivot row data
Real [:,:] temp_array; // Stores matrix data when switching rows
Real r; // Ratio value used for row operations
// Limit to handle floating point errors
Real float_error = 10e-10;
algorithm
// Transfer input matrix A into variable B
B := A;
while h <= m and k <= n loop
// Dealing with floating point errors
for i in 1:m loop
for j in 1:n loop
if abs(B[i,j]) <= float_error then
B[i,j] := 0;
end if;
end for;
end for;
// Finding the pivot
pivot_column := {B[i,h] for i in h:m};
// Get position index of lowest row with greatest pivot number
c:= h-1;
for element in pivot_column loop
c := c+1;
if abs(element) == max(abs(pivot_column)) then
max_row := c;
end if;
end for;
// No pivot in this column, move on to next column
if B[max_row, k] == 0 then
k := k+1;
else
// Swap Rows h <-> max_row
temp_array := B;
temp_array[h] := B[max_row];
temp_array[max_row] := B[h];
B:= temp_array;
// Divide pivot row by pivot number
B[h] := B[h]/B[h,k];
// For all rows below the pivot
for i in (h+1):m loop
// Store the ratio of the row to the pivot
r := B[i,k] / B[h,k];
// Set lower part of pivot column to zero
B[i,k] := 0;
// Operations on the remaining row elements
for j in (k+1):n loop
B[i,j] := B[i,j] - B[h,j] * r;
end for;
end for;
// Move on to next pivot row and column
h := h+1;
k := k+1;
end if;
end while;
// The matrix is now in row echelon form
// Set values of (h,k) to (m,n)
h := m;
k := n;
while h >= 1 and k >=1 loop
// Dealing with floating point errors
for i in 1:m loop
for j in 1:n loop
if abs(B[i,j]) <= float_error then
B[i,j] := 0;
end if;
end for;
end for;
// Finding the pivot
pivot_row := {B[h,i] for i in 1:k};
// Get position index k of pivot
c := 0;
for element in pivot_row loop
c := c+1;
if element <> 0 then
break;
end if;
end for;
k := c;
// No pivot in this row, move on to next row
if B[h, k] == 0 then
h := h-1;
else
// Perform row operations
for i in 1:(h-1) loop
r := B[i,k];
B[i] := B[i] - B[h] * r;
end for;
// Move on to next pivot row and column
h := h-1;
k := k-1;
end if;
end while;
// The matrix is now in reduced row echelon form
end GaussJordan;
Tugas 3 Metode Numerik
Di akhir kelas, Pak Dai memberikan kami tugas yaitu menyelesaikan persoalan berikut:
Berikut adalah penyelesaiannya: