1、计算方法大作业共21页实验(shyn)一 牛顿(ni dn)下山法一、实验(shyn)目的:1、掌握牛顿下山法求解(qi ji)方程根的推导原理。2、理解牛顿下山法的具体算法与相应程序的编写。二、实验内容:采用牛顿下山法求方程2x3-5x-17=0在2附近的一个根。三、实验实现:1、算法:下山因子从开始,逐次将减半进行试算,直到能使下降条件成立为止。再将得到的循环求得方程根近似值。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)
2、/df(p0); while (abs(f(p1) = abs(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、运行(ynxng)结果:四、实验(shyn)体会:牛顿下山法可以较快求的方程(fngchng)结果,对于该题,只需要5步。运用计算机的数值迭代法可以很快求得满足精度要求的结果。实验(shyn)二 矩阵的列主元三角分解一、实验目的:学会矩阵的三角分解,并且能够用MATLAB编写相关程序,实现矩阵的三角分解,解方程组。二、实验内
3、容:用列主元消去法求解方程组(实现PA=LU) 要求输出: (1)计算解X; (2)L,U; (3)正整型数组IP(i),(i=1,n) (记录(jl)主行信息)。三、实验(shyn)实现:1、算法(sun f): 列主元三角分解和普通Dooliitle分解不同(b tn),第k步分解时为了避免用绝对值很小的数作除数,引进量 若,则将矩阵的第t行与第k行元素互换,再进行正常的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、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);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(
5、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(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=);di
6、sp(L);disp(U=);disp(U);disp(IP=);disp(IP);y=Lb;x=Uy;disp(y=);disp(y);disp(x=);disp(x);4、运行(ynxng)结果:四、实验(shyn)体会:列主元三角分解(fnji)原理很简单,用MATLAB实现起来却不容易。对MATLAB还需要深入学习。实验(shyn)三 SOR迭代法一、实验(shyn)内容:采用SOR方法求解方程采用0向量为初始向量,迭代以结束,但迭代次数不操作1000次,松弛因子使用 其中,给出使用的各松弛因子及对应迭代次数。二、实验实现:1、程序代码如下:A= 10 -1 0; -1 10 -2;
7、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(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,
8、4); end end S= R(J,4)运行(ynxng)结果为:实验(shyn)七 数值积分实验(shyn)内容利用(lyng)改进Euler法,经典四级四阶Rung-Kutta方法求其数值解,: 分别取h=0.2,0.1,0.05时数值解。 (注:初值问题的精确解) 1、改进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+
9、y2)/2; x=x+h; Y(j,i)=y; end h Y(j,1:n) disp(_);(2)运行结果(ji gu)如下:2、经典(jngdin)R-K方法(1)程序代码如下(rxi):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/2); K3 = h*f(x+h/2,y+K2/2); K4 = h*f(x+h,y+K3); y = y+(K1+2*K2+
10、2*K3+K4)/6; x = x+h; Y(j,i)=y;end h Y(j,1:n) disp(_);end(2)运行结果(ji gu)如下:内容摘要(1)实验一 牛顿下山法实验目的:掌握牛顿下山法求解方程根的推导原理(2)IP(p,:)=C2(1:n)(3)I=ones(size(Q)(4)x=AlnQa=exp(x(1)alpha=x(2)beta=x(3)(2)运行结果如下:实验六 数值积分实验目的掌握复化梯形公式和复化Simpson公式,会编写程序,上机进行实例计算(5)程序代码如下:format longf=inline(log(1+x)./x)(6)sum=sum+f(x(i)(7)程序代码如下:f=inline(log(1+x)./x)(8)x=x+h(9)Y(j,i)=y
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1