フォートランと電気回路に詳しい方急募！
■ このスレッドは過去ログ倉庫に格納されています
0001higantes ◆VxtSPywA
NGNGシミュレーション板にも書いたのですが、急いでますのでこちらにも書かせていただきます。
とある大学の工学部四年生です。
もう卒業研究をしなくてはならないのにその前の課題が終わりません。
ある回路における電流の値を求めるプログラムなのですが、
自分ではどこが間違ってるのかさっぱりわかりません。
http://www.geocities.co.jp/SiliconValley-SanJose/4145/
上にフォートランのソースと回路図を置きましたので詳しい方が知識を貸していただけると小生興奮の至りであります。
0002デフォルトの名無しさん
NGNG0004デフォルトの名無しさん
NGNG坊やだからさ
どっちも間違ってるよ
0005デフォルトの名無しさん
NGNG0006higantes ◆VxtSPywA
NGNGフォートランとかはやっぱシミュレーション板で正しいのでしょうか。
0007名無しさん
NGNG0009デフォルトの名無しさん
NGNG0011名無しさん
NGNGもう手遅れだと思われ。
0012デフォルトの名無しさん
NGNG0015名無しさん
NGNG0016デフォルトの名無しさん
NGNGそれっぽい講義を担当してる他の教授をあたれ。
0017デフォルトの名無しさん
NGNGヘタレな貴方のためになぜ私たちが時間と労を割かねばならないのですか？
このようなスレッドを立てて本当に誰か協力してくれると思っているのですか？
以上、現役留年生のたわごとでした。
0020デフォルトの名無しさん
NGNGどっちか言うとフォートランは他人が見て判り辛いからね。
ルンゲクッタの部分はいいから、
回路をどうプログラムにしたかを書いてみなさい。
0021デフォルトの名無しさん
NGNG0022higantes ◆VxtSPywA
NGNGはい。とりあえず書いてみます。
今後それぞれの素子についている数字を、その枝の数字とします。
R1,R2,R3のところに流れてる電流をX,Y,Zと置きます。
枝１，２，３、４，５，６，７，８，にかかる電圧をV1,V2,V3・・・V8とすると、
V1=V4+V6
V2=-V4+V5+V7
V3=-V5+V8
という式がでます。
続く
0023higantes ◆VxtSPywA
NGNGすると
V1=V4+V6
↓
L1(dX/dT)+R1*X=-R4*X+R4*Y-L6(dX/dT)-R6*X-L6(d It/dT)-R6*It
V2=-V4+V5+V7
↓
L2(dY/dT)+R2*Y=R4*X-R4*Y-R5*Y+R5*Z-L7(dY/dT)-R7*Y-L7(d It/dT)-R7*It
V3=-V4+V5+V7
↓
L3(dZ/dT)+R3*Z=R5*Y-R5*Z-L8(dZ/dT)-R8*Z-L8(d It/dT)-R8*It
という式がでます。括弧の中は微分です。
0024デフォルトの名無しさん
NGNG0025higantes ◆VxtSPywA
NGNG(dX/dT)=(-(R1+R4+R6)*X+R4*Y+IT*W*COS(W*T+THETA)*(-L6)-R6*IT*SIN(W*T+THETA))/(L1+L6)
(dY/dT)=(-(R2+R4+R5+R7)*Y+R4*X+R5*Z+IT*W*COS(W*T+THETA)*(-L7)-R7*IT*SIN(W*T+THETA))/(L2+L7)
(dZ/dT)=(-(R3+R5+R8)*Z+R5*Y+IT*W*COS(W*T+THETA)*(-L8)-R8*IT*SIN(W*T+THETA))/(L3+L8)
という上の式がでます。
さっき言い忘れましたけどＴというのが時間です。
電源の式が
IT*sin(W*T+THETA)で与えられているので、
(d It/dT)=IT*W*COS(W*T*THETA) になります。
Ｗは角速度です。ＴＨＥＴＡはＴ＝０の時の位相です。
0027higantes ◆VxtSPywA
NGNGhttp://www.geocities.co.jp/SiliconValley-SanJose/4145/circuit.jpg
0028higantes ◆VxtSPywA
NGNG俺だけど。
>>25の続き
この式をルンゲクッタの連立常微分方程式解法を使ってプログラムを組んでみました。
一応これがメインの話なんですけども。
002920
NGNGつまり
１）電流、電圧に注目して微分方程式をたてた
２）それをルンゲクッタで解いてみた
どちらが正しいか切り分けるために
○転がってるルンゲクッタのライブラリに置き換えてみて答えが出るかみる
○ルンゲクッタの部分にLCだけのパルス波過渡応答を与えて答えが出るか見る
いっそ、どうせ遅れついでだから行列にして汎用的に解いてみたら？
0030higantes ◆VxtSPywA
NGNG・ルンゲクッタの部分はまるまる本に書いてあるのを写しただけなので問題なく動作するはずです。
・ＬＣだけの回路と言うことは積分を使うと言うことですか？
・行列式は使わずに、ルンゲクッタを使って解け、という課題なのでそれはできなさそうです。
複素数使ってやればすぐにできる気がするんですけど・・・
003120
NGNGそれでルンゲクッタが正しく組めているなら、
最初に立てた方程式が正しくない訳だ
003220
NGNG003420
NGNGR4,R5が消えて、全部が同じ電流になれば 合ってると思うよ
003620
NGNG微分方程式の解き方が判らないないなら、
転がってるmathmaticaで
転がってなければ mathcadの評価版でもダウンロードしてね
003820
NGNGとにかく それで方程式が正しいようなら
ルンゲクッタで解いてる部分が正しくないんだから
で、それもどっかから拾ったなら、その拾った時の状態に出来るだけ近いように戻してみなさいな
全部正しい筈という事なら、単純に与えた条件に対してステップが荒すぎるんだろ
0039外野
NGNG訳わからんもん。
で、電流の初期値はどうやって決めるの？
WRITE(6,*) '----SHOKICHI----'
WRITE(6,*) 'I1,I2,I3'
READ(5,*) X,Y,Z
適当に与えると収束していくって話？
0040デフォルトの名無しさん
NGNG１よ！ お前は問題が解決したらスレ放り出して、それで満足か？
自分さえ満足すればOKなのか？
0041higantes ◆VxtSPywA
NGNGずっと徹夜なもので書き込む気力が無かったですごめんなさい。
なんか方程式うんぬん以前のバグが存在するらしいのでそれ探してるんですが全然みつかりません。
0042デフォルトの名無しさん
NGNG出来合いのシミュレータで解ける問題といても誰も誉めてくれんぞ
0043デフォルトの名無しさん
NGNGあのプログラムをそのままBASICに直して電流値をプロットしたら
位相と周波数が同じで振幅だけ異なる３つのサインカーブになったんだけど。
もしこれでいいんだったらロジックはあってるんじゃない？
■ このスレッドは過去ログ倉庫に格納されています