3.数值计算介绍.ppt

4.3 概率分布与统计分析 4.3.2 随机数组和统计分析 例题：随机数据的统计量 A=randn(1000,4); Amax=max(A); Amin=min(A); CM=mean(A); MA=mean(mean(A)); C=cov(A); P=corrcoef(A); 4.4 多项式运算 4.4.1 多项式运算函数 对于多项式:a(x)=a1xn+a2xn-1+……+anx+an+1,用行列式表示:a=[a1,a2,…, an, an+1] c=conv(a,b):a(x) ×b(x)乘积多项式的系数向量； [q,r]=deconv(b,a):b(x) /a(x) 多项式的商的系数向量q,余数的系数向量r [r,p,k]=residue(b,a):因式分解b(x) /a(x) =r1/(x-p1) +r2/(x-p2)+……rn/(x-pn)+k(x)，r-为分子系数向量，p-极点系数向量，k为余数系数向量 r=roots(a):多项式根 a=ploy(r):当r为一维数组，r为多项式根，a为多项式系数向量，当r为方阵，a表示矩阵r所对应的特征多项式的系数向量 Ploy2str(q,’x’)：将多项式系数向量q写成容易读的多项式 4.4 多项式运算 4.4.1 多项式运算函数 例：求（s2+2)(s+4)(s+1)/(s3+s+1)的商及余式 p1=conv([1,0,2],conv([1,4],[1,1])); p2=[1,0,1,1]; [q,r]=deconv(p1,p2) poly2str(q,’s’) poly2str(r,’s’) 4.4 多项式运算 4.4.2 多项式拟合 已知变量x,y之间的函数关系y(x)=a1xn+a2xn-1+……+anx+an+1,通过一组已知数据{xi,yi|i=1,2,…,m},确定系数[a1,a2,…, an, an+1]，指令为： p=polyfit(x,y,n):求x，y数组所给数据的n阶拟合多项式系数向量p 例：x0=0:0.1:1; y0=[-0.447,1.978,3.11,5.25,5.02,4.66,4.01,4.58,3.45,5.35,9.22]; n=3; p=polyfit(x0,y0,n) 4.4 多项式运算 4.4.2 最小二乘法 根据试验数据，给出一组方程： y1= a1x1n +a2x1n-1+……+anx 1+an+1 y2= a1x2n +a2x2n-1+……+anx2+an+1 … ym= a1xmn +a2xmn-1+……+anxm+an+1 上述方程简写：Y=XaT，已知Y,X,求a aT=X\Y 4.4 多项式运算 4.4.2 最小二乘法 例： x0=0:0.1:1; y0=[-0.447,1.978,3.11,5.25,5.02,4.66,4.01,4.58,3.45,5.35,9.22]’; n=3; m=length(x0); n=3; X=zeros(m,n+1); for k=1:n X(:,n-k+1)=(x0.^k); end X(:,n+1)=ones(m,1); a=(X\y0)’ 谢谢！ * * 授课教师：李军求 联系方式Email:lijunqiu@bit.edu.cn 2015年9月 工程软件训练 北京理工大学机械与车辆学院 上机安排 5位老师具体负责的班级 翟丽老师：1、2班；李军求老师：3、4班；吴维老师：5、6班；张虹老师：7、8班；张志强老师：9班 机房1(120人)：1班（26人）、2、3、4班 机房4(60人）：5、6班 机房5(120人)：7、8、9、1班（11人） 第一次上机时间：9月7日13:20~16:40 第一次上机内容：符号计算和数值计算 参考教材：Matlab教程2012a,张志涌等，北航出 内容章节 数值阵列及向量化运算 4 3 数值计算 4 4 数值阵列及向量化运算 主要内容： 矩阵、数组及向量化编程 数值数组的创建和寻访 关系操作和逻辑操作 3.1 数值计算的特点和地位 符号运算得到解析