西安交大并行計(jì)算作業(yè)

1、.word格式,并行計(jì)算與程序設(shè)計(jì)作業(yè)班級(jí)：姓名：學(xué)號(hào)：1.1至1.3節(jié)作業(yè)1 .調(diào)試課件中的所有程序,并完成作業(yè),同時(shí)在程序結(jié)果中需要輸出個(gè)人信息；代碼：Programmainwrite(*,*)班級(jí)：,write(*,*)姓名：,print*,學(xué)號(hào)：end2.編寫一個(gè)數(shù)值求解方程的程序,方程為為(,),誤差小于1e-6.代碼：a=1b=5if(f(a)*f(b).LT.0)thenwrite(*,*)inter:(,a,b,)Loop1:dowhile(abs(f(a)-f(b).gt.10e-6).and.專業(yè).專注.word格式,$(abs(a-b).gt.10e-6)c=(a+b)/

2、2if(f(a)*f(c).le.0)thenb=celsea=cendifenddoLoop1write(*,*)x=,celsewrite(*,*)Pleaseinputrealintervalendifendrealfunctionf(x)f=4.1*x*3-5.3*x-11.8end結(jié)果：專業(yè).專注.word格式,1.4節(jié)作業(yè)1 .采樣簡(jiǎn)單離散求和法求下面積分值代碼：read(*,*)a,b,nh=(b-a)/(2.0*n)s=0.0 x=a+hf2=0.0f4=0.0loop1:doi=1,n-1x=x+hf2=f2+f(x)x=x+hf4=f4+f(x)enddoloop1專業(yè).專

3、注.word格式,s=h/3.0*(f(a)+f(b)+4.0*f4+2.0*f2)write(*,150)s100format(1x,a=,f8.2,2x,b=,f8.2,$2x,n=,i4)150format(1x,s=,f16.7)endrealfunctionf(x)f=x*2/(sin(x)+1)enddoubleprecisioni,ai,ydoubleprecisionsum=0i=1do10while(1/i.gt.1e-5)sum=1/i+sumi=i+110continuewrite(*,*)sum=,sumend結(jié)果：專業(yè).專注.word格式,F:Fortraritest

4、lDebugF:FortraritestlDebugJeM.e1.6節(jié)作業(yè)1 .調(diào)試課本中的所有程序；結(jié)果略2.用雙精度型數(shù)據(jù)計(jì)算：直到-代碼：doubleprecisionsum=0i=1do10while(1/i.gt.1e-5)sum=1/i+sumi=i+110continuewrite(*,*)sum=,sumend結(jié)果:專業(yè).專注.word格式,三角形三個(gè)頂點(diǎn)的坐標(biāo)分別為A(1.5,2.5),B(-2.5,1),C(1,-1),采用復(fù)型數(shù)據(jù)類型求三角形的面積和重心代碼:complexa,b,c,center,area,perimetera=(1.5,2.5)b=(-2.5,1)c=

5、(1,-1)center=(a+b+c)/3perimeter=(a+b+c)/2area=sqrt(perimeter*(perimeter-a)*(perimeter-b)*$(perimeter-c)write(*,10)center=,centerwrite(*,10)area=,area10format(1x,a,f6.3,f6.3)End結(jié)果：專業(yè).專注.word格式,3.4.嘗試編寫課件例子中對(duì)應(yīng)的破譯程序.有興趣可以嘗試修改原編碼程序中后移位數(shù) 當(dāng)前為5變化的情況,比方移動(dòng)位數(shù)為當(dāng)前位置序號(hào)代碼:programmaincharacter*80line1,line2charact

6、er11read(*,*)line1n=len(line1)do10i=1,nl1=line1(i:i)if(lge(l1,A).and.lle(l1,Z).or.$(lge(l1,a).and.lle(l1,z)thenm=ichar(l1)+iif(m.ge.ichar(a).and.m.le.ichar(z).or.專業(yè).專注.word格式,$(m.ge.ichar(A).and.m.le.ichar(Z)thenelsem=m-26endifline2(i:i尸char(m)elseline2(i:i)=l1endif10continuewrite(*,*)line2end結(jié)果：代碼1

7、.7節(jié)作業(yè)1 .調(diào)試課件中的所有程序；專業(yè).專注.word格式,(結(jié)果略)2.自由落體位移公式為：其中-,其中是初始位置,是初始速度,為重力加速度,為經(jīng)歷時(shí)間,編寫一段程序根據(jù)時(shí)間值求位移.設(shè)=1.2m,=2.5m/s(用格式輸入語句讀入t,并格式輸出結(jié)果.)代碼：programmainrealtwrite(*,*)Pleaseinputthetimeread(*,*)tcallsum(t,s)write(*,10)time=,twrite(*,10)displacement=,s10format(1x,a,f10.2)endsubroutinesum(t0,s)realt0,s,s0,v0,

8、gs0=1.2v0=2.5g=9.8s=(g*t0*2)/2+v0*t0+s0end專業(yè).專注.word格式,結(jié)果:1.8節(jié)作業(yè)1.采用三種方法矩形法、梯形法和Simpson法編程對(duì)進(jìn)行積分,積分上下限等參數(shù)需通過read語句從外部調(diào)入,并最終比擬三種方法與理論解進(jìn)行比擬；代碼：C矩陣法realerrorread(*,*)a,b,nx=ah=(b-a)/nf0=x*sin(x)s=0.0do10i=1,nsi=f0*h專業(yè).專注.word格式,s=s+six=x+hf0=x*sin(x)continuewrite(*,100)a,b,nwrite(*,200)sformat(1x,a=,f10

9、.3,3x,b=,$f10.3,3x,n=,i4)format(1x,s=,f15.8)error=sin(b)-sin(a)-b*cos(b)+a*cos(a)-swrite(*,300)errorformat(1x,error=,f16.10)endC梯形法realerrorread(*,*)a,b,nx=ah=(b-a)/ns=0.0do10i=1,nsi=(x+(i-1)*h)*sin(x+(i-1)*h)+專業(yè).專注10100200300.word格式,$(x+(i-1)*h)*sin(x+i*h)*h/2.0s=s+si10continuewrite(*,100)a,b,nwrit

10、e(*,200)s100format(1x,a=,f10.3,3x,$b=,f10.3,3x,n=,i4)200format(1x,s=,f15.8)error=sin(b)-sin(a)-b*cos(b)+a*cos(a)-swrite(*,300)error300format(1x,error=,f16.10)endCsinpson法realerrorread(*,*)a,b,nh=(b-a)/(2.0*n)s=0.0fa=fun(a)fb=fun(b)x=a+hf2=0.0專業(yè).專注.word格式,f4=fun(x)do10i=1,n-1x=x+hf2=f2+fun(x)x=x+hf4=

11、f4+fun(x)continues=h/3.0*(fa+fb+4.0*f4+2.0*f2)write(*,100)a,b,nwrite(*,150)sformat(1x,a=,f8.2,2x,b=,f8.2,$2x,n=,i4)format(1x,s=,f16.7)error=sin(b)-sin(a)-b*cos(b)+a*cos(a)-swrite(*,200)errorformat(1x,error=,f16.10)endrealfunctionfun(x)fun=x*sin(x)end矩陣法結(jié)果：專業(yè).專注10100150200結(jié)果.word格式,梯形法結(jié)果：Sinpson法結(jié)果：分

12、析：從上面得到的結(jié)果可以看出矩陣法誤差最大,sinpson法次之,梯形法最優(yōu),但是從數(shù)學(xué)上分析,Sinpson法應(yīng)得到最好的結(jié)果,在專業(yè).專注.word格式,此題中誤差較小,在可以接受的誤差范圍之內(nèi),如果遇到變化較大的函數(shù)進(jìn)行積分,可以預(yù)見的是,sinpson法最優(yōu),梯形法次之,矩形法誤差最大.3.編程對(duì)任意系數(shù)一元屢次方程進(jìn)行求根,求根范圍和方程系數(shù)通過read語句從外部調(diào)入(采用一種方法即可,注意無根的處理)；代碼：programmain5read(*,*)x1,x2,a1,a2,a3,a4f1=a1*x1*3+a2*x1*2+a3*x1+a4f2=a1*x2*3+a2*x2*2+a3*x

13、2+a4if(sign(f1,f2).eq.f1)goto5f=1.020if(abs(x1-x2).gt.1e-5).and.$abs(f).gt.1e-6)thenx=x2-(x2-x1)/(f2-f1)*f2f=a1*x*3+a2*x*2+a3*x+a4if(sign(f,f1).eq.f)thenx1=xf1=felse專業(yè).專注.word格式,x2=xf2=fendifgoto20endifif(abs(f).gt.1e-6)x=(x1+x2)/2.0write(*,100)x100format(1x,x=,f15.7)end結(jié)果：4.改編課件中求函數(shù)極值程序,實(shí)現(xiàn)最大值的求解,函數(shù)

14、為:代碼：reallow,high,x1,x2read(*,*)low,high專業(yè).專注.word格式,write(*,200)x1=low+0.618*(high-low)x2=high-0.618*(high-low)10if(high-low.gt.1e-4)thenf1=3*x1*x1-8.0*x1-7.0f2=3*x2*x2-8.0*x2-7.0write(*,202)x1,f1,x2,f2if(f1.gt.f2)thenhigh=x1x1=x2x2=high-0.618*(high-low)elselow=x2x2=x1x1=low+0.618*(high-low)endifgo

15、to10endiff1=3*x1*x1-8.0*x1-7.0f2=3*x2*x2-8.0*x2-7.0if(f1.gt.f2)thenf=f2專業(yè).專注.word格式,x=x2elsef=f1x=x1endifwrite(*,204)x,f200format(12x,x1,14x,f1,$13x,x2,13x,f2/)202format(1x,4f15.7)204format(0,x=,f10.6,5x,$f(x)=,f16.10)end結(jié)果:專業(yè).專注.word格式,flfl1.G1800091.G18000912,(J90228112,(J9022811.38280001.3828000-

16、12.3262281-12.32622811 1 382382 的-12-32fi2281-12-32fi22811.23607661.2360766-12.304564-12.3045641.472105B1.472105B-12-2755&04-12-2755&041.3821.382的回目-12-3262281-12-32622811.36200081.3620008-12.32&2281-12.32&22811.32623511.3262351-12-3331823-12-3331823i.3262391i.3262391-12-3331523-12-3

17、3315231.29181901.2918190-12-3281631-12-32816311.347550?1.347550?-12-3327265-12-33272651.32623911.3262391-12-3331823-12-33318231.32623911.3262391-12-3331&23-12-3331&231,31310861,3131086-12.3321066-12.33210661.33439391.3343939-12-3333302-12-33333021.32623911.3262391-12-3331823-12-33318231.3394

18、0981.3394098-12-3332224-12-33322241,33439391,3343939-12.3333302-12.33333021.33439391.3343939-12.3333302-12.33333021.33127031.3312703-12.3333206-12.33332061.33630051.3363005-12-3333073-12-33330731.33439391.3343939-12-3333302-12-33333021.33439391.3343939-12-3333302-12-33333021.33319191.3331919-12-3333330-12-33333301.33319191.3331919-12.3333330-12.33333301.33246351.3324635-12-3333311-12-33333111.2336541.233654-12-333333B-1