《经济数学》实训课

《经济数学》实训课 用MATLAB解决数学问题任课老师金慧萍

主要内容 • 一、Matlab基本知识复习 • 二、用Matlab求积分与求偏导数 • 三、用Matlab求行列式与矩阵运算 • 四、用Matlab解线性方程组

一、Matlab基本知识复习 1.1 启动与退出的命令启动：鼠标双击桌面上的Matlab 图标退出：在Matlab 的环境中，键入 quit(或 exit) 并回车

一、Matlab基本知识复习 1.2 使用Matlab软件时,需特别注意几点 ① 必须在英文状态下输入 ② clc并回车=清屏 ③ clear并回车=取消变量 ④ 乘法运算符号 * ，不能省略 ⑤ 分号；表示先不运算 ⑥ inf= ⑦ pi= ⑧ syms x y =定义变量 x y

一、Matlab基本知识复习 1.3 常用函数的表示函数表示开平方 sqrt(x) or x^(1/2) 幂函数 x^n 以e为底的指数 exp(x) 指数函数 a^x 自然对数 log(x) 以10为底的对数 log10(x) 正弦函数 sin(x) 反正弦函数 asin(x) 绝对值 abs(x).

一、Matlab基本知识复习 1.4 复习上学期介绍的一些命令解方程： solve (‘方程’，’变量’) 求极限： limit (函数，变量，趋近值) 求导数： diff(函数名，变量名，n)

二、用Matlab求积分与求偏导数 2.1 求不定积分命令格式：int(函数名) （表示求不定积分）例1：求 >> syms x >> int(2*x) ans = x^2 注：求不定积分得到的结果，只是被积函数的一个原函数，并没有加C。例2：求 >> int(sin(3*x-2)) ans = -1/3*cos(3*x-2)

二、用Matlab求积分与求偏导数 2.2 求定积分命令格式：int(函数名，a，b) （表示求在[a,b]区间内的定积分）例3 求 >> int(sin(x),0,pi) ans = 2 例4 求（即求四分之一个单位圆的面积） int(sqrt(1-x^2),0,1) ans = 1/4*pi

二、用Matlab求积分与求偏导数 2.3 求偏导数命令格式：diff(函数名，变量名，阶数)（与求导数格式相同）例5 设求 >> 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

二、用Matlab求积分与求偏导数 2.3 求偏导数例6 设求二阶偏导数 >> 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

作业1 1.求不定积分： 2.求定积分： 3.设：求：

三、用Matlab求行列式与矩阵运算 3.1 行列式的计算命令格式：det（A）例7 计算 >> 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

三、用Matlab求行列式与矩阵运算 3.2 矩阵的表示命令格式： (以上是三行四列的矩阵，同一行用空格分隔，行与行用分号分隔) 例8 输入矩阵 >> [1 2 3 ;3 2 1] ans = 1 2 3 3 2 1

三、用Matlab求行列式与矩阵运算 3.3 矩阵的基本运算：加法、减法、数乘、乘法例9 设求 >> 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

三、用Matlab求行列式与矩阵运算 3.3 矩阵的基本运算：加法、减法、数乘、乘法例10 设求 >> 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

三、用Matlab求行列式与矩阵运算 3.4 求逆矩阵命令格式：inv（A）求例11 设 >> 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

四、用Matlab解线性方程组 4.1 线性方程组有唯一解时，即命令格式：（注：当然这种情况也可以用后面的rref命令来求解）例12 解线性方程组 >> A=[1 -2 3 1;1 1 -1 -1;2 -1 1 0;2 2 5 -1]; >> B=[7;2;7;18]; >> A\B ans = 3.0000 1.0000 2.0000 0.0000 >> rref([1 -2 3 1 7;1 1 -1 -1 2;2 -1 1 0 7;2 2 5 -1 18]) ans = 1 0 0 0 3 0 1 0 0 1 0 0 1 0 2 0 0 0 1 0

四、用Matlab解线性方程组 4.2 化增广矩阵为最简形命令格式：rref( ) 例13 解线性方程组 >> 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 1 -2 0

作业2 1.求行列式 3.解线性方程组 2.设求：

《 经济数学 》 实训课