数值计算报告.doc

数值计算 题目名称 数值计算 学 院 机电工程学院 专业班级 13微电子制造2班 学 号 3113000453 姓 名 肖铠斌 2016-12-30 摘 要 学习数值计算，的程序。 1 1.1问题描述 1 1.2计算程序 1 1.3算例 2 2 Lagrange多项式插值 3 2.1问题描述 3 2.2计算程序 3 2.3算例 3 3 Gauss 积分方法 4 3.1.问题描述 4 3.2计算程序 4 3.3算例 4 4 Euler 方法求常微分方程初值问题 5 4.1.问题描述 5 4.2计算程序 4 4.3算例 4 5Newton 迭代法求非线性方程 6 5.1.问题描述 6 5.2计算程序 6 5.3算例 7 参考文献 8  一 、编写 Gauss 列主元消元法求解线性方程组的程序，要求附有算例。 1、问题描述：计算机中运算的时候常会碰到两个问题。? 1）一旦遇到某个主元等于0，消元过程便无法进行下去。? 2）在长期使用中还发现，即使消元过程能进行下去，但是当某个主元的绝对值很小时，求解出的结果与真实结果相差甚远。 2、程序： % 输出的量：系数矩阵和增广矩阵的秩RA,RB, 方程组中未知量的个数n和有关方程组解及其解X的信息. function [RA,RB,n,X]=gaus(A,b) B=[A b]; n=length(b); RA=rank(A); RB=rank(B); zhica=RB-RA; %求矩阵A、B的秩，并作比较 if zhica0, disp(RA~=RB，此方程组无解.) return end if RA==RB if RA==n disp(RA=RB=n，此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1); %生成全零矩阵 for p= 1:n-1 [Y,j]=max(abs(B(p:n,p))); %返回行向量Y和j，y向量记录abs(B(p:n,p)) 的每列的最大值，j向量记录每列最大值的行号 C=B(p,:); %把矩阵B中第p行所有列的值全都赋给矩阵C for k=p+1:n m= B(k,p)/B(p,p); %消元计算 B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1); end end b=B(1:n,n+1); A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 %从n-1循环到1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q); %回代计算 end else disp(RA=RBn，此方程组有无穷多解.) end End 3、算例： 求解线性方程组 在命令窗口输入： A=[5 2 1;2 8 -3;1 -3 -6]; b=[8;21;1]; [RA,RB,n,X] =gaus (A,b) RA=RB=n，此方程组有唯一解. RA = 3 RB = 3 n = 3 X = 1 2 -1 二 编写多项式插值的程序，要求附有算例 1、问题描述：Lagrange差值法满足插值条件的、次数不超过n的多项式是存在而且是唯一的。 2、程序： function yy=lagrange(x,y,x0) //定义拉格朗日差值函数 xlen=length(x); //X向量的维数 ylen=length(y); //y向