实验2数组矩阵及数值计算matlab.docVIP

2013年上学期MATLAB语言上机实验内容 姓名： 班级： 学号：08 实验2: 数组（矩阵）及数值计算 1、指令 A=ones(4,4); B=A(3,[2 3]) C=A(3,:)1 写出输出结果： B = 1 1 C = 0 0 0 0 2、 指令B=magic(3)；写出求矩阵B的特征值D和特征向量V的语句： B=magic(3); [V,D]=eig(B) V = -0.5774 -0.8131 -0.3416 -0.5774 0.4714 -0.4714 -0.5774 0.3416 0.8131 D = 15.0000 0 0 0 4.8990 0 0 0 -4.8990 3、A=rand(100,100)；用向量化运算的方法，用1条语句实现求矩阵A中大于0.4的所有元素个数：A=rand(100,100); S=sum(sum(A0.4)) S = 5977 4要求在闭区间上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。 A = 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832 A=linspace(0,2*pi,10) A = 0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832 5、由指令rand(state,0),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。 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.95012928514718 0.48598246870930 0.45646766516834 0.44470336435319 0.92181297074480 0.23113851357429 0.89129896614890 0.01850364324822 0.61543234810009 0.73820724581067 0.60684258354179 0.76209683302739 0.82140716429525 0.79193703742704 0.17626614449462 大于0.5的元素的全下标 行号 1 3 2 3 3 2 3 1 2 列号 1 1 2 2 3 4 4 5 5 大于0.5的元素的单下标 1 3 5 6 9 11 12 13 14 6、在时间区间 [0,10]中，绘制曲线。要求分别采取“标量循环运算法”和“数组运算法”编写两段程序绘图。 t=linspace(0,10,200); N=length(t); y1=zeros(size(t)); for k=1:N y1(k)=1-exp(-0.5*t(k))*cos(2*t(k)); end subplot(1,2,1),plot(t,y1),xlabel(t),ylabel(y1),grid on y2=1-exp(-0.5*t).*cos(2*t); subplot(1,2,2),plot(t,y2),xlabel(t),ylabel(y2),grid on 7、先运行指令x=-3*pi:pi/15:3*pi; y=x; [X,Y]=meshgrid(x,y); warning off; Z=sin(X).*sin(Y)./X./Y； 产生矩阵Z。（1）请问矩阵Z中有多少个“非数”数据？（2）用指令surf(X,Y,Z); shading interp