?? estimateareamc.m
字號:
function [Area] = EstimateAreaMC(xpoly,ypoly,maxsize,NbPoints,Method,VerboseOutput)
switch Method
case 'Halton'
P = haltonset(2);
RandomPoints = maxsize .* net(P,NbPoints);
k = 0;
case 'Sobol'
P = sobolset(2);
RandomPoints = maxsize .* net(P,NbPoints);
case 'Standard'
RandomPoints = maxsize .* rand(NbPoints,2);
otherwise
error('Invalid Method');
end
IN = inpolygon(RandomPoints(:,1),RandomPoints(:,2),xpoly,ypoly);
Area = maxsize .* maxsize * sum(sum(IN)) ./ NbPoints;
hold on;
if (VerboseOutput)
h = gcf;
plot(RandomPoints(IN(:,end),end-1),RandomPoints(IN(:,end),end),'g.','LineWidth',1.5);
plot(RandomPoints(~IN(:,end),end-1),RandomPoints(~IN(:,end),end),'rx','LineWidth',1.5);
set(h,'WindowStyle','Docked');
disp(['Area of the Polygon -> ' num2str(polyarea(xpoly,ypoly))]);
disp(['Estimated Area of the Polygon -> ' num2str( Area)]);
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -