このブログは、旧・はてなダイアリー「檜山正幸のキマイラ飼育記 メモ編」(http://d.hatena.ne.jp/m-hiyama-memo/)のデータを移行・保存したものであり、今後(2019年1月以降)更新の予定はありません。

今後の更新は、新しいブログ http://m-hiyama-memo.hatenablog.com/ で行います。

2017-06-01から1ヶ月間の記事一覧

3次元方向記述とその実現

銃座付きのジープ(簡易combat/fighting vehicle、TechnicalsまたはImprovised fighting vehicleと呼ぶ)。ただし、レプリカ。*1より本格的な銃座:*2天文関係なら:*3*4地球上なら:ゴー☆ジャス*5 *1:画像:http://www.gizmodo.jp/2013/09/jeep.html より *…

姿勢記述の色々

次の論文によると、 http://publications.lib.chalmers.se/records/fulltext/160346.pdf 自動車の記述座標系は、 ISO8855:1991 で決まっているらしい。 ISO8855:1991 Road vehicles -- Vehicle dynamics and road-holding ability -- Vocabularyhttps://www.…

力学用語

mechanics 力学、語源的には機械学だろうが、かなり広い範囲の意味 kinematics 運動学。運動記述学というのがふさわしそう。 statics 静力学 dynamics 動力学 分類基準はあまり分からない。静力学は時間を扱わなくて、動力学は時間を扱うのだろう、たぶん。…

maximaとの通信方法、ソケットが既に使われているときの挙動

なんとmaximaはネットワーククライアントになる。サーバーじゃなくてクライアント!ここが既にハマリ所だが、使用されているソケットにクライアントとして接続しに行くと、次のように死ぬ。 bin > .\maxima.bat -s 4000 Connecting Maxima to server on port…

scilabでは困ること

scilabでは、 --> A = [sqrt(3) -1; 1 1] A = 1.7320508 -1. 1. 1. --> A*A ans = 2. -2.7320508 2.7320508 0. これでは不便なこともある。 --> a 未定義の変数: a --> A = [a -1; 1 1 ] 未定義の変数: a

wxmxフォーマットと数学ノートブック

SageMath notebook, Jupyter Notebook などノートブック形式が主流(?)らしい。Isabelle(イザベル)もproof documentを指向していた。Maximaもノートブック形式指向らしいが、wxMaximaドキュメント形式は、まだセッション記録でしかなくて、とてもドキュ…

平方根の計算

√2 * √3 = √6 はそのままではできない。rootscontract関数に渡す必要がある。 rootscontract(sqrt(2)*sqrt(3)); しかし、√2 * √2 = 2 は自動的に計算される。定義によるからだと思うが、メカニズムはよく分かってない。

Lispとのシンボル対応

マニュアルによると: Maxima Lisp foo $foo 名詞のfoo %foo ?foo foo しかし、実際にやってみると、 :lisp (setq foo 'bar) :lisp foo /* => BAR */ Lisp側での名前が全部大文字の名前になってしまう。なんだ??いずれにしても、:lisp $名前 は探索に役に…

初期化ファイルと出力の準備

基本的なディレクトリは、maxima_userdir と maxima_tempdir 変数に入っている。ファイルの検索位置(ディレクトリ)は、file_search_maxima, file_search_lisp変数に入っている。C:/Users/hiyama/maximaディレクトリがfile_search_maxima変数に入っているの…

wxMaximaのワークシートと設定

wxMaximaがどうも使いにくいのは、コマンドラインシェルではなくて、セルの集まりのワークシートをモデルにしていたのか。コマンドラインとして入力するのではなくて、セルの内容を入力して、セルが評価されるというスプレッドシート方式だったのか。操作モ…

描画の各軸のスケール/描画出力方法

プロットオプション ['same_xy, true], ['same_xyz, true] でスケールが揃う。 wxplot2d([parametric, cos(t), sin(t)], [t, 0, 2*%pi] , [same_xy, true]) ; contour_plot(1/(1+x^2+y^2),[x,-3,3],[y,-3,3], [same_xy, true]); plot3d(10/(1+x^2+y^2),[x,-3…

初期化ファイルはtcl構文だ!

maximaのdotfile - 檜山正幸のキマイラ飼育記 メモ編 の件、あのファイルはtcl構文だろう。 http://wiki.tcl.tk/724 http://www.interq.or.jp/japan/s-imai/tcltk/array.html tcl/tkベースのmaximaのためだとすると、wxMaximaには関係ないだろう。maxima-ini…

ワケわからんバグ

以下の1番のみエラー。 block( [div], 1 + 1) block( [di], 1 + 1) block( [divv], 1 + 1) block( [div_], 1 + 1) 意味不明なバグだ。変数の名前によって事情が変わる。そもそも、名前(記号アトム)は局所変数になるはずのものなのに。ソース c:/Installed/…

ファイルのありか、ラムダ式、ブロック式

.lisp, .macファイルのありか $MAXIMA_HOME/share/maxima/5.40.0/src/ の下 $MAXIMA_HOME/share/maxima/5.40.0/share/ の下 ラムダ式 lambda([ローカル変数リスト], 処理本体) block式 Lispのprogn形式やRのlocal式のようなもの。入れ子のブロック構造を作る…

グラフの描画

1パラメータの3dカーブ - 檜山正幸のキマイラ飼育記 メモ編 の敷衍。 1変数→1変数 plot2d(関数式, [変数名, 下限, 上限]) e.g. plot2d(x^2, [x, -2, 2]) 1変数→1変数 共端である複数の関数(Set射)を同時に plot2d([式1, 式2], [変数, 下限, 上限]) e.g. pl…

構文解析と束縛力定義

$MAXIMA_HOME/share/maxima/5.40.0/src/nparse.lisp

ダメなところ

パッケージマネージャがない。http://maxima.sourceforge.net/3rdpartycode.html にサードパーティコードのリストがある。 タブ補完が使えない。とーーっっても不便。→ Ctrl+K があるが、、、 DELキーが上に遡って消す、勘弁してくれ。 ヘルプシステムがカス…

OSインターフェイスとロード

MaximaのOS、特にファイルシステム機能、 意外に検索はしにくい。「maxima ls chdir」とかやってミソ。 http://maxima.sourceforge.net/docs/manual/maxima_74.html によると次が使える。load("operatingsystem"); を忘れずに。 chdir (dir) mkdir (dir) rmd…

変数のマネジメント

values() は現在具体的な実体を割り当てられ定義された文字の一覧を返す。 functions() は現在具体的な定義を持つ関数の一覧を返す。 変数の開放 kill(変数名)は変数名に割り当てられた定義を消す。 kill(all)はすべての割り当てられた定義を消す。 変数はた…

基本操作

等号の判定には、is(式 = 式), e.g. is( 1 + 1 = 2) 方程式を解くには、solve(方程式, 未知数変数), e.g. solve(x^2 - 1 = 0, x) 微分するには、diff(式, 変数), または diff(式, 変数, 微分回数), e.g. diff(x^2, x), diff(x^2, x, 2) tex(式) でtexのテキ…

スペースファイターと銃座

ミレニアムファルコン号による空中レーザー戦 vs タイファイター - 檜山正幸のキマイラ飼育記 メモ編の続き。見るべき動画は: Star Wars The Force Awakens Millenium Falcon Scene https://www.youtube.com/watch?v=8sarFZJl3h0 STAR WARS 1-7 ~ SPACE SCE…

微分できない、暗黙のラムダ変数がない

diff(x^2, x) 変数xがバインドされていると変なことになる。xを純粋な変数にしたいなら kill(x) しないといけない。「関数表示だから、暗黙にラムダ変数だろう」という常識が効かないで、大域変数として値を探してしまうのだ。

デカルト微分圏の定義

定義を探すことが多いので、まとめておくのがよいな。クラットウェル(https://arxiv.org/pdf/1208.4070.pdf)の一般化デカルト微分圏。射の結合は、図式順で併置。射影はπ0とπ1。 A generalized Cartesian differential category consists of a Cartesian c…

なめらか代数理論

C∞代数をローヴィエル理論と見ている。面白い。 http://www.kurims.kyoto-u.ac.jp/EMIS/journals/CMUC/pdf/cmuc9603/michor.pdf

多様体が定義できる圏

参考:「包含付き圏:対象を集合っぽく扱うために」http://d.hatena.ne.jp/m-hiyama/20101102/1288670139 category with inclusionsにさらに条件を足す。 そのinclusionをopen inclusionと呼ぶ。 他に、a family of models を指定する。 category with open-…

多様体の定義

比較的標準的 Manifolds, sheaves, and cohomology Torsten Wedhorna https://www2.math.uni-paderborn.de/fileadmin/Mathematik/People/wedhorn/Lehre/SkriptMannigfaltigkeiten.pdf プレ多様体の定義 Let α ∈ N ∪ {∞}. A locally R-ringed space (X, OX) i…

球の具体的なチャートとアトラス

球を地球に例える。グリニッジ子午線(本初子午線)と赤道の交点である緯度0経度0の地点を仮に経緯原点と呼ぶ。球面上の点の名前 n : 北極 s : 南極 e0 : 経緯原点 アフリカ・ギニア湾 e1 : 東経90度 インド洋 e2 : 東経180度 太平洋・フィジー諸島あた…

ラグランジュ定式化とオイラー定式化

流体力学で、ラグランジュ流の定式化とオイラー流の定式化があるらしい。 ラグランジュ 流れに乗る粒子を、時間に沿って追跡する オイラー 特定の点の流速などの変化を見る。 オイラー定式化が「場」を使った定式化となっていると言われるが、ラグランジュだ…

仔牛肉を使ったキジ料理

ユ(Yoo)の論文 http://www.math.northwestern.edu/~bwill/conformal/Phil1.pdf からまた引用(孫引き)する。マレー・ゲル=マン(Murray Gell-Mann)の言葉だそうだ。 it is like the French cuisine that one cooks pheasant between slices of veal, di…

1パラメータの3dカーブ

2パラメータの3dサーフェイスが基本で、特にカーブ用の関数は用意されてない。2パラメータのうち片方をダミーにして使う。 plot3d([cos(t), sin(t), (1/3)*t], [t, -4*%pi, 4*%pi], [k,0, 1], [grid,500,1])ここで、kはダミー変数で無意味。tだけがパラメー…