MATLAB实习报告-1

MATLAB --实习报告班级：******姓名：**学号：***********指导老师：***日期：2011/11/24目录一．MATLAB简介 (3)二．上机实习题目 (3)1.第2章第1题 (3)2.第3章第2题 (3)3.第5章第2题 (3)4.第4章第1题 (4)三．题目实现过程 (4)1.牛顿插值以及三次样条插值（第一个实验题） (4)a.牛顿插值 (4)b.三次样条插值 (4)c.牛顿插值截图 (5)d.三次样条插值截图 (6)2.3次，4次多项式的曲线拟合 (7)a.先输入表格中对应的数据， (7)b.3次拟合 (8)c. 3次拟合截图 (9)d. 4次拟合 (9)e. 4次拟合截图 (10)3.高斯消去法解线性方程组 (11)a.高斯消去法源代码 (11)b.第一个小题 (12)c.第二个小题 (13)4.雅克比迭代法与SOR方法 (13)a. 雅克比迭代法的计算公式 (13)b. 超松弛迭代法的计算公式 (14)c.得到希尔矩阵 (14)d.得到b矩阵 (14)e.雅克比迭代法实现的函数 (14)f.SOR迭代法实现的函数 (15)g.对于雅克比迭代法，通过执行以下代码 (16)h. 对于SOR迭代法，执行相对应代码 (16)四．心得与感想 (19)一．M ATLAB简介MATLAB是一个软件，用来进行科学计算。

在实际生活或者在科学研究中，大量的计算有时候是不可避免的，在这个时候，我们可以选择很多种方式来解决我们的问题。

但是选择一个好的软件对于我们来说有时候却是非常困难的，尤其是现在各种软件层出不穷，当然里面有好的，可是有大多数的软件却是鱼目混珠，在里面滥竽充数。

如果我们不避开这些，那么我们可能面临的是用很大的工作量来完成一个极其简单的问题，或者说在付出了时间与精力之后却完不成问题。

所以必须选择一个好的计算软件来解决问题。

MATLAB 就是这样一个软件，在你能够熟悉并且流畅运用它之后，那它将是一把能够解决很多问题的钥匙。

MATLAB是Matrix Laboratory，即矩阵实验室的缩写。

它是一个有Math Work公司（由Moler，Little，Bangert在1984年在加利福利亚成立）开发的软件包。

用来实现数值计算或者符号操作。

另外MATLAB还拥有M-file这一功能，在这里面你可以编写自己的程序（当然这个程序不是windows程序），在这里你就可以实现很多计算，虽然说MATLAB自带的函数在大多数情况下已经够用，但是总会有某些时候你觉得自己编写的程序运用起来更加灵活，而且在一些具体问题的时候，这个优势将会是一大突破，可以让我们很轻松解决很多问题。

除了自己编写外，MATLAB 自带的很多优秀的程序，如plot（绘图）函数，cond（求矩阵条件数）函数，关于矩阵的一些相关函数，都是相当经典的。

总之，熟悉使用MATLAB可以让我们的工作变得很简单，经过这次的学习，我们以后的很多学习工作问题都可以很方便地解决了。

二．上机实习题目1.第2章第1题4进行插值。

用图给出{(x i,y i),x i=0.2+0.08i,i=0,1,11,10}，P4(x)及S(x)。

2.第3章第2题拟合曲线，用图示数据曲线及相应的三种拟合曲线。

3.第5章第2题用列主元高斯消去法解线性方程组Ax=b。

(1)[3.01 6.03 1.991.27 4.16−1.230.987−4.819.34][x1x2x3]=[111];(2) [3.00 6.03 1.991.27 4.16−1.230.990−4.819.34][x1x2x3]=[111].分别输出A，b，det (A),解向量x，(1)中A的条件数。

分析比较(1),(2)。

4.第4章第1题给出线性方程组H n x=b，其中系数矩阵H n为希尔伯特矩阵：H n=(h ij)∈,i,j=1,2,…,n。

R n×n,h ij=1i+j−1假设x∗=(1,1,…,1)T∈R n×n,b=H n x∗。

若取n=6,8,10，分别用雅克比迭代法及SOR迭代法(ω=1,1.25,1.5)求解。

比较计算结果。

三．题目实现过程下面开始进行计算：1.牛顿插值以及三次样条插值（第一个实验题）此题要求利用给定点，及给定点的函数值进行牛顿插值以及三次样条插值。

