仮想計算機構

IT業界と無縁な派遣社員のブログ

数学

極限としてデルタ関数が得られるもの

原島鮮『初等量子力学』10.2 より 連続関数がデルタ関数へ接近していく様子を可視化します。 ツールは SageShell を使います。 その1 sage: l = [plot(sin(a*x)/pi/x,(x,-15,15))+text("$a={}$".format(a),(15*0.7,a/pi*0.8),fontsize=20, bounding_box={'b…

楕円の速度ベクトル、加速度ベクトル

楕円 の速度ベクトル、加速度ベクトルのアニメーションを作る。可視化に必要な関数を下記で定義する。 sage: t = var('t') sage: x = 9*cos(t) sage: y = 4*sin(t) sage: vx,vy = x.derivative(),y.derivative() sage: ax,ay = vx.derivative(),vy.derivativ…

Canvasに描いた線をフーリエ級数で表現する

フーリエ級数を用いて絵を描くという試みをいろんなところで見かけます。 一番有名なのは以下の動画でしょうか。 初音ミクを描いちゃったという方もいらっしゃいます。 Bad Apple Curve!! 〜フーリエ記述子でアニメーション作ってみた〜車輪の再発明でいいじ…

直進形汎トロコイドの描画

以前、Schemeで描画するときに座標の取り方を指定できるようにしました。 【奥村本】座標の取り方を変えて線を描画する - 仮想計算機構これを利用して色々な線を描画していきます。今回は曲線!ってことで以下を参考にしました。 周期関数による幾何学模様の…

SLIB(Scheme)で集合を使ってみる

以下の環境でSLIBを使っていきます。 CPU : Intel Core i5 2.60GHz OS : WIndows 10 言語 : Gauche 0.9.9 準備 定義 基本操作 adjoin union intersection subset? member-if もう少し高度な操作 some every notany notevery おわりに 準備 今回のテーマは集…