MATLAB初级使用.ppt

(补例1) :利用随机投点方法计算圆周率 坐标原点为圆心,第一象限内单位圆部分恰被正方形[0,1]×[0,1]所围,投入1000个均匀随机点,则随机点落入1/4单位圆内的概率为 function P=mlab76(N) X=rand(N,1);Y=rand(N,1); II=find(X.^2+Y.^2=1); m=length(II);P=4*m/N; plot(X(II),Y(II),.b),hold on theta=0:.01:pi/2; x=cos(theta);y=sin(theta); plot(x,y,r) 3.1560, 3.1520, 3.1400,······ 例2. 计算下面两条曲线所围区域面积 function S=area2(N) if nargin==0,N=2000;end X=2*rand(N,1)-1; Y=2*rand(N,1); II=find(Y=1+sqrt(1-X.^2)Y=abs(X)); m=length(II);S=4*m/N; x1=0:0.01:1;y1=x1; x2=1:-0.01:-1;y2=1+sqrt(1-x2.^2); x3=-1:0.01:0;y3=-x3; fill([x1,x2,x3],[y1,y2,y3],c) S=2.5460 例3. 计算两个半径为1的直交圆柱面所围成体积 x2 + y2 = 1 , x