韩强数值分析实报告.docVIP

数值分析 实验报告 系 别 电子信息系 专 业 计算机科学与技术 班级学号 4090418 姓 名 韩强 指导教师 万 聪 2011年 6月20日 目录 实验一 11 一．实验目的： 11 1.熟使悉用matlab 11 2.编写程序，完成实题目验 11 二．实验题目： 11 三．实验步骤 11 （二）.编写拉格朗日插值函数的程序 14 实验二 19 一．实验题目 19 二实验内容 20 实验三 38 一．实验目的： 38 二．实验题目： 38 三．实验步骤： 39 实验四 43 一．实验目的 43 二．实验题目 43 三． 实验步骤 43 实验一 一．实验目的： 1.熟使悉用matlab 2.编写程序，完成实题目验 二．实验题目： 编写一个拉格朗日插值函数，对不多于9个点的插值节点都可以求出插值函数，任意给定输入x值都可以求出y值。 例如：（0,1）,(1,1) ,(4,2) ,(9,3) ,(16,4) ,(25,5) ,(36,6) ,(49,7) ,(64,8)，这实际上是一个平方根函数随意输入x在[0，64]上的值都可以输出一个y值。 进一步可以在这区间上画出图像。 三．实验步骤 例如： Y=e-t/2.5sin3t t=0:pi/50:5*pi; y=exp(-t/2.5).*sin(3*t); plot(t,y,-b) y （一）1.程序如下： 输入一个函数，该函数输入一个数字N，计算结果为到N的斐波那契数列 function [ f ] = fib( n ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here f=[1 1]; if n==1n0 f=[1]; else if n==2 f=[1 1]; else for i=3:n f(i)=f(i-2)+f(i-1); end end end g=f(n); end 2.输出函数如下： 3.输出图形如下： （二）.编写拉格朗日插值函数的程序 1. 函数： f = 1/(x^2+1) 将其保存在f.m 文件中，具体程序如下： function y = f1(x) y = 1./(x.^2+1); 2. 拉格朗日插值函数：将其保存在lagrange.m 文件中，程序代码如下如下： function y = lagrange(x0,y0,x) m = length(x); /区间长度/ n = length(x0); for i = 1:n l(i) = 1; end for i = 1:m for j = 1:n for k = 1:n if j == k continue; end l(j) = ( x(i) -x0(k))/( x0(j) - x0(k) )*l(j); end end end y = 0; for i = 1:n y = y0(i) * l(i) + y; end 3. 建立测试程序，保存在text.m文件中，实现画图： x=-5:0.001:5; y=(1+x.^2).^-1; p=polyfit(x,y,n); py=vpa(poly2sym(p),10) plot_x=-5:0.001:5; f1=polyval(p,plot_x); figure plot(x,y,r,plot_x,f1) 输入n=6,出现下面的图形： 通过上图可以看到当n=6是没有很好的模拟。 于是重新运行text.M并选择n=11 由此可见n=11时的图像是可以很好的实现模拟。 实验二 一．实验题目 用不同的方法计算积分 取不同的步长h，分别用复合梯形公式及复合辛普森公式计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h，使得精度不能再被改善。 二实验内容 （一）编写一