min117の日記

初期desireもち。趣味Mac,メインFedora,仕事xp。

数学/統計

連続複利の (1+r)**2 をChatGPTに徹底的に説明させる

youtu.be youtu.be 時点tに投資家が持っているお金がW(0)とするとき、1年後に口座に入っているのはW(0)+W(0)r なので、W(0)で囲むと(1+r)W(0)になるのはわかるけど、これが複利になったときにt年後に W(t)が(1+r)**t * W(0) と表現されるのが分からない。説…

Python数学 y=a**x と y=b**x の交わるトコ(積分の面積部分)に色を付ける

Pythonで、このグラフの斜線に色つけられるかなー、とか。 youtu.be 動画のなかで出てきた式は違うよ、ってChatGPTさんは言うのだが… ま、いいや。コードしてみる。 y = a**x と y = -x**2 + b で囲まれた部分の面積 できた。 もうね、ChatGPTあれば何でもで…

これまで書いたソースコードをjpg画像にしてAmazon Photo の容量無制限に保存

srcをjpgに ChatGPTでやりたかったことがどんどん解決してく。 Linux で例えば python3 myKINRIsrc2jpg.py 'sample.py' とかして、sample.pyのソースコードを画像に落とす処理を実装できる?パイソンじゃなくても例えば bash myKINRIsrc2jpg.sh 'sample.py' …

macOSのPython3のGUI(Qt5)で数学的幾何学模様を描いてみたい(Kritaとか使わずに)

