我有一个矩阵,A
:
A=[ x.^2 + y , 0;
0 , x.^2 + y ]
我还有一个运算符矩阵D
:
D = [d/dx , 0;
0 , d/dy ]
我希望能够相乘D*A
并得到一个如下所示的矩阵:
B = [ diff(A(1,1),x) , 0 ;
0 , diff(A(2,2),y) ]
显然我不能用 diff() 函数做到这一点,因为该函数不是可以乘以函数的运算符。那么如何使用符号运算符来解决这个问题呢?实际上,我的矩阵很大,因此在没有运算符乘法的情况下执行是不可取的。
另外,假设我找到了一种生成上述B
矩阵的方法,如下所示:
B =
[ 2*x, 0
0, 1]
B
例如,我如何评估x=2, y=1
;
我的尝试:
subs(B,x,2,y,1)
但这显然是sym.subs
函数的不正确参数
我也试过:
subs(B,2,1)
这也不起作用,所以我的另一个问题是如何在矩阵中替换x
和替换。y
B