IEEE802.15.3cの位相雑音モデルの位相雑音データ生成するMATLABソースコードを教えて下さい。

位相雑音のパワースペクトラムは以下の式です。

PSD(f)=PSD(0){1+(f/fz)^2}/{1+(f/fp)^2}
PSD(0)=-80dBc/Hz
fz=100MHz
fp=1MHz

です