MATLAB語(yǔ)言編程介紹

1、MATLAB語(yǔ)言編程介紹介紹nMATLAB中各種命令可以完成許多單一的任務(wù)，對(duì)于某些命令可以完成許多單一的任務(wù)，對(duì)于某些較為復(fù)雜的問(wèn)題，僅靠現(xiàn)有的命令或函數(shù)來(lái)解決，往往是較為復(fù)雜的問(wèn)題，僅靠現(xiàn)有的命令或函數(shù)來(lái)解決，往往是難以達(dá)到目的難以達(dá)到目的 。為此，要運(yùn)用。為此，要運(yùn)用MATLAB編程語(yǔ)言編制程編程語(yǔ)言編制程序，形成序，形成M-文件。文件。n程序是使計(jì)算機(jī)完成各項(xiàng)運(yùn)算的命令集，運(yùn)行一個(gè)編制好程序是使計(jì)算機(jī)完成各項(xiàng)運(yùn)算的命令集，運(yùn)行一個(gè)編制好的程序，計(jì)算機(jī)會(huì)從第一條命令行開始，一行接一行地執(zhí)的程序，計(jì)算機(jī)會(huì)從第一條命令行開始，一行接一行地執(zhí)行相應(yīng)的命令，直到終止。行相應(yīng)的命令，直到終止。n程

2、序一般包含：數(shù)據(jù)輸入，數(shù)據(jù)處理和結(jié)果輸出三部分，程序一般包含：數(shù)據(jù)輸入，數(shù)據(jù)處理和結(jié)果輸出三部分，其中數(shù)據(jù)處理是核心。程序編寫調(diào)試完成后，需要存盤，其中數(shù)據(jù)處理是核心。程序編寫調(diào)試完成后，需要存盤，形成永久性文件，可以隨時(shí)對(duì)它進(jìn)行調(diào)用或修改。文件形成永久性文件，可以隨時(shí)對(duì)它進(jìn)行調(diào)用或修改。文件 名以字母開頭，但不能用專用變量名，如名以字母開頭，但不能用專用變量名，如pi等。等。 nMATLAB中每一個(gè)命令都是一個(gè)命令都是一個(gè)M-文件。文件。M-文件中輸入、輸出命令命令 n直接賦值輸入： 提示對(duì)話輸入(input命令) x=2; x=input(請(qǐng)輸入?yún)?shù) x=); a=1,2;3,4; a=i

3、nput(請(qǐng)輸入矩陣 a=); s= any string!; s=input(Please input s=); n直接輸出： 格式控制輸出(fprintf命令) x fprintf(x=%.0f, y=%.5fn,pi,pi); disp(a,b); fprintf(x=%5g, y=%10.5gn,N,pi);M-文件中循環(huán)控制命令命令(for命令命令)n格式： for i=n1:(step):n2 commands; end n作用：重復(fù)執(zhí)行命令集作用：重復(fù)執(zhí)行命令集commands.n例例2：求奇數(shù)和：求奇數(shù)和：s=1+3+5+(2k-1) n=input(請(qǐng)輸入 n=); s=0;

4、 for i=1:2:n s=s+i; fprintf(i=%.0f, s=%.0fn,i,s); endn例例1：求和：求和：s=1+2+n n=input(請(qǐng)輸入 n=); s=0; for i=1:n s=s+i; fprintf(i=%.0f, s=%.0fn,i,s); endn例例3：求階乘：求階乘：p=12 3 n=n! n=input(請(qǐng)輸入 n= ); p=1; for i=1:n p=p*i; fprintf( i=%.0f, p=%.0fn ,i,p); endn例例4：求：求e：e=1+1+1/2!+1/3!+1/n! n=input(請(qǐng)輸入 n= ); p=1;e=1

