关于机器人独立关节PD控制的一个问题

在动力学代码模块有一行代码
S=inv(D0)(tol-C0dq)
想请教一下这个是怎么计算的？
我能理解这是求解ddq，但不知道这个计算的意义，或者作用。