数学运算

dot(x,y) 点乘
cross(x,y) 叉乘
diag(V,K) 生成对角矩阵，或者取矩阵主对角线元素。
inv(A) 矩阵求逆，等价于A^-1。
rank(A) 获取矩阵的秩
det 获取矩阵对应行列式的值
roots 求多项式的根
poly 求矩阵的特征多项式

1
2
3
4
5
6
7
8
9
10
11
12
%求r对应的特征多项式以及其根
>> r=[1 3 5];

>> p=poly(r)
p =
     1    -9    23   -15
 
>> roots(p)
ans =
    5.0000
    3.0000
    1.0000
1
2
3
4
5
6
7
%求函数的一重积分与二重积分。
>> f2=@(x)sin(x)+cos(x);
ans2=quad(f2,1,10)
f4=@(x)x./(x.^4+4);
ans4=quad(f4,1,10)
f5=@(x,y)sin(y).*(x+y)./(x.^2+4);
ans5=dblquad(f5,1,10,1,10)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
% 解下列线性方程组
        4 x1 +2 x3 +3 x4 = 9
        3 x2 +2 x3 + 3x4 = 8
        2 x1 + x4 = 10
        x3 -2 x2 + 3x4 =7
答案：A=[4 0 2 3; 0 3 2 3; 2 0 0 1; 0 -2 1 3;];    %3分
      b=[9 8 10 7]';                         %2分
      x=A\b; 或x=inv(A)*b;               %3分

 
>> f = 2*x1 + 3*x2 + x3 + 2*x4 -8;
>> g = x1 + 3*x2 + x4 - 6;
>> h = x1 - x2 + x3 + 8*x4 - 1;
>> i = 7*x1 + x2 - 2*x3 + 2*x4 -5;
>> [x1,x2,x3,x4] = solve(f,g,h,i)

方程ax=b(a为非奇异)两种求解方法： x=inv(a)*b — 采用求逆运算解方程
x=a\b — 采用左除运算解方程
solve(f1,f2, …fn) —— 求n个方程的解
quad 一元函数的积分，dblquad 二重函数的积分。

函数运算

  • 嵌套函数
    funtion [x,y]=fun(a,b,c) %x、y为返回值，a、b、c为输入参数。
  • 匿名函数
    fun=@(x,y)x+y

函数文件的第一行必须以function开头，脚本文件可调用函数文件。