《經济数学》实训课.pptVIP

Matlab中绘制带网格的曲面图使用mesh函数。 *例23： 绘制 的图像。 syms x y z r x=linspace(-10,10,50); y=linspace(-10,10,50); [xx,yy]=meshgrid(x,y); r=sqrt(xx.^2+yy.^2); z=sin(r)./r; mesh(z) 命令格式： int(函数名) 求不定积分 int(函数名，a，b) 求在[a,b]区间内的定积分 例24：求 。 int(2*x) ans = x^2 注：求不定积分得到的结果，只是被积函数的一个原函数，并没有加C。 例25： 求 int(sin(3*x-2)) ans = -1/3*cos(3*x-2) 10.求不定积分和定积分 例26 求 int(sin(x),0,pi) ans = 2 例27 求 （即求四分之一个单位圆的面积） int(sqrt(1-x^2),0,1) ans = pi/4 11.求偏导数 命令格式：（与求导数格式相同） diff(函数，变量，阶数) 例28 设 求 syms x y diff(x^3+3*x^2*y-2*y^2,x,1) ans = 3*x^2+6*x*y diff(x^3+3*x^2*y-2*y^2,y,1) ans = 3*x^2-4*y 例29 设 求二阶偏导数 diff(x^3+x^2*y^2-y^3,x,2) ans = 6*x+2*y^2 diff(x^3+x^2*y^2-y^3,y,2) ans = 2*x^2-6*y diff(diff(x^3+x^2*y^2-y^3,x,1),y,1) ans = 4*x*y diff(diff(x^3+x^2*y^2-y^3,y,1),x,1) ans = 4*x*y 12. 矩阵及其运算 12.1 矩阵的表示 命令格式：[1 2 3 4;2 3 4 5;3 4 5 6] 同一行用空格或逗号分隔，行与行用分号或回车符分隔。 例30 产生三行四列的矩阵 [1 2 3 4;5 6 7 8;9 10 11 12] ans = 1 2 3 4 5 6 7 8 9 10 11 12 12.2 矩阵的基本运算 （1）加法、减法、数乘 例31 设 求 syms A B A=[1 2 3;3 2 1];B=[-1 5 4; 0 -3 1]; A+B ans = 0 7 7 3 -1 2 2*A-3*B ans = 5 -11 -6 6 13 -1 （2）矩阵乘法 例32 设 求 A=[1 2 3;3 2 1];B=[1 2;0 1;-1 3]; A*B ans = -2 13 2 11 B*A ans = 7 6 5 3 2 1 8 4 0 12.3 行列式的计算 命令格式：det（A） 例33 计算 det([1 2 3;3 1 2;2 3 1]) ans = 18 det([1 1 1 1;1 3 5 7;1 9 25 49;1 27 125 343]) ans = 768 12.4 逆矩阵 命令格式：inv（A） inv([1 1 1 1;0 1 1 0;0 0 1 1;0 0 0 1]) ans = 1 -1 0 -1 0 1 -1 1 0 0 1 -1 0 0 0 1 求 例34 设 12.5 解线性方程组 （1） 化增广矩阵 为最简形 命令格式：rref( ) 例35 解线性方程组 rref([1 -2 3 1 1 7;1 1 -1 -1 -2 2;2 -1 1 0 -2 7;2 2 5 -1 1 18]) ans = 1 0 0 0 -2 3 0 1 0 0 -1 1 0 0 1 0 1 2 0 0 0