matlab第二次实验报告,matlab第二次上机实验报告.doc

matlab第二次上机实验报告.doc

电子科技大学电子工程学院标准实验报告(实验)课程名称MATLAB与数值分析学生姓名学号指导教师1、实验名称实验二线性方程组求解和函数的数值逼近2、实验目的通过上机实验，使学生对病态问题、线性方程组求解和函数的数值逼近方法有一个初步的理解。实验涉及的核心知识点病态方程求解、矩阵分解和方程组求解、LAGRANGE插值。实验重点与难点算法设计和MATLAB编程3、实验内容1对高阶多项式20112KPXXX编程求下面方程的解190PX并绘图演示方程的解与扰动量的关系。2对，生成对应的HILBERT矩阵，计算矩阵的条件数；通过先确定解获得常向量20NB的方法，确定方程组NHXB最后，用矩阵分解方法求解方程组，并分析计算结果。3对函数215FX1,X的CHEBYSHEV点，COS21KXN,2KN编程进行LAGRANGE插值，并分析插值结果。4、实验数据及结果分析1对高阶多项式20112KPXXX编程求下面方程的解190PX并绘图演示方程的解与扰动量的关系。P1,1FORI220N1,IPCONVP,N求多项式乘积ENDMZEROS1,21M的最高次幂为20，有21项HOLDONX120D1,0,01,05,1FORI15DELTDIM2DELTYROOTSPM 求多项式的根PLOTX,Y, O , COLOR ,I/5,I/20,I/10ENDTITLE 方程PX0的解与扰动量DELT的关系 LEGEND DELT1 , DELT0 , DELT01 , DELT05 , DELT1 024681012141618200102030405060与与PX0与与与与与与DELT与与与DELT1DELT0DELT01DELT05DELT12对，生成对应的HILBERT矩阵，计算矩阵的条件数；通过先确定解获20N得常向量B的方法，确定方程组NHXB最后，用矩阵分解方法求解方程组，并分析计算结果。FORN220HHILBNFPRINTF \N\NN10DCONDHN∞D ,N,CONDH, INF COND求矩阵范数X1NBHX LULUHLU分解XU\L\B利用LU分解求线性其次方程组的根XX FPRINTF \NX’ FORI1NFPRINTF 82F ,XIENDFPRINTF \NX’ FORI1NFPRINTF 82F ,XIENDEND输出结果如下N2CONDHN∞2700000E001X’100200X’100200N3CONDHN∞7480000E002X’100200300X’100200300N4CONDHN∞2837500E004X’100200300400X’100200300400N5CONDHN∞9436560E005X’100200300400500X’100200300400500N6CONDHN∞2907028E007X’100200300400500600X’100200300400500600N7CONDHN∞9851949E008X’100200300400500600700X’100200300400500600700N8CONDHN∞3387279E010X’100200300400500600700800X’100200300400500600700800N9CONDHN∞1099652E012X’100200300400500600700800900X’100200300400500600700800900N10CONDHN∞3535369E013X’1002003004005006007008009001000X’1002003004005006007008009001000N11CONDHN∞1229476E015X’10020030040050060070080090010001100X’10020030040050060070179990110001100WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND2692153E017INCONDAT48INUNTITLED7AT3N12CONDHN∞3714499E016X’100200300400500600700800900100011001200X’10020030040050358773473897294811221196WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND2739612E018INCONDAT48INUNTITLED7AT3N13CONDHN∞3650152E017X’1002003004005006007008009001000110012001300X’100200300406450873256302825844613128421061149WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND2448199E019INCONDAT48INUNTITLED7AT3N14CONDHN∞4084635E018WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND4455948E017INUNTITLED7AT7X’10020030040050060070080090010001100120013001400X’1002002974442031585101513226598505386200542165185671562WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND1024999E018INCONDAT48INUNTITLED7AT3N15CONDHN∞9756108E017X’100200300400500600700800900100011001200130014001500X’100200300406446892248251204424721053910503954117822026WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND9721674E019INCONDAT48INUNTITLED7AT3N16CONDHN∞1028629E018WARNINGMATRIXISCLOSETOSINGULARORBADLYSCALEDRESULTSMAYBEINACCURATERCOND7948463E017INUNTITLED7AT7X’1002003004005006007008009001000110012001300140015001600X’1002003004034727621072131546061734398661240453