- 积分
- 0
- 注册时间
- 2010-12-15
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-5-13 13:36:55
|
显示全部楼层
来自 陕西西安
6# dongsheng01
-
- function y=mianjiju(x,y)
- plot(x,y)
- n=size(x,2);
- for i=1:n-1
- k(i)=(x(i+1)-x(i))/(y(i+1)-y(i));
- beijihanshu{i}=@(yy)(k(i)*(yy-y(i))+x(i)).*yy;
- f(i)=quad(beijihanshu{i},y(i),y(i+1));
- end
- k(n)=(x(n)-x(1))/(y(n)-y(1));
- beijihanshu{n}=@(yy)(k(n)*(yy-y(1))+x(1)).*yy;
- f(n)=quad(beijihanshu{n},y(n),y(1));
- y=sum(f);
复制代码
以上是我按照Q=xy,P=0来算的,我取的x=[1 5 10 2];y=[1 2 5 3];
结果为:
需要注意的是,x,y对应为一个坐标,并且是与你的点的连线顺序有关的。我看你贴的程序里面x,y的顺序不太对,而且程序你也没给全。你可以按照你的P,Q的取法编个程序比较一下 |
|