计算机实习报告.doc

计 算 机 实 习 报 告 姓名：桂天才班级：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=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 试题分析；高斯积分为 ；(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=(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; 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)-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）.再者，通过这次实习，我更加熟悉了wor