应用统计软件MATLAB——05积分与符号运算汇编.ppt

MATLAB求积分与符号运算 学习 Matlab 命令 计算不定积分 定积分的概念 不定积分和广义积分 二重积分的计算 曲线积分 例3 . * * 授课教师：唐跃龙 学习内容： 1. 学习Matlab命令 1). 求和命令 sum 调用格式： sum(x) 给出向量x的各元素的累加和。若x为矩阵，则是一个元素为每列列和的行向量。 例1 x=[1,2,3,4,5,6,7,8,9,10]; sum(x) ans=55 例2 x=[1,2,3;4,5,6;7,8,9] x= 1 2 3 4 5 6 7 8 9 sum(x) ans=12 15 18 symsum(s,n) symsum(s,k,m,n) 当x的元素很有规律，例如为s(k)时，可用 symsum求得x的各项和，即 2). 求和命令 symsum 调用格式： symsum(s(k),1,n)=s(1)+s(2)+….+s(n) symsum(s(k),k,m,n)=s(m)+s(m+1)+….+s(n) syms k n symsum(k,1,10) ans=55 symsum(k^2,k,1,n) ans = 1/3*(n+1)^3-1/2*(n+1)^2+1/6*n+1/6 int(f(x)) int(f(x,y),x) 3). 积分命令 int 调用格式： int(f(x),a,b) int(f(x,y),x,a,b) 2. 计算不定积分 例4 解：输入命令： syms x; y=x^2*log(x); y1=diff(y); %对函数求导 y0=int(y1); %对导数积分的原函数 y2=int(y) y2 = 1/3*x^3*log(x)-1/9*x^3 例5 解：输入命令： syms x; syms a real; y=[sqrt(a^2-x^2),(x-1)/(3*x-1)^(1/3),x^2*asin(x)]; int(y,x) ans = [ 1/2*x*(a^2-x^2)^(1/2)+1/2*a^2*asin((1/a^2)^(1/2)*x), 1/15*(3*x-1)^(5/3)-1/3*(3*x-1)^(2/3), 1/3*x^3*asin(x)+1/9*x^2*(1-x^2)^(1/2)+2/9*(1-x^2)^(1/2)] syms a x; f=simple(int(x^3*(cos(a*x))^2)); f = 1/16*(4*a^3*x^3*sin(2*a*x)+2*x^4*a^4+6*x^2*a^2*cos(2*a*x)-6*a*x*sin(2*a*x)+3-3*cos(2*a*x))/a^4 f1=x^4/8+(x^3/(4*a)-3*x/(8*a^3))*sin(2*a*x)+ (3*x^2/(8*a^2)-3/(16*a^4))*cos(2*a*x); r=simple(f-f1) r = 3/16/a^4 3 定积分的概念 x=linspace(0,1,21); 定积分为一个和式极限，取f(x)=exp(x)，积分区间为[0,1]，等距划分为20个子区间 y=exp(x); 选取每个子区间的的端点，计算端点处的函数值 y1=y(1:20);s1=sum(y1)/20; 取区间的左端点乘以区间长度全部加起来 s1=1.6757 y2=y(2:21);s2=sum(y2)/20; 取区间的右端点乘以区间长度全部加起来 s2=1.7616 plot(x,y);hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i)],[0,0,y(i),y(i)],’r’) end for i=1:20 fill([x(i),x(i+1),x(i+1),x(i)],[0,0,y(i+1),y(i+1)],’b’) end 若取右端点，则 从图上可以看出： 当点取得越来越多时，s2-s1的值会越来越小，可试取50个点计算，看结果如何。下面按等分区间计算 syms k n; s=symsum(exp(k/n)/n,k,1,n); limit(s,n,inf) ans = exp(1)-1 结果与上面一样。 例7 解： 输入指令 syms x; I=int(exp(x),0,1) 得结果: I=exp(1)-1 4 计算定积分和广义积分 例8 解： 输入指令 syms x; int(abs(x-1),0,2) ans = 1. int( ) 还可以求函数边界的定积分问题。 syms x t; f=(-2*(x^2)+1)/(2*