Matla7.0符号计算 上机实验1.docVIP

Matlab7.0符号计算 实验报告 ⒈目的 本实验旨在向学生介绍一种解决专业问题的快速有效且具有强大功能的科学与工程计算软件。在处的值，并将结果设置为以下五种精度，即分别为小数点之后1位、2位、10位、20位、50位有效数字； x=pi-1; y=x+cos(x)-sin(x) y = 0.7598 y1=vpa(y,1) y1 = .8 y2=vpa(y,2) y2 = .76 y3=vpa(y,10) y3 = .7598193629 y4=vpa(y,20) y4 = .75981936291375673509 y5=vpa(y,50) y5 = .75981936291375673508952104384661652147769927978516 二、设为符号变量，，，试进行如下运算。 （1）， （2）， （3）对进行因式分解， （4）求的反函数。 syms x f g; f=x^4+2*x^2+1; g=x^3+6*x^2+3*x+5; （1） a=f+g a = x^4+8*x^2+6+x^3+3*x （2） b=f*g b = (x^4+2*x^2+1)*(x^3+6*x^2+3*x+5) expand(b) ans = x^7+6*x^6+5*x^5+17*x^4+7*x^3+16*x^2+3*x+5 （3） factor(f) ans = (x^2+1)^2 （4） finverse(g) Warning: finverse(x^3+6*x^2+3*x+5) is not unique. In sym.finverse at 43 ans = 1/2*(-60+4*x+4*(117-30*x+x^2)^(1/2))^(1/3)+6/(-60+4*x+4*(117-30*x+x^2)^(1/2))^(1/3)-2 三、指出下面各条语句的输出结果 f=sym(2*u); subs(f,u,2) f2=subs(f,u,u+2) a=3; subs(f2,u,a+2) subs(f2,u,a+2) syms x y f3=subs(f, u ,x+y) subs(f3,[x,y],[1,2]) subs(f3,[x,y],[x+y,x+y]) f=sym(2*u); subs(f,u,2) ans = 4 f2=subs(f,u,u+2) f2 = 2*(u+2) a=3; subs(f2,u,a+2) ans = 14 subs(f2,u,a+2) ans = 2*((a+2)+2) syms x y f3=subs(f,u,x+y) f3 = 2*x+2*y subs(f3,[x,y],[1,2]) ans = 6 subs(f3,[x,y],[x+y,x+y]) ans = 6*x+4*y 四、简化 syms x; f=((1/x^3)+(6/x^2)+12/x+8)^(1/3)； simplify(f) ans = ((2*x+1)^3/x^3)^(1/3) 五、设a、b定义如下，试上机输出factor(a)和factor(b)的结果，并指出哪个结果才是12345678901234567890的因式分解，为什么？ a=sym(12345678901234567890) b=sym(12345678901234567890) y1=factor(a) y2=factor(b) y1 = (2)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541) y2 = ^11*(7)*(7324703)*(117570121) a=sym(12345678901234567890) a = 12345678901234567890 b=sym(12345678901234567890) b = 6028163525993441*2^(11) y1=factor(a) y1 = (2)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541) y2=factor(b) y2 = (2)^11*(7)*(7324703)*(117570121) fact