結論 うまくいった。 きっかけ・経緯 この記事おもしろそう。てかThinkITサイト最近お気に入り。 thinkit.co.jp ほいでmacOSでこれ試したい。 qiita.com しかし出勤の時間が来たようだ。 早出して早帰りして夜試そう。 大成功 とりあえず動いた(ChatGPTに聞…

Pythonからはじめる数学入門 で作った関数グラフで自然対数の底e(ネイピア数)を調べる

ネイピア数e をグラフで可視化したい このyoutubeに出てきたグラフをPython3で描きたいのよ www.youtube.com できた y = 2**x (y=2のx乗)とその接線 derivative(微分) ↑↓ ネイピア数はこの間にあるらしい (↑↓それぞれの直線がy軸と接したところの間) …

天空の城ラピュタを(大人になってから)見て思うこと

2024/9/3(火) ホーキング博士によれば、時間と空間は、実は同じものだそうだ。そして、時間と空間の両方を通り抜けられるものは、重力しかない。 飛行石は唯一、重力に逆らえるものだ。それは、空間の歪曲にも、時間の前後にも、影響されない。時空を越え…

Pythonからはじめる数学入門p120 ユーザーが入力した数式のグラフを描く

コレ前からやりたかった 入力した数式のグラフを描く。 $ python3 myKINRImath15.py y = 2x + 3 のグラフ y = xの4乗 -30かけるxの2乗 +10x のグラフ y = xの5乗 -30かけるxの3乗 +50x のグラフ こりゃ楽しい。Webからの入力インターフェースも作ってみるか…

コラッツ予想をPythonで

コラッツ予想 正の整数nに対して、偶数なら2で割る、奇数なら3倍して1を足す を繰り返すと 「必ず1で終わる(と思う)」というコラッツ予想(Collatz problem) 例えば 2 → 1 3 → 10 → 5 → 16 → 8 → 4 → 2 → 1 4 → 2 → 1 5 → 16 → 8 → 4 → 2 → 1 6 → 3 → 10 …

電圧電流チェッカー Makerfire USB Cメーターテスター を買ってみる

これ見たので www.youtube.com ポチってみる。 電圧電流チェッカー Makerfire USB Cメーターテスター Type-C マルチメーターテスター TFTカラーディスプレイ 電流 電圧 抵抗 温度 通電時間 容量AMP充電器負荷インピーダンスなど表示 QC2.0 3.0 usb電流計 (Bl…

高次元の関数を微分して「1つ下の次元の関数」を返す(微分電卓)Pythonからはじめる数学p200

Derivative Calculator 微分電卓 www.try-it.jp $ cat myKINRImath11.pyprint('Pythonからはじめる数学入門 20240825')print('p200 微分電卓 指定された関数yを微分して次元を下げた後の y\' を返す')print('-------------------------------') from sympy i…

Pythonからはじめる数学入門p107 級数を出力する(sympy使用)

級数とは 一定の法則に従って変化する数を、一定の順に(無限に)並べた数列の和 8個の級数 20個の級数 Pythonからはじめる数学入門 作者:Amit Saha オライリージャパン Amazon min117.hatenablog.com www.youtube.com $ cat myKINRImath09.py print('Python…

Pythonからはじめる数学入門p50 ニュートンの万有引力の法則

万有引力は距離の2乗に反比例する force = G * (m1 * m2) / (dist**2) # -*- coding: utf-8 -*- import pandas as pdimport pandas_datareader as web import matplotlib.pyplot as plt# %matplotlib inlineimport matplotlib.font_manager as fm# plt.rcPar…

Pythonからはじめる数学入門 p35 線形グラフを画像に保存する

線形グラフ Pythonからはじめる数学入門 作者:Amit Saha オライリージャパン Amazon 本のコードから改良したところ ・x軸のタイトルを自由に(変数kawacomで)変えられるようにした ・x軸に実行時間とスクリプト名を表示 ・x軸にもy軸にも、上部のタイトルに…

Python3のハッカーっぽい3D画像をgifアニメに

ランキング参加中Python ハピタス登録で1,000円分になるURL ハッカーっぽい〜♪ できた 前回のこの記事のコードを引数をとるように改善。 min117.hatenablog.com $ vim myKINRI3DArg.py デッパリの山の大きさを決める21行目を引数で指定できるようにした。 ba…

Python3でハッカーっぽい3Dの画像を保存する(多変量正規分布をplotしてpng画像に保存)

ランキング参加中プログラミング ハピタス登録で1,000円分になるURL ハッカーっぽい3Dクール画像 できた 理屈はよく分からんがとにかくこういう画像をPythonで描きたかった。 コードは写経 $ vim myKINRI3D.py 39行目だけ自分で足した。png画像に落とす処理…

オライリー本「プログラミング文体練習(Pythonで学ぶ40のコードスタイル)」

ランキング参加中プログラミング ハピタス登録で1,000円分になるURL 恋愛小説「高慢と偏見」Pride and Prejudice プロジェクト・グーテンベルクのサイトからDLできる。 www.gutenberg.org この恋愛小説の全単語を分析し、頻出単語上位25までを並べるだけのプ…

神?そんなもんいるわけないじゃん→いや、ちょっと待てよ…数学だけじゃ説明つかんことあるぞ…(カント・純粋理性批判+柄谷行人・トランスクリティーク)

イデア=4次元以上の人間が理解できない時空 なのか 神なんぞ信じないけど、数学だけじゃ説明つかんことがあるぞ… ja.wikipedia.org ja.wikipedia.org www.youtube.com www.youtube.com www.youtube.com そう。人間が認識できる知の範囲には限界あるよって言…

Pythonから始める数学入門 写経してJupyterで動かす(Python3)

本のとおりだと情報が古くてmatplotlibでグラフ描けないので、自分なりにコードいじったけど、いろいろ勉強になった。記録に保存。 ついでに英単語も覚える。 逆比例 inversely proportional 逆 invesely 比例 proportional 線形 linear P22非線形 nonlinear…

素数か判定するExcelプログラムを書く(VBA)

入力された数字が素数かどうかを判定したい。関数とか公理公式を使わないで独自で実装する。 素数とは… 1またはその数自身でしか割れない数。 例えば13。 1または13でしか割れない → 素数です。 例えば19。 1または19でしか割れない → 素数です。 じゃあ21は…

python3 カンマ区切りCSVファイル内の「金額列のカンマ区切り」のみ消し去る

例えばこんな財務データファイルがある。 irbank.net ファイル名はSheet1a.csvで保存。 売上高列に続く数字は "116,861,000,000円 "。 つまり、金額がカンマで区切られてしまっている。これだとCSVの各項目(フィールド)の区切りのカンマなのか金額のカンマ…

python3 のpandasで「上場企業名,当期利益,自己資本」のカンマ区切りcsvファイルを読み込む

前回記事に引き続き上場企業のデータを使ってcsvをいじる。 min117.hatenablog.com カンマ区切りのcsvファイルを読み込んで 「上場企業名,当期利益,自己資本」だけを取り出して カンマ区切りで出力(ファイル名はoutt5.csv)する。 $ cat japan-all-stock-fi…

awk (入力じゃなくて)出力の区切り文字を指定する(カンマ区切りで各フィールドを出力したい)

bash、awk、python3でcsvファイルを自由自在にしたい。都度調べるのがダルいので手に覚えさせる。 サンプルcsvはこれ。上場企業の当期利益と資本金を一覧にしてみる。 hesonogoma.com とりあえずアタマの50件だけダウンロード。 $ cat japan-all-stock-finan…

bash カンマ区切りのcsvを1列目でソートする sort -t, -k1n

bashでカンマ区切りのCSVを1列目でソートしたい。こんなcsv。 こうやる。 $ cat outt3.txt | sort -t, -k1n -t, 区切り文字はカンマ,だよ -k1 キーが1列目である -k1n キーが1列目である かつ それを数字(n)としてソートする genzouw.com こうなる。 末尾…

Python3 正規表現のグラフを描く

偏差が平均との差であること、 その差を二乗したもの(正負反転)が分散であること、 二乗を元に戻すためにルートを取ったものが標準偏差であること、 は分かった。 あとはPythonに実データを食わせて標準偏差をとって正規分布になるか確認してみたい。 まず…

sed ヘッダ行だけをタテに表示する (Excelの行列を入れ替えて貼り付け)

こんな感じのカンマ区切りのcsvファイルがあったとして $ vim test.csv 1列目(ヘッダ)だけをタテに表示したい。 sedでできた。こうなる。 $ cat test.csv | head -n 1 | sed -e 's/,/\'$'\n/g' 's/,/\'$'\n/g' ではsとgの前後にシングルクオートが要るらし…

csvの任意の列2つをタテヨコにとった散布図scatterを表示するだけのpythonコード

csvの任意の列2つをタテヨコにとったグラフ(散布図)を作るだけのコード。こういうシンプルな機能を揃えたい。 サンプルcsvはこんな感じ。銀行の残高をcsvに出力したもの。 $ cat dpwithdraw.csv コードはこれ。python3。 11〜15行目 csvの列名を取り出して…

matplotlibのグラフ凡例の豆腐文字化けをIPAフォントで解消する → 失敗

jupyterでcsv読み込んでグラフ化したら豆腐文字化けした。 このサイト参考に直す。 qiita.com 調査に使ったpythonコードは3つ。 import matplotlib matplotlib.matplotlib_fname() matplotlib.get_configdir() matplotlib.matplotlib_fname() で表示されるフ…

vim慣れてるとjupyterのpython3環境が楽しい

自分のiMacにはjupyter環境がインストールされてた。たぶん去年このへんの記事を作ったんだと思う(うろ覚え)。 qiita.com 前回使ったPaizaのjupyter環境からダウンロードしたipynbファイルを取り込んだら即使えた。 min117.hatenablog.com 取り込みは簡単…

jupyter notebook(Python3)で アンスコムの4つ組(Anscombe's quartet)を試す

前回記事で、WebさえあればPython3でグラフ描ける環境ができた。 min117.hatenablog.com Shellのコンソールだけだとグラフ描けなくて困ってたけど、これなら色々と試せる。 以前からいじりたかったコレ。 さっそく試してみる。 P87の散布図。 「アンスコムの…

PaizaCloudのjupyter notebook環境でPython3を利用する

統計をかじり始めた。「何よりまず手を動かす」したいので、すぐに動かせる環境に飛びつく。 PaizaCloudにアカウントを作ってみた。 paiza.cloud 「Webの開発環境構築を挫折したことはありませんか?」 そう。環境作るの時間かかってそれで満足しちゃうのが…