實(shí)驗(yàn)二MATLAB程序設(shè)計(jì)(含實(shí)驗(yàn)報(bào)告)

1、實(shí)驗(yàn)二 MATLAB程序設(shè)計(jì)一、 實(shí)驗(yàn)?zāi)康?.掌握利用if語(yǔ)句實(shí)現(xiàn)選擇結(jié)構(gòu)的方法。2.掌握利用switch語(yǔ)句實(shí)現(xiàn)多分支選擇結(jié)構(gòu)的方法。3.掌握利用for語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。4.掌握利用while語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。5.掌握MATLAB函數(shù)的編寫(xiě)及調(diào)試方法。二、 實(shí)驗(yàn)的設(shè)備及條件計(jì)算機(jī)一臺(tái)（帶有MATLAB7.0以上的軟件環(huán)境）。 M文件的編寫(xiě)：?jiǎn)?dòng)MATLAB后，點(diǎn)擊File|New|M-File，啟動(dòng)MATLAB的程序編輯及調(diào)試器（Editor/Debugger），編輯以下程序，點(diǎn)擊File|Save保存程序，注意文件名最好用英文字符。點(diǎn)擊Debug|Run運(yùn)行程序，在命令窗口查看

2、運(yùn)行結(jié)果，程序如有錯(cuò)誤則改正三、 實(shí)驗(yàn)內(nèi)容1.編寫(xiě)求解方程的根的函數(shù)（這個(gè)方程不一定為一元二次方程，因的不同取值而定），這里應(yīng)根據(jù)的不同取值分別處理，有輸入?yún)?shù)提示，當(dāng)時(shí)應(yīng)提示“為恒不等式!”。并輸入幾組典型值加以檢驗(yàn)。（提示：提示輸入使用input函數(shù)）2.輸入一個(gè)百分制成績(jī)，要求輸出成績(jī)等級(jí)A+、A、B、C、D、E。其中100分為A+，90分99分為A，80分89分為B，70分79分為C，60分69分為D，60分以下為E。要求：（1）用switch語(yǔ)句實(shí)現(xiàn)。 （2）輸入百分制成績(jī)后要判斷該成績(jī)的合理性，對(duì)不合理的成績(jī)應(yīng)輸出出錯(cuò)信息。（提示：注意單元矩陣的用法）3.數(shù)論中一個(gè)有趣的題目：任意

3、一個(gè)正整數(shù)，若為偶數(shù)，則用2除之，若為奇數(shù)，則與3相乘再加上1。重復(fù)此過(guò)程，最終得到的結(jié)果為1。如：2®13®10®5®16®8®4®2®16®3®10®5®16®8®4®2®1運(yùn)行下面的程序，按程序提示輸入n=1,2,3,5,7等數(shù)來(lái)驗(yàn)證這一結(jié)論。請(qǐng)為關(guān)鍵的Matlab語(yǔ)句填寫(xiě)上相關(guān)注釋，說(shuō)明其含義或功能。%classic "3n+1" problem from number theory.while 1 n=in

4、put('Enter n,negative quits:'); if n<=0 breakend a=n; while n>1 if rem(n,2)=0 n=n/2; else n=3*n+1; end a=a,n; end aend4. 編寫(xiě)一個(gè)函數(shù)，計(jì)算下面函數(shù)的值，給出標(biāo)量x的值，調(diào)用該函數(shù)后，返回y的值。function y=myfun1(x)選擇一些數(shù)據(jù)測(cè)試你編寫(xiě)的函數(shù)。5. 編寫(xiě)一個(gè)函數(shù)求向量x中元素的平均值、最大值、最小值、均方根值。function m_x,max_x,min_x,rms_x=myfun2(x)方均根值(Root Mean Squa

5、re)的計(jì)算公式為：用下面數(shù)據(jù)測(cè)試你寫(xiě)的函數(shù)：（1）x=sin(0:0.01:6*pi)（2）x=rand(1,200)，得到的x為200個(gè)（0，1）之間均勻分布的隨機(jī)數(shù)。6.根據(jù)，求的近似值。當(dāng)n分別取100、1000、10000時(shí)，結(jié)果是多少？思考題：有一分?jǐn)?shù)序列： 編寫(xiě)一段程序，求前16項(xiàng)的和。四、實(shí)驗(yàn)報(bào)告要求（包含預(yù)習(xí)報(bào)告要求和最終報(bào)告要求）最終報(bào)告要求預(yù)習(xí)報(bào)告要求 1.實(shí)驗(yàn)名稱 2.實(shí)驗(yàn)?zāi)康?3.實(shí)驗(yàn)設(shè)備及條件 4.實(shí)驗(yàn)內(nèi)容及要求 5.實(shí)驗(yàn)程序設(shè)計(jì) 指程序代碼。 6.實(shí)驗(yàn)結(jié)果及結(jié)果分析 實(shí)驗(yàn)結(jié)果要求必須客觀,有數(shù)據(jù)的可以記錄數(shù)據(jù)，沒(méi)有數(shù)據(jù)的簡(jiǎn)單描述實(shí)驗(yàn)現(xiàn)象。結(jié)果分析是對(duì)實(shí)驗(yàn)結(jié)果的理論

