北京理工大学徐特立学院数值大作业上机实验预案.docx

北京理工大学徐特立学院数值分析课后上机实验选做 教材：数值计算方法（2011第一版）.丁丽娟,程杞元.高等教育出版社 ^以下代码作者原创^ 超链接： 1.2（ HYPERLINK \l 题目1_2 题目， HYPERLINK \l 原理1_2 原理， HYPERLINK \l 截图1_2 截图， HYPERLINK \l 程序1_2 代码） 2.2（ HYPERLINK \l 题目2_2 题目， HYPERLINK \l 原理2_2 原理， HYPERLINK \l 截图2_2 截图， HYPERLINK \l 程序2_2 代码） 3.1（ HYPERLINK \l 题目3_1 题目， HYPERLINK \l 原理3_1 原理， HYPERLINK \l 截图3_1 截图， HYPERLINK \l c程序3_1 c代码， HYPERLINK \l m程序3_1 m代码） 5.1（ HYPERLINK \l 题目5_1 题目， HYPERLINK \l 原理5_1 原理， HYPERLINK \l 截图5_1 截图， HYPERLINK \l 程序5_1 代码） 5.3（ HYPERLINK \l 题目5_3 题目， HYPERLINK \l 原理5_3 原理， HYPERLINK \l 截图5_3 截图， HYPERLINK \l 程序5_3 代码） 第一章：数值计算中的误差 2、题目简介：利用pi/4=1-1/3+1/5-1/7。。。级数计算pi的近似值。 输入：误差值 输出：求和项数，并输出pi值 工具：C语言 运行环境：VC-6.0 计算公式及原理：利用pi/4=1-1/3+1/5-1/7。。。级数计算pi的近似值，由数学原理可知误差会小于首次舍弃的项,可以编写循环实现。 程序运行结果截图： 程序代码：（c语言） #includestdio.h void main() { printf(第一章第2题求pi,欢迎使用，请按提示操作。\n); int i=1,n=0,k=1; double e,pi,er; printf(请输入误差(例如1e-4):); scanf(%lf,e); printf(请稍候。。。\n); er=e; pi=0; while(er=e) { pi+=k*1.0/i;k=-k;er=1.0/i; i+=2;n++; } pi*=4; printf(%d项求和后可以达到%.10lf精度，这时pi=%.10lf\n,n,e,pi); getchar(); getchar(); } 第二章： 解线性方程组的直接方法 2、题目简介：用MATLAB软件编程实现追赶法求解三对角方程组的算法，并考虑梯形电阻电路问题，电路如下： 其中电路中的各个电流{，，…，}须满足下列线性方程组： 设，，运用求各段电路的电流量。 输入：三列对角线元素向量，右侧常数元素向量 输出：三对角方程的解 工具：m语言 运行环境：MATLAB R2012.b 计算公式与原理： 上述方程组可用矩阵表示为： 根据三对角方程追赶法原理以及公式可解。 程序运行结果截图： 程序代码：（matlab） %说明：追赶法是适用于三对角矩阵的线性方程组求解的方法，并不适用于其他类型矩阵。 %定义三对角矩阵A的各组成单元。方程为Ax=d function x=zhuiganfa a=input(输入左下角对角元素a（n）(例如[1 2 3]):); b=input(输入对角元素b（n）(例如[1 2 3 4])：); c=input(输入左上角对角元素c（n）(例如[1 2 3])：); d=input(输入列矩阵d（n）例如[1 2 3 4]：); n=length(b); u0=0;y0=0; %“追”的过程 L(1)=b(1); y(1)=(d(1))/L(1); u(1)=c(1)/L(1); for i=2:(n-1) L(i)=b(i)-a(i-1)*u(i-1); y(i)=(d(i)-y(i-1)*a(i-1))/L(i); u(i)=c(i)/L(i); end L(n)=b(n)-a(n-1)*u(n-1); y(n)=(d(n)-y(n-1)*a(n-1))/L(n); %“赶”的过程 x(n)=y(n); for i=(n-1):-1:1 x(i)=y(i)-u(i)*x(i+1); end 第三章：解线性方程组的迭代法 1、题目简介：试分别用（1）Jacobi迭代法；（2）Gauss-Seidel迭代法；（3）共轭梯度法解线性方程组 迭代初始向量取=0，0，0，