资源描述:
《第06章_matlab数值计算_例题源程序汇总》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、技术资料第6章MATLAB数值计算例6.1求矩阵A的每行及每列的最大和最小元素,并求整个矩阵的最大和最小元素。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,
2、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]%产生一个二维矩阵xy1=std(x,0,1)y2=std(x,1,1)y3=std(x,0,2)y4=std(x,1,2)例6.5生成满足正态分布的10000×5随机矩阵,然后求各列元素的均值和标准方差,再求这5列随机数据的相关系数矩阵。X=randn(10000,5);M=mean(X)D=s
3、td(X)R=corrcoef(X)共享知识技术资料例6.6对下列矩阵做各种排序。A=[1,-8,5;4,12,6;13,7,-13];sort(A)%对A的每列按升序排序-sort(-A,2)%对A的每行按降序排序[X,I]=sort(A)%对A按列排序,并将每个元素所在行号送矩阵I例6.7给出概率积分的数据表如表6.1所示,用不同的插值方法计算f(0.472)。表6.1概率积分数据表x0.460.470.480.49f(x)0.48465550.49375420.50274980.5116683x=0.46:0.01:0.49;%给出x,f(x)f=[0.4846
4、555,0.4937542,0.5027498,0.5116683];formatlonginterp1(x,f,0.472)%用默认方法,即线性插值方法计算f(x)interp1(x,f,0.472,'nearest')%用最近点插值方法计算f(x)interp1(x,f,0.472,'spline')%用3次样条插值方法计算f(x)interp1(x,f,0.472,'cubic')%用3次多项式插值方法计算f(x)formatshort例6.8某检测参数f随时间t的采样结果如表6.2,用数据插值法计算t=2,7,12,17,22,17,32,37,42,47,5
5、2,57时的f值。表6.2检测参数f随时间t的采样结果t051015202530f3.10252.256879.51835.92968.84136.25237.9t35404550556065f6152.76725.36848.36403.56824.77328.57857.6T=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=interp1(T,F,X)%
6、用线性插值方法插值F1=interp1(T,F,X,'nearest')%用最近点插值方法插值F1=interp1(T,F,X,'spline')%用3次样条插值方法插值F1=interp1(T,F,X,'cubic')%用3次多项式插值方法插值例6.9设z=x2+y2,对z函数在[0,1]×[0,2]区域内进行插值。x=0:0.1:1;y=0:0.2:2;[X,Y]=meshgrid(x,y);%产生自变量网格坐标Z=X.^2+Y.^2;%求对应的函数值interp2(x,y,Z,0.5,0.5)%在(0.5,0.5)点插值interp2(x,y,Z,[0.50.6
7、],0.4)%在(0.5,0.4)点和(0.6,0.4)点插值interp2(x,y,Z,[0.50.6],[0.40.5])%在(0.5,0.4)点和(0.6,0.5)点插值%下一命令在(0.5,0.4),(0.6,0.4),(0.5,0.5)和(0.6,0.5)各点插值interp2(x,y,Z,[0.50.6]',[0.40.5])例6.10某实验对一根长10米的钢轨进行热源的温度传播测试。用x表示测量点(米),用h表示测量时间(秒),用T表示测得各点的温度(℃),测量结果如表6.2所示。表6.3钢轨各点温度测量值Txh02.557.51009