a.牛顿插值要实现牛顿插值，要用到以下代码% 调用格式：yi=Lagran_(x,y,xi)% x,y 数组形式的数据表function fi=Lagran_(x,f,xi)fi=zeros(size(xi));np1=length(f);for i=1:np1z=ones(size(xi));for j=1:np1if i~=j,z=z.*(xi-x(j))/(x(i)-x(j));endendfi=fi+z*f(i);endb.三次样条插值要实现三次样条插值，要用到以下代码function S=csfit(X,Y,dx0,dxn)N=length(X)-1;H=diff(X);D=diff(Y)./H;A=H(2:N-1);B=2*(H(1:N-1)+H(2:N));C=H(2:N);U=6*diff(D);B(1)=B(1)-H(1)/2;U(1)=U(1)-3*(D(1));B(N-1)=B(N-1)-H(N)/2;U(N-1)=U(N-1)-3*(-D(N));for k=2:N-1temp=A(k-1)/B(k-1);B(k)=B(k)-temp*C(k-1);U(k)=U(k)-temp*U(k-1);endM(N)=U(N-1)/B(N-1);for k=N-2:-1:1M(k+1)=(U(k)-C(k)*M(k+2))/B(k);endM(1)=3*(D(1)-dx0)/H(1)-M(2)/2;M(N+1)=3*(dxn-D(N))/H(N)-M(N)/2;for k=0:N-1S(k+1,1)=(M(k+2)-M(k+1))/(6*H(k+1));S(k+1,2)=M(k+1)/2;S(k+1,3)=D(k+1)-H(k+1)*(2*M(k+1)+M(k+2))/6;S(k+1,4)=Y(k+1);endc.牛顿插值截图在输入对应的两个数组后，首先进行4次牛顿插值，下面是运行截图然后根据结果绘图d.三次样条插值截图首先得到的是三次样条插值后的多项式系数矩阵，然后利用polyval()函数计算在对应点的值，并且通过绘图函数将结果表示在图像上面，通过与上面的牛顿插值所得到的图像进行比较，可以发现，利用三次样条插值得到的结果更加的精确。

2.3次，4次多项式的曲线拟合a.先输入表格中对应的数据，然后首先通过调用plot()函数，绘图b.3次拟合通过以下代码>> xx =0 0.1000 0.2000 0.3000 0.5000 0.8000 1.0000>> yy =1.0000 0.4100 0.5000 0.6100 0.91002.02002.4600>> s1 = polyfit(x, y, 3);>> temp = 0:0.001:1;>> y_temp = zeros(size(temp));>> y_temp = polyval(s1, temp);>> hold on;>> plot(temp, y_temp)c. 3次拟合截图可以得到下面的截图，这是看到3次拟合后的多项式与原来的图形比较d. 4次拟合通过下面代码>> xx =0 0.1000 0.2000 0.3000 0.5000 0.80001.0000>> yy =1.0000 0.4100 0.5000 0.6100 0.91002.02002.4600>> s2 = polyfit(x, y, 4);>> y_temp = polyval(s2, temp);>> hold on;>> plot(temp, y_temp, 'r')e. 4次拟合截图通过进行4次拟合，我们会感觉不同，从图像上可以看出（不同次数的拟合在图像上用不同的颜色标示出来以便观察，其中3次拟合是蓝色，4次的是红色）以上都是用多项式进行拟合，下面代码用一个线性函数来进行拟合>> xmean = mean(x);>> ymean = mean(y);>> sumx2 = (x - xmean) * (x - xmean)';>> sumxy = (y - ymean) * (y - ymean)';>> A = sumxy / sumx2;>> B = ymean - A * xmean;>> temp = 0:0.01:1;>> y_temp = A * temp + B;>> plot(temp, y_temp, 'green')的到下图其中绿色的线是这次拟合的结果，很明显，误差很大在进行解决这个题目时候用的拟合函数是MATLAB自带的拟合函数polyfit()，这里是多项式拟合。

polyfit()函数得到的是一个矩阵，然后运用另一个MATLAB自带的函数polyval()来进行计算给定点的函数值。

到这里这一题完成。

3.高斯消去法解线性方程组对于这道题目，要求用高斯消去法解两个线性方程组，所以利用MATLAB 的M-file编写一个文件便可以解决两个小题。

而高斯消去法的核心是逐渐利用消去法是原来的系数矩阵变成一个上三角矩阵。

然后一个迭代便可以解出方程的每一个解。

a.高斯消去法源代码function X = p178_2_GAUSS(A, b)[n, m] = size(A);X = zeros(n, 1);temp = zeros(1, m);temp_b = 0;i = 1;for j = 1: (m - 1)if (A(i, j) ~= 0)for k = (i + 1):nif (A(k, j) ~= 0)temp = A(k, :) + A(i, :) * (-A(k, j) / A(i, j));temp_b = b(k) + b(i) * (-A(k, j) / A(i, j));A(k, :) = temp;b(k) = temp_b;endendendi = i + 1;end;Abdisp('det(A) is ...');x = det(A);disp(x);disp('cond(A) is ...');x = cond(A);disp(x);X(n) = b(n) / A(n, n);for i = (n - 1):-1:1temp_b = 0;for j = (i + 1):ntemp_b = temp_b + A(i, j) * X(j);endX(i) = (b(i) - temp_b) / A(i, i);endend这个程序要求输入两个参数，一个是系数矩阵，另外一个是矩阵b。