メインカテゴリーを選択しなおす
MATLAB,Python,Scilab,Julia比較 第5章【バックナンバー】
MATLAB,Python,Scilab,Julia比較するシリーズの第4章。第4章では分類問題で最終的にはニューラルネットワークや最適化アルゴリズムの話だった。第5章はフーリエ解析学から高速フーリエの話がメインとなる。
偶関数、奇関数を駆使する数学パズルを実施。 細かいことは置いておいて、雰囲気のみでざっくり解説。 奇関数が確定すれば0にできる。 偶関数が確定すれば線対称を利用して積分範囲を半分にした上で2倍にすればOK。
前回の数学パズルを真面目に解いてみる。 まずは平方根の関数の正体を探る。 偶関数、奇関数の特性を利用しまくって定積分を最適化しまくる。 ほとんどが0に消えて、半円の方程式だけが残る。 さらに偶関数の特性を利用して四分円にする。 半径2の円を四等分すれば答えが出る。
複雑な関数も無限次元ベクトルと見なすと力業で解くことが可能。 複雑な定積分を無限次元ベクトルとして表現。 これをプログラムとして解いていく。
MATLAB,Python,Scilab,Julia比較 第5章 その15【複雑な定積分⑤】
複雑な定積分をMATLABで求めた。 同様に円周率が答えとして算出。 小数点第6位まで一緒。 Nを増やせばもっと精度は上がる。
算額（その1040）八十三藤沢町保呂羽保呂羽山保呂羽神社明治26年(1893)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html直線の上に甲円，乙円が互いに接して載っている。その2つの円の上に丙円が載っている。丙円のてっぺんまでの高さはいかほどか。甲円の半径と中心座標をr1,(x1,r1)乙円の半径と中心座標をr2,(0,r2)丙円の半径と中心座標をr3,(x3,y3)とおき，以下の連立方程式を解く。高さはy3+r3である。include("julia-source.txt");usingSymPy@symsr1::positive,x1::positive,r2::positive,r3::positiv...算額（その1040）
算額（その1039）八十二藤沢町保呂羽保呂羽山保呂羽神社明治7年(1874)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html大円の中に中円2個，小円5個を容れる。小円の直径が1寸のとき，中円の直径はいかほどか。大円の半径と中心座標をR,(0,0)中円の半径と中心座標を2r3,(x32,0);x32=x3-3r3<0小円の半径と中心座標をr3,(x3,r3),x(3-3r3,0)とおき以下の連立方程式を解く。include("julia-source.txt");usingSymPy@symsR::positive,r3::positive,x3::positivex32=x3-3r3eq1=x32^2+r...算額（その1039）
算額（その1041）八十三藤沢町保呂羽保呂羽山保呂羽神社明治26年(1893)山村善夫：現存岩手の算額，昭和52年1月30日，熊谷印刷，盛岡市.http://www.wasan.jp/yamamura/yamamura.html直線の上に甲円，丙円，丁円が隣同士互いに接して載っている。さらに，丙円と丁円の上に乙円が載っている。乙円のてっぺんまでの高さはいかほどか。乙円の半径がわかれば，算額1040と同じ問題になる。甲円の半径と中心座標をr1,(x1,r1)乙円の半径と中心座標をr2,(x2,y2)丙円の半径と中心座標をr3,(0,r3)丁円の半径と中心座標をr4,(x4,r4)とおき，以下の連立方程式を解く。高さはy2+r2である。r1,r3,r4が既知なので，x1とx4は簡単に計算できる。ちなみに，「算...算額（その1041）
算額（その401）埼玉県所沢市下新井熊野神社明治6年(1873)山口正義：毛呂周辺の算額https://yamabukiwasan.sakura.ne.jp/22moroshuuhen.pdf大円内に正7角形と小円が入っている。大円と小円の直径を与えたときに正7角形の一辺の長さを求めよ。大円と小円の半径をr0,r1，正7角径が内接する円の半径をr2として，以下の方程式を解く。当時sin(π/7)とかsin(π/7)どうやって計算していたのだろうか。include("julia-source.txt");usingSymPy@symsr0::positive,r1::positive,r2::positive,x2::positive,y2::positive,x3::positive,y3::positi...算額（その401）
算額（その402）愛知県岡崎市六所神社(1779)http://www.st.nanzan-u.ac.jp/info/ma-thesis/2015/ss/m14ss006.pdf甲，乙，丙の正方形がある。その面積の和は61歩，甲乙の一辺の差は2寸，乙丙の一辺の差は1住んである時，丙の一辺の長さを求めよ。甲，乙，丙の一辺の長さをそれぞれ「甲」，「乙」，「丙」とし，以下の方程式を解く。usingSymPy@syms甲::positive,乙::positive,丙::positive;eq1=甲^2+乙^2+丙^2~61eq2=甲-乙~2eq3=乙-丙~1res=solve((eq1,eq2,eq3),(甲,乙,丙))1-elementVector{Tuple{Sym,Sym,Sym}}:(6,4,3)正方形...算額（その402）
算額（その403）東京都あきる野市二宮二宮神社寛政6年(1794)https://yamabukiwasan.sakura.ne.jp/page3.html鈎股弦の中に内接する全円と5個の等円がある。全円の直径が44寸，弦が110寸のとき，等円の直径はいかほどか。全円の半径と中心座標をr1,(r1,r1)等円の半径をr2として以下の連立方程式を解く。include("julia-source.txt");usingSymPy@syms鈎::positive,股::positive,弦::positive,r1::positive,r2::positive,x::positive,y::positive;eq1=鈎+股-弦~2r1eq2=鈎^2+股^2~弦^2eq3=弦~股-x+鈎-y+8r2eq4=r1/...算額（その403）
以前にRustでSpigotアルゴリズムをプログラミングして、円周率を計算しました。今回は、JuliaでSpigotアルゴリズム実装して、円周率を計算していきます。Rustでの実装の様子はこちらの記事をご覧ください。今回実装するソースについ
番外編でRustで作ったライフゲームをJuliaに移植してみました。まぁ、Juliaらしい書き方ができてるかどうかは別として、まぁ、Juliaで書いたって感じですrust版は過去の記事をご覧くださいJuliaで作るライフゲームでは、やってい
以前、 Julia を勉強しようと買った本（珍しく、印刷された本を買いました） [1] が積読（つんどく）状態になっていたのですが、先日、JupyterLab で R を使えるようにしたのを機に、Julia も使えるようにしました。Julia をインストールした後、P...
Juliaでニューラルネットワークを使ってタイタニック生存者予測をする
長いタイトルになりました。今回は、JuliaのFluxを使ってニューラルネットワークを計算して、タイタニック号の生存者予測をしていきます。Fluxを使ってみたかった、、ぐらいの話なので、精度はとか。。パラメーターは。。。とかはあまり突っ込ま
juliaでグラフを描画する方法を解説してきます。今回は応用編1ですCSVやDataFramesなどを使ってタイタニックのデータを加工して、グラフに描画していきましょうPlotsのパッケージのインストールやインポートの方法はJuliaでグラ
新しい言語や環境ができると、2バイト文字の日本語の壁にぶち当たるのは世の常。けど、最近はよくよく対応できるようになってるのもまた現実。というわけで、今回は、julia Plotsで日本語を扱う方法を2通り解説していきます。実行環境場合win
