2矩阵数组符号运算.ppt

单变量非线性方程求解格式：fzero(filename,x0,tol,trace)求解x0最近的解。函数名，初值，精度，是否显示迭代信息例题：求f(x)=x-1/x+5在x0=-5和1作为迭代初值的解。functionf=aa(x)f=x-1/x+5x1=-5.1926x2=0.1926x1=fzero(aa,-5,0)x2=fzero(‘aa’,1,1)非线性方程(组)求解格式：fsolve(filename,x0,option)x=fsolve(@(x)sin(3*x),[14],optimset(Display,off))x=1.04724.1888Matlab的符号计算基本含义：运算对象是没赋值的符号变量，运算结果以标准的符号形式表达，可以获得任意精度的解。建立符号变量1.定义单个符号量sym函数：符号变量名=sym(符号字符串)。；2.定义多个符号量syms函数：syms=arg1arg2arg3arg4…建立符号表达式1.利用单引号生成：y=‘1/sqrt(2*x)’;2.利用函数sym生成：y=sym(‘1/sqrt(2*x)’)y=sym(‘[a,b;c,d]’)用sym函数生成符号表达式aex+bsin(x).f=sym(‘a*e^x+b*sin(x)’)f=a*e^x+b*sin(x)通过findsym函数确定表达式中的符号变量.symsabcxf=a*x^2+b*x+c;f,findsym(f)f=a*x^2+b*x+cans=a,b,c,x用subs函数可以将符号表达式中的符号变量用数值代替.f=sym(‘x^2+y^2’);g=subs(f,x,3);f,g,subs(g,4)f=x^2+y^2g=9+y^2ans=25基本符号运算使用数值运算符：+，-，*，/，^2.符号表达式symadd,symsub,symmul,symdiv提取分子和分母[n,d]=numden(s),将符号表达式有理化后，分子分母存放为n和d，符号表达式为符号矩阵，n和d分别是分子矩阵和分母矩阵。a11=sym(‘1/(x-1)’);%建立符号变量a12=sym(1/(x+1));a21=a11+a12;a22=a11*a12;a=[a11,a12;a21,a22];%符号矩阵a11,a12,a21,a22,a,[n,d]=numden(a)%numden(a)提取分子和分母a=[1/(x-1),1/(x+1)][1/(x-1)+1/(x+1),1/(x-1)/(x+1)]n=[1,1][2*x,1]d=[x-1,x+1][(x-1)*(x+1),(x-1)*(x+1)]a11=1/(x-1)a12=1/(x+1)a21=1/(x-1)+1/(x+1)a22=1/(x-1)/(x+1)符号表达式的因式分解与展开factor(f)对符号表达式f分解因式；expand(f)对符号表达式进行展开；collect(f)合并同类项；collect(f,v)对符号表达式按v合并同类项。symsabxyA=a^4-b^4;A1=factor(A);%对A分解因式B=(x+y)^5;A2=expand(B);%对B展开A3=collect(A2,x);%对A2按x展开A4=collect(A2,y);%对A2按y展开A5=factor(sym(‘120’));%对符号整数分解因式A1,A2,A3,A4,A5A1=(a-b)*(a+b)*(a^2+b^2)A2=x^5+5*x^4*y+10*x^3*y^2+10*x^2*y^3+5*x*y^4+y^5A3=x^5+5*x^4*y+10*x^3*y^2+10*x^2*y^3+5*x*y^4+y^5A4=x^5+5*x^4*y+10*x^3*y^2+10*x^2*y^3+5*x*y^4+y^5A5=(2)^3*(3)*(5)符号表达式的化简符号表达式与数值表达式之间的转换符号表达式之间的转换simplify(f)应用函数规则对f进行化简；s