我必须在计算过程中解决一些方程式。
我实现了这个目标W1 = solve(zero);
其中zero
就是这样的
[ w1_1*w2_1 + w1_2*w2_2, w2_1^2 + w2_2^2 - 1, w2_1*w3_1 + w2_2*w3_2]
[ w1_1*(w1_1 + w3_1) + w1_2*(w1_2 + w3_2) - 1, w2_1*(w1_1 + w3_1) + w2_2*(w1_2 + w3_2), w3_1*(w1_1 + w3_1) + w3_2*(w1_2 + w3_2) - 1]
我得到的是
W1 =
w1_1: [4x1 sym]
w1_2: [4x1 sym]
w2_1: [4x1 sym]
w2_2: [4x1 sym]
w3_1: [4x1 sym]
w3_2: [4x1 sym]
我可以访问这个结构的元素并构造一个符号矩阵:
>> W = [W1.w1_1(1) W1.w1_2(1); W1.w2_1(1) W1.w2_2(1); W1.w3_1(1) W1.w3_2(1)]
W =
[ 2^(1/2)/2, 0]
[ 0, -1]
[ 2^(1/2)/2, 0]
但问题是:矩阵大小取决于数据,因运行而异。如果没有明确地解释所有元素,我怎么能以更通用的方式做到这一点呢?