よろしくお願いいたします。
Scilabを用いて二重振り子の数値解析を行なっています。
この先scilabが必要になるのでいろいろ使ってみているのですが、
うまくいきません。
オイラー法を用いて二重振り子の運動方程式を動かそうとしているのですが、
動きません。
ソースは
//Matrix
ccc1=zeros(20000,1);
ccc2=zeros(20000,1);
s1=zeros(20000,1);
s2=zeros(20000,1);
c1=zeros(20000,1);
c2=zeros(20000,1);
m1=0.10; //一重目質量
m2=0.08; //二重目質量
l1=0.4; //一本目長さ
l2=0.3; //二本目長さ
dt=0.05; //時間刻み幅
g=9.80665; //Gravitational acceleration
c1(1)=%pi;
c2(1)=%pi*0.5;
s1(1)=0; //1.
s2(1)=0; //2.
//Eular method
for t=0:1:100
z1=-m2*l1*s1(t+1)*s1(t+1)*sin(c1(t+1)-c2(t+1))*cos(c1(t+1)-c2(t+1));
z2=m2*g*cos(c1(t+1)-c2(t+1))*sin(c2(t+1));
z3=-m2*l2*s2(t+1)*s2(t+1)*sin(c1(t+1)-c2(t+1));
z4=(m1+m2)*g*sin(c1(t+1));
z5=l1*m1+l1*m2*sin(c1(t+1)-c2(t+1))*sin(c1(t+1)-c2(t+1));
ccc1=(z1+z2+z3+z4)/z5;
ccc2(t+2)=ccc2(t+1)+dt*ccc2(t+1);
s1(t+2)=s1(t+1)+dt*ccc1(t+1);
c1(t+2)=c1(t+1)+dt*s1(t+1);
z6=(sin(c1(t+1)-c2(t+1))*(m2*l2*s2(t+1)*s2(t+1)*cos(c1(t+1)-c2(t+1))));
z7=(sin(c1(t+1)-c2(t+1))*((m1+m2)*l1*s1(t+1)*s1(t+1)));
z8=(m1+m2)*g*cos(c1(t+1)-c2(t+1))*sin(c1(t+1));
z9=(m1+m2)*g*sin(c2(t+1));
z10=(l2*m1+l2*sin(c1(t+1)-c2(t+1))*sin(c1(t+1)-c2(t+1)));
ccc1=(z6+z7+z8+z9)/z10;
ccc2(t+2)=ccc2(t+1)+dt*ccc2(t+1);
s2(t+2)=s2(t+1)+dt*ccc2(t+1);
c2(t+2)=c2(t+1)+dt*s2(t+1);
end
なのですが、
ccc2(t+2)=ccc2(t+1)+dt*ccc2(t+1);
s1(t+2)=s1(t+1)+dt*ccc1(t+1);
c1(t+2)=c1(t+1)+dt*s1(t+1);
のあたりで、
「添字に誤りがあります」とエラーが出てしまいます。
20000行一列の行列を用意して、C言語のように配列を用いているのですが、、、。
ご教授よろしくお願いいたします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 数学 単振り子とルンゲ・タック法 1 2022/07/15 00:05
- Visual Basic(VBA) エクセルVBAで教えて頂きたいのですが? 2 2022/12/31 20:28
- C言語・C++・C# C言語で再起関数とポインタを用いて文字列反転をする方法がわかりません。 4 2023/04/29 20:32
- 経済学 今期の消費量をC1、来期の消費量をC2、消費者の効用関数をu=5 (C1)^6(C2)^4とする。 2 2023/01/29 10:07
- Excel(エクセル) エクセル関数教えてください 3 2022/06/21 10:22
- 数学 a1,a2, a3をベクトル空間Vのベクトルとする。a1+a2,a2+a3,a3+a1が一次独立のと 2 2022/10/02 15:55
- Excel(エクセル) エクセルの関数式を教えてください。 2 2022/11/29 21:09
- その他(Microsoft Office) IF関数について教えてください 2 2022/05/10 13:31
- 物理学 写真の問題について質問なのですが、解説を見る限りc1とc2の電気量は等しくないように思えます。 しか 4 2023/07/01 21:41
- Excel(エクセル) エクセルVBAでセルに表示されているとおりの数値を取得したい(時間の計算結果) 1 2022/03/30 17:52
関連するカテゴリからQ&Aを探す
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
既婚で現役のAV女優さんは居ま...
-
VBA ソートすると、1、11、...
-
相対参照から絶対参照に変換す...
-
別ブックの空白行に転記
-
ExcelVBAで指定文字(この場合...
-
CDレコの曲の消し方を教えてく...
-
YmobileからSoftbankに乗り換え...
-
射精したあとの匂いって他人に...
-
精液のにおいがほとんど無いの...
-
おっぱいを舐める
-
夫にセックスがないのなら他人...
-
1日3回セックスって多いですか...
-
彼女をオカズにして抜くのって...
-
先日彼氏とラブホに行ったら電...
-
彼とのエッチで、彼がイクのが...
-
男の精子ってどんな匂いですか、
-
手マンしたあと それとなく指の...
-
手マンした手って臭いですか?
-
手マンしたあと彼氏の指に私の...
-
彼氏の精液の匂いが無味無臭(...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
既婚で現役のAV女優さんは居ま...
-
VBA ソートすると、1、11、...
-
CDレコの曲の消し方を教えてく...
-
VBAが止まります。
-
EXCELで3行を一組にして結合す...
-
大昔から、クンニ、フェラって...
-
シンナーの夏型と冬型の違いは?
-
私は今年で60歳で孤独です。40...
-
別ブックの空白行に転記
-
女性が頼まれなくてもフェラす...
-
直線コネクタの中央にコネクタ...
-
エクセル最終行の下に貼り付け
-
相対参照から絶対参照に変換す...
-
データの平均を1分値にまとめる...
-
Word 黒塗り部分の文字のみ削除...
-
4次元について
-
ウォークマンa30についてです。...
-
Excelで抽出・連続印刷したいです
-
最適な組み合わせの自動計算
-
ExcelVBAで指定文字(この場合...
おすすめ情報