matlab数值数组及其运算.ppt

《matlab数值数组及其运算.ppt》由会员分享，可在线阅读，更多相关《matlab数值数组及其运算.ppt（40页珍藏版）》请在三一办公上搜索。

1、第3章数值数组及其运算,一维数组二维数组高维数组特殊数组关系操作和逻辑操作,3.1 简单数组,x=0.1*pi.2*pi.3*pi.4*pi.5*pi.6*pi.7*pi.8*pi.9*pi pix=Columns 1 through 7 0 0.3142 0.6283 0.9425 1.2566 1.5708 1.8850 Columns 8 through 11 2.1991 2.5133 2.8274 3.1416 y=sin(x)y=Columns 1 through 7 0 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511 Columns 8 thr

2、ough 11 0.8090 0.5878 0.3090 0.0000,3.2 一维数组的创建和寻访,3.2.1 一维数组的创建逐个元素输入法冒号生成法 x=a:inc:b若(b-a)是inc的整数倍，则所生成数组的最后一个元素等于 b，否则小于 b。inc省略时，默认其取值为 1。inc可以取正数或负数。inc取正时，要保证ab。定数线性采样法 x=linspace(a,b,n),x=(0:0.1:1)*pix=Columns 1 through 7 0 0.3142 0.6283 0.9425 1.2566 1.5708 1.8850 Columns 8 through 11 2.1991

3、 2.5133 2.8274 3.1416,x=linspace(0,pi,11)x=Columns 1 through 7 0 0.3142 0.6283 0.9425 1.2566 1.5708 1.8850 Columns 8 through 11 2.1991 2.5133 2.8274 3.1416,3.2.2一维数组的子数组寻访和赋值,x(3)%寻访数组x的第三个数组ans=0.6283x(1 2 5)%寻访数组x的第1,2,5个元素组成的子数组ans=0 0.3142 1.2566 x(1:3)%寻访数组x前三个元素组成的子数组ans=0 0.3142 0.6283,x(3:en

4、d)%数组x从第三个元素开始到最后一个元素ans=Columns 1 through 7 0.6283 0.9425 1.2566 1.5708 1.8850 2.1991 2.5133 Columns 8 through 9 2.8274 3.1416 x(3:-1:1)%寻访数组x前三个元素倒排构成的子数组ans=0.6283 0.3142 0 x(1 2 3 4 4 3 2 1)%对元素可进行重复寻访ans=Columns 1 through 7 0 0.3142 0.6283 0.9425 0.9425 0.6283 0.3142 Column 8 0,x(3)=0 x=Columns

5、 1 through 7 0 0.3142 0 0.9425 1.2566 1.5708 1.8850 Columns 8 through 11 2.1991 2.5133 2.8274 3.1416 x(1 4)=1 1x=Columns 1 through 7 1.0000 0.3142 0 1.0000 1.2566 1.5708 1.8850 Columns 8 through 11 2.1991 2.5133 2.8274 3.141,3.3二维数据的创建,整个输入数组必须以方括号“”为其首尾；数组的行与行之间必须用分号“;”或回车键【Enter】隔开；数组元素必须由逗号“,”或空格

6、分隔。,a=2.7358;b=33/79;C=1,2*a+i*b,b*sqrt(a);sin(pi/4),a+5*b,3.5+i C=1.0000 5.4716+0.4177i 0.6909 0.7071 4.8244 3.5000+1.0000i,3.4二维数组元素的标识,“全下标”标识两个下标组成：行下标，列下标。如A(3,5)“单下标”标识只用一个下标来指明元素在数组中的位置“逻辑1”标识诸如寻找数组中所有大于某值的元素的问题,“单下标”与“全下标”的转换关系：以（mn）的二维数组A为例，若“全下标”元素位置是“第r行，第 c列”，那么相应的“单下标”为 l=(c-1)*m+r。sub2

7、ind根据全下标算出单下标ind2sub根据单下标算出全下标,A=1 2 3;4 5 6;7 8 9A=1 2 3 4 5 6 7 8 9,A(2,3)ans=6sub2ind(3 3,2,3)ans=8,A(7)ans=3 i,j=ind2sub(3 3,7)i=1j=3,A=-4-2 0 2 4-3-1 1 3 5L=abs(A)3L=1 0 0 0 1 0 0 0 0 1islogical(L)ans=1X=A(L)X=-4 4 5,3.5二维数组的子数组寻访和赋值,A=zeros(2,4)A=0 0 0 0 0 0 0 0 A(:)=1:8%全元素赋值方式A=1 3 5 7 2 4 6

8、 8 s=2 3 5%产生单下标数组行数组s=2 3 5A(s)ans=2 3 5,Sa=10 20 30%Sa是长度为3的“列数组”Sa=10 20 30A(s)=Sa%单下标方式赋值A=1 20 30 7 10 4 6 8 A(:,2 3)=ones(2)把A的第2、3列元素全赋值为1A=1 1 1 7 10 1 1 8,3.6数组运算和矩阵运算,3.6.1标量数组运算,g=1 2 3 4 5 6 7 8 9 10 11 12 g-2ans=-1 0 1 2 3 4 5 6 7 8 9 10,2*g-1ans=1 3 5 7 9 11 13 15 17 19 21 23 2*g/5+1an

9、s=1.4000 1.8000 2.2000 2.6000 3.0000 3.4000 3.8000 4.2000 4.6000 5.0000 5.4000 5.8000,3.6.2数组数组运算,g=1 2 3 4 5 6 7 8 9 10 11 12h=1 1 1 1 2 2 2 2 3 3 3 3,g+hans=2 3 4 5 7 8 9 10 12 13 14 15 2*g-hans=1 3 5 7 8 10 12 14 15 17 19 21,g.*hans=1 2 3 4 10 12 14 16 27 30 33 36,g*h?Error using=*Inner matrix di

