1、精選優(yōu)質(zhì)文檔-傾情為你奉上計算物理結(jié)課報告姓 名：_學(xué) 號：_班 級：_指導(dǎo)老師：_所用工具：Matlab8.0。3(A).用龍貝格法編程計算積分： I=0.10.60.02792（2-x）1.449x+10.8 1-x1.2 xdx要求誤差不超過10-5。程序代碼：function y=question3(a,b,f,w)a1=b-a;f1=subs(f,a);f2=subs(f,b);Y=(f1+f2)*a1*0.5;T(1,1)=Y;Y1=f1+f2;k=2;w0=1;while abs(w0)>w h=a1./(2(k-1);u=0; for j=1:2(k-1)+1; z=a+

2、h*(j-1); f4=subs(f,z); u=u+f4; end T(k,1)=(2*u-Y1)*h*0.5; w0=T(k,1)-T(k-1); k=k+1;endfor i=2:k-1 for j=i:k-1 T(j,i)=(4(i-1)*T(j,i-1)-T(j-1,i-1)./(4(i-1)-1); endendy=T(k-1,k-1);y=vpa(y,7);運行結(jié)果：5.(B)試用二級歐拉近似法計算取步長h=0.1。程序代碼：function Y0=question5(f,x0,y0,b,h)syms t y;n=floor(b-x0)./h)+1;x(1)=x0;Y(1)=y0

3、;f1=subs(f,'t',x(1); f1=subs(f1,'y',Y(1);Y(2)=Y(1)+h*f1;for i=3:n x(i-1)=x(i-2)+h; f2=subs(f,'t',x(i-1); f2=subs(f2,'y',Y(i-1); Y(i)=Y(i-2)+2*h*f2; endY0=Y(n);Y0=vpa(Y0,6);運行結(jié)果：>> syms t y;f=1+exp(-t)*sin(y);Y0=question5(f,0,0,1,0.1) Y0 = 1.284857. (B)設(shè)fx=x3-1.8

4、x2+0.15x+0.65=0 ，試用二分法求該方程在區(qū)間0.5,1.5內(nèi)根的近似值。 程序代碼：function x0=Dichotomy(a,b,z,w)s=z;while (b-a)>w z1=subs(z,a); c=(a+b)/2; z2=subs(s,c); if z1*z2>0 a=c; else b=c; end endx0=c;x0=vpa(x0,6);運行結(jié)果：8. (B)利用牛頓迭代法求解 e2x+x-4=0在x>0.5的根的近似值。 程序代碼：function x1=Newton(x0,f,w1) z=diff(sym(f);w2=0.1;i=0;wh

5、ile w2>w1 xa=subs(f,x0); xb=subs(z,x0); x1=x0-xa./xb; w2=abs(x1-x0); x0=x1; i=i+1; if i>=200 break; endendx1=vpa(x1,8);運行結(jié)果：10. (B)在某次實驗中，需要觀察水分的滲透速度，測得時間t與水的重量W的數(shù)據(jù)如下： t/s 1 2 4 8 16 32 64 W/g 4.22 4.02 3.85 4.59 3.44 3.02 2.59 設(shè)已知t和W的關(guān)系為W=Ats ，試用最小二乘法確定參數(shù)A，s.先需要對W=Ats求對數(shù)化為線性方程：lnW=lnA+slnt.程序

6、代碼：T=1 2 4 8 16 32 64;W=4.22 4.02 3.85 4.59 3.44 3.02 2.59;for i=1:7 T1(i)=log(T(i); W1(i)=log(W(i);endf,a,b=Niehe2(T1,W1);A=exp(a);s=b;plot(T,W,'*');hold ont=0:0.01:65;w=A*(t.s);plot(t,w);hold onlegend ('節(jié)點','圖像');function f,a,b=Nihe(x,y)m=length(x);z1=0;z2=0;z3=0;z4=0;z5=0;z6=0;for i=1:mz2=z2+x(i);z3=z3+y(i); z5=z5+x(i)*x(i);z6=z6+y(i)*x(i);endz1=z1+m;z4=z4+z2; a=(z3*z5-z2*