MATLAB 用for语句实现循环结构 例题

MATLAB 用for语句实现循环结构 例题

for语句：

说明：

1. for语句针对向量的每一个元素执行一次循环体。

for k=[1,3,2,5]
    k
end

2. 退出循环之后，循环变量的值就是向量中最后的元素值。

for k=1:2:10
end
k

3. 当for语句中向量为空时，循环体一次也不执行。

for k=1:-2:10
    k
end

例题：计算圆周率Π

方法一：

y=0;
g=-1;
n=input('n=?');
for i=1:n
    g=-g;
    y=y+g*1/(2*i-1);
end
pai=4*y

方法二：向量求和的方法实现程序

%向量求和的方法实现程序
n=input('n=?');
x=1:2:(2*n-1);
y=(-1).^(2:n+1)./x;
pai=sum(y)*4

方法三：利用定积分的近似值求Π的近似值

程序如下：

a=0;
b=1;
n=input('n=?');
h=(b-a)/n;
x=a:h:b;
f=sqrt(1-x.*x);
s=[];
for k=1:n
    s1=(f(k)+f(k+1))*h/2;
    s=[s,s1];
end
pai=4*sum(s)

方法四：利用蒙特卡洛法求Π的近似值

概率P=落在圆内的点数/所投点的总数

s=0;
n=input('n=?');
for i=1:n
    x=rand(1);
    y=rand(1);
    if x*x+y*y