Takuya Miyashita

Julia演習問題：最大値・最小値の検出 †

はじめに †

MATLABのpeaksっぽいfunctionを作成しておく．

function peaks(N=49::Int)
    function formula(x::T, y::T) where T<:Float64
        z = 3(1-x)^2*exp(-(x^2)-(y+1)^2)-10(x/5 - x^3 - y^5)*exp(-x^2-y^2)-(1/3)exp(-(x+1)^2-y^2)
    end
    vec = collect(Float64, linspace(-3.,3.,N))
    x = repmat(vec',N,1)
    y = repmat(vec,1,N)
    z = formula.(x,y)
    return (x,y,z)
end

Himmelblau's functionとか，自前でデモ用の2次元データを作成できる人はそれでも可．

問題 †

1. 上記のpeaksをpeaks.jlというファイルに保存して，functionをincludeしなさい．

   julia> include("peaks.jl")

2. function peaks()に適当な値を入力し，出力値を確認しなさい．引数の型に注意すること．

   julia> X,Y,Z = peaks() # 以下 julia>は省略

3. Zの最大値と最小値，及びそれらのインデックス(行番号，列番号)を求めて図に示しなさい．
   

   findmin(Z)
   findmax(Z)

   #以下はPlotsパッケージの使用を前提としたヒント
   using Plots
   Plots.pyplot() # or gr(), plotlyjs() など
   Plots.surface(x,y,Z) # 3D
   Plots.scatter!(xx,yy,maxval) # !マークは上書き等の破壊的なコマンドを表す

解答例（参考 †