ngspiceはこちら(https://ngspice.sourceforge.io/index.html)。Linuxでは最新版はソースコードからビルドするしかないが、マニュアル通りやればいい。この辺の情報はググれば見つかるんだけど、いざ細かいことをしようとすると、PSpiceやらLTSpiceやらに駆逐されていて情報が見つからない。
PySpiceはこちら(https://pyspice.fabrice-salvaire.fr/releases/v1.4/index.html#)。まぁpipでいける。
今回何をしたいのかというと、パラメータを変えながらシミュレーションして結果の重ね書きのグラフを作成するってこと。そんなのngspiceの回路記述のalterparamでできるじゃん(?)って思うんだけど、これがうまくいかない。alterparamした後はresetせんといかんのだけど、前のシミュレーション結果も忘れてしまうっぽくて、後で重ね書きしようとすると、そんなベクタは知らんって怒られる。
では、使ってみる。
こちら(https://pyspice.fabrice-salvaire.fr/releases/v1.4/examples/filter/low-pass-rc-filter.html)の例を試してみる。が、単位の書き方がうざい。"@u_"を付けることで、技術者が見慣れた単位を付けることができるんだけど、いや、むしろうざいから。そして、このまま実行しても色々怒られる。"in"はPythonの予約語だから使うんじゃねーっとか、その他意味わからんものがいろいろ。
で、色々いじくった結果、これでいけた。
-
import math -
import numpy as np -
import matplotlib.pyplot as plt -
import PySpice.Logging.Logging as Logging -
logger = Logging.setup_logging() -
from PySpice.Plot.BodeDiagram import bode_diagram -
from PySpice.Spice.Netlist import Circuit -
circuit = Circuit('Low-Pass RC Filter') -
circuit.SinusoidalVoltageSource('input', 'vin', circuit.gnd, amplitude=1) -
R1 = circuit.R(1, 'vin', 'vout', 1e3) -
C1 = circuit.C(1, 'vout', circuit.gnd, 1e-6) -
break_frequency = 1 / (2 * math.pi * float(R1.resistance * C1.capacitance)) -
print("Break frequency = {:.1f} Hz".format(break_frequency)) -
simulator = circuit.simulator(temperature=25, nominal_temperature=25) -
analysis = simulator.ac(start_frequency=1, stop_frequency=1e6, number_of_points=10, variation='dec') -
freqs=np.array(analysis.frequency) -
vout=np.array(analysis.vout) -
fig = plt.figure() -
ax1=plt.subplot(2,1,1) -
ax1.grid(True) -
ax1.plot(freqs,20*np.log10(np.absolute(vout)), '-') -
ax1.set_ylabel('Amplitude [dB]') -
plt.xscale('log') -
ax2=plt.subplot(2,1,2) -
ax2.grid(True) -
ax2.plot(freqs,np.angle(vout), 'r-') -
ax2.set_xlabel('frequency [Hz]') -
ax2.set_ylabel('Phase [rad.]') -
plt.xscale('log') -
fig.tight_layout() -
plt.show()
低賃金長時間労働への抗議のためパジャマで通勤するってのが隣の超大国で起こっているらしい。労働法によって労働時間や賃金についてそんなに無茶なことはできんはずだけど、抗議したくなるほどってことなんやろうなー、、、ちなみに日本も25年前は○○ムカー○”打刻してから○フ。○に戻って深夜まで○事するってのは普通だった。そのころの日本(または今の隣の超大国)と今の日本を表向きの生産性で比較したらいかんよね、、、
コメントをお書きください