異方性の大きいBi系超伝導体(Bi-2212)で最近観測されてきたジョセフソン・プラズマ共鳴吸収の温度依存性を説明するため、時間依存のギンツブルク・ランダウ方程式を用いた数値シミュレーションを行った。実験条件を再現するため、シミュレーションでは熱揺らぎ、点及び柱状欠陥、磁場勾配等を含めて行った。その結果、点欠陥及び磁場勾配が共存する系ではプラズマ振動数が温度減少と共に減少するような振る舞いが観測される一方、柱状欠陥の導入は、磁束融解線以上でプラズマ振動数を増加させる働きがあることが分かった。