matlab实验函数编写与程序设计,实验6matlab程序设计.doc

实验6matlab程序设计.doc

实验报告6学生姓名学号指导老师实验时间47地点综合实验大楼实验题目实验6MATLAB程序设计2实验目的掌握函数的编写与调用。掌握程序调式方法。实验内容1、编写程序实现函数M,R,CMYMAXA，求给定矩阵A的最大值M及最大值所在行列坐标R、C。2、三对角矩阵完全由向量A、B、C确定。编写函数12100NNNAACBTCTTRIDIAGONALMATRIXA,B,C，用于在给定A、B、C时构造矩阵三对角矩阵T，并调用TRIDIAGONALMATRIX构造如下矩阵，并计算Z51，观察Z51具有的特点(对称且中心对称541054Z)3、给定M和N，构造分块矩阵，G的大MMMZIIGNZI块小为，其中IM为M维的单位矩阵，ZM为与2题中的Z5类似M维矩阵。要N求函数头为GHHIM,N，并给出G33HHI3,3的结果，计算G33的逆，观察其特征。G33的逆参考如下4、课后探索题搜索资料，研究计算矩阵G的的快速算法，并编写程序实现之。实验指导代码示例1、3题中矩阵G的逆的计算程序FUNCTIONGINVERSEOFHHIM,NCOMPUTETHEINVERSIONOFHHI,DENOTEDBYGHHIZII000IZI000IZI000IZI000IZIZTMTM 3IMTM101111\\1110FOREXAMPLET41,0,0,0,01,1,0,0,00,1,1,0,00,0,1,1,00,0,0,1,0TMTM 1100012100012100012100011WROTEBY郑成勇ISPEYEM单位矩阵构造ZZDIAG4ZEROS1,M254R1M1IDRRMZID1ZIDM11YCELLN,1Y{1}0Y{2}INVZIFORR3NY{R}INVZY{R1}END实际上，G是中心对称的，因而只需计算主、副对角线上方的元素(含线上元素)GCELLNG{1,1}INVZIY{N}IDCFIXN1/2FORR2IDCG{R,R}INVZY{NR1}Y{R}ENDFORC2NFORR1IDCIFRCR2N坐标R,CR1位于副对角线下方RCR1N1BREAKENDIFRRC1N1G{R,CR1}G{R1,CR1}ELSEG{R,CR1}G{NCR2,NR}ENDIFR1G{R,CR1}ZIY{R}\G{R,CR1}ELSEG{R,CR1}ZY{R}\G{R,CR1}ENDENDEND2、根据输入成绩输出相应等级IF语句实现FUNCTIONFUNC1A 请输入你的成绩 IFA100|A90MSGBOX 你的成绩等级为A ELSEIFA80MSGBOX 你的成绩等级为B ELSEIFA70MSGBOX 你的成绩等级为C ELSEIFA60MSGBOX 你的成绩等级为D ELSEIFA0MSGBOX 你的成绩等级为E ENDSWICTH语句实现FUNCTIONFUNC11NDLG 请输入你的成绩 NSTR2NUMCHARNSWITCHFIXN/10CASE{10,9}MSGBOX 你的成绩等级为A CASE8MSGBOX 你的成绩等级为B CASE7MSGBOX 你的成绩等级为C CASE6MSGBOX 你的成绩等级为D CASENUM2CELL05MSGBOX 你的成绩等级为E OTHERWISEMSGBOX 请输入正确的分数 , ERROR , ERROR END3、函数递归调用计算斐波拉契数列FUNCTIONFFIBONACCIN用于求FIBONACCI数列的函数文件IFN2FFFIBN1FFIBN2ELSEF1END实验结果与分析1程序FUNCTIONM,R,CMYMAXAROW,COLSIZEAMA1,1R1C1FORI1ROWFORJ1COLIFAI,JMMAI,JRICJENDENDEND运行结果A12132356759167284965MYMAXAROW3COL5ANS28矩阵A的最大值为28，坐标为第三行第五列。2程序FUNCTIONTTRIDIAGONALMATRIXA,B,CNLENGTHATZEROSNFORR1NTR,RARENDFORR1N1TR,R1CRTR1,RBREND运行结果A45554B11111CBZ5TRIDIAGONALMATRIXA,B,CZ54100015100015100015100014INVZ5ANS02638005510011500024000060055102202004600009700024001150046002184004600011500024000970046002202005510000600024001150055102638可以看出Z5的逆对称且中心对称。3程序FUNCTIONGHHIM,NIMEYEMZMDIAG4ONES1,M244DIAGONESM1,1,1DIAGONESM1,1,1MNMNGZEROSMNFORR1MMNIFR1||RN1M1GRRM1,RRM1ZMIMELSEGRRM1,RRM1ZMENDENDFORR1MN1MGRRM1,RMR2M1IMGRMR2M1,RRM1IMEND运行结果G33HHI3,3G33310100000141010000013001000100410100010151010001014001000100310000010141000001013INVG33ANSCOLUMNS1THROUGH6042620139300595013930071400393013930346401393007140107100714005950139304262003930071401393013930071400393034640107100464007140107100714010710285701071003930071401393004640107103464005950039300262013930071400393003930046400393007140107100714002620039300595003930071401393COLUMNS7THROUGH9005950039300262003930046400393002620039300595013930071400393007140107100714003930071401393042620139300595013930346401393005950139304262观察G33的逆，可以将G33的逆按顺序分为9个33的矩阵，第一个33的矩阵关于03464中心对称，第二个33的矩阵关于01071中心对称，第三个33的矩阵关于00464中心对称，第四个33的矩阵关于01071中心对称，第五一个33的矩阵关于02857中心对称，第六个33的矩阵关于01071中心对称，第七个33的矩阵关于00464中心对称，第八个33的矩阵关于01071中心对称，第九个33的矩阵关于03464中心对称。