第4MAT六B的数值计算教学提纲.pptVIP

第4章 MATLAB 的数值计算 ;4.1数值微积分;在MATLAB数值计算中，既没有专门的求极限指令，也没有专门的求导指令。但MATLAB提供了与“求导”概念有关的“求差分”指令。 dx=diff(X) %计算向量X的前向差分 FX=gradient(F) %求一元(函数)梯度 [FX, FY] =gradient(F) %求二元（函数）梯度;dx=diff(X) %计算向量X的前向差分 FX=gradient(F) %求一元(函数)梯度 [FX, FY] =gradient(F) %求二元（函数）梯度;数值极限和导数的应用应十分谨慎;数值极限和导数的应用应十分谨慎;d=pi/100; t=0:d:2*pi; x=sin(t); dxdt_diff=diff(x)/d; dxdt_grad=gradient(x)/d; ;Sx=sum(X) 　　 % sum按列向求和得(1×n)数组 Scs=cumsum(X) %沿X列向求累计和, 仍是数组, 第(i, k)个元素是X数组第k列前i个元素的和。最后一行等于Sx;clear; d=pi/8; t=0:d:pi/2; y=0.2+sin(t); s=sum(y); s_sa=d*s; % s_sa=sum(t, y), 近似矩形法积分 s_ta=trapz(t,y); %梯形法积分;4.1.3 计算精度可控的数值积分;（1） syms x Isym=vpa(int(exp(-x^2),x,0,1)) % x.^2数组乘方亦可 Isym = 0.74682413281242702539946743613185 ;（1）符号计算法 syms x y s=vpa(int(int(x^y,x,0,1),y,1,2)) Warning: Explicit integral could not be found. s = 0.40546510810816438197801311546435 ;4.1.4 函数极值的数值求解 ;【例4.1-7 】已知 ，在-10≦x≦10区间，求函数的最小值。;【例4.1-7 】已知 ，在-10≦x≦10区间，求函数的最小值。;【例4.1-7 】已知 ，在-10≦x≦10区间，求函数的最小值。;例4.1-8: f(x,y)=100(y-x2)2+(1-x)2在区间[-5,5]的极小值;4.1.5 常微分方程Ordinary Differential Equation的数值解;例4.1-9求解：;4.2矩阵和代数方程;4.2.1 矩阵运算和特征参数;format rat %有理格式显示 A=magic(2) + j*pascal(2) A = 1 + 1i 3 + 1i 4 + 1i 2 + 2i A1=A A1=1 - 1i 4 - 1i %共轭转置 3 - 1i 2 - 2i A2=A. A2=1 + 1i 4 + 1i %非共轭转置，数组运算操作 3 + 1i 2 + 2i;计算矩阵标量特征参数秩,迹,行列式的指令; A=reshape(1:9,3,3); r=rank(A) %求秩 d3=det(A) %非满秩矩阵的行列式一定为零 d2=det(A(1:2,1:2)) %求子式的行列式 t=trace(A);【例4.2-4】矩阵标量特征参数的性质。;dx=diff(X) %计算向量X的前向差分 FX=gradient(F) %求一元(函数)梯度 [FX, FY] =gradient(F) %求二元（函数）梯度;4.1.3 计算精度可控的数值积分;4.1.4 函数极值的数值求解 ;rank(A) %求秩（Rank） det(A) %求行列式（Determinant） trace(A)