02matlab符号计算分析.pptVIP

MATLAB的符号运算;符号运算的功能;一、符号运算的基本操作;特点： ? 运算对象可以是没赋值的符号变量 ? 可以获得任意精度的解 　　 ;2. 符号变量与符号表达式; ;3.符号矩阵的创建 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 ?用matlab函数sym创建矩阵（symbolic 的缩写) 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用 标识;例如：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方 括号，这是与 matlab数值矩阵 的一个重要区别。;?用字符串直接创建矩阵;? 符号矩阵的修改;;例如：A =[ a, 2*b] [3*a, 0] A1=sym(A,2,2, 4*b) A1 =[ a, 2*b] [3*a, 4*b] ;【例2.2-4】 （1）syms a b x; f=a*sin(x)+b f1=subs( f, sin(x), log(y)) class(f1) f1 =b + a*log(y) ans =sym ;符号矩阵运算 数值运算中，所有矩阵运算操作指 令都比较直观、简单。例如：a=b+c; a=a*b ；A=2*a^2+3*a-5等。 而符号运算就不同了，所有涉及符 号运算的操作都有专用函数来进行;例1：f= 2*x^2+3*x-5; g= x^2+x-7; h= symadd(f,g) h= 3*x^2+4*x-12 ;例2：f=cos(x);g= sin(2*x); syms x f=cos(x); g=sin(2*x); f/g+f*g ans = cos(x)/sin(x)+cos(x)*sin(x) ;diff(f) — 对缺省变量求微分(62) diff(f,v) — 对指定变量v求微分 diff(f,v,n) —对指定变量v求n阶微分 int(f) — 对f表达式的缺省变量求积分(68) int(f,v) — 对f表达式的v变量求积分 int(f,v,a,b) — 对f表达式的v变量在（a,b) 区间求定积分 ;【例2.3-2】;Ans:;【例2.3-5】;例2.3-9求不定积分;【例2.3-10】;f2 = [ a*x, b*x^2] [ 1/x, sin(x)] The integral of f is +- -+ | 2 3 | | a x b x | | ----, ---- | | 2 3 | | | | log(x), -cos(x) | +- -+ ;【例2.3-11】; ;例2.3-12 求阿基米德螺线r=a*θ(a0)在θ=0到 间的曲线长度函数，并求出a=1, 时的曲线长度。;【例2.3-12】;（2）a=1时的曲线长度 L_2pi=subs(L,[a,phi],sym([1,2*pi])) L_2pi_vpa=vpa(L_2pi) L_2pi = log(2*pi + (4*pi^2 + 1)^(1/2))/2 + pi*(4*pi^2 + 1)^(1/2) L_2pi_vpa = 21.2562941482090988007025122725661088234709310476 ;（3）图形绘制;se