Tutorial kalkulus tentang optimisasi

download di: https://drive.google.com/file/d/0B1irLqfPwjq0bW9GeFZrbWRTT0k/edit?usp=sharing

Convert List Generic ke list tipe tertentu di java

Berikut ini adalah contoh bagaimana mengconvert list generic ke list tipe tertentu di java

Belajar optimisasi di MATLAB

Sebenarnya untuk kreatif g perlu mahal, asalkan kita banyak baca. Berikut ini saya akan membahas sebuah contoh dasar pada buku kalkulus James stewart edisi 5 pada bab 4.7 yakni soal latihan nomor 57.

Yang hasilnya saya gambar ulang menjadi

Pada gambar di atas, terlihat jelas bahwa sudut \(\theta\) dapat dinyatakan sebagai fungsi \(x\) yakni $$ \begin{equation} \theta = 180^o -tan^{-1} \left ( \frac{2}{x} \right ) - \tan^{-1} \left( \frac{5}{3-x} \right ) \end{equation} $$ yang nilai maksimumnya dapat dinyatakan sebagai $$ \begin{eqnarray} \frac{d \theta}{dx } &= & 0 + \frac{2}{x^2} \frac{1}{1 + \frac{4}{x^2} } - \frac{5}{(3 - x)^2} \cdot \frac{1}{1 + \frac{25}{(3 - x)^2}} \\ & = & \frac{2}{x^2 + 4} - \frac{5}{(3 - x)^2 + 25} = 0 \end{eqnarray} $$ atau \begin{eqnarray} 3 x^2 + 12 x - 12 = 0 \end{eqnarray} yang akarnya adalah \(x = - 6 \pm 6\sqrt{2} \). Di mana yang memenuhi adalah \(x = -6 + 6 \sqrt{2} = 6 (\sqrt{2} - 1) \approx 2.4\) Untuk mengilustrasikan jawaban tersebut, saya sudah membuat sebuah program sederhana dalam bahasa MATLAB yang sourcenya pembaca bisa coba di rumah, ditulis dengan menggunakan MATLAB 2009. Adapun sourcenya adalah:

% jawaban soal nomor 57 BAB 4.7 dari buku James Stewart kalkulus
% edisi 5
% author: Mohammad Fajar
function soal57
clc; 
f = figure('menubar', 'none', 'resize', 'off', ... 
    'position', [200, 100, 800, 450]); 
ax1 = axes('units', 'pix', 'position', [50 50 340 270]... 
   ,  'xlim', [0 2*pi], 'ylim', [-1.5, 1.5]); 

ax2 = axes('units', 'pix', 'position', [430 50 340 270]... 
   ,  'xlim', [0 2*pi], 'ylim', [-1.5, 1.5]); 

slider =  uicontrol('style', 'slider', 'position', [50 380 90 20]); 
text_ = uicontrol('style', 'text', 'position', [150 380 40 20 ]); 

set(slider, 'value', 0.1);
set(slider, 'sliderstep', [0.05 0.05]);
ss = 0:0.01:3; 
theta = 180 - atand(2./ss) - atand(5./(3- ss));
maksimum = max(theta); % jawabannya
while ishandle(f)
    xx = get(slider , 'value');
    xx = xx * 3;
    set(text_, 'string', num2str(xx)); 
    trigx = [0 0;5 2;0 0]; 
    trigy = [0 3-xx;0 3;3-xx 3];
    z = [1 1;1 1;1 1];
    p = patch(trigx, trigy, z, 'parent', ax2); 
    axis([-1 ,6, -1, 4]);
    thetax = 180 - atand(2./xx) - atand(5./(3 - xx)); 
    p2 = plot(ss, theta, 'parent', ax1);
    xlabel(ax1, 'x');
    ylabel(ax1, '\theta'); 
    hold on; 
    p3 = line([xx xx], [0 thetax], 'linewidth', 4, 'parent', ax1);
    tt = text(xx+.1,thetax-3,['\theta = ', num2str(thetax)],'parent', ax1 ...
        , 'fontsize', 9);
    tt2 = text(0.8,3-xx,['\leftarrow \theta = ', num2str(thetax), '^o'],'parent', ax2 ...
        , 'fontsize', 12);
    % ini hanya untuk mengakali karena step slider dan step linspace pada 
    % MATLAB tidak sama :)
    if floor(thetax) == floor(maksimum)
        set(p3, 'color', 'r');
        set(p, 'facecolor', 'y'); 
    end
    drawnow;
    if ishandle(p), delete(p);end
    if ishandle(p2), delete(p2);end
    if ishandle(p3), delete(p3);end
    if ishandle(tt), delete(tt);end
    if ishandle(tt2), delete(tt2);end
