第二讲MATLAB的数值计算幻灯片.pptVIP

逻辑关系函数 any 向量的任意元素不为0则返回真 all 向量的所有元素不为0则返回真 xor 逻辑或非 isempty 判断空矩阵 isequal 判断相等矩阵 isnumeric 判断数值矩阵 islogical 判断逻辑数组 logical 转换数值为逻辑型 isnan 判断不定数 isinf 判断无穷大元素 isfinite 判断有限大元素 find 寻找非零元素坐标 a=magic(5); a(:,3)=zeros(5,1) a = 17 24 0 8 15 23 5 0 14 16 4 6 0 20 22 10 12 0 21 3 11 18 0 2 9 a1=all(a(:,1)10) a1 = 0 a2=all(a3) a2 = 1 1 0 0 0 a11=any(a(:,1)10) a11 = 1 a22=any(a10) a22 = 1 1 0 1 1 a=[1:5]; a=1./a a = 1.0000 0.5000 0.3333 0.2500 0.2000 f1=find(a) f1 = 1 2 3 4 5 f1=find(abs(a)0.4|abs(a)0.23) f1 = 1 2 5 积分与微分 Newton-Cotes系列数值求积 将积分区间[a,b]划分为n等分，步长h=(b-a)/n，选取等距节点xk=a+kh，则求积分公式： Ck(n)由Cotes系数表给出。 n=0 时为矩形公式； n=1时为梯形公式; n=2时为Simpson公式； n=3时为Cotes公式； … 矩形求积公式 Cumsum 元素累积和函数 Cumsum(x) 对于向量x，函数返回一个向量，向量的第n个元素为向量x的前n个元素的和；对于矩阵X，函数返回一个矩阵，矩阵的列对应X的列的累积和返回值 Cumsum(x,DIM) 参数DIM指明求和从第DIM维开始。 例：x=[0 1 2;3 4 5]; cumsum(x) ans = 0 1 2 3 5 7 梯形求积分公式 Z=trapz(Y) 用梯形方法计算Y的积分近似值。对向量Y，函数返回Y的积分；对矩阵Y，函数返回一个向量，向量各元素为矩阵Y对应列向量的积分值。 Z=trapz(X,Y) Z=trapz(X,Y,DIM) Z=trapz(Y,DIM) 例：求积分 d=pi/1000; t=0:d:3*pi; nt=length(t); y=fun(t); sc=cumsum(y)*d; scf=sc(nt) scf = 0.901618619310013 z=trapz(y)*d z = 0.900840276606884 建立被积函数： Fun.m function y=fun(t) y=exp(-0.5*t).*sin(t+pi/6); 自适应法-simpson法求积 Quad Q=quad(‘F’,a,b) 求函数F(x)从a到b的相对误差为1e-3的积分近似值； Q=quad(‘F’,a,b,tol) 向量tol给出相对误差和绝对误差； Q=quad(‘F’,a,b,tol,TRACE) TRACE不为零时给出图形； 例：求积分 建立被积函数： fun1.m function f=fun1(x) f=x./(4+x.^2); 在MATLAB命令窗口中： quad(fun1,0,1) ans = 0.111572382538912 自适应法-cotes法求积 Quad8 Q=quad8(‘F’,a,b) 求函数F(x)从a到b的相对误差为1e-3的积分近似值； Q=quad8(‘F’,a,b,tol) 向量tol给出相对误差和绝对误差； Q=quad8(‘F’,a,b,tol,TRACE) TRACE不为零时给出图形； 例：求积分 建立被积函数：fun2.m function f=fun2(x) f=exp(-x/2); 在M