MATLAB数学实验_12343.pptVIP

数据类型转换 syms x y; s= diff(x^2*exp(-y),x,2) t=diff(x^2*exp(-y),x); t=diff(t,y); t=subs(t,x,1); t=subs(t,y,2) syms x; g=log(x)*sin(x); gt=taylor(g, x, 1,Order,6); %展开到5次式, Order=6 gt=vpa(gt,4) s=dsolve(‘Dy==a*y+b’) %字符串方式 syms a b y(t); %符号函数方式 s=dsolve(diff(y)==a*y+b) s=simple(dsolve(‘D2y==cos(2*x)-y’,‘y(0)==1’,‘Dy(0)==0’,‘x’)) %初值问题 s= simple(dsolve(‘x*D2y-3*Dy==x^2’,‘y(1)==0’,‘y(5)==0’,‘x’)) %边值问题 MATLAB数学实验 第七章 MATLAB符号计算 第七章 MATLAB符号计算 7.1 符号对象 7.2 符号矩阵和符号函数 7.3 符号微积分 7.4 符号方程和符号微分方程 7.5 便捷函数作图 7.6 符号计算局限性和Mupad调用 １、符号对象的定义 s=sym(str) 将数值或字符串str转化为 符号对象s, 数值为有理数表示; s= sym(num,’d’) 将数值表达式转化为符号 表达式，数值用小数表示； syms var1 var2 ? 定义var1, var2, ?为符号变量 subs(s, old, new) 将符号表达式s中的符号 变量 old用new代替 7.1 符号对象 比较: a=pi^2, b=pi^2, c=sym(pi^2),d= sym(pi^2,d), e=sym(pi^2) 2 .计算精度和数据类型转换 digits(n) 将数值计算精度设为n位； x=vpa(s) 求s的数值结果 x=vpa(s,n) 用n位精度求s的数值结果 double(s) 将符号对象转化为双精度数值 char(s) 将符号对象转化为字符串 例 pi与2^10000 vpa(pi),vpa(pi,100) 2^10000,sym(2)^10000 1.矩阵 MATLAB大部分矩阵和数组运算符及指令都可以用于符号如：inv，./，.\，eig等 syms a b;A=[1,2;a b] det(A),inv(A) 7.2 符号矩阵和符号函数 2.符号函数计算 factor(expr) 对expr作因式分解 expand(expr) 将expr展开 collect(expr,v) 将expr按变量v合并同类项 simple(expr) 将expr化简 g=finverse(f,v) 求函数f(v)的反函数g(v) fg=compose(f, g) 求函数f(v)和g(v)的 复合函数f(g(v)) [n,d]=numden(expr) 分式通分，n分子d分母 funtool 函数分析图形界面 2.符号函数计算 symfun(expr, arg) 定义符号函数 syms fun(var1, var2, …) 定义符号函数 latex(expr) 数学公式的Latex输出 ccode(expr) 数学公式的C语言代码 matlabFunction(expr) 数学公式的Matlab函数, 注意这里F是大写 funtool 函数分析图形界面 例1. 已知 f (x,y) = (x-y)3, g(x,y) = (x+y)3 求 h= f*g 并展开。 解 ?syms x y;f=(x-y)^3;g=(x+y)^3; h=f*g,hs=expand(h) hf=factor(hs) s=subs(h,y,x^2+x+1) fun=symfun(f*g,[x,y]) s=fun(x,x^2+x+1) scol=collect(s,x) ssim=simple(scol) funtool 1．极限和级数 limit(s,x,a) 返回符号表达式s当x-a时的极限 limit(s,x,a,’right’) 返回s当x-a时的右极限 例2 l