初めまして。fdtd法のサンプルコードが表示されません。

[小久保雄太]

初めまして。小久保と申します。

「はじめての音響数値シミュレーション　プログラミングガイド」を購入しまして、fdtd法のサンプルプログラムの中のsampleFDTD1.pyは普通に表示されたのですが、sampleFDTD2.pyとsampleFDTD3.pyが表示されません。何か他に入れなければならないものがありましたら、教えていただきたいです。

[ホシカズマ]

On Tuesday, September 24, 2019 at 11:30:03 AM UTC+9, ホシカズマ wrote:

初めまして，星と申します。

書籍通りのプログラムですと，matplotlibとOSの関係から，画像が表示されないという問題が発生してしまうことを確認しております。

sampleFDTD1.pyが動作するようでしたら，他のFDTDプログラムを動作させるためのインストールに不足はないと思います。

描画部分に若干の変更を加えたプログラムを

http://www.openacoustics.org/?page_id=127

にアップロードしましたので，こちらを試してもらえますでしょうか。

よろしくお願いいたします。

[小久保雄太]

星　様

ご連絡ありがとうございます。無事表示されました。

2019年9月24日火曜日 11時30分44秒 UTC+9 ホシカズマ:

[小久保雄太]

追加で質問したいのですが、サンプルプログラム中の75~85行目にある

A=dt/(Ro*dx)

Ux2[0,:]=(1./(1.+Z*A))*( (1.-Z*A)*Ux1[0,:]-2*A*P1[0,:])

Ux2[-1,:]=(1./(1.+Z*A))*( (1.-Z*A)*Ux1[-1,:]+2*A*P1[-1,:])

Uy2[:,0]=(1./(1.+Z*A))*( (1.-Z*A)*Uy1[:,0]-2*A*P1[:,0])

Uy2[:,-1]=(1./(1.+Z*A))*( (1.-Z*A)*Uy1[:,-1]+2*A*P1[:,-1])

はどういった意味があるのでしょうか？

2019年9月24日火曜日 11時30分44秒 UTC+9 ホシカズマ:

[小久保雄太]

FDTD3のプログラムです。

2019年12月2日月曜日 15時12分26秒 UTC+9 小久保雄太:

[Sato]

星さま

初めまして、佐藤と申します。

「はじめての音響数値シミュレーション」の書籍を購入したのですが、サンプルコードがうまく動かずにこちらに辿り着きました。

FTDT(program5-23)のアニメーション描画ができず、困っております。

この描画に変更を加えたというプログラムを再度ご提供いただけませんでしょうか？

なお、当方の環境はWindows10で、Python3.8を使用しています。

書籍通りのプログラムだと、figwindowが出て、すぐに閉じてしまいます。

下記のように変更を加えてみたところ、静止画にはなりましたが、アニメーションになりません。

    oldsurf = surf

    draw()

    #fig.show()

plt.show()

よろしくお願いいたします。

2019年9月24日火曜日 11時30分44秒 UTC+9 ホシカズマ:

[ホシカズマ]

佐藤様

星です。

お手数おかけしております。

以下はPython2のコードですが，描画部分を

surf = ax.plot_surface(xx,yy,P2,rstride=1, cstride=1,cmap=cm.jet)

ax.set_zlim3d(-1,1)

ax.collections.remove(oldsurf)

oldsurf=surf

pause(.01)

としてパラパラ漫画のように動作できることを確認しております。

参考いただけましたら幸いです。

現在，サーバの調子が悪く，最新版のファイルを以下に置いております。

http://sydney.arch.jcn.nihon-u.ac.jp/20200320.zip

合わせてご利用ください。

よろしくお願いいたします。

2020年5月8日金曜日 10時40分11秒 UTC+9 Sato: