Matlab第3章数值数组及向量化运算.ppt

关系操作符： 4：关系操作和逻辑操作 = = == ~= 注意“==”和“=”的区别。a=5,b=6 a==b ans = 0 a=b a=6 ! 【例3.5-1】 A=1:9,B=10-A r0=(A4) r1=(A==B) A = 1 2 3 4 5 6 7 8 9 B = 9 8 7 6 5 4 3 2 1 r0 = 1 1 1 0 0 0 0 0 0 r1 = 0 0 0 0 1 0 0 0 0 4：关系操作和逻辑操作 【例3.5-2】 t=-3*pi:pi/10:3*pi; y=sin(t)./t; tt=t+(t==0)*eps; yy=sin(tt)./tt; subplot(1,2,1),plot(t,y),axis([-9,9,-0.5,1.2]), xlabel(t),ylabel(y),title(残缺图形) subplot(1,2,2),plot(tt,yy),axis([-9,9,-0.5,1.2]) xlabel(tt),ylabel(yy),title(正确图形) 4：关系操作和逻辑操作 逻辑操作符： 4：关系操作和逻辑操作 | ~ xor xor用法:两个元素有且只有一个为0,给出1,否则给出0 xor(5,6)=0 xor(5,0)=1 【例3.5-3】 （1）逻辑、关系操作的组合 A=[-2,-1,0,0,1,2,3] L1=~(A1) % L2=(A0)(A2) % A = -2 -1 0 0 1 2 3 L1 = 1 1 1 1 1 0 0 L2 = 0 0 0 0 1 0 0 4：关系操作和逻辑操作 （2）xor的功能 A,B=[0,-1,1,0,1,-2,-3] C=xor(A,B) % A = -2 -1 0 0 1 2 3 B = 0 -1 1 0 1 -2 -3 C = 1 0 1 0 0 0 0 4：关系操作和逻辑操作 【例3.5-4】见教材140页 t=linspace(0,3*pi,500);y=sin(t); %产生正弦波，第一张图 z1=((tpi)|(t2*pi)).*y; %第二张图 w=(tpi/3t2*pi/3)+(t7*pi/3t8*pi/3); wn=~w; % 第三张图 z2=w*sin(pi/3)+wn.*z1; % 第四张图 subplot(4,1,1),plot(t,y,:r),axis([0,10,-1.5,1.5]) ylabel(y),grid on subplot(4,1,2),plot(t,z1,:r),axis([0,10,-0.2,1.5]),ylabel(z1) subplot(4,1,3),plot(t,wn,:r),axis([0,10,-0.2,1.5]),ylabel(wn) subplot(4,1,4),plot(t,z2,-b),axis([0,10,-0.2,1.5]),ylabel(z2) xlabel(t) 4：关系操作和逻辑操作 第3堂课总结 数组运算的应用 “非数”和“空”数组，isnan，isempty 关系操作符和逻辑操作符 find，ind2sub,sub2ind,mean,length,size等函数的使用方法 第二次作业 已知矩阵 ，运行指令B1=A.^(0.5), B2=A^(0.5), 可以观察到不同运算方法所得结果不同（1）请分别写出根据B1, B2恢复原矩阵A的程序（2）用指令检验所得的两个恢复矩阵是否相等。 在时间区间 [0,10]中，绘制 曲线。要求分别采取“标量循环运算法”和“数组运算法”编写两段程序绘图。（注意：体验数组运算的简捷。） 先运行指令x=-3*p