使用程式Matlab,割圓術求圓周率

2021-03-03 21:51:49 字數 2988 閱讀 6885

1樓:我行我素

function calpi(n) %定義函式bai,n是計算次du數,正整數

syms a; %定義符zhi號變數

for i=1:n %開始迴圈

a=sqrt(2-sqrt(4-a^dao2)); %計算end

a=subs(a,'a','1');%替代計算y=3*2^n*vpa(a,n+5);%求最終內結果容pi=y %顯示最終結果

利用內切圓求圓周率的matlab程式

2樓:神話夢碎

^n=1;

pi4=1/n;

delta=abs(pi-pi4*4);

jingdu=0.001;

ticwhile delta>jingdun=n+1;

pi4=pi4+(-1)^源(n-1)*1/(2*n-1);

delta=abs(pi-pi4*4);

ends1=sprintf('精度達到%f時的項數%d',jingdu,n);

s2=sprintf('pi的近似值是%f', pi4*4);

disp(s1),disp(s2)

h=toc;

s3=sprintf('使用迴圈語句的時間是%f秒',h);

disp(s3)

tick=1:1:1000;

pi_4=sum((-1).^(k-1).*1./(2*k-1));

cha=abs(pi-4*pi_4);

h1=toc;

st1=sprintf('精度達到%f時的項數%d',cha(end),k(end));

st2=sprintf('pi的近似值是%f', pi_4*4);

st3=sprintf('使用向量的時間是%f秒',h1);

disp(st1);

disp(st2);

按上面程式試試

matlab計算圓周率

3樓:匿名使用者

^第一bai問哈du~

zhia=1;

n=0;

my_pi=0;

while(a>0.01)

n=n+1;

my_pi=my_pi-1/(2*n-1)*(-1)^n;

a=1/(2*n-1)*4;

endmy_pi=my_pi*4

第二問dao哈專~屬

a=1;

n=0;

my_pi=0;

while(a>0.001)

n=n+1;

my_pi=my_pi-1/(2*n-1)*(-1)^n;

a=abs(pi-4*my_pi);

endn

利用matlab求圓周率

4樓:我行我素

n=1;

pi4=1/n;

delta=abs(pi-pi4*4);

jingdu=0.001;

ticwhile delta>jingdun=n+1;

pi4=pi4+(-1)^(n-1)*1/(2*n-1);

delta=abs(pi-pi4*4);

ends1=sprintf('精度達到%f時的項數%d',jingdu,n);

s2=sprintf('pi的近似值是%f', pi4*4);

disp(s1),disp(s2)

h=toc;

s3=sprintf('使用迴圈語句的時間是%f秒',h);

disp(s3)

tick=1:1:1000;

pi_4=sum((-1).^(k-1).*1./(2*k-1));

cha=abs(pi-4*pi_4);

h1=toc;

st1=sprintf('精度達到%f時的項數%d',cha(end),k(end));

st2=sprintf('pi的近似值是%f', pi_4*4);

st3=sprintf('使用向量的時間是%f秒',h1);

disp(st1);

disp(st2);

按上面程式試試

用matlab寫程式求圓周率(精度e=10^-7),求大神幫助,萬分感謝 5

5樓:我行我素

圓周率來的計算

方法很多,我提供自乙個:bai

n=1;

absol=0.1;

while absol>1e-7

syms a

for i=1:n

a=sqrt(2-sqrt(4-a^2));

enda=subs(a,'a','1');

y=3*2^n*vpa(a,n+5);

absol=abs(double(pi-y));

disp('絕對誤差大於du1e-7,再計算zhi')

n=n+1;

enddisp('絕對誤差小於dao1e-7,不再計算')s=sprintf('計算次數n=%d',n);

disp(s)

pi=y

6樓:匿名使用者

(1)不要顯bai示的在語句後面加上

du分號zhi

(2)不知道你有沒有儲存在dao

正確的路徑上

(3)在命令符內下先clear一下

我修改容了一點點,你看一下:

function calpi(n)

syms a;

for i=1:n

a=sqrt(2-sqrt(4-a^2));

enda=subs(a,'a','1');

y=3*2^n*vpa(a,n+5);

pi=y

儲存在work資料夾下,calpi.m

>>clear

>>calpi(100)

>>pi =

3.1415926535897932384626433832795你自己再弄一下,應該可以的!

編寫matlab程式分別使用if語句和switch語句實

function fun a input a b input b c input c x input x if x 0.5 x 1.5 y a x zhi2 b x c elseif x 1.5 x 3.5 y a sin c b x elseif x 3.5 x 5.5 y log abs b x...

求MATLAB簡單潮流計算程式,求MATLAB簡單潮流計算程式。。。。 20

function lianxuchaoliu clear clc n 9 節點數 nl 9 支路數 isb 1 平衡節點號 pr 0.00001 誤差精度 b1 1 4 0.0576i 0 1.05 1 4 5 0.017 0.092i 0.158i 1 0 5 6 0.039 0.17i 0.35...

matlab中,使用什麼程式能增加影象的對比度

matlab中,使用什麼程式能增加影象的對比度的方法。如下參考 1 介紹影象在matlab中的運算。首先,介紹影象閱讀。複製使用imread命令 如aa png檔案,我們要讀取,有兩種方式 一是直接雙擊,下面出現對話方塊,點選完成。二是使用命令 aa imread aa png 2 影象生成後,以a...