当前位置:文档之家› 雅可比迭代实验报告记录

雅可比迭代实验报告记录

雅可比迭代实验报告记录
雅可比迭代实验报告记录

雅可比迭代实验报告记录

————————————————————————————————作者:————————————————————————————————日期:

雅可比迭代法求解线性方程组的实验报告

一、实验题目

分别利用雅可比迭代法和高斯-塞德尔迭代法求解以下线性方程组:

使得误差不超过 0.00001。

二、实验引言

1.实验目的

①掌握用迭代法求解线性方程组的基本思想和步骤,熟悉计算机fortran 语言;

②了解雅可比迭代法在求解方程组过程中的优缺点。

2.实验意义

雅克比迭代法就是众多迭代法中比较早且较简单的一种,求解方便实用。

三、算法设计

1.雅可比迭代法原理:

设有线性方程组Ax=b 满足0≠ii a , 将方程组变形为: x=Bx+f, 则雅可比(Jacobi)迭代法是指f Bx X k k +=+)1(,即 由初始解逐步迭代即可得到方程组的解。 算法步骤如下:

步骤1.给定初始值)0()0(2)0(1,,,n x x x ?,精度e,最大容许迭代次数M ,令k=1。

步骤2.对i=1,2,…,n 依次计算

)0()1()0()1(11|

|)