5、; for i=1:n p=p*i; p1=1/p; e=e+p1; fprintf( i=%.0f, p=%.0f, e=%.8f n ,i,p,e); endn例例5：求：求pi：pi/4=1-1/3+1/5-1/7+1/9-1/11+ n=input(請(qǐng)輸入 n= ); pi=0;k=-1; for i=1:2:n k=(-1)*k; pi=pi+k*1/i; fprintf( i=%.0f, pi=%.8f n ,i,4*pi); endM-文件中條件循環(huán)命令命令(while命令命令)n格式： while (condition is true) commands; end n作用：重復(fù)

6、執(zhí)行命令集作用：重復(fù)執(zhí)行命令集commands.n例例6：求和：求和：s=1+2+n3000 s=0;k=1; while s3000 s=s+k; fprintf( k=%.0f, s=%.0fn ,k,s); k=k+1; endn例例7：求：求e：e=1+1+1/2!+1/3!+ (1/n!=1.0e-8 p=p*i; r=1/p; e=e+r; fprintf( i=%.0f, r=%.9f, e=%.9f n ,i,r,e); i=i+1; endM-文件中選擇控制命令命令(if命令命令)n單項(xiàng)選擇控制單項(xiàng)選擇控制n格式： if (condition is true) command

7、s; end n作用：若條件成立，則執(zhí)行命令集作用：若條件成立，則執(zhí)行命令集 commands. 否則，不執(zhí)行。否則，不執(zhí)行。n例例8：求：求n個(gè)實(shí)數(shù)中最大的數(shù)個(gè)實(shí)數(shù)中最大的數(shù)M. a=input(請(qǐng)輸入數(shù)組 an= ); m,n=size(a); M=a(1); for i=2:n if M0 sn=1;elseif x=0 sn=0;else sn=-1;endfprintf(x=%.5f, sn=%.0fn,x,sn);數(shù)值計(jì)算問(wèn)題數(shù)值計(jì)算問(wèn)題n問(wèn)題問(wèn)題1：求無(wú)理數(shù)的近似值求無(wú)理數(shù)的近似值 先求先求 的近似值，再設(shè)計(jì)通用程序的近似值，再設(shè)計(jì)通用程序.n問(wèn)題問(wèn)題2：用兩種方法求用兩種方法求

8、Pi的近似值的近似值)2, 1( nAAn2n例10：用二分法求函數(shù)x2-2=0的正實(shí)根.0)()(,2 , 1 , , 2)(2bfafbaxxf :2/ )() 1bac; ),| )(|0)( ifcgrcfcf);()(, 0)()( elseifcfbfcbafcf);()(, elsecfafca1). gotoor ;stop, ,| )(| if)2cgrcf輸入初值：Kerbaxf , , ),(2/ )(baccr ?0)(rf?0)()(rfaf2/ )(bac?|ercrrb ra 輸出結(jié)果：STOPNoNoNoYesYesYesn例11：用切線法求函數(shù)x2-2=0的正

9、實(shí)根.0)( ,2 , 1 , , 2)(2 xfbaxxf,.2 , 1 , 0,)()( 2) 2,) 110kxfxfxxxkkkk2). gotoor ;stop, ,| )(| if) 3kkxgrxfn例12：用泰勒公式求pi的近似值。)() 12() 1(.753 arctan)(12121753kkkxokxxxxxxxf) 12(1) 1(.715131141kk特點(diǎn)：收斂速度慢，循環(huán)次數(shù)多特點(diǎn)：收斂速度慢，循環(huán)次數(shù)多.n例13：用梯形求積公式求pi的近似值。10102|arctan414) 1xdxxnkkhaxnabhxfxfxfhTdxxfknnkknba, 1 , 0 , , )()(2)(2)(1102)梯形公式梯形公式nkkhaxnabhxxfhTxfxxfxfxfhTdxxfknkkknnnkkknkknba, 1 , 0 , , ).)2(21 )()2(2)(2)(4)(111111023)復(fù)化梯形公式復(fù)化梯形公式n例14：用Simpson求積公式求pi的近似值。1) 辛普生公式辛普生公式nkkhaxnabhxxfxfxfxfhSdxxfknkkknnkknba,