实验二--MATLAB的符号计算与可视化.doc

《实验二--MATLAB的符号计算与可视化.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验二MATLAB的符号计算与可视化1、目的和要求（1）熟练掌握MATLAB符号表达式的创建、代数运算及化简。（2）熟悉符号方程的求解。（3）熟练掌握MATLAB二维曲线、三维图形的绘制。（4）熟练掌握各种特殊图形的绘制。2:实验内容：（1）完成教材实验三第1节“1.创建符号表达式和符号表达式的操作”中（1）-（5）部分的内容，分别用sym和syms创建符号表达式f和g，并对它们进行相关操作，思考每一条命令的作用是什么，并提交命令行和结果；3:实验内容：（1）完成教材实验三第1节“1.创建符号表达式和符号表达式的操作”中（1）-（5）部分的内容，分别用sym和syms创建符号表达式f和g，

2、并对它们进行相关操作，思考每一条命令的作用是什么，并提交命令行和结果；（1）创建符号表达式：使用sym命令创建符号表达式：f=sym('sin(x)')f=sin(x)>>g=sym('y/exp(-2*t)')g=y*exp(2*t)使用syms命令创建符号表达式：>>symsxyt>>f=sym(sin(x))f=sin(x)>>g=sym(y/exp(-2*t))g=y*exp(2*t)（2）：自变量的确定：>>symvar(g)ans=[t,y]>>symvar(g,1)ans=y>>findsym(g,2)ans=y,t（3）：用常数替换符号变量：>>x=0:10;>>y=sub

3、s(f,x)y=Columns1through800.84150.90930.1411-0.7568-0.9589-0.27940.6570Columns9through110.98940.4121-0.5440（4）：符号对象与数值的转换和任意精度控制：>>f1=subs(f,'5')f1=sin(5)>>y1=double(f1)y1=-0.9589>>y2=eval(f1)y2=-0.9589采用digits和vpa实现任意精度控制：>>digitsDigits=32>>vpa(f1)ans=-0.95892427466313846889315440615599>>vpa(f1,10

4、)ans=-0.9589242747（5）：求反函数和复合函数用finverse函数求f,g的反函数>>f=sym('sin(x)');>>g=sym('y/exp(-2*t)')g=y*exp(2*t)>>finverse(f)Warning:finverse(sin(x))isnotunique.ans=asin(x)>>finverse(g)ans=y/exp(2*t)>>finverse(g,'t')ans=log(t/y)/2用compose函数求f,g的复合函数>>compose(f,g)ans=sin(y*exp(2*t))>>compose(f,g,'z')ans=sin(

5、z*exp(2*t))（2）自建两个一元四次符号表达式，分别进行其符号表达式的加、减、乘等运算，并提交命令行和结果；>>symsx>>f=x^4+2*x^3-5*x^2+6*x+8f=x^4+2*x^3-5*x^2+6*x+8>>class(f)ans=sym>>g=2*x^4-5*x^3+8*x^2+7*x-2g=2*x^4-5*x^3+8*x^2+7*x-2>>f+gans=3*x^4-3*x^3+3*x^2+13*x+6>>f-gans=-x^4+7*x^3-13*x^2-x+10>>f*gans=(x^4+2*x^3-5*x^2+6*x+8)*(2*x^4-5*x^3+8*x^2+

6、7*x-2)（3）自建一个可化简一元五次多项式和一个三角函数符号表达式，依次使用pretty,horner,factor,simplify和simple等函数对该表达式进行化简，并提交命令行和结果；symsxy>>f=x^5+x^4+2*x+2f=x^5+x^4+2*x+2>>g=cos(y)^2-sin(y)^2g=cos(y)^2-sin(y)^2>>class(f)ans=sym>>class(g)ans=sym>>pretty(f)54x+x+2x+2>>horner(f)ans=x*(x^3*(x+1)+2)+2>>factor(f)ans=(x+1)*(x^4+2)>>simp

7、lify(f)ans=(x^4+2)*(x+1)>>simple(f)simplify:x^5+x^4+2*x+2radsimp:x^5+x^4+2*x+2simplify(100):(x^4+2)*(x+1)combine(sincos):x^5+x^4+2*x+2combine(sinhcosh):x^5+x^4+2*x+2combine(ln):x^5+x^4+2*x+2factor:(x+1)*(x^4+2)ex