工程应用数学8-5.ppt

§8.5 MATLAB与积分运算 §8.5 MATLAB与积分运算 数学表达式 MATLAB命令格式 不定积分 int(f(x)) 或int(f(x),x) 注意：积分表达式只是函数f(x)的一个原函数，后面没有带任意常数 C 定积分 int(f(x),a,b) 或int(f(x),x,a,b) 对函数f(x)中指定的符号变量x计算从a到b的定积分 * 例8-60 求不定积分 解：在MATLAB中输入以下命令： syms x int(x^4/(1+x^2)) 按回车键，显示结果为： ans =1/3*x^3-x+atan(x) §8.5 MATLAB与积分运算 * 例8-61 求不定积分 解：在MATLAB中输入以下命令： syms x int(1/((sin(x)^2)*(cos(x)^2))) ans = 1/sin(x)/cos(x)-2/sin(x)*cos(x) §8.5 MATLAB与积分运算 * 例8–62 求不定积分 (a0) 解：在MATLAB中输入以下命令： syms x a int(1/((x^2+a^2)^(1/2))) 按回车键，显示结果为： ans = log(x+(x^2+a^2)^(1/2)) §8.5 MATLAB与积分运算 * 例8–63 计算定积分解： 在MATLAB中输入以下命令： syms x k=int(sqrt(sin(x)^3-sin(x)^5),0,pi)； simplify(k) 按回车键，显示结果为： ans = 4/5 §8.5 MATLAB与积分运算 * 例8–64 计算定积分 解：在MATLAB中输入以下命令： syms x int(x^2*log((1-x)/(1+x)),-1/2,1/2) 输出结果为： ans =0 §8.5 MATLAB与积分运算 * 例8–65 计算由心形线r=a(1+cosθ),(a 0) 所围成的面积. 解：在MATLAB中输入以下命令： syms theta a int(a^2*(1+cos(theta))^2,theta,0,pi) 按回车键，显示结果为： ans = 3/2*pi*a^2 §8.5 MATLAB与积分运算 * 例8–66 一圆柱形水池，直径是20米，高10米，内盛满了水，现将水抽尽，求克服重力所作的功． 解：在MATLAB中输入以下命令： syms x g int(10^5*pi*g*x,0,10) 按回车键，显示结果为： ans = 5000000*pi*g §8.5 MATLAB与积分运算 * 例8–67 设长为l，质量均匀分布的杆，在杆的一端的延长线上距端点为a的位置有一质量为m的质点．求杆与质点之间的引力． 解：引力微元 因此引力 §8.5 MATLAB与积分运算 * 在MATLAB中输入以下命令： syms mu k m x a l ; int((mu*k*m)/(a+l-x)^2,x,0,l) 按回车键，显示结果为： ans = mu*k*m*l/a/(a+l) 再输入以下命令，可令结果显示得更清晰 pretty(ans) §8.5 MATLAB与积分运算 *