《计算机实习报告.doc》由会员分享，可在线阅读，更多相关《计算机实习报告.doc（8页珍藏版）》请在三一办公上搜索。
1、 计 算 机 实 习 报 告 姓名：桂天才班级：01010904学号: 20093015192011年12月20日（说明：完成任务一的四道编程题，全部使用matlab编程来求解比较简单，matlab更接近于数学语言。）1. 第一题：用牛顿法求解附近的实根1.1. 计算方法：牛顿迭代法又称为切线法，试题分析；牛顿迭代法其迭代关系为 根据这一迭代关系式，利用循环就可以得出结果。其中。1.2. 实验源程序：.M文件function x,k,f1=funNewton(f)format long k=1f1=0x(1)=2x(2)=x(1)-(x(1)3-3*x(1)-1)/(3*x(1)2-3)f1=
2、x(2)-x(1)while abs(f1)=0.0005 k=k+1 x(k+1)=x(k)-(x(k)3-3*x(k)-1)/(3*x(k)2-3) f1=x(k+1)-x(k)endendk = 3x = 2.000000000000000 1.888888888888889 1.879451566951567 1.879385244836671f1 = -6.632211489598916e-0051.3. 结果分析：由其结果可以得到 x=1.879385244836671保留4位有效数字可得：x=1.879.2. 用三点高斯积分公式计算下列积分（1） （2）2.1 试题分析；高斯积分
3、为 ；(1).(2).通过建立.m文件以及matlab编程即可得到结果。2.2 实验源程序：.M文件为function fj=gausjf(f,a,b) % feval(f,x)为f 关于x函数。syms a b fx1=(a+b)/2-sqrt(3/5)*(b-a)/2;x2=(a+b)/2+sqrt(3/5)*(b-a)/2;fj=(b-a)*(5*feval(f,x1)/9+8*feval(f,(b+a)/2)/9+5*feval(f,x2)/9)/2;f1=sqrt(x+1.5);f2=1/x;fj1=gausjf(f1,-1,1) fj2=gausjf(f2,1,3) end fj=
4、(b-a)/2)*(5/9)*f(b+a)/2-sqrt(3/5)*(b-a)/2)+(8/9)*f(b+a)/2)+(5/9)*f(b+a)/2+sqrt(3/5)*(b-a)/2)fj1 = 1.9220fj2 = 1.09862.3 结果分析：由matlab运算所得结果（1）.积分结果为： 1.9220（2）.积分结果为： 1.09863. 用高斯消去法求解下列方程。3.1. 计算方法：高斯消去法其基本思想为，给线性方程组通过消元手续加工成单位上三角，在进行直下往上的回代运算。3.2. 实验源程序：.M文件为：clear;function RA,RB,n,X=gaus(A,b)B=A b
5、; n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA;if zhica0,disp(RA=RB.)returnendif RA=RB if RA=ndisp(RA=RB=n.) X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-
6、sum(A(q,q+1:n)*X(q+1:n)/A(q,q); endelse disp(RA=RB run(G:360data重要数据我的文档MATLABLanguage.m)f =0.121414.3 结果分析：故其在1.1300节点的插值函数值为：f =0.121415. 计算机实习总结：（1）.通过这次实习，我更加熟悉了matlab 这个软件的使用，同时将编程与计算方法结合起来，我的实践能力得到了提高，（2）.再者，通过这次实习，我更加熟悉了word排版。（3）。通过这次实习我对算法设计也有了进一步的了解，将数学与编程联系起来，增加了学习的趣味性，也增强了我对于数学的理解和matlab的了解。总之，这次实习我还是学到一些知识。
