虚数単位を含む数式の実部と虚部を一つのグラフに描きたい

Question

matplotlibを使って虚数単位を含む数式の結果のプロットをしたいです。

###発生している問題・エラーメッセージ
x and y must have same first dimensionのエラーメッセージが出ます

###該当のソースコード
import numpy as np, scipy.integrate as spi, scipy as sp, scipy.special as sps
import matplotlib as mpl, matplotlib.pyplot as plt
import sys

mpl.rcParams['font.family'] = 'serif'
mpl.rcParams['font.serif'] = 'Times New Roman'
mpl.rcParams['font.size'] = 15

def DelanyBazley(f, sigma):
    fDivSigma = f / sigma
    return 1. + 9.08 * (fDivSigma ** -0.75) + 11.9j * (fDivSigma ** -0.73)

fig, ax = plt.subplots(1, 1, figsize=(120/25.4, 100/25.4))

sigma = 300
fs = np.arange(125, 5000, 20)

for fI in range(len(fs)):
    f = fs[fI]
    z = np.real(DelanyBazley(f, sigma))

ax.plot(fs, z)

print fs, z

ax.grid()
ax.set_xscale('log')
ax.set_xlim(125, 5000)
ax.set_xticks([125, 250, 500, 1000, 2000, 4000])
ax.set_ylim(0, 10)
ax.minorticks_off()
ax.xaxis.set_major_formatter(mpl.ticker.FormatStrFormatter('%d'))

fig.savefig('delany-bazley.pdf')
plt.show()


###補足情報(言語/FW/ツール等のバージョンなど)
pythonとmatplotlibを使って言います

Comments

[KSwordOfHaste]

pythonのコードの場合字下げがないと閲覧者が混乱します。質問を編集し、コード部分の行範囲を選択した後<code>ボタンを押してください。

Answer 1

for fI in range(len(fs)):
f = fs[fI]
z = np.real(DelanyBazley(f, sigma))

の部分を

z = np.real(DelanyBazley(fs, sigma))

にできませんか。

エラーが出るとしたら、それを追記していただけると回答できるかもしれません。