shingoushori's dialy

音信号処理を専ら研究していた元博士後期課程の学生によるメモ

Mathematica 9 での ImplicitPlot の代替策 → フォルマント分布図

陰関数を用いてプロットできる関数"ImplicitPlot"が,
Mathematica 9で使えなくなっていることに気づいた.

↓こちらの記事によると,「関数"ContourPlot"を用いよ」とのこと.
アップグレード情報:Graphics`ImplicitPlot`

とりあえず,ImplicitPlotをContourPlotに書き換えれば良さそう.

 

気づいたきっかけは,拓殖大学幹教授が公開なさっている,
フォルマント分布図の再描画を試みたこと.
私が書き換えを試したところ,ひと工夫必要であった.

 oval[cx0_, cy0_, ca_, cb_, c_] := (((cx - cx0) Cos[c] - (cy - cy0) Sin[c])/
      ca)^2 + (((cx - cx0) Sin[c] + (cy - cy0) Cos[c])/cb)^2;
formant = ContourPlot[{
   oval[500, 900, 260, 100, -Pi .37] == 1,
   oval[335, 1250, 340, 105, -Pi .47] == 1,
   oval[290, 2500, 700, 85, -Pi .505] == 1,
   oval[525, 2250, 600, 105, -Pi .47] == 1,
   oval[900, 1400, 540, 160, -Pi .33] == 1},
  {cx, 0, 1400}, {cy, 0, 3500}, AspectRatio -> 1, PlotPoints -> 100]