资源描述：
《数值分析上机报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、数值分析上机报告思想。Neiu=2;%%%初始值x0x0=input('inputinitialvaluex0>>');k=0;%迭代次数max=100;%最大迭代次数R=eval(subs(f,'x0','x'));%求解f(x0)，以确定初值x0时否就是解iu*eval(subs(f,'x0','x'))/eval(subs(df,'x0','x'));R=x1-x0;x0=x1;k=k+1;if(eval(sub
2、s(f,'x0','x'))<1e-10);breakendifk>max;%如果迭代次数大于给定值，认为迭代不收敛，重新输入初值ss=input('mayberesultiserror,chooseanep(ss,'y')x0=input('inputinitialvaluex0>>');k=0;elsebreakendendendk;%给出迭代次数x=x0；%给出解结果分析和讨论：x2?0在[1，2]内的根。(??5*10?6,下同)1.用二
3、分法计算方程sinx?2计算结果为x=1.40441513061523；f(x)=-3.797205105904311e-007；k=18；由f(x)知结果满足要求，但迭代次数比较多，方法收敛速度比较慢。2.用二分法计算方程x3?x?1?0在[1，1.5]内的根。计算结果为x=1.32471847534180；f(x)=2.209494846194815e-006；k=17；由f(x)知结果满足要求，但迭代次数还是比较多。3.用Neaxaiki?k(k?1)并将第r行和第k行的元素进行交换，以使得当前的akk的数值比0要大的多。这种列主元的
4、消去法的主要步骤如下：1.消元过程对k=1,2,…,n-1,进行如下步骤。1)选主元，记
5、ark
6、?aiki?k若
7、ark
8、很小，这说明方程的系数矩阵严重病态，给出警告，提示结果可能不对。2)交换增广阵A的r，k两行的元素。arj?akj(j=k,…,n+1)3)计算消元aij?aij?aikakj/akk(i=k+1,…,n;j=k+1,……,n+1)2.回代过程对k=n,n-1,…,1,进行如下计算篇三：数值分析上机实验报告-李宝君《数值分析实验报告》指导老师：代新敏姓名：李宝君学号：2011020917专业：机械制造及其自动化院系：机
9、械工程学院贵州大学2011级研究生第一题一程序说明：1Householder其基本思想是：利用初等反射阵H?E?2uuT，将矩阵的每一行向量变换成所需方向的行向量，从而最终得到想要的三对角阵。它的乘法运算次数仅是Givenr方法的一半，且只需要作n?2次开方运算。归纳起来，对换矩阵三对角化的算法步骤为：（1）令A0?A，aij?aij，已知Ar?1，即Ar?1?(aij)。（2）sr?((1)(r)i?r?1?(an(r)21/2ir))T(r)(r)(r)(r))??（3）?r?sr2?ar(r?1,rsr.，ur??0,?0,ar?1,
10、r?sign(ar?1,r)sr,ar?2,r?an,r?.（4）yr?Ar?1ur/?r。（5）kr?1Turyr/?r。2（6）qr?yr?krurTT（7）Ar?Ar?1?(qrur?urqr)，r?1,2,?n?2。2松弛法简称SOR法，其基本思想是在GS方法已求出x(m)，x(m?1)的基础上，经过重新组合而得到新的序列，而此新序列使收敛速度加快。其算法如下：i?1nx(m)i?(1??)x(m?1)i??(?bxj?1(m)ijj??bijx(jm?1)?gi)j?i?1当?>1时称为超松弛法。3用列主元素消去法求解BX=
11、b：n阶方程组的系数矩阵为：a11a12…a1nb1a21a22…a2nb2an1an2?annbGauss消去法的算法为：ij=aij/ajj(ajj!=0)j=1,2?n,i=j+1,j+2?n⑴ij=aij-lik-1ak-1i,j=k,k+1?n,k=2,3?ni=bi-lik-1bk-1i=k,k+1?n,k=2,3?n⑵xi=(bi-∑aijxj)/aiii=n,n-1?1,j=i+1,i+2,?n列主元素消去法是在Gauss消去法的基础上选主元，选取绝对值最大（或尽量大）的元素为主元，使lij绝对值很小。二计算程序:#incl
12、udestdafx.h#includestdio.h#includemath.h#defineMAX9voidSOR(double(*p)[MAX],intn);//定义一个函数
