【正文】 “非數(shù)”和“空”數(shù)組 非數(shù)NaN【】非數(shù)的產生和性質演示。（1）a=0/0,b=0*log(0),c=infinf Warning: Divide by zero.a = NaNWarning: Log of zero.b = NaNc = NaN （2）0*a,sin(a) ans = NaNans = NaN （3）a==nan ans = 0 （4）a~=nana==bbc ans = 1ans = 0ans = 0 （5）class(a)isnan(a) ans =doubleans = 1 【】非數(shù)元素的尋訪rand(39。state39。,0)R=rand(2,5)。R(1,5)=NaN。R(2,3)=NaN R = NaN NaN isnan(R) ans = 0 0 0 0 1 0 0 1 0 0 Linear_index=find(isnan(R))[r_index,c_index]=ind2sub(size(R),Linear_index)。disp(39。r_index c_index39。),disp([r_index c_index]) Linear_index = 6 9r_index c_index 2 3 1 5 “空”數(shù)組【】關于“空”數(shù)組的算例。（1）a=[],b=ones(2,0),c=zeros(2,0),d=eye(2,0),f=rand(2,3,0,4) a = []b = Empty matrix: 2by0c = Empty matrix: 2by0d = Empty matrix: 2by0f = Empty array: 2by3by0by4 （2）class(a)isnumeric(a)isempty(a) ans =doubleans = 1ans = 1 which a ndims(a) size(a) a is a variable.ans = 2ans = 0 0 （3）b_c1=b.*c b_c2=b39。*c b_c3=b*c39。 b_c1 = Empty matrix: 2by0b_c2 = []b_c3 = 0 0 0 0 （4）a==b ans = 0 b==c ans = Empty matrix: 2by0 cd ans = Empty matrix: 2by0 a==0 Warning: Future versions will return empty for empty == scalar parisons.ans = 0 a~=0 Warning: Future versions will return empty for empty ~= scalar parisons.ans = 1 （5）A=reshape(4:5,2,5) A = 4 2 0 2 4 3 1 1 3 5 L2=A10find(L2) L2 = 0 0 0 0 0 0 0 0 0 0ans = [] （6）A(:,[2,4])=[] A = 4 0 4 3 1 5 關系操作和邏輯操作 關系操作【】關系運算示例。A=1:9,B=10A,r0=(A4),r1=(A==B) A = 1 2 3 4 5 6 7 8 9B = 9 8 7 6 5 4 3 2 1r0 = 1 1 1 0 0 0 0 0 0r1 = 0 0 0 0 1 0 0 0 0 【】關系運算運用之一：求近似極限，修補圖形缺口。t=2*pi:pi/10:2*pi。y=sin(t)./t。tt=t+(t==0)*eps。yy=sin(tt)./tt。subplot(1,2,1),plot(t,y),axis([7,7,]),xlabel(39。t39。),ylabel(39。y39。),title(39。殘缺圖形39。)subplot(1,2,2),plot(tt,yy),axis([7,7,])xlabel(39。t39。),ylabel(39。yy39。),title(39。正確圖形39。) Warning: Divide by zero. 極限處理前后的圖形對照 邏輯操作【】邏輯操作示例。注意邏輯運算和關系運算之間的優(yōu)先級次序。（）A=3:3。L1=~(A0) L2=~A0 L3=~A L4=A2amp。A1 L1 = 1 1 1 1 0 0 0L2 = 0 0 0 1 0 0 0L3 = 0 0 0 1 0 0 0L4 = 0 0 1 1 0 0 0 【】邏輯操作應用之一：逐段解析函數(shù)的計算和表現(xiàn)。本例演示削頂整流正弦半波的計算和圖形繪制。t=linspace(0,3*pi,500)。y=sin(t)。 %處理方法一：z1=((tpi)|(t2*pi)).*y。 % 3w=(tpi/3amp。t2*pi/3)+(t7*pi/3amp。t8*pi/3)。% 4w_n=~w。 % 5z2=w*sin(pi/3)+w_n.*z1。 % 6subplot(1,3,1),plot(t,y,39。:r39。),ylabel(39。y39。)subplot(1,3,2),plot(t,z1,39。:r39。),axis([0 10 1 1])subplot(1,3,3),plot(t,z2,39。b39。),axis([0 10 1 1]) %處理方法二：z=(y=0).*y。 % 11a=sin(pi/3)。z=(y=a)*a+(ya).*z。 % 13plot(t,y,39。:r39。)。hold on。plot(t,z,39。b39。)xlabel(39。t39。),ylabel(39。z=f(t)39。),title(39。逐段解析函數(shù)39。)legend(39。y=sin(t)39。,39。z=f(t)39。),hold off 圖 關系、邏輯函數(shù)20 / 2