我是一位机械工程师,对数学的背景知识最少。我很想使用CVXPY编写有关转子平衡问题的简单代码。如此简单而强大的功能令人惊讶。我的问题如下:
在获得用户输入的矩阵ALPHA(M,N)
和A(M,1)
后
W=cvxpy.Variable((N,1),complex=True)
objective2=cvxpy.Minimize(cvxpy.norm((ALPHA*W+A),"inf"))
prob2.solve()
W很复杂,因为它返回平衡转子的每个N
平面的权重及其角度。
我的问题是如何对W
施加约束。例如,如何确保每个W
平面的N
小于特定值?
在获得约束质量矩阵wc(N,1)
的用户输入后,以下代码似乎可以解决问题:
const=[]
const += [cp.norm(W2[i])<=wc[i]for i in range (N) ]
prob3=cp.Problem(objective2,const)
prob3.solve()