數(shù)學實驗13-matlab編程

1、 數(shù)學實驗理學院數(shù)學學科 李換琴 信息01、02、03班上機時間安排10月15日周五晚上18：3022：0010月29日周五上午10：0012：0011月21日周日上午08：3012：0012月01日周三上午08：0010：00信息04、05、06班上機時間安排10月16日周六上午08：3012：0010月29日周五晚上18：0020：0011月21日周日下午14：3018：0012月04日周六上午10：0012：00計算機01、02、03班上機時間安排10月10日周日上午08：3012：0010月29日周五上午8：009：5011月13日周六上午08：3012：0012月01日周三上午10：

2、0012：00計算機04、05、06班上機時間安排10月10日周日下午14：3018：0010月29日周五晚上20：0022：0011月13日周六晚上18：3022：0012月04日周六上午08：0010：00MATLAB語言編程介紹MATLAB中各種命令可以完成許多單一的任務，對于某些較為復雜的問題，僅靠現(xiàn)有的命令或函數(shù)來解決，往往是難以達到目的 。為此，要運用MATLAB編程語言編制程序，形成M-文件。程序是使計算機完成各項運算的命令集，運行一個編制好的程序，計算機會從第一條命令行開始，一行接一行地執(zhí)行相應的命令，直到終止。程序編寫調(diào)試完成后，需要存盤，形成永久性文件，可以隨時對它進行調(diào)用

3、或修改。文件名以字母開頭，但不能用專用變量名，如pi,ans,eps等。 程序結構程序一般包含：a=input(請輸入數(shù)組 an= ); m,n=size(a); M=a(1); for i=2:nif M n2.省略格式： for i=n1:n2這里step=1.例1：求n個奇數(shù)和：s=1+3+5+(2n-1) n=input(please input n=); s=0; for i=1:n s=s+(2*i-1); fprintf(i=%.0f, s=%.0fn,i,s); end 程序aa1例2：求階乘：p=12 3 n=n! n=input(請輸入 n= ); p=1; for i=1

4、:n p=p*i; fprintf( i=%.0f, p=%.0fn ,i,p); end aa2.m例3：求e：e=1+1+1/2!+1/3!+1/n! n=input(請輸入 n= ); p=1;e=1; for i=1:n p=p*i; p1=1/p; e=e+p1; fprintf( i=%.0f, p=%.0f, e=%.8f n ,i,p,e); end aa3.mM-文件中條件循環(huán)命令(while命令)格式： while (condition is true) commands; end 作用：當條件成立時，執(zhí)行命令集commands,直到條件不成立.例4：求和：s=1+2+n3

5、000 s=0;k=1; while s+k3000 s=s+k; fprintf( k=%.0f, s=%.0fn ,k,s); k=k+1; end aa4.m例5：求e：e=1+1+1/2!+1/3!+ (1/n!=1.0e-8p=p*i; r=1/p; e=e+r; fprintf( i=%.0f, r=%.9f, e=%.9f n ,i,r,e); i=i+1; end aa5.mM-文件中選擇控制命令(if命令)單項選擇控制格式： if (condition is true) commands; end 作用：若條件成立，則執(zhí)行命令集 commands. 否則，不執(zhí)行。例6：求n個

6、實數(shù)中最大的數(shù)M. a=input(請輸入數(shù)組 an= ); m,n=size(a); M=a(1); for i=2:n if M0 z=1; elseif x=0 z=0; else z=-1; end y=x2*z aa7.m用M-文件定義函數(shù)格式 function y1,y2=funname(x,y,a,n)其中： function 為關鍵字； y1,y2,為輸出列表； funname為自定義的函數(shù)名； (x,y,a,n,)為函數(shù)參數(shù)，用來傳遞相 關數(shù)值. 存盤時，要求用funname作為函數(shù)名. 例7：建立符號函數(shù)sgn(x) function sn=sgn(x) if x0 sn=

7、1; elseif x=0 sn=0; else sn=-1; end以sgn作為文件名存盤，即建立了函數(shù)。調(diào)用： 在命令區(qū)執(zhí)行 : sn=sgn(10)或sn=sgn(-2)例8：建立計算(a+b)n,(a-b)n的函數(shù) function p,q=fun1(a,b,n) p=(a+b).n; q=(a-b).n; fun1.m 例編制程序，從鍵盤輸入a,b,n,計算(a+b)n,(a-b)n aa8.ma=input(a=);b=input(b=);n=input(n=)；p,q=fun1(a,b,n);fprintf(a+b)n=%.4f,(a-b)n=%.4fn,p,q)例10：用二分法

