T_NAKAの阿房ブログ

「自力で正規確率分布する乱数を作ってみる（２）、（３）」で掲題の "von Neumannの棄却法"を使ってますが、解りづらいので、自分自身のための備忘録として書いておきます。 簡単な例を考えましょう。 まず、0 から ｎ までの一様分布で、この分布の乱数を考えます。 ここでこの乱数で m (\(0\leq m\leq n\…

前記事のプログラムをちょっと直して標題のプログラムを求めます。 矢印で示した行の右辺を２乗に直しただけです。 結果：

「調和振動子モデルの固有関数を求めるプログラム」のプログラム部分を Python3 で書き直してみます。 結果：

「エルミート多項式を再帰的プログラムで実装してみる」のプログラム部分を Python3 で書き直してみます。 結果： （\(n=0\) の場合は基準点が \(+1\) であることに注意。） \(H_{0}(s)=1\) \(H_{1}(s)=2s\) \(H_{2}(s)=4s^{2}-2\) \(H…

前記事は実際のデータを使ったのですが、sklearn.datasets にあるものを利用しました。 実際の状況でが外部からデータを取り込むのが普通でしょう。それを考えて、ここでは同じデータ内容ですが、ファイルからデータを読み取って対処しています。忘れてしまうので備忘録として書いておきました。 まず、ネットからデータを貰って次のよ…

やはり、実際のデータを使って線形回帰をやってみました。 と言っても、ここでもアヤメのデータですが。。 import sklearn.datasets as datasets import matplotlib.pyplot as plt iris=datasets.load_iris() X=iris.data[:,[…

先週はFisdomで１件、AsukaAcademyで１件修了証をいただきました。 １．Python入門 2019 私は[MIT] Python入門と[MIT] コンピュータサイエンスとプログラミング入門でPythonを勉強しましたが、それでも知らないことがあることをこの講座で気づかされました。なので、受講して良かった…

先週は Fisdom の講座の修了証を頂きました。 １．未経験者のためのPHPとMySQLを用いたWebデータベースシステム入門 2019 これは、ある在庫管理用のWebを一から作っていくというものでした。 これ自体はエクセルを上手く使えばオフラインで実装できる程度のものですが、これを（ApacheとMySQLを搭…

先週は Fisdom の講座の修了証を頂きました。 １．Webで学ぶ JavaScript2019 　 Visual Studio Code を使ってJavaScriptを勉強する講座でした。このVisual Studio Code を今一つ理解していないので、アラームが出たりして躓きました。そういうときは、js を直…

先週は Fisdom の講座の修了証を頂きました。 １．Pythonで学ぶビジネスデータ分析入門 2019 具体的には、Deep Learningまでは行かない機械学習の勉強で、jupyter Notebook を使って Python というよりも、pandas や scikit-learn などの応用ですね。ちょっ…

今回は残りの「線形重回帰　演習２」をやってみます。 [演習問題]----------------------------------- 線形重回帰　演習２ 3. "03_iris.csv" を DataFrame として読み込み、品種名が "Iris-setosa"のもののみ取り出す。 4. 品種名が "Iris-s…

「JMOOCの新規理工系基礎講座」で紹介した「Pythonで学ぶビジネスデータ分析入門」を受講しています。７単元までありまして、現在６単元を勉強中です。この６単元と７単元のヴォリュームが多いのでなかなか進みませんが、その中で演習問題があったので、一応正解を見ないでやってみました。 [演習問題]-------------------…

"Introduction to Computer Science and Programming in Python" の Problem Set 1　には問題がＡ,Ｂ,Ｃ に分かれています。今回は Part C についてちょっと訳してみます。 [問題集1 パートC]---------------------------- 配…

"Introduction to Computer Science and Programming in Python" の Problem Set 1　には問題がＡ,Ｂ,Ｃ に分かれています。今回は Part B についてちょっと訳してみます。 [問題集1 パートB]---------------------------- 配…

"Introduction to Computer Science and Programming in Python" の Problem Set 1　には問題がＡ,Ｂ,Ｃ に分かれています。今回は Part A についてちょっと訳してみます。 [問題集1 パートA]---------------------------- 配…

Asuka Academy の 「[MIT] コンピュータサイエンスとPython入門 Part 1」 の修了証を頂きました。 これは 6.0001 Introduction to Computer Science and Programming in Python の前半を日本語化したもので、１コマ50分×６回で全部で5…

「[MIT] コンピュータサイエンスとPython入門 Part 1 」が始まってますが、「10問の修了認定テストにて（リトライなしで）80％以上の正解率」という修了条件がちょっと厳しく、修了証を貰えないかもしれません。まあ目的はプログラミング能力向上なのでそれでも良いと思っています。この講座は "Introduction to Com…

「PCアプリ環境が古くて授業についていけない、、、」で Python3 が上手く動きませんでした。ブルーバックス「入門者の Python」を入手したので、この本にしたがって Anaconda を再ロードして Spyder 環境で Python3 を少しずつ動かしています。 さて、以前に Python2 で書いていたプログラムがそのまま…

ちょっと簡単な練習をしてみます。参考データは「天皇の一覧」の神武天皇から孝明天皇の享年を使いました。つまり初代から121代の天皇のデータです。 まず、区間数などのパラメータを与えずにヒストグラムを描いてみました。 享年が100歳を上回る方々がいらっしゃるのは不自然ですね。これが正規分布とどれだけフィットするかを QQ…

まあせっかく『統計ソフト「Ｒ」超入門 』 を読んだので、「ダーウィンのデータ（３）」を使って分散分析をしてみました。結果がこの「ダーウィンのデータ（３）」と同じにならないと拙いでしょう。。それを確かめました。 まず、エクセルのデータを書き直しました。 これをＲにインポートします。Excel ファイルからインポートは簡…

先月、『統計ソフト「Ｒ」超入門 』というブルーバックスの新刊が出ていました。統計解析をするためにいろいろなソフトウェアツールがありますが、無償で一般性があるのはやはり「Ｒ」なんでしょうね。何年か前にちょっとやってみようと思っていたのですが、ダウンロードしただけで中々始められませんでした。コマンド入力とか、ちょっとやってられないと感じたか…

暗黒通信団の『"距離"のノート』という小冊子を眺めていますが、そこで掲題のp-ノルムというのが出てきます。 この本で、a1=1.4、a2=2.1、a3=5.7 の場合の p=1～7 の値が表になっていました。 手慰みのために、これを計算するためのプログラムを作ってみました。 まず、p-ノルムの定義です。 　 プログ…

これも数検１級の問題ですが、こういうのはどうやって解いていいものやら、、 問題： 正の整数 に対して 　 とする。このとき、 の中に異なる整数は何個含まれるか？ ただし、 は実数 に対して、 以下の最大の整数を表す。 考察： ちょっと計算すると、 　、　　、　 　、　　、　 …… ………

これも数検１級の問題ですが、これはそんなに難しくないのでプログラム化も考えてみました。 問題： ３つのベクトル 　 について、 を求めよ。 解： このレベルだと定義通り順番に計算した方が良いのか？ という公式を使った方が良いのか？迷うところです。とりあえずこの公式を使うと、 　 　　　 順番に…

先週は Fisdom の講座1つの修了証を頂きました。 １．誰でも教えられる！ プログラミング授業 政府が、小中学校で2020年からプログラミングを必修化しようとしていますが、これを教える立場の先生への教育というスタンスのようですね。 プログラミング教育の意義・目的の説明 → 授業の準備の方法 → プログラミングの…

光子０個、１個、２個、３個、…、に平均αの２乗個のポアソン分布の係数を掛けて線形結合式で分布状態を描けるプログラムを考えてみました。 これは、調和振動子モデルの固有関数を求めるプログラムと分布を描くプログラムを組み合わせたものです。αは実数を手入力することになります（例えば"1"ではなく"1.0"と入力）。 αが増えるにつれ波形がく…

Python でグラフを描く方法を忘れてしまったので、「自由落下時間における一般相対論とNewton力学の比較」を例にしておさらいをしてみました。 プログラム 結果グラフ

Lab 4: Control Flow with While　の　Problem 1 - Login security を読んでみます。 Problem 1 - Login security コンピュータサイエンスのセキュリティの重要な側面の1つは、ハッシュという概念である。テキストを取って、何とかそれを数値に変換する。 これは…

Lab 3: Control Flow with If-Else-Elif　の　Problem 2 - Cafe menu を読んでみます。 Problem 2 - Cafe menu 単純なカフェメニューを最初に表示し（下記の例を参照）、選択肢の番号を入力するようにユーザーに求め、適切な処置を印刷するか、または選択肢が有効では…