matlabm文件和函數(shù)句柄學(xué)習(xí)教案

1、會計學(xué)1matlabm文件和函數(shù)句柄文件和函數(shù)句柄第一頁，編輯于星期六：十九點 二十分。第1頁/共37頁第二頁，編輯于星期六：十九點 二十分。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頁第三頁，編輯于星期六：十九點 二十分。for循環(huán)的重要說明第3頁/共37頁第四頁，編輯于星期六：十九點 二十分。for n=1:10 x(n)=sin(

2、n*pi/10); n=10;endxx = Columns 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頁第五頁，編輯于星期六：十九點 二十分。 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頁第六頁，編輯于星期六：十九點 二十分。for n=1:5 fo

3、r m=5:-1:1 A(n,m)=n2+m2;end disp(n)end12345AA =2 5 10 17 265 8 13 20 2910 13 18 25 3417 20 25 32 4126 29 34 41 50 第6頁/共37頁第七頁，編輯于星期六：十九點 二十分。 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頁第八頁，編輯于星期六：十九點 二

4、十分。x=zeros(1,10); % preallocated memory for xfor n=1:10 x(n)=sin(n*pi/10);end; 第8頁/共37頁第九頁，編輯于星期六：十九點 二十分。 num=0;EPS=1; while (1+EPS)1EPS=EPS/2;num=num+1;end numnum =53 EPS=2*EPSEPS =2.2204e-016第9頁/共37頁第十頁，編輯于星期六：十九點 二十分。if expression commandsend if expression commands evaluated if True else command

5、s evaluated if Falseendif expression1 commands evaluated if expression1 is True elseif expression2 commands evaluated if expression2 is True elseif else commands evaluated if no other expression is True end 第10頁/共37頁第十一頁，編輯于星期六：十九點 二十分。 EPS=1; for num=1:1000 EPS=EPS/2; if (1+EPS) numnum = 53 EPS=1;

6、for num=1:1000 EPS=EPS/2; if (1+EPS)1continue end EPS=EPS*2 breakendEPS =2.2204e-016 numnum = 53 第11頁/共37頁第十二頁，編輯于星期六：十九點 二十分。第12頁/共37頁第十三頁，編輯于星期六：十九點 二十分。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,

7、Mark,Rank,Rank);for i=1:5 switch S(i).Marks case 100 S(i).Rank=滿分; case a S(i).Rank= 優(yōu)秀; case b S(i).Rank= 良好; case c S(i).Rank= 及格; otherwise S(i).Rank=不及格; endend 第13頁/共37頁第十四頁，編輯于星期六：十九點 二十分。disp(學(xué)生姓名 , 得分 , 等級);disp( )for i=1:5;disp(S(i).Name,blanks(6),num2str(S(i).Marks),blanks(6),S(i).Rank);en

8、d;學(xué)生姓名 得分 等級 Jack 72 及格Marry 83 良好Peter 56 不及格Rose 94 優(yōu)秀Tom 100 滿分 第14頁/共37頁第十五頁，編輯于星期六：十九點 二十分。第15頁/共37頁第十六頁，編輯于星期六：十九點 二十分。N=4;A=magic(3);try A_N=A(N,:),catch A_end=A(end,:),endA_end = 4 9 2lasterrans =Index exceeds matrix dimensions. 第16頁/共37頁第十七頁，編輯于星期六：十九點 二十分。第17頁/共37頁第十八頁，編輯于星期六：十九點 二十分。第18頁/

9、共37頁第十九頁，編輯于星期六：十九點 二十分。第19頁/共37頁第二十頁，編輯于星期六：十九點 二十分。circle.mfunction sa = circle(r,s)%CIRCLEplot a circle of radii r in the line specified by s.% r指定半徑的數(shù)值% s指定線色的字符串%sa圓面積% circle(r)利用藍實線畫半徑為 r 的圓周線.% circle(r,s)利用串 s 指定的線色畫半徑為 r 的圓周線.% sa=circle(r)計算圓面積，并畫半徑為 r 的藍色圓面.% sa=circle(r,s)計算圓面積，并畫半徑為 r

10、的 s 色圓面.% 編寫于1999年4月7日，修改于1999年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頁第二十一頁，編輯于星期六：十九點 二十分。第21頁/共37頁第二十二頁，編輯于星期六：十九點 二十分。第22頁/共37頁第二十三頁，編輯于星期六：十九點 二十分。指令含義nargin在函數(shù)體內(nèi)，用于獲取實際輸入?yún)⒘縩ar

11、gout在函數(shù)體內(nèi)，用戶獲取實際輸出參量nargin(fun)獲取fun指定函數(shù)的標稱輸入?yún)⒘繑?shù)nargout(fun)獲取fun指定函數(shù)的標稱輸出參量數(shù)inputname(n) 在函數(shù)體內(nèi)使用，給出第 n個輸入?yún)⒘康膶嶋H調(diào)用變量名第23頁/共37頁第二十四頁，編輯于星期六：十九點 二十分。第24頁/共37頁第二十五頁，編輯于星期六：十九點 二十分。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*pi;ss=

12、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頁第二十六頁，編輯于星期六：十九點 二十分。第26頁/共37頁第二十七頁，編輯于星期六：十九點 二十分。assigninxyq.mfunction y=assigninxyq(x)y=sqrt(x);

13、t=x2;assignin(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頁第二十八頁，編輯于星期六：十九點 二十分。第28頁/共37頁第二十九頁，編輯于星期六：十九點 二十分。第29頁/共37頁第三十頁，編輯于星期六：十九點 二十分。第30頁/共37頁第三十一頁，編輯于星期六：十九點 二十分。y=eval(CEM) 執(zhí)行CEM指定的計算y1,y2,=eval(CEM) 執(zhí)行CEM代表的函數(shù)文件，并輸出結(jié)果 第3

14、1頁/共37頁第三十二頁，編輯于星期六：十九點 二十分。t=pi;eval(theta=t/2,y=sin(theta)theta = 1.5708y = 1第32頁/共37頁第三十三頁，編輯于星期六：十九點 二十分。y1,y2,=feval(FH,arg1,arg2,)執(zhí)行函數(shù)句柄FH指定的計算y1,y2,=feval(FN,arg1,arg2,)執(zhí)行函數(shù)名字符串FN指定的計算y1,y2,=feval(FIL,arg1,arg2,)執(zhí)行內(nèi)聯(lián)函數(shù)FIL指定的計算 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). rand(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頁第