MATLAB第3章MATLAB数值计算-复习.pptVIP

* ZHU Jun College of Physics Science Technology YANGZHOU UNIVERSITY CHINA 第3 章 MATLAB 数值运算 多项式 多项式的表达和创建 polyval(p, x)：计算x点中多项式值 roots(a)：求多项式的根 插值和拟合 interp1(x,y,xi,spline) p=polyfit(x,y,n) 请分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估测 X=9.5 时 Y 的值。 x=1:10 y=[16 32 70 142 260 436 682 1010 1432 1960] p2=polyfit(x,y,2) yp2=polyval(p2,9.5) p3=polyfit(x,y,3) yp3=polyval(p3,9.5) yi1=interp1(x,y,9.5) yi2=interp1(x,y,9.5,spline) 数值微分 diff函数：计算两个相邻点的差值 dy(x)/dx的数值微分为dy=diff(y)./diff(x)。 50 35 27 20 15 12 8 5 3 0 v 20 15 12 9 7 5 3 2 1 0 t t=[0 1 2 3 5 7 9 12 15 20]; v=[0 3 5 8 12 15 20 27 35 50]; a=diff(v)./diff(t) plot(t,v,-*,t(1:length(t)-1),a,-o) 50 35 27 20 15 12 8 5 3 0 v 20 15 12 9 7 5 3 2 1 0 t t=[0 1 2 3 5 7 9 12 15 20]; v=[0 3 5 8 12 15 20 27 35 50]; ti=0:20; vi=interp1(t,v,ti,spline); a=diff(vi)./diff(ti) plot(ti,vi,-*,ti(1:length(ti)-1),a,-o) 3.3.2 数值积分 cumsum(矩形积分)：cumsum(x)*h trapz(梯形积分)：z=trapz(x,y) quad(辛普森积分)：q=quad(f,a,b) quadl(科茨积分，也称高精度数值积分)：q=quadl(f,a,b) 4种近似方法的精度由低而高，和 trapz比较，quad、quadl不同之处在于这两者类似解析式的积分式，只需设定上下限及定义要积分的函数；而 trapz 是针对离散点数据做积分。 50 35 27 20 15 12 8 5 3 0 v 20 15 12 9 7 5 3 2 1 0 t t=[0 1 2 3 5 7 9 12 15 20]; v=[0 3 5 8 12 15 20 27 35 50]; s=trapz(x,y) x=linspace(0,1,100); y=x./(x.^2+4); t=cumsum(y*(1/99)); t(99) x=linspace(0,1,100); y=x./(x.^2+4); t=trapz(x,y) t=quad(x./(x.^2+4),0,1) t=quadl(x./(x.^2+4),0,1) a=[1/2 1/3 1;1 5/3 3;2 4/3 5]; b=[1;3;2]; c=a\b c = 4 3 -2 3.4 线性方程组 3.5 稀疏矩阵 S=sparse(A)：将一个满矩阵A 转化为一个稀疏矩阵S。若S本身就是一个稀疏矩阵，则sparse(S)返回S。 S=sparse(i,j,s,m,n)：在第i 行、第j 列输入数值s，矩阵共m 行n 列，输出S 为一个稀疏矩阵，给出(i,j)及s。 S=sparse(i,j,s)：比较简单的格式，只输入非零元的数据s 以及各非零元的行下标i和列下标j。 S=sparse(m,n)：是sparse([],[],[],m,n,0)的省略形式，用来产生一个m×n 的全零矩阵。