- 締切済み
scilabについて質問
scilabについて質問です。 ランダムウォークの二次元plotをするにあたって、 初期位置を原点として、ステップ幅（S）は原点からの距離dに対して、 S=e(-d/μ) に比例する。つまり、原点からの距離が遠ざかるとステップ幅が小さくなるようにdを決めたいのですが、 原点をgrandで決めたとき、dをどのようにscilabで書けばいいのかがわかりません。 どなたかご教示願います。
- none99jp
- お礼率62% (5/8)
- フリーウェア・フリーソフト
- 回答数1
- ありがとう数0
- みんなの回答 （1）
- 専門家の回答
みんなの回答
- hot-IC
- ベストアンサー率87% (7/8)
関連するQ&A
- ランダムウォークのような問題です。教えて頂けませんでしょうか。
ランダムウォークのような問題です。教えて頂けませんでしょうか。 1次元ランダムウォークの少し変形バージョンのようなものだと思うのですが、 以下の試行を十分大きい回数行ったあとの、粒子の確率分布について どなたか教えて頂けませんでしょうか。 (解法、物理などで似たような命題(ブラウン運動だよ、のような)、もしくは参考書など) 試行:============================================================ 初期位置x=x0にいる1次元粒子を考える。 粒子は位置xにいるとき、ガウス分布f(x)=a*exp(-(x-b)^2/c)の最大ステップ幅で、しかも φを0<φ< 2πの一様分布の乱数だとしたとき、ステップ幅 f(x)*sinφ で移動する。 ================================================================ これを十分な回数繰り返します。 x0, a, b, cの値をいろいろ変えたときの挙動が知りたいと考えています。 aが小さいときにはあまり初期位置から動きませんので、x0を中心としたunipolarな分布で 一方aが大きいときにはbipolarな分布になるのかな...とまではプログラムや直感的方法で 定性的には考えたのですが、それ以上が詰められません。 どなたか分かる方、よろしくお願い致します。
- ベストアンサー
- 物理学
- 中心極限定理の高次元化
一次元シンプルランダムウォークSnはnが十分大きいとき、その確率分布が一次元の正規分布に収束します。 では一般のD次元シンプルランダムウォークSnもnが 十分大きいときにD次元の正規分布に収束するのでしょうか？ ともあれ一般次元のシンプルランダムウォークSnのnが十分大きいときに収束する確率分布をご存知の方は教えてください。
- ベストアンサー
- 数学・算数
- 単振動の運動方程式 Scilabについて
摩擦のない床の上にある物体がばねで壁につながれているとします。このとき物体に働く力はばねの弾性力だけということになります。 単振動の運動方程式について考えています。 これについて、以下のようにScilabを用いて考えます。 いかのScilabでの考察が分からなくて困っています。 // 重りの質量と、ばねのばね定数、最初に引っ張る長さ m = 1; k = 1; L = 1; // 初期値 x0 = [L, 0]' ; // [ ]' は転置 t0 = 0; // 連立微分方程式の係数行列 w = k/m; M = [ [0 1] ; [-w 0] ] ; // 導関数の定義 deff("xdot = df(t,x)", "xdot=M*x"); // x=[x1, x2]'となっている // 数値解の計算 t = t0: 0.01: 5; // t0-5秒までを0.01区切りで y = ode(x0, t0, t, df); // yも行列なので、物体の位置は1行目の要素をとればよい plot2d(t, y(1,:), 5); // 5 は色を指定する数 xgrid(2); 以上がScilabでの考察なのですが、「初期値」とは何なのでしょうか?? 「連立微分方程式」はどの関係を式であらわしているのでしょうか?? 「導関数」はどんな式ですか??読み取れません…また、どのようにその式が導かれるのかも知りたいです。（これが一番知りたいです!!） 回答よろしくお願いします。
- ベストアンサー
- 物理学
- 一次元のランダムウォーク
下の問が解けなくて困っています。解説をして頂けると非常にありがたいです。お願いしますm(__)m （問）一次元のランダムウォークする粒子を考える。離散的な時間（ｔ＝0,1,2,...）を考え、時刻ｔ＝0では原点にいるとする。この粒子の各時刻での変位は確率pで１、確率p-1で-1とする。時刻tでの粒子の位置ｘ（ｔ）の期待値、分散の値を一般のpで計算しなさい。
- 締切済み
- 数学・算数
- ランダムウォーク
ランダムウォークの問題 ある参考書の問題をやっていた所、わからなくて解答解説を見ようと思ったのですが載っていませんでした。まわりに聞ける人もいないので教えて頂けないでしょうか。よろしくお願いします。 問題 対称なランダムウォークで (1) Ｐ｛Ｓ1≧０、Ｓ２≧０、・・・、Ｓ２ｎ≧０ ｜ Ｓ０＝０｝＝（１/2）*U2n (2) Ｐ｛Ｓ1＞０、Ｓ２＞０、・・・、Ｓ２ｎ＞０ ｜ Ｓ０＝０｝＝U2n ちなみにSnはn歩目の状態でU2nは2n歩目に原点に戻る確率つまりＰ｛S2n=0|S0=0｝のことです。 そしてその確率の値はU2n=2nCn*1/2^(2n)です。（Cはコンビネーションです）
- 締切済み
- 数学・算数
- 2次元単純ランダムウォークは原点にいつ帰ってくるか？
原点スタートの2次元単純ランダムウォークを考えます。すなわちX_iを独立同分布でそれぞれ確率1/4で(±1,0),(0,±1)のいずれかのベクトルを取る確率ベクトルとし、S_n=X_1+X_2+…+X_nとします。S_nを2次元単純ランダムウォークと呼ぶことにします。n>0のとき初めてS_n=0となったとするとき、H=nとおきます。Hを原点への最初再帰時刻（あるいは最初到達時刻）と呼ぶことにします。P(H<∞)=1、E(H)=∞がよく知られています。つまりいつかは必ず帰っては来るものの、その期待値は∞だということです。これらのことは大抵の確率論の本には書いています。そこで最初再帰時刻Hの分布を知りたいと思いました。時刻2nにしか戻らないので、P(H=2n)はいくらか？という問題といってもいいです。 このことは組み合わせ論の問題だと考えられます。つまりZ^2格子上で原点からスタートした長さ2nの道で、最初と最後以外に原点を通らないものは何通りあるか？という問題と同値です。もし途中で原点を通ってもよいのであれば、その道の本数は4項分布を用いて簡単に書き下せます。つまりi回右に、i回左に、j回上に、j回下に動いたとして、i+j=nであればよいので、Σ_{i+j=n}(2n;i,i,j,j)(1/4)^{2n}とかけるわけです。ただし(2n;i,i,j,j)=(2n)!/{i!i!j!j!}は4項係数とします。2n回目で初めて原点に戻ってくるというのは、2n回目に原点にいるような道から、2k回目に初めて原点に戻ってきて、その後2(n-k)回で原点に戻ってくる（この間は何回戻ってもよい）ような道の数をk=1～n-1まで引いてやればいいですから、p_n:=P(H=2n)の満たす漸化式を導くことは容易です。しかしそれが4項係数の和を含み、またp_nを決めるのにp_1～p_{n-1}の情報が全部いるという難解なものです。うまい数え方を見つけてp_nを簡単に表すよい方法はないものでしょうか？ちなみに1次元では簡単に解けましたが、同じ方法では出来そうもありませんでした。
- 締切済み
- 数学・算数
- scilabでのグラフ描画時の軸設定
scilabで2次元グラフを描くときの軸の設定 scilabのカラー描画であるMatplotを使ってグラフを描画しているのですが、縦軸と横軸の交点が思い通りに設定できません。 また、グラフと軸がずれて配置されてしまいます。 現在このようなプログラムでグラフを描画しています。 u=[1:1:5;1:1:5;1:1:5;1:1:5;1:1:5]; umin=min(u(:)); umax=max(u(:)); u=255*(u-umin)/(umax-umin); Matplot(u);//行列プロット set(gcf(),'color_map',jetcolormap(256)); colorbar(umin,umax);//color barの表示 g=gca(); g.axes_visible='off'; g.children.data=g.children.data($:-1:1,:);//原点を左下にする xtics=0:1:4;//x軸目盛の設定 drawaxis(x=xtics,y=0.5,dir='d',tics='v'); g.children(1).tics_labels=string(-2:1:2); ytics=0:1:4;//y軸目盛りの設定 drawaxis(x=0.5,y=ytics,dir='l',tics='v'); g.children(1).tics_labels=string(1:1:5); このプログラムを実行すると、縦軸が0.5で、横軸が-1.5のところで交わってしまいます。 また、カラー表示されたグラフが縦軸・横軸共に1.5くらいずつずれてしまいます。 ちなみにこのプログラムは以下のURLのサイトのプログラムを参考にしたものです。 http://vision.kuee.kyoto-u.ac.jp/lecture/dsp/?menu=mat_figure 初心者の質問で申し訳ありませんが、かなり困っています。ご指導、アドバイスよろしくお願いいたします。
- 締切済み
- その他（プログラミング・開発）
- 解析力学のラグランジアンの対称性について
「２つの粒子（質量m1、m2）が３次元空間内にあり、原点から距離に比例する中心力（比例定数は共通のk）を受けて運動しており、重力は無視する。」 (1)このときのラグランジアンの対称性は？ (2)２つの粒子をばねでつないだ場合対称性はどう変わるか？」 という問題なのですが、とくに(2)がどうすれば解けるのかわかりません。 解き方を教えてください！よろしくお願いします。
- 締切済み
- 物理学
- 同次座標 比例 記法
同次座標というのがあります。 たとえば、３次元空間で、原点からの方向を考えたいとします。 ３次元空間の各点に対して、原点からの位置ベクトルを考えたとき、 (1,2,3)や(2,4,6)や(-3,-6,-9)などを同一視したい。 このとき、比例を用いて、 1:2:3=2:4:6=-3:-6:-9 とかけるので、一般に、 x:y:z=1:2:3 (ある実数kにおいて、x=k,y=2k,z=3k) と表せます。 これは「:」という記号のすばらしさだと思います。 ところで、それは３つの数の組に同じ数をかけたものを等しいとみなしました。 今度は、３つの数の組に同じ数を足したものを等しいとみなしたいとします。 (1,2,3)に対して、(2,3,4)や(-1,0,1)などを同一視したい。 そういった座標は一般には知られていますか? 便利な記法はありますか？ その他、なんでも情報があればいただきたいです。
- ベストアンサー
- 数学・算数
補足
ご回答ありがとうございます。 S=e^(-d/μ) の間違いでした。すいません。m(__)m eは自然対数の底です。S,dは上記のとおりで、μは定数です。 つまり、原点から遠ざかり、dが大きくなるとステップ幅（S）が小さくなるんです。それで、dをプログラム上どのように書けばいいのかわからないんです。 すごく基本的なことだと思うのですが、基本的なことが分かっていないため困っています。 宿題なので丸投げはできないのですが、いくら考えてもわかりません。ヒントだけでもいただければありがたいのですが。宿題という理由から、ソースの公開はご容赦ください。