10、mensions must agree.,g./hans=1.0000 2.0000 3.0000 4.0000 2.5000 3.0000 3.5000 4.0000 3.0000 3.3333 3.6667 4.0000,3.6.3数组运算和矩阵运算的对照,A=zeros(2,3);A(:)=1:6;A=A*(1+i)A_A=A.A_M=AA=1.0000+1.0000i 3.0000+3.0000i 5.0000+5.0000i 2.0000+2.0000i 4.0000+4.0000i 6.0000+6.0000iA_A=1.0000+1.0000i 2.0000+2.0000i 3.

11、0000+3.0000i 4.0000+4.0000i 5.0000+5.0000i 6.0000+6.0000iA_M=1.0000-1.0000i 2.0000-2.0000i 3.0000-3.0000i 4.0000-4.0000i 5.0000-5.0000i 6.0000-6.0000i,3.6.4标准数组生成函数,ones(3)ans=1 1 1 1 1 1 1 1 1 ones(1,2)ans=1 1 magic(3)ans=8 1 6 3 5 7 4 9 2,a=1 2 3 4a=1 2 3 4 diag(a)ans=1 0 0 0 0 2 0 0 0 0 3 0 0 0 0

12、 4 diag(a,2)ans=0 0 1 0 0 0 0 0 0 2 0 0 0 0 0 0 3 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0,3.6.5数组操作函数,a=-4:4a=-4-3-2-1 0 1 2 3 4 A=reshape(a,3,3)A=-4-1 2-3 0 3-2 1 4,flipud(A)ans=-2 1 4-3 0 3-4-1 2 fliplr(A)ans=2-1-4 3 0-3 4 1-2,3.7数组构作技法综合,A=1 2 3 4 5 6 7 8 9,A(3,3)=0A=1 2 3 4 5 6 7 8 0,A(2,6)=1A=1 2

13、 3 0 0 0 4 5 6 0 0 1 7 8 0 0 0 0,A(:,4)=4A=1 2 3 4 0 0 4 5 6 4 0 1 7 8 0 4 0 0,B=A(3:-1:1,1:3)或B=A(end:-1:1,1:3)B=7 8 0 4 5 6 1 2 3 B=A(3:-1:1,:)B=7 8 0 4 0 0 4 5 6 4 0 1 1 2 3 4 0 0 C=A B(:,1 3)C=1 2 3 4 0 0 7 0 4 5 6 4 0 1 4 6 7 8 0 4 0 0 1 3,3.8“非数”和“空”数组,3.8.1非数NaNNaN性质：NaN参与运算所得的结果也是NaN，即具有传递性。

14、非数没有“大小”概念，因此不能比较两个非数的大小。非数的功用：真实记述0/0,/,0运算的后果；避免可能因0/0,/,0运算而造成程序执行的中断；在数据可视化中，用来裁剪图形。,a=0/0Warning:Divide by zero.(Type warning off MATLAB:divideByZero to suppress this warning.)a=NaN isnan(a)ans=1,3.8.2“空”数组某维长度为 0或若干维长度均为 0的数组都是“空”数组。“空”数组的功用：在没有“空”数组参与的运算时，计算结果中的“空”可以合理地解释“所得结果的含义”。运用“空”数组对其他非

15、空数组赋值，可以改变数组的大小，但不能改变数组的维数。,a=a=b=ones(2,0)b=Empty matrix:2-by-0 c=zeros(2,0)c=Empty matrix:2-by-0 isempty(a)ans=1 A=reshape(-4:5,2,5)A=-4-2 0 2 4-3-1 1 3 5,3.10关系操作和逻辑操作,在所有关系表达式和逻辑表达式中，作为输入的任何非0数都被看作是“逻辑真”，而只有0才被认为是“逻辑假”。关系表达式和逻辑表达式的计算结果，即输出，是一个由0和1组成的“逻辑数组”。在此数组中的1表示“真”，0表示“假”。逻辑数组是一种特殊的数值数组。与“数值

16、类”有关的操作和函数对它也同样适用；它又自身的特殊用途，如数组寻访等。,C,B=10-A,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;%该自变量数组中，存在0值y=sin(t)./t;%在t=0处，计算将产生NaNtt=t+(t=0)*eps;%使0元素被一个“机器零”小数代替yy=sin(tt)./tt;%用sin(eps)/eps近似代替sin(0)/0极限subplot(1,2,1),plot(t,y),a

17、xis(-7,7,-0.5,1.2),xlabel(t),ylabel(y),title(残缺图形)subplot(1,2,2),plot(tt,yy),axis(-7,7,-0.5,1.2),A=1:9,L1=(A5),L2=(A3)&(A7)A=1 2 3 4 5 6 7 8 9L1=1 1 1 1 1 0 0 0 0L2=0 0 0 1 1 1 0 0 0,t=linspace(0,3*pi,500);y=sin(t);%产生正弦波z1=(t2*pi).*y;%获得整流半波w=(tpi/3%获得消顶整流半波subplot(2,2,1),plot(t,y,:r),ylabel(y)subplot(2,2,2),plot(t,z1,:r),axis(0 10-1 1)subplot(2,2,3),plot(t,z2,-b),axis(0 10-1 1),3.11数组大小,A=reshape(1:8,2,4)A=1 3 5 7 2 4 6 8 s=size(A)s=2 4 r,c=size(A)r=2c=4,numel(A)ans=8 length(A)ans=4,