syms x1 x2 ;%定義變數
a = [x1 x2];
f = [x1*x1+t*x1;x2*cos(x2)]; %定義函式，以矩陣的形式 
x = jacobian(f,a); % 求取雅可比矩陣，會發現x是sym型別的
b = [1 2]; 
y = subs(x,a,b); %賦值
%此時結果
y =
[ 5/2,               0]
[   0,   cos(2) - 2*sin(2)]
y = eval(y); %轉化為double型普通矩陣，
%結果
y =
    2.5000      0
       0          -2.2347
 

相關推薦