实验08数据处理及多项式计算第6章.docxVIP

精品文档 精品文档 PAGE PAGE21 精品文档 PAGE 附参考答案： 《MATLAB软件》课内实验 王平 实验08数据处理与多项式计算 （第6章MATLAB数值计算） 一、实验目的 掌握数据统计和分析的方法。 掌握数值插值与曲线拟合的方法及其应用。 掌握多项式的常用运算。 二、实验内容 检验一组随机数的性质 利用MATLAB提供的rand函数生成30000个符合均匀分布的随机数，然后检验随机数的性质：A=rand(30000,1); 均值和标准方差。mean(A),std(A), 最大元素和最小元素。max(A),min(A), (3)大于0.5的随机数个数占总数的百分比。 sum(A0.5)/30000, 命令及运行结果（建议在命令窗口中逐条输入命令） ： 学生成绩表处理 将100个学生5门功课的成绩存入矩阵P中，进行如下处理： P=round(45+(95-45)*rand(100,5)); 分别求每门课的最高分、最低分及相应学生序号。[M5,K5]=max(P),[m5,k5]=min(P), (2) 分别求每门课的平均分和标准方差。mean(P),std(P), (3) 5门课总分的最高分、最低分及相应学生序号。[M1,K1]=max(sum(P,2)), [m1,k1]=min(sum(P ,2)), (4) 将5门 课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh。 [zcj,xsxh]=sort(sum(P ,2),descend); 提示：上机调试时，为避免输入学生成绩的麻烦，可用取值范围在[45,95]之间的随机矩阵来表示学生成绩。 命令及运行结果（建议在命令窗口中逐条输入命令） ： P=round(45+(95-45)*rand(100,5)); [M5,K5]=max(P)%(1) M5= 94 95 94 95 93 K5= 12 18 26 59 16 [m5,k5]=min(P) m5= 46 45 46 46 45 k5= 99 23 44 33 48 disp(mean(P));% (2) 71.4300 68.3300 70.2300 68.9100 69.6500 disp(std(P)); 14.8720 14.4998 13.7340 14.5137 12.8844 S=sum(P,2);%(3) [M1,K1]=max(S) M1= 401 K1= 85 [m1,k1]=min(S) m1= 268 k1= 35 [zcj,xsxh]=sort(S,descend);%(4) disp([zcj,xsxh]); 85 68 93 15 97 26 三次样条插值应用 某气象观测得某日 6:00~18:00 之间每隔 2h的室内外温度（ 0C）如实验表 1所示。 实验表 1室内外温度观测结果（ 0C） 时间h6 8 10 12 14 16 18 室内温度t1 18.0 20.0 22.0 25.0 30.0 28.0 24.0 室外温度t2 15.0 19.0 24.0 28.0 34.0 32.0 30.0 试用三次样条插值分别求出该日室内外 6:30~17:30 之间每隔2h各点的近似温度（0C）。 程序： clc; h=6:2:18; t1=[18.0,20.0,22.0,25.0,30.0,28.0,24.0]; t2=[15.0,19.0,24.0,28.0,34.0,32.0,30.0]; H=6.5:2:17.5; T1=spline(h,t1,H); %interp1(h,t1,H, spline) T2=spline(h,t2,H); %interp1(h,t2,H, spline) disp(num2str([H;T1;T2])); 运行结果： 曲线拟合应用 已知lgx在[1,101]区间  10个整数采样点的函数值如实验表  2所示。 实验表2lgx 在10个采样点的函数值 x 1 11 21 31 41 51 61 71 81 91 101 lgx0 1.04141.32221.49141.61281.70761.78531.85131.90851.95102.0043 试求lgx  的5次拟合多项式  p(x)，并绘制出  lgx  和  p(x)在[1,101]  区间的函数曲线。 程序： clc; x=1:10:101; lgx=[0,1.0414,1.3222,1.4914,1.6128,1.7076,1.7853,1.8513,1.9085,1.9510,2.0043]; P=polyfit(x,lgx,5); %求5次拟合多项式的系数和常数项 X=linspace(1,101,100); %变量取值 lgX=log10(