matlab寻找闭合,关于 求闭合曲线的面积！

我最初抄了个程序：

x=[0.15 0.26 0.38 0.50 0.64 0.76 0.89 1.02 1.14 1.15 1.03 0.91 0.79 0.67 0.55 0.43 0.32 0.19  -0.15 -0.26 -0.39 -0.51 -0.62 -0.75 -0.89 -1.01 -1.14 -1.15 -1.03 -0.90 -0.79 -0.67 -0.54 -0.42 -0.30 -0.18 0.15];

y=[3.99 5.42 6.91 8.33 10.01 11.44 12.89 14.28 15.66 15.53 13.88 12.15 10.50 8.80 7.18 5.51 3.91 2.23 -2.45 -3.81 -5.39 -6.89 -8.32 -9.82 -11.38 -12.79 -14.20 -14.02 -12.33 -10.63 -8.91 -7.23 -5.56 -3.86 -2.23 -0.59 3.99];

plot(x,y);

n=length(x);

s=0;

for i=1:n-1

a=x(i)*y(i+1)-x(i+1)*y(i);

s=s+a;

end

S=1/2*s;  %多变形面积

S

S =

-3.6368.

但当我加入绝对号时

x=[0.15 0.26 0.38 0.50 0.64 0.76 0.89 1.02 1.14 1.15 1.03 0.91 0.79 0.67 0.55 0.43 0.32 0.19  -0.15 -0.26 -0.39 -0.51 -0.62 -0.75 -0.89 -1.01 -1.14 -1.15 -1.03 -0.90 -0.79 -0.67 -0.54 -0.42 -0.30 -0.18 0.15];

y=[3.99 5.42 6.91 8.33 10.01 11.44 12.89 14.28 15.66 15.53 13.88 12.15 10.50 8.80 7.18 5.51 3.91 2.23 -2.45 -3.81 -5.39 -6.89 -8.32 -9.82 -11.38 -12.79 -14.20 -14.02 -12.33 -10.63 -8.91 -7.23 -5.56 -3.86 -2.23 -0.59 3.99];

plot(x,y);

n=length(x);

s=0;

for i=1:n-1

a=x(i)*y(i+1)-x(i+1)*y(i);

s=abs(s+a);

end

S=1/2*s;  %多变形面积

S

S =

0.2423

这是怎么个情况，，应该相等才对呀！！

1.jpg

(17.94 KB, 下载次数: 8)

2011-12-3 23:11 上传

这是图形曲线