《东南大学《数值分析》上机题.docx》由会员分享，可在线阅读，更多相关《东南大学《数值分析》上机题.docx（10页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、数值分析上机题1,其精确值为始十七（1）编制按从大到小的顺序与=匕+舌 + +，计算S,V的通用程序。（2）编制按从小到大的顺序Sn = J；），J +白，计算Sn的通用程序。（3）按两种顺序分别计算S”，S, S1（）6,并指出有效位数。（编制程序时用单精度）（4）通过本上机题，你明白了什么？程序代码(matlab编程)： clccleara=single(1./(2:10A7 .A2-l);SI (1)=single(0);SI (2)=1/(2A2-1);for N=3:10A2SI (N)=a(l);for i=2:N-lSI (N)=S1(N)+a(i);end endS2(1)=s
2、ingle (0);S2(2)=1/(2A2-1);for N=3:10A2S2(N)=a(N-l);for i=linspace (N-2,1rN-2)S2 (N) =S2 (N) +a (i); end endSI表示按从大到小的顺序的SnS2表示按从小到大的顺序的Sn 计算结果从大到小的顺序的值从小到大的顺序的值精确值有效位数从大到小从小到大0.7400490.740050.74004965%0.7498520.74990.749944%0.7498520.7499990.74999936通过本上机题，看出按两种不同的顺序计算的结果是不相同的，按从大到小的顺序计算的值与精确习题623.(
3、上机题)常微分方程初值问题数值解(1)编制R&方法的通用程序；(2)编制AB4方法的通用程序(由RK4提供初值)；(3)编制AB4-AM4预测校正方法的通用程序(由R&提供初值)；(4)编制带改进的AB4-AM4预测校正方法的通用程序(由R&提供初值)；(5)对于初值问题=(0x 0,当/ (一夕5)时，Newton迭代序列收敛于根 石。试确定尽可能大的3。2 .试取若干初始值，观察当公 (8,7)，(T2), (-5,6)，(氏1)，(l,oo)时Newton序列是 否收敛以及收敛于哪一个根。MATLAB程序问题1clccleardx=0 . 5;x (1)=0.5;while(dxle-6
4、)i=l；error=l;while (errorle-8)x(i+1)=x(i)-(l/3*x(i)A3-x (i)/ (x (i)A2-l);error=abs(x(i + 1)-x (i);i=i+l;endif(x(i)=0)x(1)=x(1)+dx;elsedx=dx/2;x (1)=x (1)-dx;endend经计算，最大的g为0.774596问题2clcclearx2(1)=lel4;i=l；error=l;while (errorle-8)x2(i + l)=x2(i)-(l/3*x2(i)3-x2 (i)/ (x2 (i)A2-l);error=abs (x2 (i + 1
5、)-x2 (i);i=i+l;if (ile4)breakendend对于不同得初始值收敛于不同的根，/在（-8,j）内收敛于-在（-0.774, 0.774）内收敛于0,在（1, +8）内收敛于百，但在内（0.774, 1）和（-1, 0.774）均可能收敛于6和。分析：对于不同的初值，迭代序列会收敛于不同的根，所以在某个区间内求根对于初值的选取有很大的 关系。产生上述结果的原因是区间不满足大范围收敛的条件。数值分析上机题339.(上机题)列主元三角分解法对于某电路的分析，归结为求解线性方程组RI=V。 (1)编制解n阶线性方程组Ax=b的列主元三角分解法的通用程序；(2)用所编制的程序解线
6、性方程组RI=V,并打印出解向量，保留五位有效数；(3)本编程之中，你提高了哪些编程能力？程序：clcclearA= 31, -13,0,0,0, -10,0,0,0-13,35,-9,0,0,0,0,00,-9,31,-10,0,0,0,0,00,0,-10,79,-30,0,0,0,-90,0,0,-30,57,-7,0,-5,00,0,0,0,-7,47,-30,0,00,0,0,0,0,-30,41,0,00f0,0r0r-5,0r 0,27,-20r 0r0,-9,0r 0,0,-2,29;b=-工5,27,-23,0,-20,12,-7,7,10，m,n=size (A);Ap=A
7、,b;x=zeros (n,1); for i=l:m-l, maxaAmaxi=max(abs(Ap(i:end,j);maxi=maxi+i-l;if(maxa=0)mid=Ap(maxi,:);Ap(maxi, :)=Ap (i,:);Ap(i,:)=mid; for k=i:mAp(i + 1:m, :)=Ap(i + 1:m, :)-Ap(i + 1:mr j)*(Ap(i, :) ./maxa);endendendfor i=linspace(m,1,m)x(i)=(Ap(i,end)-Ap(i,1:end-1)*x)/Ap(i,i); end结果：方程的解为(保留5位有效数字)：
8、xl=-0.28923, x2= 0.34544, x3= -0.71281, x4= -0.22061, x5= -0.43040, x6= 0.15431, x7= -0.057823, x8= 0.20105, x9= 0.29023。习题4巧0123456789102.513. 304. 044. 705. 225. 545. 785.405. 575. 705. 80编制求第一型3次样条插值函数的通用程序; 已知汽车曲线型值点的数据如下：37.(上机题)3次样条插值函数端点条件为)，o=O.8,%。二0.2。用所编制程序求车门的3次样条插值函数S(x),并打印出S(i+O. 5) (
9、i=0, 1, -9)o 程序：(1) clc clear%x=0,1,2,3,4,5,6,7,8,9,10;y=2.51,3.30,4.04,4.7,5.22z 5.54,5.78,5.40,5.57,5.70,5.80;yl=0.8;yend=0.2;% n=size (x,2)-1;h=x(2:end)-x(1:end-1);miu=h(1:end-1) ./ (h (1:end-1)+h(2:end);lamda=l-miu;f1=yl,(y(2:end)-y(1:end-1)./h,yend;%fxn-1,xnf2=fl(2:end)-fl(1:end-1)./h(l),h(l:en
10、d-1)+h(2:end) , h(end);%fxn-1, xn, xn+1 A=2 , *eye(n+1);A(2:end,1:end-1)=A(2:end,1:end-1)+diag(miu,1 1);A(1:end-1z 2:end)=A(1:end-1,2:end)+diag(1,lamda);M=A(6*f21 );Sx=y(1:end-1) 1, (y(2:end)-y(1:end-1) ./h) , -(1/3*M(1:end-1)+l/6*M(2:end).*hf),1/2*M(1:end-l),1/6*(M(2:end)-M(l:end-l)./h1；% xx=input(
11、 x= );for j =2:n+lif xxx (j)S=Sx (j-1, :)* 1,xx-x (j-1)z (xx-x (j-1)A2, (xx-x (j-1)3 1;break endend clcclear%x=0,1,2,3,4,5,6,7,8,9,10;y=2.51,3.30,4.04,4.7,5.22z 5.54,5.78,5.40,5.57,5.7 0,5.80;yl = 0 . 8;yend=0.2;% n=size (xz 2)-1;h=x(2:end)-x(1:end-1);miu=h(1:end-1) ./ (h (1:end-1)+h(2:end);lamda=l-
12、miu;f1=yl,(y(2:end)-y(1:end-1)./h,yend;%fxn-1,xnf2=fl(2:end)-fl(1:end-1)./h(l),h(l:end-1)+h(2:end),h(end);%fxn-1, xn, xn+1 A=2.*eye(n+1);A(2:end,1:end-1)=A(2:end,1:end-1)+diag(miu,1 ,);A (1:end-1z 2:end)=A(1:end-1z 2:end)+diag(1,lamda);M=A(6*f21);Sx=y(1:end-1) 1z (y(2:end)-y(1:end-1) ./h) 1 -(1/3*M(
13、1:end-1)+l/6*M(2:end).)，1/2*M(1:end-l),1/6*(M(2:end)-M(l:end-1)./h1；%for i=0:9xx=i + 0 . 5;for j =2:n + 1if xx2.90856S(l.5)=3.67843S (2.5)=4.38147S(3.5)=4.98819S(4.5)=5.38328S(5.5)=5.7237S(6.5)=5.59441S(7.5)=5.42989S(8.5)=5.65976S(9.5)=5.7323习题五重积分的计算23 (上机题)重积分的计算题目：给定积分/(/) =1(，/(x,y)dx)dy。取初始步长h和
14、k,及精度。应用复化梯形公式，采用逐 次二分步长的方法，编制计算I的通用程序。计算至相邻两次近似值之差的绝对值不超过为止。1)用所编程序计算积分/(7)=(火，+p2)公)力，取 = 0.5*10-5。程序：clcclear% examplef=inline(!tan(x.人2+y.人2) ，, ，x，, ，y。;a=0;b=pi/3;c=0;d=pi/6;% define error=l; k=l;n=l;while(error0.5e-5)x,y=meshgrid(linspace(c,d,2Ak+1),linspace(a,b, 2Ak+1);h=(b-a)/2Ak;l=(d-c)/2A
15、k;z=f(x,y);zl = z (1:end-1,1:end-1);z2 = z (1:end-1,2:end);z3=z (2:end,1:end-1);z4 = z (2:end,2:end);t(k)=h*l/4*(sum(sum(zl)+sum(sum(z2)+sum(sum(z3)+sum(sum(z4);% extrapolationif(k=2)T (1, k-D=4/3*t (k) -l/3*t (k-1) ;% T (1) error=min(error,abs(t(k)-t (k-1);if(k=3)T(2, k-2)=16/15*T(1,k-1)-1/15*T(1,k
16、-2); % T (2) error=min(rror,abs(T(1,k-1)-T(1,k-2);if(k=4)T(3, k-3) =64/63*T(1,k-2)-1/63*T (1, k-3) ; % T (3) error=min(error,abs(T(2,k-2)-T (2,k-3); if (k=5)error=min(error,abs(T(3,k-3)-T (3,k-4); end endendendk=k+l;end计算结果：12345T(f) 0.5197965 0.3879734 0.3503495 0.3400887 0.3374218T(1)(f) 0.344032 0.337808 0.336668 0.336533 0.336521T (2)(f) 0.337393 0.336592 0.336524 0.336521T0.3377090.336650.33653160.3367464I(f)=0.33652二分6次