测控软件技术基础作业报告.doc

中国地质大学（武汉） 测控技术与仪器系测控软件技术基础作业报告一、题目一：1、代码：%内容 ： 生成一个1到30并且以3递进的数列，计算其长度、均值、最值和方差%作者 ： 柏明明 测控一班 20121002120%备注 ： 测控软件与技术基础第一次作业Num = 1:3:30;Length = length(Num) %长度Average = mean(Num) %均值 Max = max(Num) %最大值Min = min(Num) %最小值Var = sum(Num - mean(Num).2)/Length %方差Var1 = var(Num) %经检验用该公式计算结果与方差定义式结果不同var2 = (std(Num)2 %经检验用该公式计算结果与方差定义式结果不同2、 运行结果：Length = 10；Average = 14.5000；Max = 28；Min = 1；Var = 74.2500；Var1 = 82.5000；var2 = 82.5000。3、 结论：该向量的长度、均值、最值和方差结果如上所示，但有一点值得注意，在求方差时通过三种不同求解方法分别得到了三个结果，Var、Var1、var2。其中Var是严格用方差定义式求得的，Var1是用matlab自带的求方差函数var()求得的，Var2是用matlab自带的求标准差函数std()平方后求得的，Var1值等于Var2值而与Var值不等，说明matlab在求方差时不是按照方差定义式求解的。二、题目二：1、代码：%内容 ： 计算y = 3*x + 1 从1到4的积分syms a;x = 1:0.001:4;y = 3*x + 1;JiFen1 = sum(0.001*y) %用细分曲边梯形法求积分JiFen2 = int(3*a+1,1,4) %公式求积分2、 运行结果：JiFen1 = 25.5085；JiFen2 = 51/2。3、 结论：y = 3*x + 1 从1到4的积分结果如上所示，其中JiFen1是通过近似求曲边梯形的面积得到的，JiFen2是通过matlab自带的求定积分函数int()求得，两者结果近似相等，说明matlab在求定积分时也是按照求曲边梯形面积近似得到的，只不过它分割的更为细致，因而结果更为精确。3、 题目三：1、代码：%内容 ：编写程序，生成一份班级的学科成绩汇总表，并保存在Excel表格中。要求学生人数不少于40，每人至少6们课，成绩计分式包括百分制和五分制两种。Grade = 41 42 43 44 45 1;42 43 44 45 46 2;43 44 45 46 47 3;44 45 46 47 48 4;45 46 47 48 49 5;46 47 48 49 50 2;47 48 99 50 51 3;48 49 90 51 52 4;49 50 51 52 53 5;50 51 92 53 94 2;51 52 93 54 55 3;52 53 94 95 56 4;53 54 55 56 57 5;54 55 96 57 58 2;55 56 57 58 99 3;56 57 98 59 60 4;57 58 59 60 61 5;58 59 90 91 62 2;59 60 61 62 63 3;60 61 62 63 64 4;61 62 63 64 65 5;62 63 64 65 66 2;63 64 95 66 67 3;64 65 66 67 68 4;65 66 97 68 69 5;66 67 68 69 70 2;67 98 69 90 71 3;68 99 70 71 92 4;69 90 71 72 73 5;80 91 72 73 74 4;81 92 73 94 75 5;82 93 74 75 76 5;83 94 75 76 77 2;84 95 76 77 78 2;85 96 77 78 79 1;86 97 98 79 80 3;97 98 99 80 81 4;98 79 80 81 92 3;99 80 81 82 83 2;90 81 82 83 84 1; xlswrite(Grade.xls,Grade,A1:F40); %xlswrite(文件名，数组名,起始单元格：终止单元格)函数 2、 运行结果：见 附录。3、 结论：xlswrite(文件名，数组名,起始单元格：终止单元格)函数在使用时各参数分别如括号中所示，当没有该文件时会自动创建该文件并写入。4、 题目四：1、代码：%内容 ： 对上题存储在excel中的班级数据，统计成绩的最值和均值，以及不及格人数及大于90的人数，绘制直方图展示成绩分布特性 grade = xlsread(Grade.xls);Lesson1 = grade(:,1); %Lesson1成绩Lesson2 = grade(:,2); %Lesson2成绩Lesson3 = grade(:,3); %Lesson3成绩Lesson4 = grade(:,4); %Lesson4成绩Lesson5 = grade(:,5); %Lesson5成绩Lesson6 = grade(:,6); %Lesson6成绩 Max1 = max(Lesson1) %最大值Max2 = max(Lesson2)Max3 = max(Lesson3)Max4 = max(Lesson4)Max5 = max(Lesson5)Max6 = max(Lesson6) Min1 = min(Lesson1) %最小值Min2 = min(Lesson2)Min3 = min(Lesson3)Min4 = min(Lesson4)Min5 = min(Lesson5)Min6 = min(Lesson6) Mean1 = mean(Lesson1) %均值Mean2 = mean(Lesson2)Mean3 = mean(Lesson3)Mean4 = mean(Lesson4)Mean5 = mean(Lesson5)Mean6 = mean(Lesson6) figure(1); %直方图绘制hist(Lesson1);xlabel(X轴);ylabel(Y轴);title(Lesson1直方图展示); figure(2);hist(Lesson2);xlabel(X轴);ylabel(Y轴);title(Lesson2直方图展示); figure(3);hist(Lesson3);xlabel(X轴);ylabel(Y轴);title(Lesson3直方图展示); figure(4);hist(Lesson4);xlabel(X轴);ylabel(Y轴);title(Lesson4直方图展示); figure(5);hist(Lesson5);xlabel(X轴);ylabel(Y轴);title(Lesson5直方图展示); figure(6);hist(Lesson6);xlabel(X轴);ylabel(Y轴);title(Lesson6直方图展示); Ls1Not = length(find(Lesson160) %每门课不及格人数Ls2Not = length(find(Lesson260)Ls3Not = length(find(Lesson360)Ls4Not = length(find(Lesson460)Ls5Not = length(find(Lesson560)Ls6Not = length(find(Lesson690) %每门课优秀人数Ls2OK = length(find(Lesson290)Ls3OK = length(find(Lesson390)Ls4OK = length(find(Lesson490)Ls5OK = length(find(Lesson590)Ls6OK = length(find(Lesson690)2、运行结果：（1） 数值结果：a.各门成绩最大值：Max1 = 99；Max2 = 99；Max3 = 99；Max4 = 95；Max5 = 99；Max6 = 5。b.各门成绩最小值：Min1 = 41；Min2 = 42；Min3 = 43；Min4 = 44；Min5 = 45；Min6 = 1。c.各门成绩均值：Mean1 = 64；Mean2 = 67.5000；Mean3 =73；Mean4 =66.2500；Mean5 = 67.2500；Mean6 = 3.2750。d.各门成绩未通过人数：Ls1Not = 19；Ls2Not = 18；Ls3Not =10；Ls4Not =15；Ls5Not =13；Ls6Not =40。e.各门成绩优秀人数：Ls1OK =3；Ls2OK =10；Ls3OK =10Ls4OK =3；Ls5OK = 4；Ls6OK =0。（2）图像结果：图一 Lesson1直方图展示图二 Lesson2直方图展示图三 Lesson3直方图展示图四 Lesson4直方图展示图五 Lesson5直方图展示图六 Lesson6直方图展示3、 结论：对上题存储在excel中的班级数据，统计各门课成绩的最值和均值，以及不及格人数及大于90的人数如运行结果中数值结果所示，各门课成绩直方图展示如图像结果所示。5、 题目五：1、 代码：%内容 ： 绘制 y = x4 + 3*x + 1 x = -10:0.01:10;y = x.4 + 3*x + 1;figure(1);plot(x,y,r:);xlabel(X轴);ylabel(Y轴);legend(y=x4+3*x+1);title(y=x4+3*x+1曲线图);grid on;2、 运行结果：图七 y = x4+3*x+1 绘制结果3、 结论：y = x4 + 3*x + 1绘制结果如图七所示。要掌握绘制二维函数图像中常用的函数的用法，比如figure()、plot()、xlabel()、ylabel()、legend()、title()等等。6、 题目六：1、代码：%内容 ： fft测试：幅度为0.7频率为50HZ的信号与幅度为1频率为120HZ的信号叠加后的频率观察 Fs = 1000; %信号采样频率T = 1/Fs; %信号采样时间L = 1000; %信号长度t = (0:L-1)*T; %时间轴x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); %50HZ和120HZ信号的叠加y = x + 2*randn(size(t); %符合正太分布的噪声的引入figure(1); plot(Fs*t(1:50),y(1:50);title(引入正态分布噪声的信号的时域图像);xlabel(t(ms);ylabel(y(t); NEET = 2nextpow2(L); Y = fft(y,NEET)/L; f = Fs/2*linspace(0,1,NEET/2+1);figure(2);plot(f,2*abs(Y(1:NEET/2+1);title(引入正态分布噪声的信号的频域图像);xlabel(f（HZ）);ylabel(|Y(f)|);2、 运行结果：图八 引入正态分布噪声的信号的时域图像图九 引入正态分布噪声的信号的频域图像3、 结论：a.通过对图九的观察可以发现在50HZ和120HZ左右明显有两个峰值，且其幅度跟设置幅度相近，说明该傅里叶变换可以反映所设置的时域信号的频域特性。b.在仿真时域信号时，常引入随机数rand()或randn()做为噪声。7、 题目七：1、代码：%内容： 统计一下班上同学的来源管理，用饼状图表示出来。 filename = 7.xls;note = 01,02,03,04,05,06,07,08,09,10,. 11,12,13,14,15,16,17,18,19,20,. 21,22,23,24,25,26,27,28,29,30,. 31,32,33,34,35,36,37,38,39,40,. 41,42,43;ndata, text, alldata = xlsread(filename);A = tabulate(ndata);baifenbi=num2str(A(:,3),%1.2f);%计算百分比baifenbi=repmat(blanks(4),length(A(:,1),1),baifenbi,.repmat(%,length(A(:,1),1);baifenbi=cellstr(baifenbi);Label=strcat(note,baifenbi);pie(A(:,3),Label)2、运行结果：图十 班上同学的来源分布图3、 结论：班上同学的来源分布如图是所示。附录一：Grade.xls：414243444514243444546243444546473444546474844546474849546474849502474899505134849905152449505152535505192539425152935455352539495564535455565755455965758255565758993565798