Perbandingan Euler dan Runge-Kutta Orde Empat pada Gerak Rolling Kapal Nonlinier Amplitudo Besar

Info proyek

• Bidang studiTeknik Perkapalan dan Kelautan
• JenjangTesis (S2)
• Jenis metodeSimulasi Numerik
• Perkiraan pengerjaan± 5 hari kerja
• Perangkatloops dan fungsi anonim; polyfit, loglog (MATLAB inti)
• HargaRp 400.000 - Rp 2.400.000

Tujuan analisis

Mensimulasikan sudut rolling nonlinier dengan Perbandingan Euler dan Runge-Kutta dan menelaah sensitivitasnya terhadap koefisien redaman rolling dan tinggi metasentrik.

Variabel

Contoh hasil

• 0.25Step size h
• 0.15Max error Euler
• 0Max error RK4
• 10811.225Error ratio
• 1.049Order Euler
• 3.964Order RK4

Pada langkah h=0.25, galat maksimum Euler untuk Nonlinear roll angle adalah 0.1497 sedangkan RK4 hanya 1.4e-05; orde konvergensi terukur 1.05 (Euler) dan 3.96 (RK4) sesuai teori.

Dihitung dari data sintetis sebagai ilustrasi, bukan dari penelitian nyata.

Tentang Perbandingan Euler dan Runge-Kutta

Perbandingan Euler dan Runge-Kutta orde 4 menunjukkan bagaimana orde metode numerik memengaruhi akurasi dan laju konvergensi galat. Berguna untuk membenarkan pemilihan metode integrasi dan menunjukkan verifikasi numerik pada skripsi/tesis.

Langkah pengerjaan

1. 1Pilih persamaan dengan solusi analitik
2. 2Implementasikan metode Euler dan RK4
3. 3Jalankan dengan beberapa ukuran langkah h
4. 4Hitung galat terhadap solusi eksak
5. 5Taksir orde konvergensi dari kemiringan log-log
6. 6Plot galat terhadap h

Asumsi yang perlu diperiksa

• Solusi eksak diketahui untuk acuan
• Ukuran langkah cukup kecil agar berada pada rezim asimtotik
• Fungsi mulus sehingga orde teoretis berlaku

Membaca dan melaporkan hasil

Laporkan tabel galat tiap h dan orde konvergensi terukur (sekitar 1 untuk Euler dan 4 untuk RK4). Jelaskan trade-off antara biaya per langkah dan akurasi.

Contoh kode MATLAB

Skrip awal singkat dengan nama variabel proyek ini. Versi yang diserahkan disesuaikan dengan data Anda dan diberi komentar lengkap.

f = @(t, y) -2*y + t;
exact = @(t) 1.25*exp(-2*t) + t/2 - 0.25;
T = 2;
roll_damping_coefficient = [0.2 0.1 0.05 0.025 0.0125];             % step sizes
err = zeros(2, numel(roll_damping_coefficient));

for k = 1:numel(roll_damping_coefficient)
    h = roll_damping_coefficient(k);
    n = round(T / h);
    ye = 1; yr = 1; t = 0;
    for i = 1:n
        ye = ye + h * f(t, ye);                    % Euler
        k1 = f(t, yr);
        k2 = f(t + h/2, yr + h/2*k1);
        k3 = f(t + h/2, yr + h/2*k2);
        k4 = f(t + h, yr + h*k3);
        yr = yr + h/6 * (k1 + 2*k2 + 2*k3 + k4);   % RK4
        t = t + h;
    end
    err(1, k) = abs(ye - exact(T));
    err(2, k) = abs(yr - exact(T));
end
nonlinear_roll_angle = err;

pe = polyfit(log(roll_damping_coefficient), log(err(1,:)), 1);
pr = polyfit(log(roll_damping_coefficient), log(err(2,:)), 1);
fprintf('Euler order: %.2f, RK4 order: %.2f\n', pe(1), pr(1));
disp(array2table(err', 'VariableNames', {'Euler', 'RK4'}, 'RowNames', string(roll_damping_coefficient)));

figure;
loglog(roll_damping_coefficient, err(1,:), 'o-', roll_damping_coefficient, err(2,:), 's-', 'LineWidth', 1.2);
legend('Euler', 'RK4'); xlabel('Step size h'); ylabel('Absolute error'); grid on;

Paket dan harga

Pilih seberapa banyak bantuan yang Anda perlukan. Setiap paket yang lebih tinggi sudah mencakup semua item paket di bawahnya. Harga belum termasuk PPN 11%.

• Pembayaran: DP 50% sebelum pekerjaan dimulai, pelunasan setelah hasil diserahkan
• Pajak: Harga belum termasuk PPN 11%
• Data dan integritas akademik: Hasil dihitung dari data penelitian Anda sendiri; Anda yang menulis dan mempertanggungjawabkan naskah, kami membantu analisis dan pemahaman
• Revisi: Maksimal 2 kali revisi termasuk dalam harga (Paket 3: 3 kali)

Pertanyaan umum