1、會計學(xué)1matlabm文件文件(wnjin)和函數(shù)句柄和函數(shù)句柄第一頁，共38頁。第1頁/共37頁第二頁，共38頁。for n=1:10 x(n)=sin(n*pi/10);endxx =Columns 1 through 7 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090Columns 8 through 10 0.5878 0.3090 0.0000第2頁/共37頁第三頁，共38頁。for循環(huán)的重要(zhngyo)說明第3頁/共37頁第四頁，共38頁。for n=1:10 x(n)=sin(n*pi/10); n=10;endxx = Co

2、lumns 1 through 7 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns 8 through 10 0.5878 0.3090 0.0000 第4頁/共37頁第五頁，共38頁。 data=3 9 45 6; 7 16 -1 5data = 3 9 45 6 7 16 -1 5for n=data x=n(1)-n(2)end;x = -4x = -7x = 46x = 1 第5頁/共37頁第六頁，共38頁。for n=1:5 for m=5:-1:1 A(n,m)=n2+m2;end disp(n)end12345AA

3、 =2 5 10 17 265 8 13 20 2910 13 18 25 3417 20 25 32 4126 29 34 41 50 第6頁/共37頁第七頁，共38頁。 n=1:10;x=sin(n*pi/10)x =Columns 1 through 7 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 0.8090Columns 8 through 10 0.5878 0.3090 0.0000 第7頁/共37頁第八頁，共38頁。x=zeros(1,10); % preallocated memory for xfor n=1:10 x(n)=sin(

4、n*pi/10);end; 第8頁/共37頁第九頁，共38頁。 num=0;EPS=1; while (1+EPS)1EPS=EPS/2;num=num+1;end numnum =53 EPS=2*EPSEPS =2.2204e-016第9頁/共37頁第十頁，共38頁。if expression commandsend if expression commands evaluated if True else commands evaluated if Falseendif expression1 commands evaluated if expression1 is True elsei

5、f expression2 commands evaluated if expression2 is True elseif else commands evaluated if no other expression is True end 第10頁/共37頁第十一頁，共38頁。 EPS=1; for num=1:1000 EPS=EPS/2; if (1+EPS) numnum = 53 EPS=1; for num=1:1000 EPS=EPS/2; if (1+EPS)1continue end EPS=EPS*2 breakendEPS =2.2204e-016 numnum = 5

6、3 第11頁/共37頁第十二頁，共38頁。第12頁/共37頁第十三頁，共38頁。for i=1:10;ai=89+i;bi=79+i;ci=69+i;di=59+i;end;c=d,c;Name=Jack,Marry,Peter,Rose,Tom;Mark=72,83,56,94,100;Rank=cell(1,5);S=struct(Name,Name,Marks,Mark,Rank,Rank);for i=1:5 switch S(i).Marks case 100 S(i).Rank=滿分(mn fn); case a S(i).Rank= 優(yōu)秀; case b S(i).Rank= 良

7、好; case c S(i).Rank= 及格; otherwise S(i).Rank=不及格; endend 第13頁/共37頁第十四頁，共38頁。disp(學(xué)生姓名 , 得分 , 等級);disp( )for i=1:5;disp(S(i).Name,blanks(6),num2str(S(i).Marks),blanks(6),S(i).Rank);end;學(xué)生姓名 得分 等級 Jack 72 及格Marry 83 良好(lingho)Peter 56 不及格Rose 94 優(yōu)秀Tom 100 滿分 第14頁/共37頁第十五頁，共38頁。第15頁/共37頁第十六頁，共38頁。N=4;A

8、=magic(3);try A_N=A(N,:),catch A_end=A(end,:),endA_end = 4 9 2lasterrans =Index exceeds matrix dimensions. 第16頁/共37頁第十七頁，共38頁。第17頁/共37頁第十八頁，共38頁。第18頁/共37頁第十九頁，共38頁。5.函數(shù)函數(shù)(hnsh)體：這部分內(nèi)容由實現(xiàn)該體：這部分內(nèi)容由實現(xiàn)該 M函數(shù)函數(shù)(hnsh)文件文件功能的功能的Matlab指令指令組成。它接收輸入?yún)⒘浚M行程序流控制，得到輸出參量。組成。它接收輸入?yún)⒘浚M行程序流控制，得到輸出參量。第19頁/共37頁第二十頁，共38頁

