Matlab用不同公式计算Π值

GUI界面介绍

• 变化图：此部分将展示出各公式计算Π时，随着迭代次数的变化，Π值的变化趋势图。
• 公式：此部分将展示出各个计算公式，便于使用者理解。
• 参数设置：使用者可在此部分设置所需要展示的公式以及迭代次数（默认为100次）。
• Π值显示：此部分将展示各公式计算出的Π值，默认保留小数点后9位，便于通过数值比较各公式间的精度。

四个计算公式

GUI效果图

部分代码

% 四种方法计算Π的公式
% Formula1
function PI = Formula1(n)
PI=1;a=1;item=1;f=1;
for i=1:n
    a=a+2;
    f=-f;
    item=f*(1/a);
    PI=PI+item;
end
digits(10);  % 设置精度
PI=vpa(PI*4);
% Formula2
function PI = Formula2(n)
PI=1;a=1;item=1;
for i=1:n
    item=(a+1)*(a+1)/(a*(a+2));
    PI=PI*item;
    a=a+2;
end
digits(10);
PI=vpa(PI*2);
% Formula3
function PI = Formula3(n)
PI=1;a=0;b=1;item=1;
for i=1:n
    a=a+1;
    b=b+2;
    item=item*(a/b);
    PI=PI+item;
end
digits(10);
PI=vpa(PI*2);
% Formula4
function PI = Formula4(n)
PI=1;a=1;item=1;
for i=1:n
    item=1+1/(a*(a+2));
    PI=PI*item;
    a=a+2;
end
digits(10);
PI=vpa(PI*2);

function pbtShow_Callback(hObject, eventdata, handles)
choice=round(get(handles.FormulaChoice,'Value'));
%disp(choice)
N=round(str2num(get(handles.epoch,'String')));
switch choice
    case 1
        for i=1:N
            PI(i)=Formula1(i);
        end
        sLegend='Formula1';
    case 2
        for i=1:N
            PI(i)=Formula2(i);
        end
        sLegend='Formula2';
    case 3
        for i=1:N
            PI(i)=Formula3(i);
        end
        sLegend='Formula3';
    case 4
        for i=1:N
            PI(i)=Formula4(i);
        end
        sLegend='Formula4';
    case 5
        for i=1:N
            PI(1,i)=Formula1(i);
            PI(2,i)=Formula2(i);
            PI(3,i)=Formula3(i);
            PI(4,i)=Formula4(i);
        end
        sLegend={'Formula1','Formula2','Formula3','Formula4'};
    otherwise
        disp('Please choice correct parameter !');
        return;
end
% 显示Π变化趋势图
axes(handles.axes1);
plot(1:N,PI);% 绘图
line([1,N],[pi,pi],'color','r','linestyle',':');% 显示参考线
ylim([2.6,3.6]);% 对Y轴设定显示范围 
legend(sLegend);
xlabel('epoch');
ylabel('value');

% 显示某一公式计算出的Π值
if choice < 5
    set(handles.showpi,'String',char(PI(N))); 
else
    set(handles.showpi,'String',''); 
end