Simulasi ode45 Gerak Tanah akibat Gelombang Seismik pada Model Satu Derajat Kebebasan

Info proyek

• Bidang studiGeofisika dan Seismologi
• JenjangTesis (S2)
• Jenis metodeSimulasi Numerik
• Perkiraan pengerjaan± 5 hari kerja
• Perangkatode45, odeset (MATLAB inti)
• HargaRp 450.000 - Rp 2.800.000

Tujuan analisis

Mensimulasikan perubahan perpindahan tanah dengan ode45 dan menilai sensitivitasnya terhadap kekakuan lapisan serta redaman material.

Variabel

Contoh hasil

• 61Accepted steps
• 0.656Mean step (s)
• 4.135Max of state 1
• 10.638Oscillation period (s)
• 0Max error vs tight solution
• 367Function evaluations

ode45 menyelesaikan sistem Ground displacement dan Layer stiffness dalam 61 langkah adaptif (rata-rata 0.656 s) dengan galat maksimum 3.9e-05 terhadap solusi toleransi ketat; periode osilasi sekitar 10.64 s.

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

Tentang Persamaan Diferensial dengan ode45

ode45 menyelesaikan persamaan diferensial biasa non-kaku dengan metode Runge-Kutta adaptif orde 4 dan 5. Dipakai untuk mensimulasikan model dinamik seperti populasi, reaksi kimia, dan sistem mekanik pada skripsi/tesis teknik atau sains.

Langkah pengerjaan

1. 1Tuliskan persamaan sebagai sistem orde satu dy/dt = f(t,y)
2. 2Definisikan fungsi turunan dan kondisi awal
3. 3Atur toleransi dengan odeset (RelTol, AbsTol)
4. 4Panggil ode45 pada rentang waktu yang diinginkan
5. 5Verifikasi hasil (kekekalan besaran atau solusi analitik)
6. 6Plot solusi terhadap waktu dan potret fase

Asumsi yang perlu diperiksa

• Sistem tidak kaku (bila kaku gunakan ode15s)
• Fungsi turunan mulus pada rentang integrasi
• Toleransi cukup ketat untuk akurasi yang dibutuhkan

Membaca dan melaporkan hasil

Laporkan model, parameter, kondisi awal, toleransi solver, dan bukti verifikasi numerik. Tampilkan grafik solusi serta tafsirkan perilaku dinamik yang diperoleh.

Contoh kode MATLAB

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

% Lotka-Volterra predator-prey model
a = 1.1; b = 0.4; d = 0.1; g = 0.4;
f = @(t, y) [a*y(1) - b*y(1)*y(2); d*y(1)*y(2) - g*y(2)];
opts = odeset('RelTol', 1e-8, 'AbsTol', 1e-10);
[layer_stiffness, ground_displacement] = ode45(f, [0 60], [10; 5], opts);   % time, states

% Conserved quantity V = d*x - g*ln(x) + b*y - a*ln(y)
V = d*ground_displacement(:,1) - g*log(ground_displacement(:,1)) + b*ground_displacement(:,2) - a*log(ground_displacement(:,2));
fprintf('Steps taken: %d\n', numel(layer_stiffness));
fprintf('Max drift of conserved quantity: %.2e\n', max(abs(V - V(1))));

figure;
subplot(1,2,1);
plot(layer_stiffness, ground_displacement(:,1), layer_stiffness, ground_displacement(:,2), 'LineWidth', 1.2);
legend('Prey', 'Predator'); xlabel('Time'); ylabel('Population');
subplot(1,2,2);
plot(ground_displacement(:,1), ground_displacement(:,2)); xlabel('Prey'); ylabel('Predator'); title('Phase portrait');

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