end

adapaun hasil runningnya adalah:

Memahami garis singgung dengan MATLAB

Berikut ini saya berikan contoh program untuk memahami makna filosofis dari garis singgung.

function garis_singgung
% author: mohammad fajar :)
clc; 
f = figure('menubar', 'none', 'resize', 'off'); 
axes('units', 'pix', 'position', [50 50 450 300]... 
   ,  'xlim', [0 2*pi], 'ylim', [-1.5, 1.5]); 
slider =  uicontrol('style', 'slider', 'position', [50 380 90 20]); 
text = uicontrol('style', 'text', 'position', [150 380 40 20 ]); 
t = 0:.1:2*pi;  
y1 = sin(t);  
y2 = cos(t); 
plot(t, y1 , t, y2);
hold on ;
line([0 2*pi], [0 0]);
hold on ;
while ishandle(f)
    value = get(slider , 'value');
    value = value * 2 * pi; 
    set(text, 'string', num2str(value)); 
    hasil = cos(value); 
    yy = (t - value ).* hasil + sin(value) ;
    p = plot(t, yy, 'r' , 'linewidth', 2); 
    p1 = line([value value ] , [0 cos(value)], 'linewidth', 4);
    p2 = plot(value, sin(value),'r.' , 'markersize', 23); 
    % method drawnow tidak mengenali settingan di atas
    axis([0 , 2*pi, -1.5, 1.5 ]); 
    drawnow; 
    if ishandle(p), delete(p); end
    if ishandle(p1), delete(p1); end
    if ishandle(p2), delete(p2); end
end
end

Adapun hasilnya adalah pada gambar berikut:

Belajar monte carlo dengan MATLAB untuk menghitung luas kurva

Berikut ini contoh program yang saya buat pake MATLAB untuk menghitung luas kurva "apa saja" dengan menggunakan metode Monte Carlo.

Yang source codenya 

function test_monte_carlo
clear all ;
clc;
f = figure('menubar', 'none', 'resize', 'off' ); 

ax = axes('units', 'pix', 'position', [40 50 340 340]);

uicontrol('units', 'pix', 'position', [400 360 60 25],... 
    'style', 'text','string','f(x) :' ,  'fontsize', 12, 'horizontalAlignment', 'left', ...
    'fontsize', 12, 'fontweight', 'bold', 'backgroundcolor', get(f, 'color'));

% input fungsi
inputFungsi = uicontrol('units', 'pix', 'position', [450 360 100 25],... 
    'style', 'edit', 'fontsize', 12);

uicontrol('units', 'pix', 'position', [400 330 60 25],... 
    'style', 'text','string','N :' ,  'fontsize', 12, 'horizontalAlignment', 'left', ...
    'fontsize', 12, 'fontweight', 'bold', 'backgroundcolor', get(f, 'color'));

inputN = uicontrol('units', 'pix', 'position', [450 330 100 25],... 
    'style', 'edit', 'fontsize', 12, 'fontweight', 'bold');

% keterangan hasil
ket_hasil  = uicontrol('units', 'pix', 'position', [400 280 100 30],... 
    'style', 'text', 'fontsize', 12, 'fontweight', 'bold', ...
    'string', '0' , 'backgroundcolor', get(f, 'color'));

% hitung button
uicontrol('units', 'pix', 'position', [400 250 100 30],... 
    'style', 'pushbutton', 'fontsize', 12, 'fontweight', 'bold', ...
    'string', 'HITUNG', 'Callback', @fungsi_hitung);

% reset button
uicontrol('units', 'pix', 'position', [400 210 100 30],... 
    'style', 'pushbutton', 'fontsize', 12, 'fontweight', 'bold', ...
    'string', 'RESET', 'Callback', @reset_fungsi);

    function fungsi_hitung(varargin)
        N = get(inputN , 'string'); 
        N = str2double(N); 
        x = 0:.01:2*pi; 
        yy = get(inputFungsi, 'string'); 
        y = eval(yy, x); 
        plot(x,y, 'parent', ax, 'linewidth', 3);
        hold on;
        miny = min(y); 
        maxy = max(y);
        deltax = linspace( 0 , 2*pi ,10000); 
        deltay = linspace( miny,maxy, 10000); 
        a = 0;
        for i=1:N,
            m = randi(length(deltax),1);
            n = randi(length(deltay),1);
            x = deltax(m); 
            y = deltay(n);
            if  y  >= 0 
                if y  <= eval(yy, x)
                    a = a + 1;
                    plot(x,y, 'g*', 'parent', ax);
                    hold on ;
                else
                   plot(x,y, 'r*', 'parent', ax);
                   hold on ;
                end
            else
                if y  >= eval(yy, x)
                    a = a + 1;
                     plot(x,y, 'g*', 'parent', ax);
                     hold on ;
                else
                     plot(x,y, 'r*', 'parent', ax);
                     hold on ;
                end
            end
        end
        axis tight; 
        luas =  (a/N) * ( (maxy- miny) * (2*pi - 0)); 
        set(ket_hasil, 'string', num2str(luas)); 
    end
    
    function reset_fungsi(varargin)
        cla reset; 
        set(ket_hasil, 'string', '0'); 
    end
