[工学]5 第五次课MATLAB.ppt

§ 3.5 符号矩阵的运算 一、基本运算 1. 符号矩阵的四则运算 ? 矩阵的加（+）、减（-）法 例如：a=sym(‘[1/x,1/(x+1);1/(x+2),1/(x+3)]’); b=sym(‘[x,1;x+2,0]’); b+a ? 矩阵的乘（*）、除（/、\）法 如： a\b § 3.6 符号微积分 一、符号极限 极限的求解可由limit函数来实现，调用格式： limit(F,x,a) 计算符号表达式F在x→a条件下的极限值； limit(F,a) 计算符号表达式中由findsym(F)返回的独立变量 趋向于a的极限值； limit(F) 计算a=0时的极限； limit(F,x,a,’right’)或limit(F,x,a,’left’) rigft,left用来 指定取极限的方向。 如： syms x a t h; limit(sin(x)/x) ans=1 limit((1+2*t/x)^(3*x),x,inf) ans=exp(6*t) limit(1/x,x,0,’right’) ans=inf § 3.7 符号代数、微分方程求解 一、线性方程组的符号解法 线性方程的符号求解函数 linsolve solve 此方法可得到方程组的精确解，所得的解析解可由函数vpa转换在浮点近似数值。 例：试对数值方程组用符号求解函数来求解。 a=sym(‘[10, –1, 0;-1 ,10 ,–2;0 ,–2 ,10]’) b=(‘[9;7;6]’) linsolve(a,b) vpa(ans,5) § 3.8 符号函数的二维图 一、符号函数的简易绘图函数ezplot 函数的调用格式： ezplot(f) 绘制f(t)的函数图，这里f为代表数学表达式的包含单 个符号变量x的字符串或符号表达式。x轴的近似范围 为[-2*pi,2*pi]。 ezplot(f, xmin, xmax)或ezplot(f, [xmin, xmax]) 例：绘制误差函数的图形 ezplot(‘erf(x)’) 或 ezplot erf(x) 绘制三角函数的图形 ezplot(‘sin(x)’) * * 内容: 上次课内容回顾: 1.数组及其运算 主要介绍了基本数组运算、数组函数运算和数组逻辑运算.注意数组运算和矩阵运算的不同之处. 2.多项式运算 主要介绍了多项式的表示方法、多项式运算 3.符号表达式的生成 主要介绍了符号运算的特点,表达式的形式及创建符号表达式的方法,创建符号表达式是后继符号各种运算的基础,因而在创建过程中创建方法和要注意的点一定要搞清楚. 4. 符号和数值之间的转换 主要介绍了如何将符号运算得到的解析解转换为数值解. 5. 符号函数的运算 主要介绍了复合函数的运算和反函数的运算. § 3.4 符号矩阵的创立 在MATLAB中创建符号矩阵的方法和创建数值矩阵的形式很相似，只不过要用到符号定义函数sym。 一、使用sym函数直接创建符号矩阵 矩阵元素是任何不带等号的符号表达式，各符号表达式的长度可以不同 如：a=sym(‘[1/s+x,sin(x) cos(x)^2/(b+x);9,exp(x^2+y^2),log(tanh(y))]’) a = [ 1/s+x, sin(x), cos(x)^2/(b+x)] [ 9, exp(x^2+y^2), log(tanh(y))] 二、用创建子阵的方法创建符号矩阵 ms=[‘[1/s,sin(x)]’;’[1 ,exp(x)]’] ms =[1/s,sin(x)] [1 ,exp(x)] 注意：要保证同一列的各元素字符串具有相同的长度，为此 在较短字符串的前后可用空格符补充。