资源描述:
《实验四 稀疏矩阵及其他特殊数组具体操作》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、实验四稀疏矩阵及其他特殊数组操作一、实验目的1.掌握生成多维数组的方法及其标识2.掌握稀疏矩阵的创建及运算3.使用元胞数组和结构数组二、实验内容1.多维数组生成多维数组可以直接输入元素赋值生成,也可以由低维数组或由函数生成。>>a=1:9>>b=reshape(a,3,3)>>c=cat(3,b,b)>>c(18)=[]%删除第18个元素查看三维数组c的元素存放顺序,可以看出三维数组是把第3维视做1页,先存放第1页的元素,在1页中先存放第1列的元素,再存放第2列的元素。练习:1>通过“全下标”元素赋值方式创建一个3行2列3页的三维数组。2>由函数ones,zeros,rand和randn直接
2、创建2行3列2页的三维数组。3>已知三维数组A。A(:,:,1)=65324425A(:,:,2)=32364564(1)将三维数组A中第13个元素的重新赋值为1。用单下标及全下标两种方式赋值。A(13)=1A(1,3,2)=1(2)将三维数组A中第2行第4列所有页的元素重新赋值为9。A(2,4,:)=9(1)求数组A各维的大小以及返回数组A行数或列数的最大值。[x,y,z]=size(A)length(A)(2)将该数组A中第10个元素删除,观察数组A的变化。A(10)=[]2>创建三维数组B,第一页为,第二页为,第三页为。重排生成数组C为2行,3列,2页。B(2,2,1)=[46;78]
3、B(2,2,2)=[35;10]B(2,2,3)=[19;23]C=reshape(B,2,3,2)2、稀疏矩阵(1)创建稀疏矩阵>>s=sparse([122334],[112233],[123456])(2)将稀疏矩阵与全元素矩阵转换>>f=full(s)>>k=f+s%稀疏矩阵与全元素矩阵的运算,注意结果的显示方式。>>h=f.*s%稀疏矩阵与全元素矩阵的运算,注意结果的显示方式。>>whos%查看各变量(3)用图形观看稀疏矩阵的非零元素分布情况。>>nnz(s)%返回非零元素个数>>spy(s)(4)根据表示稀疏矩阵的矩阵A(A为具有稀疏特征的带状矩阵),产生一个稀疏存储方式矩阵C。
4、>>B=[01112;02122;03132;41420;51520]>>d=[-3;0;3]>>C=spdiags(B,d,5,6)练习:1>设,将X转化为稀疏存储方式。sparse([141335],[124455],[419210-5])2>根据表示稀疏矩阵的矩阵,用spconvert()函数产生一个稀疏存储方式矩阵B。A=[221;31–1;433;538;6612]B=spconvert(A)3>创建1个5*5的稀疏矩阵A=,一个5*5全元素随机矩阵B,计算C=A+B,D=A.*B,查看A、B、C、D的存储空间。speye(5,5)+sparse([444445],[123455]
5、,[-234060])4>已知线性方程组AX=b,其中,,运用稀疏存储矩阵的方式求其解。B=[-1–1–1–10;22222;0–1–1–1-1]d=[-1;0;1]b=[1;0;0;0;0]A=spdiags(B,d,5,5)x=Ab3、元胞数组和结构数组的使用元胞数组和结构数组的使用举例如下。(1)创建结构数组表示3个学生的成绩。结构数组中有三个域分别为:姓名(name),学号(Id),成绩(scores)。>>student(1)=struct(‘name’,’John’,’Id’,’20030115’,’scores’,[85,96,74,82,68])>>student(2)=s
6、truct(‘name’,’Rose’,’Id’,’20030102’,’scores’,[95,93,84,72,88])>>student(3)=struct(‘name’,’Billy’,’Id’,’20030117’,’scores’,[72,83,78,80,83])(2)修改学生2的第2个成绩为73。>>student(2).scores(2)=73;>>student(2)(3)使用setfield命令进行上述修改。student=setfield(student,{2},’scores’,{2},73)(4)将scores域显示并计算平均成绩。>>all_scores=cat
7、(1,student.scores)>>average_scores=mean(all_scores)%mean为平均值函数,按列计算。(5)将平均成绩放在元胞数组中,使用3种方法创建元胞数组。方法1:>>average={‘平均成绩’,average_scores}方法2:>>averageA(1)={‘平均成绩’}>>averageA(2)={average_scores}方法3:>>averageB{1