7.微积分new概论.pptVIP

Matlab实验与分析 应用微积分 把握知识结构和内容复习 ①对每章、每节的知识点、例题、难点等能出口成章地说出来 ②将试题内容再返到上述结构上 上一讲内容 主要内容 预备知识：微积分 主要内容 数值微积分MATLAB指令 数值差分示例 x=[1 3 8 7];diff(x),diff(x,2) ans=2 5 -1 ans=3 -6 %diff(x,2)= diff(diff(x)) A=[1 3;5 2;6 5;7 7];diff(A) ans =4 -11312 数值导数示例 数值求导可处理成差分的商 clear;x=[1 1.1 1.2 1.3];y=x.^3; dy=diff(y)./diff(x) dy= 3.3100 3.9700 4.6900 dy=gradient(y,x) dy= 3.3100 3.6400 4.3300 4.6900 3*x.^2 ans= 3.0000 3.6300 4.3200 5.07 解上一例 quad(inline(exp(-x.^2)),-1,1) ans = 1.4936 quadl(@(x)exp(-x.^2),-1,1) ans = 1.4936 假奇异积分? 计算数值积分，如果用梯形法计算 clear;x=-1:0.01:1;y=x.^(1/3); trapz(x,y) 结果为 ???ans =1.1241 + 0.6490i是一个复数。实际上，这个错误是由于数值方法的特点造成的，数值方法对x^(1/3)是通过对数exp(ln(x)/3)计算，就会出现复数，这种情况称为假奇异积分。这类情况在适当定义被积函数后可用quad,quadl求解。先编写M函数 function y=fun(x) y=x.^(1/3); if x0, y=-(-x).^(1/3); end 保存后,在命令窗口用MATLAB代码: ? quadl(@fun,-1,1) ans = 1.2361e-006 +7.1365e-007i fun=inline(x.*exp(x.^2+y.^2),x,y); dblquad(fun,0,2,-2,2) ans = 881.8304 fun=inline(y.*sin(t)+z.*cos(t),t,y,z); triplequad(fun,0,pi,0,1,-1,1) ans =2.0000 主要内容 计算实验：数值微积分 中心差商精度较高 Diff采用向前差商，gradient采用中心差商，从而gradient内点误差较小 高阶导数也可用差商法得到 function d=deriv(fname,a,h0,e) %用途：中心差商公式求函数f(x)在某一点x=a的导数 %格式：d=deriv(fname,a,e) d返回数值解,fname表示f(x)， %h0为初始步长，e为精度要求。 %首先取一个初始步长，然后逐次减半，直到相邻两次计算结果误差 %小于误差限 h=h0; d=(feval(fname,a+h)-feval(fname,a-h))/2/h; d0=d+2*e; while abs(d-d0)ed0=d;h0=h;h=h0/2;d=(feval(fname,a+h)-feval(fname,a-h))/2/h; end 简单应用 deriv(@(x)sin(x),pi/6,0.1,1e-4) 2.非矩形区域上的重积分 重积分的数值计算可通过单积分组合计算 function S=dblquad2(fun,a,b,clo,dhi,n) %非矩形区域重积分 if nargin6,n=100;end x=linspace(a,b,n+1); S=0; for i=1:n S=S+dblquad(fun,x(i),x(i+1),feval(clo,(x(i)+x(i+1))/2),feval(dhi,(x(i)+x(i+1))/2)); end clear; fun=inline(sqrt(1-x.^2),x,y); clo=inline(-sqrt(1-x^2)); dhi=inline(sqrt(1-x^2)); dblquad2(fun,-1,1,clo,dhi,50) 主要内容 建模实验 例 考虑函数 f(x)=x2cos(x2+3x-4) 在[-2, 2]内的图象特征。 fun=inline(x.*x.*cos(x.*x+3*x-4));fplot(fun,[-2 2]);grid on;hold on; %-0.5，1.5附近局部极小，1附近局部极大，-0.5~0.5不明朗 h=0.01;x=-2:h:2;df=diff(feval(fun,x))/h; plo