6、評(píng)判。 7.實(shí)驗(yàn)中出現(xiàn)的問(wèn)題及解決方法 8. 思考題的回答四、 實(shí)驗(yàn)報(bào)告的提交方式Word文檔，命名方式:實(shí)驗(yàn)號(hào)_你的學(xué)號(hào)_姓名!例如本次實(shí)驗(yàn)：實(shí)驗(yàn)一_000000001_張三.doc（信息101提交報(bào)告郵箱）：E_mail: matlab_xx01(網(wǎng)絡(luò)工程101提交作業(yè)郵箱)：E_mail: Matlab_wg01（注意網(wǎng)絡(luò)班的M是大寫(xiě)的)下一次課前提交,過(guò)期不收!五、 參考文獻(xiàn)參考教材和Matlab幫助文件。1. 實(shí)驗(yàn)名稱 MATLAB程序設(shè)計(jì)2. 實(shí)驗(yàn)?zāi)康?1.掌握利用if語(yǔ)句實(shí)現(xiàn)選擇結(jié)構(gòu)的方法。2.掌握利用switch語(yǔ)句實(shí)現(xiàn)多分支選擇結(jié)構(gòu)的方法。3.掌握利用for語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的

7、方法。4.掌握利用while語(yǔ)句實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的方法。5.掌握MATLAB函數(shù)的編寫(xiě)及調(diào)試方法。3. 實(shí)驗(yàn)設(shè)備及條件 計(jì)算機(jī)一臺(tái)（帶有MATLAB7.0以上的軟件環(huán)境）。 4. 實(shí)驗(yàn)內(nèi)容及要求 把實(shí)驗(yàn)內(nèi)容的應(yīng)用題，用MATLAB的語(yǔ)法編寫(xiě)出來(lái)，并運(yùn)行成功，注意題與題之間用相應(yīng)注釋分割。 5. 實(shí)驗(yàn)程序設(shè)計(jì) %1disp('一元二次方程計(jì)算器，請(qǐng)輸入下列數(shù)值進(jìn)行計(jì)算：');a=input('a=');b=input('b=');c=input('c=');if a=0 & b=0 & c=0 disp('此為恒

8、不等式');else d=b*b-4*a*c; x=(-b+sqrt(d)/(2*a),(-b-sqrt(d)/(2*a); disp('x1=',num2str(x(1),',x2=',num2str(x(2);end%2scores=input('成績(jī)等級(jí)劃分請(qǐng)輸入一個(gè)百分制成績(jī)，查詢劃分等級(jí)：');switch scores case 100 rate='A+' case num2cell(90:99) rate='A' case num2cell(80:89) rate='B' ca

9、se num2cell(70:79) rate='C' case num2cell(60:69) rate='D' case num2cell(0:59) rate='E' otherwise disp('輸出出錯(cuò)');enddisp(rate)%3%classic "3n+1" problem from number theory.while 1 n=input('Enter n,negative quits:'); %輸入一個(gè)非負(fù)整數(shù) if n<=0 %如果輸入的數(shù)為負(fù)數(shù),則不做操作

10、breakend a=n; %將n的值賦值給a while n>1 %滿足條件n>1時(shí)循環(huán)下列語(yǔ)句 if rem(n,2)=0 %當(dāng)n能被2整除則n除以2 n=n/2; else %如果不能被整除，則n乘以3加1 n=3*n+1; end a=a,n; %輸出數(shù)組a,n，并對(duì)比n是否>1,大于則繼續(xù)上訴循環(huán) end a %輸出最后的結(jié)果aend%4 myfun1.m%-clc format compactformat long gx=input('請(qǐng)輸入x的值:');y=tran1(x);disp('y=')disp(y)%4 tran.m%-

11、function y=myfun1(x)%當(dāng) x<=0 y=sin(x);%當(dāng) 3>x>0 y=x;%當(dāng) x>3 y=-x+6;if x<=0 y=sin(x);elseif 3>x>0 y=x;else x>3 y=-x+6;end%5 tran.m%-function m_x,max_x,min_x,rms_x=myfun2(x)%分別是平均值，最大值，最小值，均方根值m_x=mean(mean(x)max_x=max(max(x) min_x=min(min(x)a=sqrt(mean(x.2);rms_x=sqrt(mean(a.2)%5

12、 myfun2.m%-x=input('平均值、最大值、最小值、均方根值的計(jì)算，請(qǐng)輸入x的值，:');m_x,max_x,min_x,rms_x=tran(x);disp('平均值為')m_xdisp('最大值為')max_xdisp('最小值為')min_xdisp('均方根值為')rms_x%6 disp('(pi2)/6 = (1/12)+(1/22)+(1/32)+(1/n2),求pi的近似值') approximation=0;n=input('請(qǐng)輸出n值：'); appr

13、oximation=sqrt(sum(1./(1:n).2)*6);approximation 6. 實(shí)驗(yàn)結(jié)果及結(jié)果分析 第一題一元二次方程(ax2+bx+c=0)計(jì)算器，請(qǐng)輸入下列數(shù)值進(jìn)行計(jì)算：a=2b=1c=3x1=-0.25+1.199i,x2=-0.25-1.199i-一元二次方程(ax2+bx+c=0)計(jì)算器，請(qǐng)輸入下列數(shù)值進(jìn)行計(jì)算：a=0b=0c=>> 2此為恒不等式第二題成績(jī)等級(jí)劃分請(qǐng)輸入一個(gè)百分制成績(jī)，查詢劃分等級(jí)：111輸出出錯(cuò)成績(jī)等級(jí)劃分請(qǐng)輸入一個(gè)百分制成績(jī)，查詢劃分等級(jí)：55E_第三題%classic "3n+1" problem from

14、 number theory.while 1 n=input('Enter n,negative quits:'); %輸入一個(gè)非負(fù)整數(shù) if n<=0 %如果輸入的數(shù)為負(fù)數(shù),則不做操作 breakend a=n; %將n的值賦值給a while n>1 %滿足條件n>1時(shí)循環(huán)下列語(yǔ)句 if rem(n,2)=0 %當(dāng)n能被2整除則n除以2 n=n/2; else %如果不能被整除，則n乘以3加1 n=3*n+1; end a=a,n; %輸出數(shù)組a,n，并對(duì)比n是否>1,大于則繼續(xù)上訴循環(huán) end a %輸出最后的結(jié)果aend第四題請(qǐng)輸入x的值:-1y

15、= -0.841470984807897-請(qǐng)輸入x的值:2y= 2-請(qǐng)輸入x的值:9ans = 1y=-3第五題請(qǐng)輸入x的值:sin(0:0.01:6*pi)平均值為m_x = -1.1256e-007最大值為max_x = 1.0000最小值為min_x = -1.0000均方根值為rms_x =0.7071請(qǐng)輸入x的值:rand(1,200)平均值為m_x = 0.4969最大值為max_x = 0.9943最小值為min_x = 0.0099均方根值為rms_x = 0.5707第六題(pi2)/6 = (1/12)+(1/22)+(1/32)+(1/n2),求pi的近似值請(qǐng)輸出n值：100approximation = 3.1321(pi2)/6 = (1/12)+(1/22)+(1/32)+(1/n2),求pi的近似值請(qǐng)輸出n值：1000approximation = 3.1