PID controller implementation about bdsim HOT 4 CLOSED

Hi @zam101 and @MartinGrignard, check out the thread in Issue #18. There is a new implementation of PID on the master branch version.

from bdsim.

Not yet, sorry. The usual trick is to make the d term s/(s+a) where a is another design parameter, set somewhere toward the open loop bandwidth of the plant

from bdsim.

I've been thinking about this and it's much more complex than I initially thought. My previous answer to @zam101 was hasty and not well thought through, it's more complex than I said then. I tried to put some of this down on a wiki page. My main focus with bdsim right now is getting it ready for the forthcoming book, and after that a derivative block and PID block are top of the feature list.

from bdsim.

Hi @petercorke,
First, thanks for the work done on this wonderful package!
I'm also trying to build a PID controller with it and, following your previous answer that recommends adding a low-pass filter to the derivative term, the expression becomes

where $\alpha = \frac{1}{\omega_0}$ and $\omega_0$ is the cutoff frequency of the system.

Focusing solely on the derivative term, if I write

derivator = bd.LTI_SISO([K_d, 0], [alpha, 1])

I get the exact same error as @zam101:

---------------------------------------------------------------------------
AssertionError                            Traceback (most recent call last)
Cell In[29], line 1
----> 1 derivator = bd.LTI_SISO([k_d, 0], [alpha, 1])

File [.../bdsim/run_sim.py:806), in BDSim.blockdiagram..new_method..block_init_wrapper(self, *args, **kwargs)
    804 def block_init_wrapper(self, *args, **kwargs):
--> 806     block = cls(*args, bd=bd, **kwargs)  # call __init__ on the block
    807     return block

File [.../bdsim/blocks/transfers.py:332), in LTI_SISO.__init__(self, N, D, x0, **blockargs)
    330 if x0 is None:
    331     x0 = np.zeros((n,))
--> 332 assert nn <= n, "direct pass through is not supported"
    334 # convert to numpy arrays
    335 N = np.r_[np.zeros((len(D) - len(N),)), np.array(N)]

AssertionError: direct pass through is not supported

Do you have some hints on how I could actually make it work?

from bdsim.