計(jì)算方法大作業(yè)

1、實(shí)驗(yàn)一 牛頓下山法一、 實(shí)驗(yàn)?zāi)康模?、 掌握牛頓下山法求解方程根的推導(dǎo)原理。2、 理解牛頓下山法的具體算法與相應(yīng)程序的編寫。二、 實(shí)驗(yàn)內(nèi)容：采用牛頓下山法求方程2x3-5x-17=0在2附近的一個(gè)根。三、 實(shí)驗(yàn)實(shí)現(xiàn)：1、 算法：下山因子從開始，逐次將減半進(jìn)行試算，直到能使下降條件成立為止。再將得到的循環(huán)求得方程根近似值。2、 程序代碼如下：function p,k=NewtonDownHill(f,df,p0)N=2000;Tol=10(-5);e=10(-8);for k=1:N lamda=1; p1=p0-lamda*f(p0)/df(p0); while (abs(f(p1) = ab

2、s(f(p0) & lamdae) lamda=lamda/2; p1=p0-lamda*f(p0)/df(p0); end if abs(p1-p0)Tol break end p0=p1;endans=p13、 運(yùn)行結(jié)果：四、 實(shí)驗(yàn)體會(huì)：牛頓下山法可以較快求的方程結(jié)果，對于該題，只需要5步。運(yùn)用計(jì)算機(jī)的數(shù)值迭代法可以很快求得滿足精度要求的結(jié)果。實(shí)驗(yàn)二 矩陣的列主元三角分解1、 實(shí)驗(yàn)?zāi)康模簩W(xué)會(huì)矩陣的三角分解，并且能夠用MATLAB編寫相關(guān)程序，實(shí)現(xiàn)矩陣的三角分解，解方程組。2、 實(shí)驗(yàn)內(nèi)容：用列主元消去法求解方程組（實(shí)現(xiàn)PA=LU) 要求輸出： （1）計(jì)算解X; （2）L,U; （3）正整型數(shù)

3、組IP(i),(i=1,n) (記錄主行信息)。3、 實(shí)驗(yàn)實(shí)現(xiàn)：1、算法： 列主元三角分解和普通Dooliitle分解不同，第k步分解時(shí)為了避免用絕對值很小的數(shù)作除數(shù)，引進(jìn)量 若，則將矩陣的第t行與第k行元素互換，再進(jìn)行正常的Doolittle分解。2、程序代碼如下：clear all;clc;A=1 1 1 1 1 1 1;2 1 1 1 1 1 1;3 2 1 1 1 1 1;4 3 2 1 1 1 1;5 4 3 2 1 1 1;6 5 4 3 2 1 1;7 6 5 4 3 2 1;b=7;8;10;13;17;22;28;n=length(A);IP=eye(n);U=zeros(n

4、);L=eye(n);m,p=max(A(:,1);C1=b(1);b(1)=b(p);b(p)=C1;C2(1:n)=IP(1,:);IP(1,:)=IP(p,:);IP(p,:)=C2(1:n);C2(1:n)=A(1,:);A(1,:)=A(p,:);A(p,:)=C2(1:n);U(1,:)=A(1,:);for k=2:n L(k:n,k-1)=A(k:n,k-1)/U(k-1,k-1); j=1; for i=k:n S(j)=A(i,k)-L(i,1:k-1)*U(1:k-1,k); j=j+1; end m,p=max(abs(S(:) ; C2(1:n)=A(k,1:n);A

5、(k,1:n)=A(k+p-1,1:n);A(k+p-1,1:7)=C2(1:n); C2(1:n)=IP(k,1:n);IP(k,1:n)=IP(k+p-1,1:n);IP(k+p-1,1:n)=C2(1:n); U(k,k:n)=A(k,k:n)-L(k,1:k-1)*U(1:k-1,k:n); if kn A(k+1:n,k)=A(k+1:n,k)-L(k+1:n,1:k-1)*U(1:k-1,k); end enddisp(L=);disp(L);disp(U=);disp(U);disp(IP=);disp(IP);y=Lb;x=Uy;disp(y=);disp(y);disp(x=

6、);disp(x);4、 運(yùn)行結(jié)果：4、 實(shí)驗(yàn)體會(huì)：列主元三角分解原理很簡單，用MATLAB實(shí)現(xiàn)起來卻不容易。對MATLAB還需要深入學(xué)習(xí)。實(shí)驗(yàn)三 SOR迭代法1、 實(shí)驗(yàn)內(nèi)容：采用SOR方法求解方程采用0向量為初始向量，迭代以結(jié)束，但迭代次數(shù)不操作1000次，松弛因子使用 其中，給出使用的各松弛因子及對應(yīng)迭代次數(shù)。2、 實(shí)驗(yàn)實(shí)現(xiàn)：1、程序代碼如下：A= 10 -1 0; -1 10 -2; 0 -2 10;b=9;7;6;x0=0;0;0;e=10(-6);n=1000;x=x0;j=18;w=0.1*j; for k=1:n for i=1:length(b) x(i)=w*(b(i)-A(

7、i,:)*x)/A(i,i)+x(i); if norm(x-x0)eps) J=J+1; h=h/2; x=a+h:2*h:b-h; R(J+1,1)=R(J,1)/2+h*sum(f(x); for K=1:min(3,J) R(J+1,K+1)= R(J+1,K)+( R(J+1,K)- R(J,K)/(4K-1); end if(J3) err=abs( R(J+1,4)- R(J,4); end end S= R(J,4)運(yùn)行結(jié)果為：實(shí)驗(yàn)七 數(shù)值積分實(shí)驗(yàn)內(nèi)容利用改進(jìn)Euler法，經(jīng)典四級四階Rung-Kutta方法求其數(shù)值解，： 分別取h=0.2,0.1,0.05時(shí)數(shù)值解。 （注：初

8、值問題的精確解） 1、 改進(jìn)Euler法（1） 程序代碼如下：x0=0;y0=1;xn=2;f=inline(4*x./y-x.*y);H=0.2 0.1 0.05;for j=1:3 h=H(j); n=(xn-x0)/h; x=x0;y=y0; for i=1:n y1=y+h*f(x,y); y2=y+h*f(x+h,y1); y=(y1+y2)/2; x=x+h; Y(j,i)=y; end h Y(j,1:n) disp(_);（2） 運(yùn)行結(jié)果如下：2、 經(jīng)典R-K方法（1） 程序代碼如下：x0=0;y0=1;xn=2;f=inline(4*x./y-x.*y);H=0.2 0.1 0.05;for j=1:3 h=H(j); n=(xn-x0)/h; x=x0;y=y0; for i = 1:n K1 = h*f(x,y); K2 = h*f(x+h/2,y+K1