第4章MATLAB符号计算.pptVIP

第4章MATLAB符号计算;4.1符号函数的计算;使用sym函数定义符号变量;使用sym函数定义符号表达式;注意：

将表达式整体定义为符号表达式时，其中的变量并未被定义为符号变量。

f=sym(m*y^2+n*y)

f=

m*y^2+n*y

f1=2*m

???Undefinedfunctionorvariablem.;符号变量的基本操作;(1)findsym函数用于寻找符号变量;(2)digits设定所用数值的精度;任意精确度的符号表达式;(3)数值型变量与符号型变量的转换;t=0.75;

sym(t)%有理数形式

ans=

3/4

sym(t,f)%浮点数形式

ans=

1.8000000000000*2^(-1)

sym(t,d)

ans=

.75000000000000000000000000000000;symsx

s1=x^3-1;

s2=x-1;

s3=s1+s2

s3=

x^3-2+x

s4=s1-s2

s4=

x^3-x;化简函数;symsxt

f=x*(x*(x-8)+6)*t;

collect(f)

ans=

t*x^3-8*t*x^2+6*t*x

collect(f,t)

ans=

x*(x*(x-8)+6)*t;symsx

f=factor(x^3-1)

f=

(x-1)*(x^2+x+1)

s=sym(sin(a+b));

expand(s)

ans=

sin(a)*cos(b)+cos(a)*sin(b);symsxy

f=x/y-y/x;

[m,n]=numden(f)

m=

x^2-y^2

n=

y*x;symsx

f=sin(x)^2+cos(x)^2

f=

sin(x)^2+cos(x)^2

simplify(f)

ans=

1

simple(f);3．替换函数;symsab

t=solve((a+b)^3+(a+b)^2+6)

t=

-1/3*(82+9*83^(1/2))^(1/3)-1/3/(82+9*83^(1/2))^(1/3)-1/3-a

1/6*(82+9*83^(1/2))^(1/3)+1/6/(82+9*83^(1/2))^(1/3)-1/3+1/2*i*3^(1/2)*(-1/3*(82+9*83^(1/2))^(1/3)+1/3/(82+9*83^(1/2))^(1/3))-a

1/6*(82+9*83^(1/2))^(1/3)+1/6/(82+9*83^(1/2))^(1/3)-1/3-1/2*i*3^(1/2)*(-1/3*(82+9*83^(1/2))^(1/3)+1/3/(82+9*83^(1/2))^(1/3))-a

[r,s]=subexpr(t,s);

;用subs函数求值

y=subs(t,a,2)

y=

-4.2188

-1.3906-1.5274i

-1.3906+1.5274i

vpa(y,2)

ans=

-4.2

-1.4-1.5*i

-1.4+1.5*i;

symsx

s=(3*x^3+x^2-1)/(x^2+1);

r=subs(s,x,1)

r=

(3*(1)^3+(1)^2-1)/((1)^2+1)

r=subs(s,x,1)

r=

1.5000;三、可视化符号函数计算器;4.2符号微积分;例子：;二、符号求导;例子：;symsxy

int(x/(x+y),x)

ans=

x-y*log(x+y)

int(x/(x+y),y)

ans=

x*log(x+y)

int(x^2/(x+2),x,1,3)

ans=

4*log(5)-4*log(3)

double(ans)

ans=

2.0433;四、积分变换;symsxt

y=exp(-x^2);

Ft=fourier(y,x,t)

Ft=

pi^(1/2)*exp(-1/4*t^2)

fx=ifourier(Ft,t,x)

fx=

exp(-x^2);Lt=laplace(y,x,t)

Lt=

1/2*pi^(1/2)*exp(1/4*t^2)*erfc(1/2*t)

Lx=ilaplace(Lt,t,x)

Lx=

exp(-x^2)

Zt=ztrans(y,x,t)

Zt=

ztrans(exp(-x^2),x,t)

Zx=iztrans(Zt,t,