matlabR2012a课后习题答案第三章.doc

数值数组及其运算 习题3及解答 要求在闭区间上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。 〖目的〗 数值计算中产生自变量采样点的两个常用指令的异同。 〖解答〗 %方法一 t1=linspace(0,2*pi,10) %方法二 t2=0:2*pi/9:2*pi %要注意采样间距的选择，如这里的2*pi/9. t1 = Columns 1 through 7 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 10 4.8869 5.5851 6.2832 t2 = Columns 1 through 7 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 Columns 8 through 10 4.8869 5.5851 6.2832 由指令rng(default),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。 〖目的〗 数组下标的不同描述：全下标和单下标。 sub2ind, int2str, disp的使用。 随机发生器的状态控制：保证随机数的可复现性。 〖解答〗 rng(default) A=rand(3,5) [ri,cj]=find(A0.5); id=sub2ind(size(A),ri,cj); ri=ri;cj=cj; disp( ) disp(大于0.5的元素的全下标) disp([行号 ,int2str(ri)]) disp([列号 ,int2str(cj)]) disp( ) disp(大于0.5的元素的单下标) disp(id) A = 0.8147 0.9134 0.2785 0.9649 0.9572 0.9058 0.6324 0.5469 0.1576 0.4854 0.1270 0.0975 0.9575 0.9706 0.8003 大于0.5的元素的全下标 行号 1 2 1 2 2 3 1 3 1 3 列号 1 1 2 2 3 3 4 4 5 5 大于0.5的元素的单下标 1 2 4 5 8 9 10 12 13 15 采用默认全局随机流，写出产生长度为1000的“等概率双位（即取-1，+1）取值的随机码”程序指令，并给出 -1码的数目。 〖目的〗 两种基本随机发生器的使用。 关系运算产生逻辑数组——可用于数组的元素的标识和寻访。 逻辑数组的应用。 如何判断两个整数数组是否相等。 〖解答〗 （1）运用均匀随机数解题法——解法1 rng default %为以下结果重现而设；产生默认随机流。详见第4.3.2节 A=rand(1,1000); a=2*(A0.5)-1; Na=sum(a==-1) Na = 512 （2）运用正态随机数解题法——解法2 randn(state,123) B=randn(1,1000); b=2*(B0)-1; Nb=sum(b==-1) Nb = 462 （3）直接发生法——解法3 c=randsrc(1,1000,[-1,1]); Nc=sum(c==-1) Nc = 482 已知矩阵，运行指令B1=A.^(0.5), B2=A^(0.5), 可以观察到不同运算方法所得结果不同。（1）请分别写出根据B1, B2恢复原矩阵A的程序。（2）用指令检验所得的两个恢复矩阵是否相等。 〖目的〗 数组运算和矩阵运算的不同。 如何判断两个双精度数组是否相等。 norm指令的应用。 〖解答〗 A=[1,2;3,4]; B1=A.^0.5 B2=A^0.5 A1=B1.*B1; A2=B2*B2; norm(A1-A2,fro) % 求误差矩阵的F-范数，当接近eps量级时，就认为实际相等 B1 = 1.0000 1.4142 1.7321 2.0000 B2 = 0.5537 + 0.4644i 0.8070 - 0.2124i 1.2104 - 0.3186i 1.7641 + 0.1458i ans = 8.4961e-016