mat的lab习题.ppt

* * 八、符号求积 * * 符号积分 积分分为定积分和不定积分。运用函数int可以求得符号表达式的积分，即找出一个符号表达式F，使得diff(F)=f，也可以说是求微分的逆运算。 语法： int(f, t) %求符号变量t的不定积分 int(f, t,a,b) %求符号变量t的积分 int(f, t, m, n) %求符号变量t的积分 * * 例 求下列积分。 (2) (4) 命令如下： x=sym(x);t=sym(t); int(abs(1-x),1,2) %求(1) 1/2 f=1/(1+x^2); int(f,-inf,inf) %求(2) pi f=x^3/(x-1)^10; I=int(f,2,3) %求(3) 138535/129024 int(4*x/t,t,2,sin(x)) %求(4) 4*log(sin(x))*x-4*log(2)*x * * 【例3.16】 求 和 的积分。 f=sym(cos(x)); int(f) %求不定积分 ans = sin(x) int(f,0,pi/3) %求定积分 ans = 1/2*3^(1/2) int(f,a,b) %求定积分 ans = sin(b) ?sin(a) int(int(f)) %求多重积分 ans = ?cos(x) diff和int命令也可以直接对字符串f进行运算。 f=cos(x); * * 【例3.16续】 求符号矩阵 的积分。 syms t x g=[2*x t^2;t*sin(x) exp(x)] %创建符号矩阵 g = [ 2*x, t^2] [ t*sin(x), exp(x)] int(g) %对x求不定积分 ans = [ x^2, t^2*x] [?t*cos(x), exp(x)] int(g,t) %对t求不定积分 ans = [ 2*x*t, 1/3*t^3] [ 1/2*t^2*sin(x), exp(x)*t] int(g,sym(a),sym(b)) %对x求定积分 ans = [ b^2 - a^2, -t^2*(a - b)] [ t*(cos(a) - cos(b)), exp(b) - exp(a)] * * * 例 求解积分问题 MATLAB求解 * * 八、符号求积 1 符号求积 a)1/2 b) 1 c) 1/4 d) 1/3 （c） * * 2 符号求积 a)pi b)2*pi c)3*pi d)1/2*pi （a） * * 3符号求积 a) 3.20808 b) 3.10808 c) 3.30808 d) 3.40808 （b） * * 九、方程符号求解 * * 符号方程的求解 代数方程 当方程不存在解析解又无其他自由参数时，MATLAB可以用solve命令给出方程的数值解。 语法： solve(eq, v) %求方程关于指定变量的解 solve(eq1, eq2, v1, v2,…) %求方程组关于指定变量的解 【例】 求方程ax2+bx+c=0和sinx=0的解。 f1=sym(a*x^2+b*x+c) %无等号 f1 = a*x^2+b*x+c solve(f1) %求方程的解x ans = -(b + (b^2 - 4*a*c)^(1/2))/(2*a) -(b - (b^2 - 4*a*c)^(1/2))/(2*a) f2=sym(sin(x)) f2 = sin(x) solve(f2,x) ans = 0 * * 【例】 求三元非线性方程组 的解。 eq1=sym(x^2+2*x+1); eq2=sym(x+3*z=4); eq3=sym(y*z=?1); [x,y,z]=solve(eq1,