贵大matlab实验五.docVIP

实验五 MATLAB符号运算 一、实验目的 1．掌握符号对象的定义方法以及符号表达式的运算法则。 2．掌握微积分的符号运算方法。 3．掌握级数求和的方法。 4. 掌握代数方程和微分方程符号求解的方法。 5. 了解数值运算与符号运算的区别。 二、实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 设计提示： 注意符号运算的相关函数的帮助查询是在 sym目录下，应采用doc sym/函数名进行查询。 三、实验内容 1．运行以下命令，把各函数的功能及参数形式抄写下(不记录过程和结果)； (1)定义符号变量 syms syms x y t (2)定义符号函数(或符号表达式)sym f= sin(t)+x^2+y^3 g=sym(‘sin(t)+x^2+y^3’) (3)符号函数的变量替换 subs h=subs(f, x, sqrt(x)) s=subs(f,{x,y},{t^2,sqrt(t)}) (4)符号表达式到数值的转换 double g= sym(pi+sqrt(2)) double(g) (5)合并同类项 collect EXPR=sym((x^2+x*exp(-t)+1)*(x+exp(-t))); expr1=collect(EXPR) % 默认合并 x 同幂项系数 expr2=collect(EXPR,exp(-t)) % 合并 exp(t)同幂项系数 (6)因式分解 factor syms a x? f1=x^2-a^2? g1=factor(f1) (7)展开 expand expand(g1) 2. 已知x=6，y=5,利用符号表达式求 提示利用sym函数定义符号常数x,y。 3．用符号方法求解下列极限或导数或积分 （1） （2） （3） （4）已知，分别求，求 4．求解非线性方程组的解；（利用solve函数） 5．求微分方程 x + 2 x + x = sin t , x ( 0 ) = x ( 0 ) = 1 的解，利用函数ezplot绘制其图形。同时计算其数值解进行比较。 四、实验报告要求（包含预习报告要求和最终报告要求） 1.实验名称 2.实验目的 3.实验设备及条件 4.实验内容及要求 5.实验程序设计 指程序代码。 6.实验结果及结果分析 实验结果要求必须客观,有数据的可以记录数据，没有数据的简单描述实验现象。结果分析是对实验结果的理论评判。 7.实验中出现的问题及解决方法 8. 思考题的回答 五.实验程序设计 clc format compact %1 %(1)定义符号变量 syms syms x y t %(2)定义符号函数(或符号表达式)sym f=sin(t)+x^2+y^3 g=sym(sin(t)+x^2+y^3) %(3)符号函数的变量替换 subs h=subs(f,x,sqrt(x)) s=subs(f,{x,y},{t^2,sqrt(t)}) %(4)符号表达式到数值的转换 double g=sym(pi+sqrt(2)) double(g) %(5)合并同类项 collect EXPR=sym((x^2+x*exp(-t)+1)*(x+exp(-t))); % 默认合并 x 同幂项系数 expr1=collect(EXPR) % 合并 exp(t)同幂项系数 expr2=collect(EXPR,expr(-t)) %(6)因式分解 factor syms a x; f1=x^2-a^2; g1=factor(f1) %(7)展开 expand expand(g1) %2 syms x y z x=6; y=5; z=(x+1)/(sqrt(3+x)-sqrt(y)) %3．用符号方法求解下列极限或导数或积分 %(1) syms x; f2=((5*x+3)/(5*x+1))^(x+1); limit(f2,x,inf) %(2) f3=(tan(x)-sin(x))/(x^2*sin(x)); limit(f3,x,0) %(3) y=exp(-x)*sin(x); diff(y) diff(y,2) int(y,0,1) %(4) syms x t A=[a*x^2 t^3;t*cos(x) log(x)]; diff(A,x) diff(diff(A,t)) diff(diff(A,x),t) %4 syms x y f4=(x^2-2*x*y+y^2-3); f5=(x^2-4*x+3); [x,y]=so