8、求函數(shù)x2-2=0的正實根.輸入初值：輸出結果：迭代次數(shù)k,方程的近似根c.NoNoNoYesYesYesa=1;b=2;er=1.0e-8;f=inline(x.2-2);k=1;c=(a+b)/2;r=c+1;while abs(r-c)err=c;if abs(f(r)1.0e-10fprintf(k=%.0f,r=%.8fn,k,r);break;endif f(a)*f(r)0 b=r;else a=r;endc=(a+b)/2;fprintf(k=%.0f,r=%.8fn,k,r);k=k+1;enderfenfa.m本次實驗三人一個小組，合作完成任務，交一份實驗報告。任務：1，2

9、任選1題，3，4任選1題， 練習1第1題任選1小題。實驗報告內(nèi)容：1）任務描述； 2）問題分析及程序（包括符號說明）； 3）運行結果 ；4）結果分析與說明。實驗報告提交時間：待下次任務完成以后，將兩次內(nèi)容一并打印，封皮上寫清楚班級，學號，姓名。課代表收齊于17周，18周周五晚上7點半至9點半交到我辦公室（理科樓225）。1. 每門課程考試閱卷完畢,任課教師都要對各班的考試成績進行統(tǒng)計,統(tǒng)計內(nèi)容包括:全班人數(shù),總得分,平均得分,不及格的人數(shù)及90分(包括90分)以上的人數(shù).請編制程序解決這一問題,并自給一組數(shù)據(jù)驗證程序的正確性.要求:使用者在提示下通過鍵盤輸入學生成績,計算機自動處理后,顯示需要

10、的結果.2. 根據(jù)我國個人所得稅計算方法,編制程序,要求:使用者在系統(tǒng)提示下通過鍵盤輸入月工資薪金收入總數(shù),計算機則在屏幕上顯示個人所得稅額,界面友好,方便使用.個人所得稅計算方法: 月個人所得稅（月工資薪金收入1600）*適用稅率速算扣除數(shù)附表：個人所得稅稅率表（工資、薪金所得適用）個人所得稅稅率表（工資、薪金所得適用）級 數(shù) 全月應納稅所得額 稅率（）稅率 速算扣除數(shù)1 不超過500元的 52 超過500元至2000元的部分 10 253 超過2000元至5000元的部分 15 1254 超過5000元至20000元的部分 20 3755 超過20000元至40000元的部分 25 137

11、56 超過40000元至60000元的部分 30 33757 超過60000元至80000元的部分 35 63758 超過80000元至100000元的部分 40 103759 超過100000元的部分 45 15375首先由計算機隨機產(chǎn)生一個 1,100 之間的一個整數(shù)，然后由用戶猜測所產(chǎn)生的這個數(shù)。根據(jù)用戶猜測的情況給出不同的提示，如果猜測的數(shù)大于產(chǎn)生的數(shù)，則顯示 “High” ，小于則顯示 “ Low ” ，等于則顯示 “You won！”，同時退出游戲。用戶最多有 7 次機會。3. 編寫程序?qū)崿F(xiàn)猜數(shù)游戲4. 編寫程序完成復印資料給出的問題。a=1.;b=2.;r=1.0e-7;ya=fc1(a);yb=fc1(b);k=1;if abs(ya)r root=a; fprintf(k=%.0f, Root x=%.4f, f(x)=%12.