山东大学数值计算实验报告1.pdfVIP

山东大学 计算机科学与技术 学院 数值计算 课程实验报告 学号： 姓名： 鸿武 班级：计机 14.1 实验题目： 1、 P44 ，Compter Problem 1.1 2 、 P45，Compter Problem 1.4 3、 高斯消元法 Computer Problems P100 2.2 ，（a ）（b ） (1) 写出直接LU 分解函数，写出直接LU 分解结果 (2) 结合LU 分解，前代、回代解方程组。写前代、回代函数 (3) 直接使用软件环境提供的函数LU 分解函数、解方程组方法比较 说明：(1)、(2)，这几部分不得直接使用软件环境提供的函数完成，可以参考课本上的算 法流程实现 实验学时：4 实验日期：2016 年 10 月 15 日 实验目的： 熟悉MATLAB 编程，消元法，LU 分解 硬件环境： PC 机 软件环境： MATLAB、JAVA、C++ 实验步骤与内容： 1、Write a program to compute the absolute and relative errors in Stirling’s approximation ! = √2 (/) for n = 1,...,10. Doestheabsoluteerrorgrow or shrink as n increases? Does the relative error grow or shrink as n increases? 根据题意编写代码如下： errorj.m function [ y ] = errorj( x ) y = sqrt(2*pi*x)*(x/exp(1))^x-factorial(x); end errorx.m function [ y ] = errorx( x ) y = (sqrt(2*pi*x)*(x/exp(1))^x-factorial(x))/factorial(x); end Stirling.m n = 15; x = 1:1:n; y1 = zeros(1,n); y2 = zeros(1,n); for i = 1:n y1(i) = abs(errorj(x(i))); y2(i) = abs(errorx(x(i))); end % scatter(x,y1) % hold on scatter(x,y2,p) hold on 当n 取 1 到 15 时，绘制出绝对误差的图像如下所示： 当n 取 1 到 100 时，绘制出绝对误差的图像如下所示： 当n 取 1 到 15 时，绘制出相对误差的图像如下所示： 当n 取 1 到 100 时，绘制出相对误差的图像如下所示： 观察以上四个图我们可以发现，随着n 的增加绝对误差不但增大，并且增大地越来越快。而 相对误差越来越小，并接近于 0 。 2.Write a program to compute the mathematical constant e, the base of natural logarithms, from the definition = (1 + 1) . Specifically, compute (1 + →∞ 1/n)n for n =10 , k = 1,2,...,20. If the programming language youusedoesnothaveanoperatorforexponentiation, you may use the equivalent n formula (1+1/n) = exp(nlog(1+1/n)), where exp and log are built-in functions. Determine the error in your successive approximations by comparing them with t