第06章_MATLAB数值计算_例题源程序汇总.docxVIP

第6章 MATLAB数值计算 例6.1求矩阵A的每行及每列的最大和最小元素，并求整个矩阵的最大和最小元素 _13 -56 78| 25 63 -235 A = 78 25 563 1 1 0 -1 A=[13,-56,78;25,63,-235;78,25,563;1,0,-1]; max(A,[],2) %求每行最大元素 min (A,[],2) %求每行最小元素 max(A) %求每列最大元素 min(A) %求每列最小元素 max(max(A)) %求整个矩阵的最大元素。也可使用命令： max(A(:)) min(min( A)) %求整个矩阵的最小元素。也可使用命令：min (A(:)) 例6.2求矩阵A的每行元素的乘积和全部元素的乘积。 A=[1,2,3,4;5,6,7,8;9,10,11,12]; S=prod(A,2) prod(S) %求A的全部元素的乘积。也可以使用命令 prod(A(:)) 例 6.3 求向量 X=(1 !, 2! , 3!,…，10!)。 X=cumprod(1:10) 例6.4对二维矩阵x,从不同维方向求出其标准方差。 x=[4,5,6;1,4,8] %产生一个二维矩阵 x y1=std(x,0,1) y2=std(x,1,1) y3=std(x,0,2) y4=std(x,1,2) 例6.5生成满足正态分布的100000随机矩阵，然后求各列元素的均值和标准方差， 再求这5列随机数据的相关系数矩阵。 X=ra ndn (10000,5); M=mea n(X) D=std(X) R=corrcoef(X) 例6.6对下列矩阵做各种排序。 1 -8 5 A= 4 12 6 13 7 13 7 -13 A=[1,-8,5;4,12,6;13,7,-13];sort(A)-sort(-A,2) [X,l]=sort(A) A=[1,-8,5;4,12,6;13,7,-13]; sort(A) -sort(-A,2) [X,l]=sort(A) 例6.7给出概率积分 %对A的每列按升序排序 %对A的每行按降序排序 %对A按列排序，并将每个元素所在行号送矩阵I f (x)=2 x f (x)= 2 x n oe dx 的数据表如表6.1所示，用不同的插值方法计算f(0.472) 表6.1概率积分数据表 x 0.46 0.47 0.48 0.49 f(x) 0.4846555 0.4937542 0.5027498 0.5116683 x=0.46:0.01:0.49; %给出 x, f(x) f=[0.4846555,0.4937542,0.5027498,0.5116683]; format long in terp1(x,f,0.472)% in terp1(x,f,0.472) %用默认方法，即线性插值方法计算f(x) in terp1(x,f,0.472, nearest)% in terp1(x,f,0.472, nearest) %用最近点插值方法计算f(x) in terp1(x,f,0.472,spli ne)in terp1(x,f,0.472,cubic)%用3次样条插值方法计算 in terp1(x,f,0.472,spli ne) in terp1(x,f,0.472,cubic) %用3次多项式插值方法计算f(x) format short 例6.8某检测参数f随时间t的采样结果如表6.2,用数据插值法计算t=2, 7, 12, 17, 22, 17, 32, 37, 42, 47, 52, 57 时的 f 值。 表6.2检测参数f随时间t的采样结果 t 0 5 10 15 20 25 30 f 3.1025 2.256 879.5 1835.9 2968.8 4136.2 5237.9 t 35 40 45 50 55 60 65 f 6152.7 6725.3 6848.3 6403.5 6824.7 7328.5 7857.6 T=0:5:65; X=2:5:57; F=[3.2015,2.2560,879.5,1835.9,2968.8,4136.2,5237.9,6152.7,... 6725.3,6848.3,6403.5,6824.7,7328.5,7857.6]; F1=i nterp1(T,F,X) %用线性插值方法插值 interp2(x,y,Z,[0.5 0.6],[0.4 0.5]) % 在(0.5,0.4)点和(0.6, 0.5)点插值 F1=in terp1(T,F,X, nearest)F1=in terp1(T,F,X,spli ne) F1=in terp1(T,F,X, nearest) F1=in terp1(T,F,X,spli