第四章 数值计算.pptVIP

4.3.1随机数发生器和统计分析指令函数含义rand(m,n)产生mxn维的[0,1]区间均匀分布随机数组．randn(m,n)产生mxn维的均值为0标准差为1的正态分布随机数组min(x)对矩阵各列分别求最小值max(x)对矩阵各列分别求最大值median(x)对矩阵各列分别求中位数mean(x)对矩阵各列分别求均值std(x)对矩阵各列分别求标准差var(x)对矩阵各列分别求方差C=conv(x)给出矩阵各列的协方差阵P=corrcoef(x)给出矩阵各列间的相关系数*第30页，共50页，星期日，2025年，2月5日例4.3-1基本统计示例randn(state,0);A=randn(1000,4);AMAX=max(A),AMIN=min(A)AMED=median(A);AMEAN=mean(A);ASTD=std(A)已知A=[1,2,3;4,5,6;7,8,9],演示max,min,median,mean,std,var,conv,corrcoef各列最大,最小值：max(A),min(A)各行最大,最小值：max(A),min(A)各列中间值：median(A)各列平均值：mean(A)各列方差,标准差：var(A),std(A)A=[1,2,3;4,5,6;7,8,9]789*第31页，共50页，星期日，2025年，2月5日4.3.2差分和累计指令函数含义del2(U,hx,hy)五点Laplaciandiff(X,m,n)沿第n维求m次差分(二维：n=1;列,n=2,行)[DZx,DZy]=gradient(z,dx,dy)对Z求x,y方向梯度prod(X,n)沿第n维求乘积sum(X,n)沿第n维求和trapz(x,Y,n)用梯形法沿第n维求函数关于自变量x的积分cumprod(X,n)沿第n维求累计乘积cumsum(X,n)沿第n维求累计和cumtrapz(x,Y,n)用梯形法沿第n维求函数Y关于x的累计积分[XS,KK]=sort(X,n)沿第n维对X元素按模增大排列*第32页，共50页，星期日，2025年，2月5日diff,trapz,cumtrapz指令的演示a=[1,2,3,4,5];b=[a,a+1,a+2,a+3]diff(a)=diff(a,1)=[1,1,1,1]diff(b)=[1,1,1,1;1,1,1,1;1,1,1,1;1,1,1,1]trapz(b)=[12162024]cumtrapz(b)00001.502.503.504.504.006.008.0010.007.5010.5013.5016.5012.0016.0020.0024.0012342345345645675678*第33页，共50页，星期日，2025年，2月5日例4.3.2-1求1+2+3+…+100以及50!x=1:1:100;sum_x=sum(x);sum_x=5050a=1:1:50;prod_a=prod(a)prod_a=3.0414e+064cumsum,cumprod的用法上两个函数分别是求向量的累计和和累计乘积.如果a=1:1:n,则cumsum(a)=[1,1+2,1+2+3,...,1+2+3+…+n]cumprod(a)=[1,1*2,1*2*3,…,1*2*3*…*n]*第34页，共50页，星期日，2025年，2月5日例4.3.2-1求f(x)=3x2在区间[0,2]的积分dt=0.001;t=(0:dt:2);y=3*t.^2;s1=dt*sum(y)s2=dt*trapz(y)s=dt*cumsum(y);s3=s(end)s