9、。circle.mfunction sa = circle(r,s)%CIRCLEplot a circle of radii r in the line specified by s.% r指定(zhdng)半徑的數(shù)值% s指定(zhdng)線色的字符串% sa圓面積% circle(r)利用藍實線畫半徑為 r 的圓周線.% circle(r,s)利用串 s 指定(zhdng)的線色畫半徑為 r 的圓周線.% sa=circle(r)計算圓面積，并畫半徑為 r 的藍色圓面.% sa=circle(r,s)計算圓面積，并畫半徑為 r 的 s 色圓面.% 編寫于1999年4月7日，修改于1999

10、年8月27日。if nargin2 error(輸入?yún)⒘刻唷?;end;if nargin=1 s=b;end;t=0:pi/100:2*pi;x=r*exp(i*t);if nargout=0 plot(x,s);else sa=pi*r*r; fill(real(x),imag(x),s)end第20頁/共37頁第二十一頁，共38頁。第21頁/共37頁第二十二頁，共38頁。第22頁/共37頁第二十三頁，共38頁。指令含義nargin在函數(shù)體內(nèi)，用于獲取實際輸入?yún)⒘縩argout在函數(shù)體內(nèi)，用戶獲取實際輸出參量nargin(fun)獲取fun指定函數(shù)的標(biāo)稱輸入?yún)⒘繑?shù)nargout(fun)

11、獲取fun指定函數(shù)的標(biāo)稱輸出參量數(shù)inputname(n)在函數(shù)體內(nèi)使用，給出第 n個輸入?yún)⒘康膶嶋H調(diào)用變量名第23頁/共37頁第二十四頁，共38頁。n先叢所在函數(shù)空間獲取變量值，用eval(exp1)計算原串表達式；若該計算失敗，則再從workspace指定的工作空間中獲取變量，再通過eval(exp2)計算替代串表達式。第24頁/共37頁第二十五頁，共38頁。evallinzzy.mfunction y1=evallinzzy(a,s)t=(0:a)/a*2*pi;y1=subevallinzzy(4,s);function y2=subevallinzzy(a,s)t=(0:a)/a*2

12、*pi;ss=a*exp(i*t);switch s case base,caller y2=evalin(s,ss); case self y2=eval(ss);end a=30;t=(0:a)/a*2*pi;sss=base,caller,self;for k=1:3y0=evallinzzy(8,sssk);subplot(1,3,k)plot(real(y0),imag(y0)end 第25頁/共37頁第二十六頁，共38頁。第26頁/共37頁第二十七頁，共38頁。assigninxyq.mfunction y=assigninxyq(x)y=sqrt(x);t=x2;assignin

13、(base,yy,t) clear;x=4;y=assigninxyq(x);disp(blanks(5),x,blanks(5),y,blanks(4),yy),disp(x,y,yy) x y yy4 2 16 第27頁/共37頁第二十八頁，共38頁。第28頁/共37頁第二十九頁，共38頁。第29頁/共37頁第三十頁，共38頁。第30頁/共37頁第三十一頁，共38頁。y=eval(CEM) 執(zhí)行CEM指定的計算(j sun)y1,y2,=eval(CEM) 執(zhí)行CEM代表的函數(shù)文件，并輸出結(jié)果 第31頁/共37頁第三十二頁，共38頁。t=pi;eval(theta=t/2,y=sin(th

14、eta)theta = 1.5708y = 1第32頁/共37頁第三十三頁，共38頁。y1,y2,=feval(FH,arg1,arg2,)執(zhí)行函數(shù)句柄FH指定(zhdng)的計算y1,y2,=feval(FN,arg1,arg2,)執(zhí)行函數(shù)名字符串FN指定(zhdng)的計算y1,y2,=feval(FIL,arg1,arg2,)執(zhí)行內(nèi)聯(lián)函數(shù)FIL指定(zhdng)的計算 x=pi/4;Ve=eval(1+sin(x) Ve = 1.7071 Vf=feval(1+sin(x),x) ? Error using = fevalInvalid function name 1+sin(x). r

15、and(seed,1);A=rand(2,2);Heig=eig;d=feval(Heig,A) d = 0.7568 -0.1488 de=eval(Heig,A) ? Error using = evalFunction eval not defined for variables of class function_handle. 第33頁/共37頁第三十四頁，共38頁。第34頁/共37頁第三十五頁，共38頁。第35頁/共37頁第三十六頁，共38頁。argout1,argout2,argoutn=FunName(argin1,argin2,arginn)hfun=FunNameargout1,argout2,argoutn=feval(hfun,argin1,argin2,arginn) 第36頁/共37頁