迭代matlab符号求解函数

Question

我正在尝试求解数组 d 的以下方程。我使用了下面的代码片段：

channel_size = 9e-3;
d = [11e-3, 12e-3]; 
sigma = 0.49;
ee = 727/9806.65; 
alpha = d-channel_size;
sym('p',[1 2])

for i = 1:2
    eqn = alpha == (4.3^(1/7))*(p^(3/11))*(((1-(sigma^2))/ee)^(7/5))/(d.^(1/6))
    S = solve(eqn, p)*0.015;
    vpa(S/13e-12)
end

事实上，我应该得到与 d(1) 和 d(2) 相对应的两个数字，但它不起作用并且出现此错误:

Error using mupadmex
Error in MuPAD command: Operands are invalid. [linalg::matlinsolve]
Error in sym/privBinaryOp (line 1693)
            Csym = mupadmex(op,args{1}.s, args{2}.s, varargin{:});
Error in sym/mrdivide (line 232)
        X = privBinaryOp(A, B, 'symobj::mrdivide');

Answer 1

声明 (d.^(1/6)) 是错误的，应该使用 (d(i)^(1/6))。另外，由于 alpha = d-channel_size，在等式中，我还应该声明 alpha(i) 而不是简单的 alpha。