Matlab语言应用.docVIP

Matlab语言应用 符号数学计算 1 符号表达式： 在Matlab中，命令sym定义符号变梁和符号表达式，例如： 例题1 x=sym(x) x = x a=sym(alpha) a = alpha rho=sym((1+sqrt(5))/2) rho = (1+sqrt(5))/2 f=rho^2-rho-1 f = (1/2+1/2*5^(1/2))^2-3/2-1/2*5^(1/2) simplify(f) ans = 0 命令syms可以用于简化多个符号变量的定义。对于，可以用以下的命令实现 syms a b c x f=a*x^2+b*x+c f = a*x^2+b*x+c 2符号表达式的计算 2．1 Matlab具有对符号表达式进行高级运算的能力。函数compose把f(x)和g(y)复合成f(g(x)),该函数有几种不同的使用方法： compose(f,s)返回f(g(y)),这里f=f(x),g=g(y); compose(f,g,z)返回f(g(z)),这里f=f(x),g=g(y); compose(f,g,x,z)返回f(g(z)),并指定x为f的独立变量； compose(f,g,x,y,z)返回f(g(z)),并指定x为f的独立变量，y为g的独立变量。 例题2 syms x y z t u; f=1/(1+x^2);g=sin(y);h=x^t;p=exp(-y/u); compose(f,g) ans = 1/(sin(y)^2+1) compose(f,g,t) ans = 1/(sin(t)^2+1) compose(h,g,x,z) ans = sin(z)^t compose(h,g,t,z) ans = x^sin(z) compose(h,p,x,y,z) ans = exp(-z/u)^t compose(h,p,t,u,z) ans = x^exp(-y/z) 2．2 函数finverse求表达式的逆函数；函数symsum对表达式求和 例题： syms x y z finverse(1/tan(x)) ans = atan(1/x) f=x^2+y; finverse(f,y) %指定y为x的独立变量 ans = -x^2+y finverse(f) %缺省独立变量为x Warning: finverse(x^2+y) is not unique. %由于逆函数不唯一，告警 In C:\MATLAB6P1\toolbox\symbolic\@sym\finverse.m at line 43 ans = (-y+x)^(1/2) 2．3 symsum函数求表达式的符号和有4种形式： symsum(f)返回； symsum(f,s)返回； symsum(f,a,b)返回； symsum(f,s,a,b)返回。 例题：求 syms n symsum((2*n-1)^2,1,n) ans = 11/3*n+8/3-4*(n+1)^2+4/3*(n+1)^3 factor(ans) ans = 1/3*n*(2*n-1)*(2*n+1) pretty(ans) 1/3 n (2 n - 1) (2 n + 1) % 结果是： 例题：求： syms n symsum(1/(2*n-1)^2,1,inf) ans = 1/8*pi^2 %结果是： 3．高等数学运算 3．1微分 符号表达式的微分以4种形式利用函数diff syms a b c x f=a*x^3+x^2-b*x-c %定义符号表达式 f = a*x^3+x^2-b*x-c diff(f) %对缺省变量微分 ans = 3*a*x^2+2*x-b diff(f,a) %对变量a微分 ans = x^3 diff(f,2) %对变量x进行二次微分 ans = 6*a*x+2 diff(f,a,2) %对变量a进行二次微分 ans = 0 3．2 积分 积分函数int(f)，其中f是一符号表达式，它力图求出符号表达式F, 使得diff(F)=f, 即返回f的不定积分。 int(f,v)