matlab實(shí)驗(yàn)九 π的近似計(jì)算

實(shí)驗(yàn)九n的近似計(jì)算【實(shí)驗(yàn)?zāi)康摹苛私鈭A周率兀的計(jì)算歷程。了解計(jì)算n的割圓術(shù)、韋達(dá)公式、級(jí)數(shù)法、拉馬努金公式、迭代法。學(xué)習(xí)、掌握MATLAB軟件有關(guān)的命令?！緦?shí)驗(yàn)內(nèi)容】利用韋達(dá)（VieTa）公式2 \:22+\/2\：2+\：2+<2y2+\：2+\2+司2...兀2 2 2 2計(jì)算n的近似值?！緦?shí)驗(yàn)準(zhǔn)備】割圓術(shù)韋達(dá)（VieTa）公式利用級(jí)數(shù)計(jì)算n拉馬努金（Ranmaunujan）公式迭代方法n的兩百位近似值【實(shí)驗(yàn)重點(diǎn)】1.圓周率的計(jì)算歷程【實(shí)驗(yàn)難點(diǎn)】1.圓周率的各種計(jì)算公式【實(shí)驗(yàn)方法與步驟】練習(xí)1用劉徽的迭代公式七=J^4=^,s.=32心.，X6,2“+1， & 62n62n+i 620計(jì)算丸的近似值。相應(yīng)的MATLAB代碼為>>clear;>>x=1;>>fori=1:30>>x=vpa(sqrt(2-sqrt(4-xA2)),15)%計(jì)算精度為15位有效數(shù)字>>S=vpa(3*2Ai*x,10)>>end計(jì)算可得x=.517638090205042S=3.105828541x=.261052384440103S=3.132628613x=.130806258460286S=3.139350203x=.654381656435522e-1S=3.141031951x=.327234632529735e-1S=3.141452472x=.163622792078742e-1S=3.141557608x=.818120805246955e-2S=3.141583892x=.409061258232818e-2S=3.141590463x=.204530736067660e-2S=3.141592106x=.102265381402739e-2S=3.141592517x=.511326923724832e-3S=3.141592619x=.255663463951308e-3S=3.1415926453.1415926513.1415926513.1415926533.1415926533.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.1415926543.141592654TOC\o"1-5"\h\zx=.127831732236766e-3 Sx=.639158661510219e-4 Sx=.319579330795908e-4 Sx=.159789665403054e-4 Sx=.798948327021645e-5 Sx=.399474163511619e-5 Sx=.199737081755909e-5 Sx=.998685408779670e-6 Sx=.499342704389851e-6 Sx=.249671352194927e-6 Sx=.124835676097464e-6 Sx=.624178380487320e-7 Sx=.312089190243660e-7 Sx=.156044595121830e-7 Sx=.780222975609150e-8 Sx=.390111487804575e-8 Sx=.195055743902288e-8 Sx=.975278719511453e-9 S練習(xí)2用韋達(dá)公式

2_v2<2+克。2+偵2+J2x2+\：2+2+顯???兀2 2 2 2計(jì)算n的近似值。相應(yīng)的MATLAB代碼為>>clear;>>x=1;>>fori=1:20;>>i>>x=vpa(x*cos(pi/2A(i+1)),30);%計(jì)算精度為30位有效數(shù)字>>pai=vpa(2/x,22)>>error=vpa(pai-2/x,22)>>end計(jì)算可得i=1pai=2.828427124746190097603i=2pai=3.061467458920718232298i=3pai=3.121445152258052404216i=4pai=3.136548490545939249125i=5pai=3.140331156954752858963i=6pai=3.141277250932772720892i=7pai=3.141513801144301048318i=8pai=3.141572940367091461583i=9pai=3.141587725277159716287error=-.3774484188e-21error=-.2223605969e-21error=-.579142540e-22error=.2408548301e-21error=-.3651343370e-21error=-.1464196722e-21error=.4482597912e-21error=.596014783e-22error=-.1189029953e-21i=10pai=3.141591421511200086037error=.3229610904e-21error=-.3774484188e-21error=-.2223605969e-21error=-.579142540e-22error=.2408548301e-21error=-.3651343370e-21error=-.1464196722e-21error=.4482597912e-21error=.596014783e-22error=-.1189029953e-21i=11pai=3.141592345570117830926error=.4967305463e-21i=12pai=3.141592576584872625535error=-.4344678192e-21i=13pai=3.141592634338562821891error=.4862357044e-21i=14pai=3.141592648776985437337error=.2757983366e-21i=15pai=3.141592652386591008149error=.956629086e-22i=16pai=3.141592653288992575505error=.4169549797e-21i=17pai=3.141592653514592792966error=-.2783657489e-21i=18pai=3.141592653570993021726error=-.1396314172e-21i=19pai=3.141592653585093078916error=-.1682315366e-21i=20