end

Menggunakan slider di MATLAB

Contoh penggunaan slider di MATLAB

Adapun source codenya:

function test_slider

figure('menubar', 'none', 'position', [300 300 200 100]); 

slider = uicontrol('style', 'slider', 'position',  [10 10, 100 30]... 
    , 'callback', @call);

ket = uicontrol('style', 'text', 'position',  [10 55, 100 40], ...
    'fontsize', 14, 'fontweight', 'bold' );

set(ket, 'string', get(slider, 'value')); 
        
    function call(varargin)
        set(ket, 'string', get(slider, 'value')); 
    end
end 

Menggunakan input field GUI di MATLAB

Berikut adalah contoh penggunaan input field di MATLAB.



adapun sourcenya adalah:

function tes_button
clc; 
figure('menubar', 'none', 'position', [300 300 200 100]); 

 uicontrol('style', 'pushbutton','position', [10 10, 100 50], 'string', 'TEST'...
 ,  'Callback', @call ); 

edit = uicontrol('style', 'edit','position', [10 65, 60 30], 'string', '0'.... 
    ,'fontsize', 12, 'fontweight','bold');

tampil = uicontrol('style', 'text', 'position', [80 65 60 25], 'string' , '0' ...
      ,'fontsize', 12, 'fontweight','bold'); 

        function call(varargin)
            a = get(edit, 'string'); 
            if ~isnan(str2double(a))
                set(tampil, 'string', a); 
            else
                set(tampil, 'string', 'SALAH');
            end
        end
end

Plot animasi sinusoidal di MATLAB

Berikut ini adalah contoh plot animasi sinusoidal di MATLAB

Adapun source codenya adalah

clc ;
awal = 0; 
akhir = 2*pi ; 
f = figure; 
temp = awal + .1;
step =0.05;
while ishandle(f) && temp <= akhir
 x = awal:step:temp;
 plot(x,sin(2.*x),'-ro'); 
 axis([0  akhir  -1 1]);
 drawnow;
 temp = temp+ step;
 pause(0.02);
end

while ishandle(f)
 awal = awal + step; 
 akhir = akhir + step; 
 x = awal:.025:akhir;
 plot(x,sin(2.*x),'-ro'); 
 axis([awal  akhir  -1 1]);
 drawnow;
 pause(0.02);
end
disp('program selesai'); 

Membuat 2 plot dalam satu figure di MATLAB

Berikut adalah source code membuat dua plot dalam satu figure di MATLAB.

function main
f = figure('menubar', 'none', 'position', [200 100 780 500]); 
panel_plot = uipanel('title', 'plot area', ...
            'units', 'pix', ...
            'positioN', [10 50 750 400] , ...
            'fontsize', 12);

g0 =axes('units','pixels','position',[50  50  300 300], ...
    'xlim', [0 (2*pi)], ... 
    'parent', panel_plot);

g1 =axes('units','pixels','position',[400 50   300 300], ...
    'xlim', [0 (2*pi)], ... 
    'parent', panel_plot);

x = 0:.1:2*pi; 
plot(x, sin(x), 'parent', g1); 
plot(x,sin(x.^2), 'parent', g0); 

Dengan contoh plot

medan vektor konservatif dan simply-connected

Salah satu hal yang menarik dari kalkulus vektor adalah persoalan bagaimana menentukan apakah suatu medan vektor bersifat konservatif atau bukan. Suatu medan vektor dikatakan konservatif jika integral yang dilakukan terhadap suatu lintasan tertutup yang berada di dalam domain medan vektor tersebut tidak bergantung pada lintasan yang diambil. Salah satu ciri dari medan vektor konservatif adalah curl yang diambil terhadap medan vektor tersebut sama dengan nol, atau bisa dikatakan medan vektor tersebut irotasional. Dengan demikian medan vektor tersebut dapat dinyatakan sebagai gradient suatu medan potensial skalar karena curl dari gradient potensial pasti nol. Contoh medan vektor konservatif adalah medan listrik dan medan gravitasi bumi.

