06.08.2013 Views

内容简介作者简介 - 科学与工程计算系

内容简介作者简介 - 科学与工程计算系

内容简介作者简介 - 科学与工程计算系

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

exp(i*x)<br />

5.2.2 符号函数的求反和复合<br />

2<br />

【例 5.2.2-1】求 f = x 的反函数<br />

syms x;f=x^2;g=finverse(f)<br />

Warning: finverse(x^2) is not unique.<br />

> In D:\MATLAB6P5\toolbox\symbolic\@sym\finverse.m at line 43<br />

g =<br />

x^(1/2)<br />

fg=simple(compose(g,f)) %验算 g(f(x))是否等于 x<br />

fg =<br />

x<br />

x<br />

【例 5.2.2-2】求 f = , g = cos( y + fai)<br />

的复合函数<br />

2<br />

1+<br />

u<br />

(1)<br />

syms x y u fai t;f=x/(1+u^2);g=cos(y+fai);fg1=compose(f,g)<br />

fg1 =<br />

cos(y+fai)/(1+u^2)<br />

(2)<br />

fg2=compose(f,g,u,fai,t)<br />

fg2 =<br />

x/(cos(y+t)^2+1)<br />

5.2.3 置换及其应用<br />

5.2.3.1 自动执行的子表达式置换指令<br />

【例 5.2.3.1-1】演示子表达式的置换表示。<br />

clear all,syms a b c d W;[V,D]=eig([a b;c d]);<br />

[RVD,W]=subexpr([V;D],W) %<br />

RVD =<br />

[ -(1/2*d-1/2*a-1/2*W)/c, -(1/2*d-1/2*a+1/2*W)/c]<br />

[ 1, 1]<br />

[ 1/2*d+1/2*a+1/2*W, 0]<br />

[ 0, 1/2*d+1/2*a-1/2*W]<br />

W =<br />

(d^2-2*a*d+a^2+4*b*c)^(1/2)<br />

5.2.3.2 通用置换指令<br />

【例 5.2.3.2-1】用简单算例演示 subs 的置换规则。<br />

(1)<br />

syms a x;f=a*sin(x)+5;<br />

(2)<br />

f1=subs(f,'sin(x)',sym('y')) %<br />

f1 =<br />

a*y+5<br />

6

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!