プログラミング覚書（Julia）

基本操作（VSCode使用を前提）

packageをインストール "]" + "add xxx" 

Alt + Enter (全実行)

Shift + Enter (ラインを実行して次のラインにカーソルを移す。RのCtrl +Enter)

Ctrl + Enter (ラインを実行)

comment/uncomment Ctrl + /  あるいは #=  ...  =#で囲む

REPL(Matlabのcommand windowに相当)に出力

println()

背景テーマはDark Visualstudioが自分は見やすい

別ファイルに作った自作のfunctionを呼ぶ：

include(file name)

スクリプト内に簡単なfunctionをつくる：

function eslowe(tc)

    es = 6.107799961+tc*(4.436518521*(10^-1)+tc*(1.428945805*(10^-2)

    +tc*(2.650648471*(10^-4)+tc*(3.031240396*(10^-6)+tc*(2.034080948*(10^-8)+tc

    *6.136820929*(10^-11))))))

    return es

end

for loop内で新しく定義された変数は、loopの外（globalかlocal変数）で定義しておかないと内部変数として外部からは認識されない。

コードを印刷

PrintCodeをインストール >  F1

出力をREPLかスクリプト内に出力させるかなどを変更

File > Preferences >  Settings > Extensions > Julia > Result Type

Arrayの計算

配列の数が同じ場合

A = rand(5)

B = rand(5)

C = A + B

配列の数が異なる場合(Broadcasting)や、配列の除算、乗算のときは .　を使用

A = rand(5)

B = 3

C = A .+ B

配列の抜き出しに際する注意点

上記の変数Gについて、G[1:2]、G[1 : 2]、G[1 :  2]ではエラーがでないが、G[1  : 2]ではエラーとなる。

配列の条件検索Matlabと同じで .& あるいは .| (スカラーの比較 には&& や||を用いてもよい e.g., i > 1 && j > 1）

データ操作

NaNを含む平均

mean(filter(!isnan, A))

最大最小

Aの中の最大最小

maximum(A), minimum(A)　

大きい方、小さい方

max(A,  B), min(A, B) 

条件指定は各条件を括弧で

findA =

    (ttimenum .> T - 15 * 60 * 1e3) .&

    (ttimenum .<= T + 15 * 60 * 1e3)

検索には

T[.findA]

missing データを扱う場合

missingの配列  f = ismissing.(A)

missingでない配列 (上のコマンドに続けて)  .!f

配列を抽出するとき（Matlabのfindに対応）findall(ismissing,  A) 

平均操作 mean(skipmissing(A))

単調増加（減少)のベクトルを作成

A = collect(1 : 2 : 10)

内挿（例は変数'albedo'の欠損値を内挿で補完）

using Interpolations

     itp = LinearInterpolation(timenum[isfinite.(albedo)],albedo[isfinite.(albedo)],extrapolation_bc=Linear())

     albedo[isnan.(albedo)]=itp(timenum[isnan.(albedo)])

ループ

for i = 0:dt:tmax

end

時間データ

1秒以下のデータを読むには

df = DateFormat("y-m-d H:M:S.s")

timenum = DateTime.(TTIME, df)

Matlabのdatenumのように*時間データを数値に変換する（平均操作や時間の丸め込みなどに利用）

timenum_value = Dates.value.(timenum)

*ただしMatlabの1が24時間に対し、下記のttimenumの1は1ミリセカンド

上記のttimenum_valueを一般時間データに変換しなおす

st = DateTime("0/12/31 0:00:00")

ttimenum = st .+ Millisecond.(ttimenum_value)

特殊な操作

デバッグ

   (1) 一度全実行 (ここで基本的なエラーを取り除く必要があるのが手間)

   (2) Breakpointを設定

   (3)  @run function name

for loop内で新しく定義された変数は、loopの外（globalかlocal変数）で定義しておかないと内部変数として外部からは認識されない（のでデバッグが難しい）。

変数のcopyについて

A=[1,2,3] 

B=A

A[1]=2

println(B[1])とすると2が出力される。これを防ぐにはA = copy(B)とする。しかし、左辺で特定の領域を参照しない場合はB=Aでもよい。上記の例で、A=2としてもBは変わらない。