ch2 - 符号计算2010a(3)

syms a b x; f=a*sin(x)+b f =

b + a*sin(x)

（2）

f1=subs(f,sin(x),'log(y)') class(f1)

f1 =

b + a*log(y) ans = sym

（3）

f2=subs(f,a,3.11) class(f2)

f2 =

b + (311*sin(x))/100 ans = sym

（4）

f3=subs(f,{a,b,x},{2,5,sym('pi/3')})

class(f3) f3 =

3^(1/2) + 5 ans = sym

（5）

format format compact f4=subs(f,{a,b,x},{2,5,pi/3})

class(f4) f4 =

6.7321 ans = double

（6）

f5=subs(f,x,0:pi/2:pi) class(f5) f5 =

[ b, a + b, b] ans = sym

（7）

t=0:pi/10:2*pi;

f6=subs(f,{a,b,x},{2,3,t}) plot(t,f6) f6 =

Columns 1 through 8

3.0000 3.6180 4.1756 4.6180 4.9021 5.0000 4.9021 4.6180

Columns 9 through 16

4.1756 3.6180 3.0000 2.3820 1.8244 1.3820 1.0979 1.0000

Columns 17 through 21

1.0979 1.3820 1.8244 2.3820 3.0000

11

54.543.532.521.5101234567图2.2-1

（8）

k=(0.5:0.1:1)';

f6=subs(subs(f,{a,b},{k,2}),x,t); size(f6) plot(t,f6) ans =