matlab中的ichol,matlab实验报告.doc

matlab实验报告.doc

学号姓名班级指导教师王朋飞MATLAB编程技术实验报告实验一熟悉MATLAB及其工具箱实验目的和要求实验目的是熟悉MATLAB界面和HELP功能。打开MATLAB软件，查看窗口布局，命令窗口的使用方法，HELP命令的使用方法，查看HELP文档。上机实验代码打开MATLAB软件，使用HELP命令查看ABS的说明，进入HELP浏览器查看工具箱。有很多TOOLBOX上机实验讨论和总结点击REFERENCEPAGEINHELPBROWSERDOCABS可以进一步进入MATLAB的HELP浏览器，可以进一步搜索和链接相关的函数。列出MATLAB常用的基本数学函数ABSX纯量的绝对值或向量的长度ANGLEZ复数Z的相角PHASEANGLESQRTX开平方REALZ复数Z的实部IMAGZ复数Z的虚部CONJZ复数Z的共轭复数ROUNDX四舍五入至最近整数FIXX无论正负，舍去小数至最近整数FLOORX地板函数，即舍去正小数至最近整数CEILX天花板函数，即加入正小数至最近整数RATX将实数X化为分数表示RATSX将实数X化为多项分数展开SIGNX符号函数SIGNUMFUNCTION。当X0时，SIGNX1。REMX,Y求X除以Y的馀数GCDX,Y整数X和Y的最大公因数LCMX,Y整数X和Y的最小公倍数EXPX自然指数POW2X2的指数LOGX以E为底的对数，即自然对数或LOG2X以2为底的对数LOG10X以10为底的对数实验二M文件的写法实验目的和要求实验目的是掌握如何编写M文件，有哪几种M文件。实际编写2个M文件实现函数调用，完成找1000以内的素数的功能。上机实验代码子函数FUNCTIONRESULTMYPRIMENRESULT1FORI2SQRTNIFMODN,I0RESULT0BREAKELSERESULT1ENDENDIFN2RESULT1END主函数调用子函数CLEARALLCLOSEALLCLCMARKZEROS1,1000FORJ21000MARKJMYPRIMEJENDNUMSUMMARKN11000INDEXFINDMARK1DISP 1000以内有 ,NUM2STRNUM, 个素数。 , 如下 DISPNINDEX结果1000以内有168个素数。上机实验讨论和总结通过编写MYPRIME子函数和主函数MAINPRIME。了解了M文件有两种M脚本文件和M函数文件。M脚本文件没有参数传递功能，但M函数文件有此功能。M函数文件的命名必须是其函数名，不可改变。M脚本文件则为完成一固定功能的模块，其运行时产生的变量均为全局变量，区别于M函数的局部变量，并且没有参数传递。M函数文件的格式有严格规定，它必须以“FUNCTION”开头，其格式如下FUNCTION输出变量函数名称(输入变量)语句；END；通过HELPPRIMES命令，知道MATLAB自带了找素数功能的函数PRIMES。实际应用中，直接调用就可以了。实验三矩阵分解编程实验目的和要求矩阵分解有很多种方法，以DOOLITTLE分解为例，理解DOOLITTLE分解的方法，实现DOOLITTLE分解的编程实现，掌握矩阵的基本操作方法。上机实验代码FUNCTIONL,UDOOLITTLEA矩阵DOOLITTLE分解M,NSIZEAIFMNERROR MATRIXSHOULDBEASQUAREMATRIX ENDIFMRANKAERROR MATRIXSHOULDBEFULLRANK ENDLEYEN,NUZEROSN,NFORK1NFORJKNTEMPSUM0FORT1K1TEMPSUMTEMPSUMLK,TUT,JENDUK,JAK,JTEMPSUMU的第K行ENDFORIK1NTEMPSUM0FORT1K1TEMPSUMTEMPSUMLI,TUT,KENDLI,KAI,KTEMPSUM/UK,KL的第K列ENDENDEND上机实验讨论和总结DOOLITTLE分解是从高斯消元法演化而来，认识其本质是高斯消元法的引申。MATLAB中自带了矩阵分解的函数。MATLAB中LU分解，L1,U1LUX的结果L是下三角的置换矩阵即L1PL2，U是上三角阵；书上用的是紧凑格式公式，可以看出方法不同结果是不唯一的。MATLAB中常用的矩阵分解函数还有LU,CHOL,ILU,QR,ICHOL,LDL,LUINC等。实验四多项式运算与绘图实验目的和要求学会用多项式对给定数据做曲线拟合，对下列数据进行曲线拟合(5次多项式)并绘制图形显示拟合效果。上机实验代码CLEARALLCLOSEALLCLCX07Y005,048,085,1,09,06,015,035PPOLYFITX,Y,5C0017DPOLYVALP,CPLOTX,Y, RO ,C,D, K 结果上机实验讨论和总结在MATLAB中，用多项式的系数来表示高次多项式；POLYVAL用来计算多项式在某点的值，用POLYFIT做多项式拟合。多项式拟合时要适当选择多项式的次数，如果选择的次数过高会出现过拟合的现象，如果选择次数过低，会出现拟合误差大的问题。其实在物理课本上我们很少看到高于4次X01234567Y00504808510906015035方的公式。通常，曲线拟合首先在5次方内尝试。