n ,2,1,0(/)(i i i i i ii ii j n

i j j ij j x x x x e i a a x a b x →-=?=≠-=∑≠=,

步骤3.求出}{max 1i n

i e e ≤≤=,若ε

?????=+--=-+-=--2.453.82102.7210321

321321x x x x x x x x x

步骤4.若,

+

<转步骤2继续迭代。若,

M

k≥表明迭代失败,停止计k→

1

M

,k

k

算。

2.算法流程图

四、程序设计

program jacobi

implicit none

integer::i,j

integer::k

save k

real,parameter::e=0.001

integer,parameter::n=3

real::x(n),y(n),b(n)

data b/7.2,8.3,4.2/

real::D

real::a(n,n)

open (unit=10,file='1.txt')

data a/10,-1,-1,-1,10,-1,-2,-2,5/

write(10,*)"**********矩阵A的形式为**********"

write(10,"(1x,3f6.2,/)")a

forall(i=1:n)

x(i)=0

end forall

k=0

100 D=0

do i=1,n

y(i)=b(i)

do j=1,n

if(i/=j) y(i)=y(i)-a(i,j)*x(j)

end do

y(i)=y(i)/a(i,i)

end do

do j=1,n

D=abs(x(j)-y(j))

end do

forall(i=1:n)

x(i)=y(i)

end forall

if(D>=e) then

k=k+1

write(10,*)"迭代次数为:",k

goto 100

else

goto 200

end if

200 write(10,*)"****************************************"

write(10,*)"用jacobi方法解得的结果X[t]为:"

write(10,"(1x,3f6.2,/)")x(:)

stop

end program

五、结果及讨论

1.实验结果

**********矩阵A的形式为**********

10.00 -1.00 -1.00

-1.00 10.00 -1.00

-2.00 -2.00 5.00

迭代次数为: 1

迭代次数为: 2

迭代次数为: 3

迭代次数为: 4

迭代次数为: 5

迭代次数为: 6

迭代次数为:7

****************************************

用jacobi方法解得的结果X[t]为:

1.10 1.20 1.30

2.讨论分析

(1)误差

从上述输出结果中可以看出,当迭代次数k增大时,迭代值x1,y1,z1

会越来越逼近方程组的精确解x=1.0,y=1.2,z=1.3。

(2)收敛性

在本题目中, 用雅可比迭代法和高斯-塞德尔迭代法分别求解该线性方程组,得到的近似根是收敛的

六、算法评价

优点:迭代法算法简单,编制程序比较容易。

缺点:迭代法要求方程组的系数矩阵有某种特殊性质(譬如是所谓对角占优阵)以保证过程的收敛性。高斯—塞德尔迭代法比雅可比迭代法收敛快(达到同样的精度所需迭代次数少),但这个结论,在一定条件下才是对的,甚至有这样的方程组,雅可比方法收敛,而高斯—塞德尔迭代法却是发散的。在雅可比迭代法求解线性方程组时,只要误差截断设计的合理,原则上可以得到很正确的解。而通常我们选取设计误差限或设计最大迭代次数的方法来控制。由于它的准确性,故在

实际应用中比较常见,对于解一般线性方程组非常有效准确。通过该算法以及编程对求解的过程,我们不难发现,雅克比迭代法的优点明显,计算公式简单,每迭代一次只需计算一次矩阵和向量的乘法,且计算过程中原始矩阵A始终不变,比较容易并行计算。然而这种迭代方式收敛速度较慢,而且占据的存储空间较大,所以工程中一般不直接用雅克比迭代法,而用其改进方法。

附:

高斯—赛德尔程序

program G-S

implicit none

integer::i,j

integer::k

save k

real,parameter::e=0.001

integer,parameter::n=3

real::x(n),y(n),b(n)

data b/7.2,8.3,4.2/

real::D

real::a(n,n)

open (unit=10,file='1.txt')

data a/10,-1,-1,-1,10,-1,-2,-2,5/

write(10,*)"**********矩阵A的形式为**********"

write(10,"(1x,3f6.2,/)")a

forall(i=1:n)

x(i)=0

end forall

k=0

100 D=0

do i=1,n

y(i)=b(i)

do j=1,n

if(i

if(i>j) y(i)=y(i)-a(i,j)*y(j)

end do

y(i)=y(i)/a(i,i)

end do

do j=1,n

D=abs(x(j)-y(j))

end do

forall(i=1:n)

x(i)=y(i)

end forall

if(D>=e) then

k=k+1

write(10,*)"迭代次数为:",k

goto 100

else

goto 200

end if

200 write(10,*)"****************************************"

write(10,*)"用Gauss-seidel方法解得的结果X[t]为:"

write(10,"(1x,3f6.2,/)")x(:)

stop

end program

**********矩阵A的形式为**********

10.00 -1.00 -1.00

-1.00 10.00 -1.00

-2.00 -2.00 5.00

迭代次数为: 1

迭代次数为: 2

迭代次数为: 3

迭代次数为: 4

****************************************

用Gauss-seidel方法解得的结果X[t]为:

1.10 1.20 1.30

数值实验报告

数值实验报告五 班级:2017级学号:**** 姓名:*** 2017.12.5 1.数值实验问题 试用雅可比迭代法,高斯-赛德尔迭代法,超松驰迭代计算线性方程组: 取=(0,0,0,松弛因子分别取w=0.1t,1要求达到精度 。试通过数值计算得出不同的松弛因子所需要的迭代次数和收敛最快的松弛因子,并指出哪些松弛因子使得迭代发散。 2.数值方法 A=, L=-, U=-, D=diag() (1)雅可比迭代公式:

D. (2)高斯-赛德尔迭代法公式: (3)超松驰迭代方法公式: 其中w为松弛因子。 3.数值结果 如下表

最后四组,测得其在前10次内迭代所产生的结果,其中每一列为一

次迭代结果,分别如图: SOR-1.6 SOR-1.7 SOR-1.8 SOR-1.9 由于计算数据限制,其前五十列数据基本为空,所以取51到60列

由此看出,最后四组数据是发散的,数据结果不稳定,不收敛。所以最后四组得不到所需数据。 4.讨论 本次实验,分别用雅可比迭代公式,高斯-赛德尔迭代公式,超松驰迭代公式计算了此线性方程组。其中,雅可比和高斯迭代能够很好的进行运算,而超松驰迭代方法中,若松弛因子取得不够恰当,则会导致整个运算失败,得不到所需的结果,迭代不收敛,发散。此外,在进行初始值的赋值中,我是对每个矩阵都进行了赋值操作,而更简便的是,调用matlab中存在的函数,对矩阵进行运算,从而简化操作和代码,也使程序适用性更广。 程序代码: 1.雅可比迭代 function [x]=yakebi(D,L,U,b,j) format long B=D\(L+U);

雅可比迭代法

2013-2014(1)专业课程实践论文 题目:雅可比迭代法

一、算法理论 设有方程组),...,2,1(1 n i b x a i j n j ij ==∑= 记作,b Ax = (1) A 为非奇异阵且),,...,2,1(0n i a ij =≠将A 分裂为U L D A --=,其中 D =????????????????nn a a a 22 11,L =-??? ????? ???? ????-00001,21323121n n n n a a a a a a U =-?? ? ?? ? ? ? ????????-0000,122311312n n n n a a a a a a 将式(1)第)....2,1(n i i =个方程用ii a 去除再移项,得到等价方程组 (),,...,2,111n i x a b a x n i j j j ij i ii i =??? ? ? ?? -=∑≠= (2) 简记作 ,0f x B x += 其中 ().,111 0b D f U L D A D I B ---=+=-= 对方程组(2)应用迭代法,得到解式(1)的雅可比迭代公式 () () ()()()()()????????? ?? ? ??- ==∑≠=+,1,...,11002010n i j i k j ij i ii k i t n x a b a x x x x x , 初始向量 (3)

其中()()()()()T k n k k k x x x x ,,...,21=为第k 次迭代向量。设()k x 已经算出,由式(3)可计算下一次迭代向量()(),,...,2,1,...;2,1,01n i k x k ==+ 显然迭代公式(3)的矩阵形式为 ()()()()???+=+,010f x B x x k k ,初始向量 其中0B 称为雅可比方法迭代矩阵。

MATLAB样例之雅克比迭代法

要求: 下面分别使用雅克比迭代法和高斯-赛德尔迭代法求一个方程组的近似解用的线性方程组是按实验要求给的: 7*x1+x2+2*x3=10 x1+8*x2+2*x3=8 2*x1+2*x2+9*x3=6 雅克比迭代法的matlab代码:(老师写的) A=[7,1,2;1,8,2;2,2,9]; b=[10;8;6]; if(any(diag(A))==0) error('error,pause') end eps=input('误差限eps='); N=input('迭代次数N='); D=diag(diag(A)); B=inv(D)*(D-A); f=inv(D)*b; K=0; x0=zeros(size(b)); while 1 x1=B*x0+f K=K+1; fprintf('第-次迭代的近似解为',K) disp(x1'); if norm(x1-x0,inf)N fprintf('迭代超限') end x0=x1; end 高斯-赛德尔迭代法matlab代码:(自己改的)

A=[7,1,2;1,8,2;2,2,9]; b=[10;8;6]; if(all(diag(A))==0) error('error,pause') end eps=input('误差限eps='); N=input('迭代次数N='); D=diag(diag(A)); B=inv(D)*(D-A); f=inv(D)*b; K=0; x0=zeros(size(b)); x00=x0; while 1 x11=B*x0+f; x00(1,1)=x11(1,1); x12=B*x00+f; x00(2,1)=x12(2,1); x13=B*x00+f; x00(3,1)=x13(3,1); x1=x00 K=K+1; fprintf('第-次迭代的近似解为',K) disp(x1'); if norm(x1-x0,inf)N fprintf('迭代超限') end x0=x1; end

数学实验“线性方程组的J-迭代,GS-迭代,SOR-迭代解法”实验报告(内含matlab程序代码)

西京学院数学软件实验任务书 课程名称数学软件实验班级数0901 学号0912020107 姓名李亚强 实验课题线性方程组的J-迭代,GS-迭代,SOR-迭代方法。 实验目的 熟悉线性方程组的J-迭代,GS-迭代,SOR-迭代方法。 实验要求运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成。 实验内容线性方程组的J-迭代;线性方程组的GS-迭代;线性方程组的SOR-迭代。 成绩教师

实验四实验报告 一、实验名称:线性方程组的J-迭代,GS-迭代,SOR-迭代。 二、实验目的:熟悉线性方程组的J-迭代,GS-迭代,SOR-迭代,SSOR-迭代方法,编程实现雅可比方法和高斯-赛德尔方法求解非线 性方程组121231 235210 64182514 x x x x x x x x +=?? ++=??++=-?的根,提高matlab 编程能力。 三、实验要求:已知线性方程矩阵,利用迭代思想编程求解线性方程组的解。 四、实验原理: 1、雅可比迭代法(J-迭代法): 线性方程组b X A =*,可以转变为: 迭代公式(0)(1)() k 0,1,2,....k k J X X B X f +???=+=?? 其中b M f U L M A M I B J 111),(---=+=-=,称J B 为求解 b X A =*的雅可比迭代法的迭代矩阵。以下给出雅可比迭代的 分量计算公式,令),....,() ()(2)(1)(k n k k k X X X X =,由雅可比迭代公式 有 b X U L MX k k ++=+) () 1()(,既有i n i j k i ij i j k i ij k i ij b X a X a X a +- -=∑∑+=-=+1 )(1 1 )() 1(, 于

Jacobi迭代法求解线性方程组实验报告

仿真平台与工具应用实践 Jacobi迭代法求解线性方程组 ^ 实验报告 : 院系: 专业班级: 姓名: 学号: 指导老师: }

一、 ; 二、 实验目的 熟悉Jacobi 迭代法原理; 学习使用Jacobi 迭代法求解线性方程组; { 编程实现该方法; 三、 实验内容 应用Jacobi 迭代法解如下线性方程组: , ?? ???=++--=+-=+-1552218474321321321x x x x x x x x x ,要求计算精度为710- 四、 ^ 五、 实验过程 (1)、算法理论 Jacobi 迭代格式的引出是依据迭代法的基本思想:构造一个向量系列(){}n X ,使其收敛至某个极限*X ,则*X 就是要求的方程组的准确解。 Jacobi 迭代

将方程组: ???????=+++=+++=+++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 22112222212111212111 )1( ~ 在假设0≠ii a ,改写成()??? ????++++=++++=++++=--n n n n n n n n n n n g x b x b x b x g x b x b x b x g x b x b x b x 112211223231212113132121 )2( 如果引用系数矩阵 ??????????=nn n n a a a a A 1111,???? ??????=0011 n n b b B 及向量??????????=n x x X 1,??????????=n b b b 1,??????????=n g g g 1, 方程组(1)和(2)分别可写为:b AX =及g BX X +=,这样就得到了jacobi 迭代格式01g BX X k k +=+用jacobi 迭代解方程组b AX =时,就可任意取初值0X 带入迭代可知式g BX X k k +=+1,然后求k k X ∞ →lim 。但是,n 比较大的时候,写方程组)1(和)2(是很麻烦的,如果直接由A ,b 能直接得到B ,g 就是矩阵与向量的运算了,那么如何得到B ,g 呢实际上,如果引进非奇异对角矩阵 ()0≠ii a ???? ??????=nn a a D 00011 将A 分解成:,D D A A +-=要求b AX =的解,实质上就有,)(DX X D A AX +-=而D 是非奇异的,所以1-D 存在,,)(X A D AX DX -+=从而有,11b D AX D X --+=我们在这里不妨令,1A D I B --=b D g 1-=就得到jacobi 迭代格式:g BX X k k +=+1 】

jacobi迭代法和Gauss-Seidel迭代法

数值计算方法实验报告(五) 班级:地信10801 序号:姓名: 一、实验题目:jacobi迭代法和Gauss-Seidel迭代法 二、实验学时: 2学时 三、实验目的和要求: 1.掌握迭代法的基础原理。 2.掌握jacobi迭代法和Gauss-Seidel迭代法的步骤。 3.能用程序语言对jacobi迭代法和Gauss-Seidel迭代法进行编程实现。 四、实验过程代码及结果 1、代码: #include #include float x[100],xk[100]; float e; int N,M=1000; float a[100][101]; void initdata() { cout<<"输入方程阶数:"; cin>>N; cout<<"输入误差限e:"; cin>>e; cout<<"输入方程系数:"<>a[i][j]; cout<<"输入初始解向量x0:"<>xk[i]; } void jocobi() { int Nx=0,times=0; while(Nx=M){cout<<"发散"<

for(int i=1;i<=N;i++) { float sum=0; for(int j=1;j<=N;j++) if(i!=j)sum+=xk[j]*a[i][j]; x[i]=(a[i][N+1]-sum)/a[i][i]; if(fabs(x[i]-xk[i])

科学学生实验报告单

四年级科学实验报告单

五年级科学实验报告单 1、唾液能消化淀粉的验证实验: 实验仪器:碘酒,滴管,试管,淀粉液、馒头等。 实验过程:取两个试管,分别加入等量的淀粉液,在其中一个试管中加入少量唾液,并摇晃,使其均匀混合。将两个试管放入温度为40摄氏度左右的温水中。过一会儿,分别往两个试管中放入一滴碘酒,观察现象。

实验现象:加入唾液的淀粉液没有变化,没有加入唾液的淀粉变蓝了。 实验结论:淀粉遇到碘酒会变成蓝色. 2、吸进的气体与呼出的气体是否相同的实验 实验仪器:水槽、玻璃吸管、集气瓶、烧杯、蜡烛、澄清的石灰水、火柴等。 实验一步骤: 1、用排水法收集呼出的气体,在水中用玻璃片将瓶口盖严,然后将瓶子从水中取出; 2 把瓶盖声上的玻璃片打开一个小口,将燃烧着的火柴慢慢放入瓶,看到什么现象?这说明什么? 实验一现象:燃烧的火柴熄灭了。 实验一结论:呼出的气体是不支持燃烧的气体。 实验二步骤: 1、按课本中的装置,经过弯玻璃管吸气,让瓶外空气经石灰水进入人体,石灰水有变化吗?(没有变化) 2经过直玻璃管向石灰水吹气,石灰水有变化吗?(有变化)这说明什么? 实验二结论:呼出的气体能使澄清的石灰水变浑浊。 概括出呼出的气体中含氧气少、二氧化碳多。推想出人体需要氧气,排出二氧化碳。 3、凸透镜成像 实验仪器:凸透镜、纸屏、蜡烛、火柴等。 实验步骤: 1、将点燃的蜡烛放于凸透镜和纸屏中间,立在桌上,使它们在一条直线上,并使火焰、镜面、纸屏的中心高度大体相同。 2、适当调整凸透镜与纸屏的距离,在纸屏上可以看到蜡烛的像吗?像是什么样的? 3、研究像的大小与成像的规律是怎样的? 实验结论:利用凸透镜形成的像都是倒立的。 1、当凸透镜距纸屏近,距蜡烛远时,形成的是缩小的像。 2、当凸透镜距纸屏远,距蜡烛近时,形成的是放大的像。 3、当凸透镜距纸屏和距蜡烛相等时,形成的是相等的像。

雅可比迭代法与矩阵的特征值

实验五 矩阵的lu分解法,雅可比迭代法 班级: 学号: 姓名:

实验五 矩阵的LU 分解法,雅可比迭代 一、目的与要求: 熟悉求解线性方程组的有关理论和方法; 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序; 通过实际计算,进一步了解各种方法的优缺点,选择合适的数值方法。 二、实验内容: 会编制列主元消去法、LU 分解法、雅可比及高斯—塞德尔迭代法德程序,进一步了解 各种方法的优缺点。 三、程序与实例 列主元高斯消去法 算法:将方程用增广矩阵[A ∣b ]=(ij a )1n (n )+?表示 1) 消元过程 对k=1,2,…,n-1 ①选主元,找{}n ,,1k ,k i k +∈使得 k ,i k a = ik a n i k max ≤≤ ②如果0a k ,i k =,则矩阵A 奇异,程序结束;否则执行③。 ③如果k i k ≠,则交换第k 行与第k i 行对应元素位置, j i kj k a a ? j=k,┅,n+1 ④消元,对i=k+1, ┅,n 计算 kk ik ik a a l /= 对j=l+1, ┅,n+1计算 kj ik ij ij a l a a -= 2) 回代过程 ①若0=nn a ,则矩阵A 奇异,程序结束;否则执行②。 ②nn n n n a a x /1,+=;对i=n-1, ┅,2,1,计算 ii n i j j ij n i i a x a a x /11,??? ? ? ?- =∑+=+ 程序与实例 程序设计如下:

#include #include using namespace std; void disp(double** p,int row,int col){ for(int i=0;i>p[i][j]; } } int findMax(double** p,int start,int end){ int max=start; for(int i=start;iabs(p[max][start])) max=i; } return max; } void swapRow(double** p,int one,int other,int col){ double temp=0; for(int i=0;i

实验记录及实验报告撰写要求1

环境工程实验报告撰写要求 环境工程实验报告应包括实验预习报告、实验原始记录和实验报告三部分,其中实验预习报告和实验原始记录需指导老师签字。实验报告采用河北科技大学实验报告标准纸手写。 实验预习报告、实验原始记录和实验报告要求如下。 一、实验预习报告要求 实验预习报告包括实验目的、实验原理、实验材料及装置、实验内容及实验步骤等,具体内容如下。 1、实验目的 实验目的要明确,在理论上验证定理、公式、算法,并使实验者获得深刻和系统的理解,在实践上,掌握使用实验设备的技能技巧和程序的调试方法。 2、实验原理 实验原理是指自然科学中具有普遍意义的基本规律,实验原理的表述的内容是实验设计的整体思路,即通过何种手段达到何种实验目的,还包括实验现象与结果出现的原因以及重要实验步骤设计的根据等。 3、实验材料与装置 实验所用的设备和材料。 4、实验步骤 从理论和实验两个方面考虑,要写明依据何种原理、定律算法、或操作方法进行实验,并详细写出理论计算过程。 其他内容参见实验报告专用纸。 二、实验原始记录 将实验现象和数据仔细地记录在实验原始记录中,做到原始记录准确、简练、详尽、清楚。如称量试材样品的重量、滴定管的读数、分光光度计的读数等,都应设计一定的表格准确记下正确的读数,并根据仪器的精确度准确记录有效数字。每一个结果至少要重复观测两次以上,符合实验要求并确知仪器工作正常后再写在实验报告上。 实验中使用仪器的类型、编号以及试剂的规格、化学式、分子量、准确的浓度等,都应记录清楚,以便总结实验完成报告时进行核对和作为查找成败原因的

参考依据。如果发现记录的结果有怀疑、遗漏、丢失等,都必须重做实验。三、实验报告要求 实验结束后,应及时整理和总结实验结果,在预习报告的基础上完成实验报告中的结果与讨论部分,包括: 1、数据处理和结果 实验数据处理和结果包括实验现象的描述,实验数据的处理等。对于实验结果的表述,一般有三种方法: (1)文字叙述 根据实验目的将原始资料系统化、条理化,用准确的专业术语客观地描述实验现象和结果,要有时间顺序以及各项指标在时间上的关系。 (2)图表和计算公式 用表格或坐标图或计算公式的方式使实验结果突出、清晰,便于相互比较,尤其适合于分组较多,且各组观察指标一致的实验,使组间异同一目了然。每一图表应有表头和计量单位,能说明一定的中心问题。 (3)曲线图 绘制曲线图,使变化趋势形象生动、直观明了。 在实验报告中,可任选其中一种或几种方法并用,以获得最佳效果。 2、问题和讨论 根据相关的实验结果及理论知识对所得到的实验结果进行解释和分析。如果所得到的实验结果和预期的结果一致,那么它可以验证什么理论?实验结果有什么意义?说明了什么问题?另外,也可以写一些本次实验的心得以及提出一些问题或建议等。 3、结论 针对这一实验所能验证的概念、原则或理论的简明总结,是从实验结果中归纳出的一般性、概括性的判断,要简练、准确、严谨、客观。 4、参考文献 本实验开展所需的文献。 注:经实验指导老师签字的实验原始记录表放在实验报告最后一起上交,同组人不能多于3人。

最新三年级科学下册实验报告单

实验一、温度和温度计 活动1:感受1号杯和2号杯里水的冷热 1号杯水() 2号杯水() 活动2:观察温度计 .观察常用液体温度计的主 要构造。 你观察温度计上有摄氏度 (℃)的标记吗? 你观察温度计上每一小格表 示多少? 最高()最低() 你观察温度计的最高温度和 最低温度是多少? 实验现象温度计里面的液柱热了就会上升,冷了就会下降。 活动3:下面的温度你会读和写吗? 28摄氏度写作: 20摄氏度写作: 零下5摄氏度写作: -21℃读作: 31℃读作: 实验要求:用温度计测量水的温度。 实验用品:400ml烧杯一个一支温度计适量冷水和一暖壶热水吸水纸废物瓶。 步骤操作要求评分标准满分得分1 清点仪器用品按材料清单清点材料用品是否齐全(5分)。 5

2 观察温度计的 零刻线、分度值 和量程。 A、观察温度计的零刻线。(10分) B、观察温度计的分度值和量程 。(10分) 20 3 用手感知水温。将手指伸入烧杯中(冷水)或将手放在烧杯 外壁(热水),手的感觉 (10分),估测水的温度(10分)。 20 4 将温度计测量 水的温度。 A、手拿温度计上端,将其竖直放入水中。(10 分) B、温度计的玻璃泡要完全浸没在水中,玻璃 泡不要碰烧杯的侧壁和底部。(10分) C、等示数稳定时再读数。读数时,要让玻璃 泡继续停留在水中。(10分) D、视线要和温度计的示数保持相平。连续三 次测水的温度分别为、、 ,平均水温为。(15分) 45 5 整理仪器,擦拭 桌面。 A、将温度计擦干放回原处。(5分) B、擦拭桌面。(5分) 10 实验三、水结冰了 一、实验名称:水结冰了 二、实验目的:观察水在不同温度下温度计的读数 三、实验步骤: 1、在试管里加入一半的纯净水,用温度计测量并记录试管里水的温度 2、拿一只保温杯(或在普通塑料杯外包裹一块干毛巾)在杯内装满碎冰, 把试管插入碎冰中,用温度计观测试管里水温的变化 3、在碎冰里加入较多的食盐,保持几分钟持续观测试管里的水温 4、观测试管里的水开始结冰时的温度 四、实验器材:试管、保温杯、温度计、碎冰块、食盐、纯净水。 水结冰了的实验记录表

雅克比迭代法和高斯-赛德尔法解线性方程组(C++)

作业:① 分别用J 法和G-S 法求解下列方程,并讨论结果。 123122*********x x x -?????? ??? ?= ??? ? ??? ??????? #include using namespace std; //J 法解线性方程 int main(){ int m,n,i,j,times=0,mtimes; double s,sum,max; cout<<"请输入系数矩阵行数m 、列数n :"<>m>>n; if(m>A[i][j]; cout<<"请输入常数向量B :"<>B[i]; cout<<"请输入最大允许误差s:"<>s; cout<<"请输入最大迭代次数:"<>mtimes; cout<<"请输入一零级向量X:"<>X[i]; T[i]=X[i];//T[]存放上一次迭代结果 }

金蝶实验报告记录

金蝶实验报告记录

————————————————————————————————作者:————————————————————————————————日期: 2

3 实验报告 一.实验目的 通过进行金蝶kis 个人全程实验,从建账到日常业务处理,报表练习,购销存业务处理等相关操作,了解该会计软件的主要功能及特点,更深入的掌握功初始化能的使用要领、日常业务处理的操作技巧及期末处理的使用方法,了解从编制记账凭证,审核记账凭证、过账、期末结账以及编制财务报表等一系列会计处理的全过程。通过个人全程实验,全面的了解和熟悉会计实务,掌握会计基本概念,对会计核算的基本程序和基本方法有了更深一步的掌握。 二.实验原理 借贷记账法规则:有借必有贷,借贷必相等。 资产=负债+所有者权益 三.实验设备 装有金蝶kis 专业版软件的电脑,金蝶实验教材 四.结果预测 根据实验教材以及掌握的知识进行实验,得到实验结果。将辉煌工业有限公司2012年发生的经济业务完整正确的体现在账套中,为该公司经营决策者及相关利益者,提供有用的财务信息,有利于做出正确的经济决策。 五、实验步骤 (一)建账 1.建账套。 (1)金蝶kis 专业版→工具→账套管理→新建。录入相关数据。 (2) “金蝶kis 专业版”,以 “manager”,密码为空,登录到“ao11\辉煌工业有限公司”,进入金蝶kis 专业版主界面。(3)“基础设置”模块, “系统参数” →“会计期间”→“设置会计期间”,启用会计年度为“2012” →“财务参数”,启用会计年度为“2012”,启用会计期间为“1”。 2.增加用户。 基础设置→用户管理→新建用户→录入相关用户及其权限。 (二)初始化 1.设置财务系统账套选项 (1)增加电汇结算方式:基础设置→结算方式→新增→代码:js02,名称:电汇 (2)选择凭证保存后立即新增:账务处理→凭证录入→查看→选项→选定“凭证保存后立即新增” →确定 2.设置业务系统账套选项 (1)增加存放地点。基础设置→核算项目→仓库→新增。依次录入。 3.基础资料 (1)增加美元币别 基础设置→币别→新增→币别代码:usd ,币别名称:美元,记账汇率: 8.7 (2)增加客户 基础设置→核算项目→客户→新增。 (3)增加供应商 基础设置→核算项目→供应商→新增。 (4)增加部门 基础设置→核算项目→部门→新增。 (5)增加职员 基础设置→核算项目→职员→新增。 (6)增加或修改会计科目

matlab Jacobi迭代法Gauss-seidel和SOR迭代

1.Jacobi迭代法 例1 用jacobi迭代法求解代数线性代数方程组,保留四位有效数字(err=1e-4) 其中A=[8 -1 1;2 10 -1;1 1 -5];b=[1 ;4; 3]。 解:编写jacobi迭代法的函数文件,保存为jacobi.m function [x,k]=jacobi(A,b,x0,eps,N) % 求解Ax=b;x0为初始列向量;eps为误差容限;N为最大迭代次数 % 输出x为近似解;k为迭代次数 n=length(A); x=zeros(n,1); for k=1:N for i=1:n ――――――― end if norm(x-x0,inf)

end x0=x; end 编写主程序如下 format long clear A=[8 -1 1;2 10 -1;1 1 -5]; b=[1 ;4; 3]; x0=[0.125; 0.4 ;-0.6 ]; % x0为初始列向量N为最大迭代次数err=1e-4; % err为误差容限 N=25; % N为最大迭代次数 [x,k]=jacobi(A,b,x0,err,N) 得到结果如下 x = 0.22492315625000 0.30561995000000 -0.49388680000000

k = 6 2.Gauss-seidel迭代法 例2 用Gauss-seidel迭代法求解代数线性代数方程组,保留四位有效数字(err=1e-4) 其中A=[8 -1 1;2 10 -1;1 1 -5];b=[1 ;4; 3]。 解:编写Gauss-seidel迭代法的函数文件,保存为gaus.m function [x,k]=gaus(A,b,x0,eps,N) % 求解Ax=b;x0为初始列向量;eps为误差容限;N为最大迭代次数% 输出x为近似解;k为迭代次数 n=length(A); x=zeros(n,1); for k=1:N for i=1:n ―――――― end if norm(x-x0,inf)

相关主题
文本预览
相关文档 最新文档