Tag: Integral

  • Integral Numerik Metode Trapezoida Dengan Matlab

    Integral Numerik Metode Trapezoida Dengan Matlab

    Integral Numerik Metode Trapezoida dilakukan dengan cara membuat bangun-bangun trapesium hayal diantara luas daerah yang dibatasi oleh sebuah fungsi.

    Integral Numerik Metode Trapezoida

    Metode Trapezoida merupakan salah satu metode numerik yang digunakan untuk menghitung aproksimasi integral dari suatu fungsi. Metode ini didasarkan pada pembagian interval integrasi menjadi sejumlah kecil sub-interval dan menghitung luas di bawah kurva dengan mengaproksimasi kurva sebagai segmen-segmen garis lurus yang menghubungkan titik-titik pada fungsi.

    Integral Metode Trapezoida

    A. Dasar Teori

    Secara matematis, jika kita ingin menghitung integral dari fungsi f(x) di antara batas a dan b:

    \int^b_af(x)dx

    Metode Trapezoida membagi interval [a,b] menjadi n sub-interval dengan panjang yang sama,

    Δx=\frac{b-a}{n}

    Titik-titik pembagi ini diberi label x0, x1, x2,…, xn, di mana:

    xi​=a+iΔx untuk i = 1, 2, 3, …, n

    Luas di bawah kurva dihitung dengan mengaproksimasi setiap sub-interval sebagai sebuah trapesium. Luas setiap trapesium adalah:

    L=\frac{1}{2}(f(x_i)+f(x_{i+1}))Δx

    Total luas semua trapesium, yang merupakan aproksimasi integral, diberikan oleh:

    \int^b_af(x)dx≈\frac{Δx}{2}[f(x_0)+2∑^{n-1}_{i=1}f(x_i)+f(x_n)]

    B. Langkah-Langkah Menggunakan Metode Trapezoida

    Tentukan Batas Integrasi dan Jumlah Sub-Interval (n):

    Misalkan kita ingin mengintegrasikan fungsi f(x) dari a ke b dan kita pilih n sub-interval.

    Hitung Lebar Sub-Interval (Δx):

    Δx=\frac{b-a}{n}

    Hitung Nilai Fungsi di Titik-Titik Pembagi:

    • Untuk setiap titik xi​=a+iΔx , hitung f(xi) untuk i=0,1,2,…,n

    Gunakan Rumus Trapezoida untuk Mengaproksimasi Integral:

    \int^b_af(x)dx≈\frac{Δx}{2}[f(x_0)+2∑^{n-1}_{i=1}f(x_i)+f(x_n)]

    C. Contoh Kasus

    Misalkan kita ingin menghitung aproksimasi integral dari f(x)=sin⁡(x) dari 0 sampai π dengan n=4:

    A. Solusi Numerik

    1. Tentukan Batas dan Sub-Interval:

    a=0, b=\pi , n = 4
    \Delta x=\frac{\pi-0}{4}=\frac{\pi}{4}

    2. Hitung Titik-Titik Pembagi dan Nilai Fungsi setiap kelipatan π/4, dimulai dari 0, π/4, π/2, 3π/4 dan π.

    f(x_0)=\sin 0 = 0
    f(x_1)=\sin \frac{\pi}{4}=\frac{\sqrt{2}}{2}
    f(x_2)=\sin \frac{\pi}{2}=1
    f(x_3)=\sin \frac{3\pi}{4}=\frac{\sqrt{2}}{2}
    f(x_4)=\sin π = 0

    3. Aproksimasi Integral:

    \int^π_0\sin x \ dx≈\frac{\frac{\pi}{4}}{2}[0+2(\frac{\sqrt{2}}{2}+1+\frac{\sqrt{2}}{2})+0]
    \int^π_0\sin x \ dx≈\frac{\frac{\pi}{4}}{2}[0+2(2.4142)+0]
    \int^π_0\sin x \ dx≈\frac{\pi}{4}(2.4142)
    \int^π_0\sin x \ dx≈ 0.6 \pi≈1.8961

    B. Solusi Analitik

    \int^π_0\sin x \ dx = - cos |^π_0=-(-1-1)=2

    Nilai dari solusi numerik ini mendekati nilai solusi analitik.

    D. Trapezioda dengan Matlab

    Buat Fungsi untuk Metode Trapezoida: Kita akan membuat fungsi yang menerima parameter fungsi yang akan diintegrasikan, batas bawah dan atas integrasi, serta jumlah sub-interval.

    function result = trapezoidal_rule(f, a, b, n)
        % f: fungsi yang akan diintegrasikan
        % a: batas bawah
        % b: batas atas
        % n: jumlah sub-interval
    
        % Hitung lebar tiap sub-interval
        h = (b - a) / n;
        
        % Hitung nilai fungsi di titik-titik pembagi
        x = a:h:b;
        y = f(x);
        
        % Terapkan rumus metode Trapezoida
        result = (h / 2) * (y(1) + 2 * sum(y(2:end-1)) + y(end));
    end

    Gunakan Fungsi untuk Mengaproksimasi Integral: Panggil fungsi trapezoidal_rule dengan parameter yang sesuai.

    % Definisikan fungsi yang akan diintegrasikan
    f = @(x) sin(x);
    
    % Batas integrasi
    a = 0;
    b = pi;
    
    % Jumlah sub-interval
    n = 4;
    
    % Hitung integral menggunakan metode Trapezoida
    approx_integral = trapezoidal_rule(f, a, b, n);
    
    % Tampilkan hasil
    disp(['Aproksimasi integral: ', num2str(approx_integral)]);
    

    Penjelasan Kode

    1. Fungsi trapezoidal_rule:
      • f: Fungsi yang akan diintegrasikan, didefinisikan sebagai fungsi anonim.
      • a dan b: Batas bawah dan atas dari integral.
      • n: Jumlah sub-interval yang digunakan dalam metode Trapezoida.
      • h: Lebar tiap sub-interval, dihitung sebagai (b−a)/n.
      • x: Vektor yang berisi titik-titik pembagi dari aaa ke b dengan jarak h.
      • y: Nilai fungsi f di titik-titik pembagi.
      • result: Aproksimasi integral yang dihitung menggunakan rumus metode Trapezoida.
    2. Menghitung Integral:
      • Fungsi anonim f didefinisikan sebagai @(x) sin(x).
      • Batas integrasi adalah a = 0 dan b = pi.
      • Jumlah sub-interval adalah n = 4.
      • Hasil aproksimasi integral ditampilkan menggunakan disp.

    E. Tugas

    Buatlah sebuah solusi integral numerik metode trapzoida untuk fungsi berikut

    \int^b_a(3x^3-5) dx

    dan

    \int^b_a(\cos x +2)dx

    keterangan

    1. ganti nilai b dengan tanggal lahir anda masing-masing dan nilai a dengan bulan lahir.
    2. Jumlah sub-interval yang digunakan n = 4
  • Belajar Matlab – Solusi Integral Metode Simpson

    Belajar Matlab – Solusi Integral Metode Simpson

    AhmadDahlan.NET – Metode Simpson adalah metode integral numerik yang digunakan untuk menghitung luas daerah yang dibatasi oleh persamaan garis f(x). Pendekatan yang dilakukan lebih detail dari Pendekatan Trapezoid dimana Daerah di bagi ke dalam dua bangun trapesium.

    Solusi Integral Metode Simpson

    Metode Indtegral dengan Pendekatan Simpson

    Gambar pada sisi kiri menunjukkan metode integral analitik untuk menghitung luas wilayah yang dibatasi oleh garis f(x) yang mulai dari a sampai b. Pada gambar pada sisi kanan adalah metode Simpson yang digunakan untuk menghitung luas daerah yang diarsis menggunakan metode numerik Simpson.

    Metode ini Simpson sama dengan metode Trapezoida namun luas wilayah di bagi ke dalam dua trapseium sehingga hasil perhitungan jauh lebih teliti dibandingkan dengan metode Simpson. Lebar trapesium di bagi menjadi dua bagian dengan lebar h.

    Metode Simpson dibagi ke dalam dua kelompok yakni :

    1. Metode simpson 1/3
    2. Metode simson 3/8

    A. Metode simpson 1/3

    Metode ini mengaproksimasi integral dengan menggunakan polinomial kuadrat. Rumus integral Simpson 1/3 untuk fungsi f(x) pada interval [a,b] adalah sebagai berikut:

    \int^b_a f(x).dx≈\frac{b-a}{6}\left[f(a)+4f\left(\frac{a+b}{2}\right)+f(b)\right]

    Untuk membagi ke lebih banyak sub-interval dapat digunakan:

    \int^b_a f(x).dx≈\frac{h}{3}\left[f(x_0)+4\Sigma^{n-1}_{i=1,3,5, ...}f(x_i)+2\Sigma^{n-2}_{i=2,4,6, ...}f(x_i)+f(x_n)\right]

    Di mana:

    • n adalah jumlah sub-interval (harus genap),
    • h=(b−a)/n
    • xi=a+ih

    Contoh Implementasi dengan Matlab

    function integral = simpson_1_3(f, a, b, n)
        if mod(n, 2) == 1
            error('n harus genap');
        end
        
        h = (b - a) / n;
        x = a:h:b;
        fx = arrayfun(f, x);
        
        result = fx(1) + fx(end);
        result = result + 4 * sum(fx(2:2:n));
        result = result + 2 * sum(fx(3:2:n-1));
        
        integral = result * h / 3;
    end
    
    % Contoh penggunaan
    f = @(x) sin(x);
    a = 0;
    b = pi;
    n = 100;  % harus genap
    
    integral = simpson_1_3(f, a, b, n);
    fprintf('Nilai integralnya adalah: %.6f\n', integral);

    B. Metode Simpson 3/8

    Metode ini mengaproksimasi integral dengan menggunakan polinomial kubik. Rumus integral Simpson 3/8 untuk fungsi f(x) pada interval [a,b] adalah sebagai berikut:

    \int^b_af(x).dx≈\frac{3h}{8}\left[ f(a)+3f(a)\left(a+\frac{h}{3}\right) +3f(a)\left(a+\frac{2h}{3}\right)+f(b)\right]

    Untuk lebih banyak sub-interval, rumus umum Simpson 3/8 adalah:

    \int^b_af(x).dx≈\frac{3h}{8}\left[ f(x_0)+3\Sigma^{n-1}_{i=1,2,4,5 ...}f(x_i)+2\Sigma^{n-3}_{i=3,6,9, ...}f(x_i)+f(x_n)\right]

    dimana n kelipatan bilangan 3.

    Contoh Implementasi dengan Matlab

    function integral = simpson_3_8(f, a, b, n)
        if mod(n, 3) ~= 0
            error('n harus kelipatan 3');
        end
        
        h = (b - a) / n;
        x = a:h:b;
        fx = arrayfun(f, x);
        
        result = fx(1) + fx(end);
        result = result + 3 * sum(fx(2:3:n));
        result = result + 3 * sum(fx(3:3:n));
        result = result + 2 * sum(fx(4:3:n-3));
        
        integral = result * 3 * h / 8;
    end
    
    % Contoh penggunaan
    f = @(x) sin(x);
    a = 0;
    b = pi;
    n = 99;  % harus kelipatan 3
    
    integral = simpson_3_8(f, a, b, n);
    fprintf('Nilai integralnya adalah: %.6f\n', integral);
    

    E. Tugas

    Buatlah sebuah solusi integral numerik metode trapzoida untuk fungsi berikut

    ∫_𝑎^𝑏(3𝑥3−5)𝑑𝑥

    dan

    ∫^b_a(cos⁡𝑥+2)𝑑𝑥

    keterangan

    1. ganti nilai b dengan tanggal lahir anda masing-masing dan nilai a dengan bulan lahir.
    2. Kerjakan masing-masing soal dengan metode 1/3 dan 3/8.
  • Belajar Matlab -Penyelesaian Integral dengan Metode Trapezoida

    Belajar Matlab -Penyelesaian Integral dengan Metode Trapezoida

    AhmadDahlan.NET – Penyelesaian Integral dengan Metode Trapezoida dilakukan dengan persamaan berikut :

    \int^b_af(x)dx=\frac{h}{2}[f(x_0)+f(x_1)]-\frac{h^3}{12}f''(ξ)

    dimana x0=a dan x1=b dan h = b – a. Pada umumnya turunan suku ke-2 f”, biasanya dapat diabaikan sehingga persamaan ini dapat disederhanakan menjadi :

    \int^b_af(x)dx=\frac{h}{2}[f(x_0)+f(x_1)]

    Pendekatan Trapezoida hanya bisa digunakan untuk persamaan yang turunan keduannya nol. Grafik dari pendekatan Trapezoida seperti berikut ini :

    Solusi persamaan integral dengan pendekatan Pendekatan Trapezoida

    Grafik pada bagian kiri adalah grafik yang menunjukkan persamaan f(x). Metode Trapezodia digunakan untuk menghitung luas daerah yang menyerupai bentuk trapesium di bawah garid f(x) dengan batas dari a sampai b.

    Perhatikan daerah antara garis f(x) dengan garis lurus antara f(x1) dan f(x0). Metode trapezodia tidak mempu menghitung luad daerah tersebut sehingga jika nilainya terlalu besar, Metode ini tidak menunjukkan hasil yang teliti.

    A. Studi Kasus Solusi Integral

    Misalkan ada sebuah persamaan f(x) = x2 dengan batas atas b = 6 dan a = 2. Bandingkan hasil keduanya dengan metode Trapezodia dengan Metode Analitik!

    a. Metode Trapezodia

    Nilai h = b – a = 6 – 2 = 4

    f(a) = a2 = 4
    f(b) = b2 = 36

    Integral Metode Trapezodia

    \int^6_2x^2dx=\frac{4}{2}[36+4]=80

    b. Metode Analitik

    \int^6_2x^2dx= \frac{1}{3}x^3|^6_2 = \frac{1}{3}(216-8) = 69,33

    Latihan

    Bandinkan hasil Metode Trapezodia dengan Metode Analitik untuk f(x) = x1 dan f(x) = x3!

    B. Membuat Script Metode Trapezodia

    Membuat script di Matlab

    clear all 
    clc
    
    a = ...
    b = ...
    
    x0=a;
    x1=b;
    h = b-a;
    
    Int_trapez = (h/2)*(f(x0)+f(x1))
    

    Fungsi external

    function y = f(x)
    y = ... 

    Tugas !

    Tuliskan script Metode Trapezodia dengan Matlab untuk menyelesaikan persamaan :

    f(x)=\sqrt{1-x}

    Bandingkan hasilnya dengan metode analitik, jika batas awal a = 1 dan batas akhir b = 3!