matlab-pi的近似值代码..docVIP

利用基于的级数来算的近似数： =16—4=16 程序代码： clear; digits(160); syms x x=sym(0); for k=1:15 k x=x+sym(-1)^sym((k)-sym(1))/(sym(2)*sym(k)-sym(1))*(sym(16)/sym(5)^(sym(2)*sym(k)-sym(1))-sym(4)/sym(239)^(sym(2)*sym(k)-sym(1))); vpa(x,40) vpa(vpa(pi,60)-x,5) end 运算结果： k =1 ans =3.183263598326359832635983263598326359833 ans =-0.041671 k =2 ans =3.14059702932606031430453110657922889815 ans =0k =3 ans =3.141621029325034425046832517116408069706 ans =-0.000028376 k =4 ans =3.141591772182177295018212291112329795027 ans =8.8141e-7 k =5 ans =3.14159268240439951724025983607357586049 ans =-2.8815e-8 k =6 ans =3.141592652615308608149350747666502755367 ans =9.7448e-10 k =7 ans =3.141592653623554761995504593820311845927 ans =-3.3762e-11 k =8 ans =3.141592653588602228662171260486978513156 ans =1.191e-12 k =9 ans =3.141592653589835847485700672251684395509 ans =-4.2609e-14 k =10 ans =3.141592653589791696917279619620105448141 ans =1.5415e-15 k =11 ans =3.141592653589793294747374857715343543379 ans =-5.6285e-17 k =12 ans =3.141592653589793236391840944671865282509 ans =2.0708e-18 k =13 ans =3.141592653589793238539324592671865282509 ans =-7.6681e-20 k =14 ans =3.141592653589793238459788161264457875102 ans =2.8552e-21 k =15 ans =3.14159265358979323846275020767549235786 ans =-1.0682e-22 从结果看取15项部分和就可以精确到第21位小数。 利用高斯公式： =48+32—20计算的近似值。 程序代码： clear; digits(160); syms x x=sym(0); for k=1:20 k x=x+sym(-1)^(sym(k)-sym(1))/(sym(2)*sym(k)-sym(1))*(sym(48)/sym(18)^(sym(2)*sym(k)-sym(1))+sym(32)/sym(57)^(sym(2)*sym(k)-sym(1))-sym(20)/sym(239)^(sym(2)*sym(k)-sym(1))); vpa(x,40) vpa(vpa(pi,60)-x,5) end 运算结果： k =1 ans =3.144388167070395654407986493430228290391 ans =-0.0027955 k =2 ans =3.141587573607883117581411853379420096518 ans =5.08e-6 k =3 ans =3.141592664765766446467954462521640720131 ans =-1.1176e-8 k =4 ans =3.141592653562973161607826508581675751712 ans =2.682e-11 k =5 ans =3.141592653589860958821758446223180726536 ans =-6.772e-14 k =6 ans =3.141592653589793061616