Matlab作业习题与答案详解（附程序）.pdfVIP

Matlab作业习题与答案详解（附程序）

clearall;clc;closeall;

x=-10:0.01:20;

y=4*sin(x)./x;

ymin=min(y)

二、蒙特卡罗算法的数值计算

当前的油位高度是2.3米，见图1。模拟油流进储油罐的过程（图

维数任选），请计算罐内油量。三维的效果图参见图2。储油罐由两部

分组成，中间是圆柱体，两端是球罐体。（本题简化自2011年

UCMCMA题《储油罐的变位识别与罐容表标定》，细节参见原题原

题附件2cumcm2010A.doc。）

图1

图2

主程序：

clc;

clearall;

closeall;

center1=[-3.375,0,1.5];%左球罐中心center2=[3.375,0,1.5];%

右球罐中心

n=10000;%每次的撒点数

delta=0.02;%层高

h=3;

en=h/delta;

Show;%画出油罐

fori=0:en-1

x=(rand(1,n)-0.5)*10;%随机生成n个点

y=(rand(1,n)-0.5)*h;

z=(rand(1,n)*delta+i*delta);

Z=[x;y;z];

[dis1

dis2]=juli(center1,center2,Z);%算出各点对应的距离

index=find(((x-4x4)dis21.5)|(x-

4|x4)dis11.625);%找出在罐内的点

plot3(x(index),y(index),z(index),.k);%画出在罐内的点

drawnow

end

子程序1：

function[dis1dis2]=juli(a,b,q)

d11=q(1,:)-a(1);

d12=q(2,:)-a(2);

d13=q(3,:)-a(3);

d1=sqrt(d11.^2+d12.^2+d13.^2);

d21=q(1,:)-b(1);

d22=q(2,:)-b(2);

d23=q(3,:)-b(3);

d2=sqrt(d21.^2+d22.^2+d23.^2);

d1(d1d2)=d2(d1d2);

dis1=d1;

dis2=sqrt(d12.^2+d13.^2);

子程序2：

functiontu=Show

%===圆柱部分==

figure(color,w)

h=3;

y=3/2*cos(s);

z=3/2*sin(s)+1.5;

x=t;

ezmesh(x,y,z,[0pi*2-44]);

axisequal

hiddenoff

holdon

%===左罐部分==

z=cos(s)*(2.6406-t.^2).^(1/2)+1.5;y=sin(s)*(2.6406-

t.^2).^(1/2);

x=t-3.375;

ezmesh(x,y,z,[0pi*2-1.625-0.625]);

axisequal

hiddenoff

holdon

%===右罐部分===

z=cos(s)*(2.6406-t.^2).^(1/2)+1.5;

y=sin(s)*(2.6406-t.^2).^(1/2);

x=t+3.375;

ezmesh(x,y,z,[0pi*20.6251.625]);

colormap(gray)

axisequal