もっと高速な方法はあるでしょうか? 島村1243 2007/12/13 14:21:55 (修正1回) └続き 島村1243 2007/12/13 14:23:05 (修正1回) ├(1)について(処理手順の改良) 山中和義 2007/12/14 08:55:07 │├!もったいない、整理すれば、4倍速。 SECOND 2007/12/14 21:34:04 ││└確かに速いです。 島村1243 2007/12/14 22:33:52 ││ └どっちでしょう? SECOND 2007/12/14 22:50:11 (修正1回) │└20秒速くなりました。 島村1243 2007/12/14 22:06:46 └(3)について各点での流線 山中和義 2007/12/15 09:43:05 ├!等ポテンシャルを色の濃淡で表現する 山中和義 2007/12/15 11:03:55
Re: (3)について各点での流線 |
返事を書く ノートメニュー |
山中和義 <drdlxujciw> 2007/12/15 11:03:55 | |
!等ポテンシャルを色の濃淡で表現する
OPTION ARITHMETIC complex LET r=1 !導体球の半径 LET a1=4 !点電荷のx座標 LET a2=r^2/a1 !影像電荷のx座標 LET q1=1 !点電荷の大きさ LET q2=-r/a1*q1 !影像電荷の大きさ DEF f(z)=q1*LOG(z-a1)+q2*LOG(z-a2) !w=f(z)=φ+iψ 吹き出し(吸い込み) SET WINDOW -6,6,-6,6 !描画範囲 SET POINT STYLE 1 FOR x=-6 TO 6 STEP 0.02 !複素平面を走査する FOR y=-6 TO 6 STEP 0.02 WHEN EXCEPTION IN LET w=f(complex(x,y)) !LET c=MOD(Re(w)*1000,256)/256 !電位 ※1000 調整要 LET c=MOD(Im(w)*1000,256)/256 !電気力線 SET COLOR MIX(1) 0,0,c !青の濃淡 PLOT POINTS: x,y USE END WHEN NEXT y NEXT x DRAW grid !座標 SET AREA COLOR "red" !点電荷の位置に赤丸印を表示 DRAW disk WITH SCALE(0.1)*SHIFT(a1,0) SET AREA COLOR "green" !接地導体球の位置に緑丸印を表示 DRAW disk WITH SCALE(0.1)*SHIFT(a2,0) END |
└これを利用すると流れの方向がわかるので、... 山中和義 2007/12/15 20:11:13 └とても数理的で良いですね。 島村1243 2007/12/15 23:15:50 (修正1回) └接地導体球の作画に問題が有る様 島村1243 2007/12/16 07:37:08 └対策案できました。 島村1243 2007/12/16 10:44:04 (修正1回) └ポテンシャル関数f(x,y)なら、電位分布は等... 山中和義 2007/12/16 13:56:21 └完璧に目的達成です! 島村1243 2007/12/16 19:04:03 (修正1回) └プログラムの続き 島村1243 2007/12/16 19:05:45