数值计算报告教程.doc

数值计算 题目名称 数值计算 学 院 机电工程学院 专业班级 13微电子制造2班 学 号 3113000453 姓 名 肖铠斌 2016-12-30 摘 要 学习数值计算引论，应用MATLAB编写Gauss 列主元消元法求解线性方程组、多项式插值、Gauss 积分方法、Euler 方法求常微分方程初值问题、Newton 迭代法求非线性方程的程序。 【关键词】：Gauss 列主元消元法求解线性方程组、多项式插值、Gauss 积分方法、Euler 方法求常微分方程初值问题、Newton 迭代法求非线性方程 目录 TOC \o 1-2 \h \u  HYPERLINK \l _Toc310 1 Gauss 列主元消元法求解线性方程组  PAGEREF _Toc310 1  HYPERLINK \l _Toc18208 1.1问题描述  PAGEREF _Toc18208 1  HYPERLINK \l _Toc31616 1.2计算程序  PAGEREF _Toc31616 1  HYPERLINK \l _Toc28090 1.3算例 2  HYPERLINK \l _Toc10582 2 Lagrange多项式插值 3  HYPERLINK \l _Toc18208 2.1问题描述 3  HYPERLINK \l _Toc31616 2.2计算程序 3  HYPERLINK \l _Toc28090 2.3算例 3  HYPERLINK \l _Toc18514 3 Gauss 积分方法 4  HYPERLINK \l _Toc18208 3.1.问题描述 4  HYPERLINK \l _Toc31616 3.2计算程序 4  HYPERLINK \l _Toc28090 3.3算例 4  HYPERLINK \l _Toc15283 4 Euler 方法求常微分方程初值问题 5  HYPERLINK \l _Toc18208 4.1.问题描述 5  HYPERLINK \l _Toc31616 4.2计算程序 4  HYPERLINK \l _Toc28090 4.3算例 4  HYPERLINK \l _Toc20548 5Newton 迭代法求非线性方程 6  HYPERLINK \l _Toc18208 5.1.问题描述 6  HYPERLINK \l _Toc31616 5.2计算程序 6  HYPERLINK \l _Toc28090 5.3算例 7  HYPERLINK \l _Toc20548 参考文献 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向量记录每列最大值的行号