中國(guó)地質(zhì)大學(xué)（武漢） 測(cè)控技術(shù)與儀器系測(cè)控軟件技術(shù)基礎(chǔ)作業(yè)報(bào)告一、題目一：1、代碼：%內(nèi)容 ： 生成一個(gè)1到30并且以3遞進(jìn)的數(shù)列，計(jì)算其長(zhǎng)度、均值、最值和方差%作者 ： 柏明明 測(cè)控一班 20121002120%備注 ： 測(cè)控軟件與技術(shù)基礎(chǔ)第一次作業(yè)Num = 1:3:30;Length = length(Num) %長(zhǎng)度Average = mean(Num) %均值 Max = max(Num) %最大值Min = min(Num) %最小值Var = sum(Num - mean(Num).2)/Length %方差Var1 = var(Num) %經(jīng)檢驗(yàn)用該公式計(jì)算結(jié)果與方差定義式結(jié)果不同var2 = (std(Num)2 %經(jīng)檢驗(yàn)用該公式計(jì)算結(jié)果與方差定義式結(jié)果不同2、 運(yùn)行結(jié)果：Length = 10；Average = 14.5000；Max = 28；Min = 1；Var = 74.2500；Var1 = 82.5000；var2 = 82.5000。3、 結(jié)論：該向量的長(zhǎng)度、均值、最值和方差結(jié)果如上所示，但有一點(diǎn)值得注意，在求方差時(shí)通過(guò)三種不同求解方法分別得到了三個(gè)結(jié)果，Var、Var1、var2。其中Var是嚴(yán)格用方差定義式求得的，Var1是用matlab自帶的求方差函數(shù)var()求得的，Var2是用matlab自帶的求標(biāo)準(zhǔn)差函數(shù)std()平方后求得的，Var1值等于Var2值而與Var值不等，說(shuō)明matlab在求方差時(shí)不是按照方差定義式求解的。二、題目二：1、代碼：%內(nèi)容 ： 計(jì)算y = 3*x + 1 從1到4的積分syms a;x = 1:0.001:4;y = 3*x + 1;JiFen1 = sum(0.001*y) %用細(xì)分曲邊梯形法求積分JiFen2 = int(3*a+1,1,4) %公式求積分2、 運(yùn)行結(jié)果：JiFen1 = 25.5085；JiFen2 = 51/2。3、 結(jié)論：y = 3*x + 1 從1到4的積分結(jié)果如上所示，其中JiFen1是通過(guò)近似求曲邊梯形的面積得到的，JiFen2是通過(guò)matlab自帶的求定積分函數(shù)int()求得，兩者結(jié)果近似相等，說(shuō)明matlab在求定積分時(shí)也是按照求曲邊梯形面積近似得到的，只不過(guò)它分割的更為細(xì)致，因而結(jié)果更為精確。3、 題目三：1、代碼：%內(nèi)容 ：編寫程序，生成一份班級(jí)的學(xué)科成績(jī)匯總表，并保存在Excel表格中。要求學(xué)生人數(shù)不少于40，每人至少6們課，成績(jī)計(jì)分式包括百分制和五分制兩種。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(文件名，數(shù)組名,起始單元格：終止單元格)函數(shù) 2、 運(yùn)行結(jié)果：見 附錄。3、 結(jié)論：xlswrite(文件名，數(shù)組名,起始單元格：終止單元格)函數(shù)在使用時(shí)各參數(shù)分別如括號(hào)中所示，當(dāng)沒(méi)有該文件時(shí)會(huì)自動(dòng)創(chuàng)建該文件并寫入。4、 題目四：1、代碼：%內(nèi)容 ： 對(duì)上題存儲(chǔ)在excel中的班級(jí)數(shù)據(jù)，統(tǒng)計(jì)成績(jī)的最值和均值，以及不及格人數(shù)及大于90的人數(shù)，繪制直方圖展示成績(jī)分布特性 grade = xlsread(Grade.xls);Lesson1 = grade(:,1); %Lesson1成績(jī)Lesson2 = grade(:,2); %Lesson2成績(jī)Lesson3 = grade(:,3); %Lesson3成績(jī)Lesson4 = grade(:,4); %Lesson4成績(jī)Lesson5 = grade(:,5); %Lesson5成績(jī)Lesson6 = grade(:,6); %Lesson6成績(jī) 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(Lesson1<60) %每門課不及格人數(shù)Ls2Not = length(find(Lesson2<60)Ls3Not = length(find(Lesson3<60)Ls4Not = length(find(Lesson4<60)Ls5Not = length(find(Lesson5<60)Ls6Not = length(find(Lesson6<60) Ls1OK = length(find(Lesson1>90) %每門課優(yōu)秀人數(shù)Ls2OK = length(find(Lesson2>90)Ls3OK = length(find(Lesson3>90)Ls4OK = length(find(Lesson4>90)Ls5OK = length(find(Lesson5>90)Ls6OK = length(find(Lesson6>90)2、運(yùn)行結(jié)果：（1） 數(shù)值結(jié)果：a.各門成績(jī)最大值：Max1 = 99；Max2 = 99；Max3 = 99；Max4 = 95；Max5 = 99；Max6 = 5。b.各門成績(jī)最小值：Min1 = 41；Min2 = 42；Min3 = 43；Min4 = 44；Min5 = 45；Min6 = 1。c.各門成績(jī)均值：Mean1 = 64；Mean2 = 67.5000；Mean3 =73；Mean4 =66.2500；Mean5 = 67.2500；Mean6 = 3.2750。d.各門成績(jī)未通過(guò)人數(shù)：Ls1Not = 19；Ls2Not = 18；Ls3Not =10；Ls4Not =15；Ls5Not =13；Ls6Not =40。e.各門成績(jī)優(yōu)秀人數(shù)：Ls1OK =3；Ls2OK =10；Ls3OK =10Ls4OK =3；Ls5OK = 4；Ls6OK =0。（2）圖像結(jié)果：圖一 Lesson1直方圖展示圖二 Lesson2直方圖展示圖三 Lesson3直方圖展示圖四 Lesson4直方圖展示圖五 Lesson5直方圖展示圖六 Lesson6直方圖展示3、 結(jié)論：對(duì)上題存儲(chǔ)在excel中的班級(jí)數(shù)據(jù)，統(tǒng)計(jì)各門課成績(jī)的最值和均值，以及不及格人數(shù)及大于90的人數(shù)如運(yùn)行結(jié)果中數(shù)值結(jié)果所示，各門課成績(jī)直方圖展示如圖像結(jié)果所示。5、 題目五：1、 代碼：%內(nèi)容 ： 繪制 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ùn)行結(jié)果：圖七 y = x4+3*x+1 繪制結(jié)果3、 結(jié)論：y = x4 + 3*x + 1繪制結(jié)果如圖七所示。要掌握繪制二維函數(shù)圖像中常用的函數(shù)的用法，比如figure()、plot()、xlabel()、ylabel()、legend()、title()等等。6、 題目六：1、代碼：%內(nèi)容 ： fft測(cè)試：幅度為0.7頻率為50HZ的信號(hào)與幅度為1頻率為120HZ的信號(hào)疊加后的頻率觀察 Fs = 1000; %信號(hào)采樣頻率T = 1/Fs; %信號(hào)采樣時(shí)間L = 1000; %信號(hào)長(zhǎng)度t = (0:L-1)*T; %時(shí)間軸x = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t); %50HZ和120HZ信號(hào)的疊加y = x + 2*randn(size(t); %符合正太分布的噪聲的引入figure(1); plot(Fs*t(1:50),y(1:50);title(引入正態(tài)分布噪聲的信號(hào)的時(shí)域圖像);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(引入正態(tài)分布噪聲的信號(hào)的頻域圖像);xlabel(f（HZ）);ylabel(|Y(f)|);2、 運(yùn)行結(jié)果：圖八 引入正態(tài)分布噪聲的信號(hào)的時(shí)域圖像圖九 引入正態(tài)分布噪聲的信號(hào)的頻域圖像3、 結(jié)論：a.通過(guò)對(duì)圖九的觀察可以發(fā)現(xiàn)在50HZ和120HZ左右明顯有兩個(gè)峰值，且其幅度跟設(shè)置幅度相近，說(shuō)明該傅里葉變換可以反映所設(shè)置的時(shí)域信號(hào)的頻域特性。b.在仿真時(shí)域信號(hào)時(shí)，常引入隨機(jī)數(shù)rand()或randn()做為噪聲。7、 題目七：1、代碼：%內(nèi)容： 統(tǒng)計(jì)一下班上同學(xué)的來(lái)源管理，用餅狀圖表示出來(lái)。 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);%計(jì)算百分比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、運(yùn)行結(jié)果：圖十 班上同學(xué)的來(lái)源分布圖3、 結(jié)論：班上同學(xué)的來(lái)源分布如圖是所示。附錄一：Grade.xls：41424344451424344454624344454647344454647484454647484954647484950247489950513484990515244950515253550519253942515293545535253949556453545556575545596575825556575899356579859604575859606155859909162259606162633606162636446162636465562636465662636495666736465666768465669768695666768697026798699071368997071924699071727358091727374481927394755829374757658394757677284957677782859677787918697987980397989980814987980819239980818283290818283841附錄二：班上同學(xué)的來(lái)源分布表：第 20 頁(yè)