Akan Tetapi pernyataan ini tidak bisa dibalik. Kita tidak bisa memastikan suatu medan vektor bersifat konservatif hanya dengan melihat apakah curl-nya nol atau bukan. Terdapat beberapa kasus di mana suatu medan vektor memiliki curl nol akan tetapi bukanlah medan vektor konservatif. Tinjau medan vektor dua dimensi berikut: \begin{equation} \vec{F} (x,y) = \left( \frac{-y}{x^2 + y^2}\right) \hat{i} + \left( \frac{x}{x^2 + y^2}\right)\hat{j} \end{equation} Yang curl-nya adalah $$ \begin{eqnarray} \nabla \times \vec{F} &=&\left[ \frac{\partial}{\partial x} \left( \frac{x}{x^2 + y^2}\right) + \frac{\partial }{\partial y} \left( \frac{-y}{x^2 + y^2}\right) \right] \hat{j} \\ &=& \left[ \frac{y^2 - x^2 }{\left( x^2 + y^2 \right)^2} + \frac{y^2 - x^2 }{\left( x^2 + y^2 \right)^2} \right] \hat{j} \\ &=& 0 \end{eqnarray} $$ Kita dapat memilih suatu kurva tertutup yang masih dalam domain medan vektor \(\vec{F}\) di mana integral yang diambil sepanjang kurva ini tidak sama dengan nol. Tinjau suatu kurva lingkaran \(C\) yang berpusat di titik \((0, 0)\) dan berjari-jari satu satuan yang dapat diparameterisasi sebagai $$ \begin{eqnarray} \vec{F }(c(t)) &= & \frac{-\sin t }{\cos^2 t + \sin^2 t} \, \hat{i} + \frac{\cos t }{\cos^2 t + \sin^2 t } \, \hat{j} \\ &=& -\sin t \, \hat{i} + \cos t \, \hat{j} \end{eqnarray} $$ dengan \( 0 <= t <= 2 \pi \). Integral yang dilakukan sepanjang kurva ini tidak sama dengan nol yakni $$ \begin{eqnarray} \int_C \vec{F} \cdot ds & = & \int_{a}^{b} \vec{F} (c(t)) \cdot c' (t) \, dt \\ &=& \int_{0}^{2 \pi} \left( -\sin t \, \hat{i} + \cos t \, \hat{j} \right) \cdot \left( - \sin t \, \hat{i} + \cos t \, \hat{j} \right) \\ &=& \int_{0}^{2\pi } 1 \, dt = 2 \pi \end{eqnarray} $$ Yang artinya terdapat suatu kurva yang mana integral \( \int_C \vec{F} \cdot ds \) tidak sama dengan nol yang mengindikasikan bahwa vektor \( \vec{F}\) bukanlah medan vektor konservatif.

Mengapa demikian adalah karena vektor \( \vec{F} \) memiliki titik singularitas yakni pada titik \((0,0)\). Artinya adanya titik singularitas ini mengakibatkan Teorema Stokes yang menghubungkan integral lintasan dengan curl vektor tidak bisa diterapkan pada medan vektor yang bersangkutan. Curl yang sama dengan nol baru bisa dipakai untuk menandai vektor konservatif jika vektor tersebut didefenisikan pada domain yang simply-connected. Adanya titik singularitas pada contoh di atas mengakibatkan domain tersebut tidak lagi simply-connected. Maksud simply connected sendiri adalah kita dapat mereduksi sembarang kurva dalam domain 2D atau permukaan dalam domain 3D menjadi sebuah titik dan masih tetap berada dalam domain tersebut. Untuk domain 2 dimensi simply connected mengharuskan bahwa seluruh domain terdefenisi secara utuh tanpa ada hole di dalam domain. Sementara untuk domain tiga dimensi sebuah hole pada domain bisa ditoleransi asalkan hole tersebut tidak menembus sampa ke luar domain. Syarat simply-connected ini merupakan syarat mutlak bagi teorema Stokes . Pada kasus dua dimensi adanya hole mengakibatkan kita tidak lagi bisa menghubungkan integral lintasan yang mengelilingi hole tersebut dengan integral permukaan pada keseluruhan hole, karena domainnya tidak terdefenisi pada hole. Sementara pada kasus tiga dimensi, walaupun terdapat suatu hole dalam daerah tertentu kita masih bisa mencari alternatif permukaan yang tidak melalui hole akan tetapi batasannya masih pada kurva yang sama dengan permukaan yang melalui hole tersebut.

Jadi kesimpulannya: suatu medan vektor bersifat konservatif jika curl yang dilakukan pada medan vektor tersebut sama dengan nol serta domain di mana medan vektor tersebut didefenisikan harus simply-connected.