数学实验报告(基于matlab).docx

西安交通大学 实验报告 课程：高数实验 班级：自动化32 姓名：张文杰 学号：2130504049 日期：2014年4月20日计算值其它公式：请按照上述公式分别编程计算并与 （1）就精度与叠加次数进行比较，能得出什么结论？？解：（1）按照以下的程序得到n取不同的值（1）式计算的的近似值format longx=1:6;%首先计算原函数的由10的1到6次方展开时的值for n=10.^x%利用循环一次性将公式展开为不同的次数 s=0; digits(25); for k=1:n s=s+4*(-1)^(k+1)/(2*k-1);%计算出的近似值 endvpa(s,25) %定义精度为25位end （2）分别计算a,b,c,d四个公式的 的近似值x=1:6;for m=10.^x %计算上面四个公式的不同展开次数的值fa=0;fb=0;fc=0;fd=0 %分别表示四个等式的值digits(25); for k=1:m fa=fa+1/(2*k-1)^2; fb=fb+(-1)^(k-1)/k^2; fc=fc+(-1)^(k-1)/(2*k-1)^3; fd=fd+sin(2*k-1)/(2*k-1)^3; endfa=vpa(sqrt(8*fa),25)fb=vpa(sqrt(12*fb),25)fc=vpa((32*fc)^(1/3),25)fd=vpa((1+sqrt(1+32*fd))/2,25) %end%（3）将体现四个公式收敛情况的图形画到一个直角坐标系上pfa=[];pfb=[];pfc=[];pfd=[];ps=[]; %绘制出体现收敛速度的图像t=100; %选取的项数为150step=5; %步长选为5for j=10:step:t% fa=0;fb=0;fc=0;fd=0;s=0; for k=1:j fa=fa+1/(2*k-1)^2; fb=fb+(-1)^(k-1)/k^2; fc=fc+(-1)^(k-1)/(2*k-1)^3; fd=fd+sin(2*k-1)/(2*k-1)^3; s=s+4*(-1)^(k+1)/(2*k-1); end pfa=[pfa,sqrt(8*fa)]; pfb=[pfb,sqrt(12*fb)]; pfc=[pfc,(32*fc)^(1/3)]; pfd=[pfd,(1+sqrt(1+32*fd))/2];% ps=[ps,s];endj=10:step:t;plot(j,pfa,b,j,pfb,g,j,pfc,r,j,pfd,y,j,ps,k) %绘制收敛的图像程序的运行结果：叠加次数函数 结果函数 结果原函数的不同展开次数的值ans=3.04183961892940324389659ans=3.131592903558553686593768ans=3.1405926538397941349956ans=3.141492653590034489496929ans=3.141582653589719775766298ans=3.141591653589774324473183N=10fa=3.10962545798864775647985fb= 3.132977195469482278866735fc=3.141526087929505717255552fd=3.141538602696021698079676N=100fa=3.138407967067091242086008fb=3.141498114035649713571274fc=3.141592586052465385648702fd=3.141592599756393155985279N=1000fa=3.141274327602740168430273fb=3.141591699614915800253812fc=3.141592653522246259001349fd=3.141592653631021470062024N=10000fa=3.141560822439948719164704fb=3.141592644041437232260705fc=3.141592653589793238462643fd=3.141592653589710071315722N=100000fa=3.141589470489344115833319fb=3.141592653494265530156326fc=3.141592653589793238462643fd=3.141592653589793238462643N=1000000fa=3.141592335279969727679372fb=3.141592653588781924867135 fc= 3.141592653589793