- 积分
- 1
- 注册时间
- 2011-8-21
- 仿真币
-
- 最后登录
- 1970-1-1
|
发表于 2011-11-17 00:57:51
|
显示全部楼层
来自 江苏扬州
本帖最后由 stats01 于 2011-11-17 00:58 编辑
clear,clc
x=[...];y=[...];z=[...];z=z-min(z);
b=[-72.8789307972411,1006.01113795262,794.748796191598,-5030.05569361460,-5030.05567734056];
fx=@(b,x,y)b(1)+b(2)*x+b(3)*y+b(4)*x.^2+b(5)*y.^2;
figure(1),clf
scatter3(x,y,z,'o','markerfacecolor','k')
stem3(x,y,z,'filled')
[x1,y1]=meshgrid(min(x):range(x)/50:max(x),min(y):range(y)/50:max(y));
z1=fx(b,x1,y1);
pause(1)
hold on
surf(x1,y1,z1);
pause(1)
shading interp
pause(1)
alpha(.8)
axis tight
xlabel('X'),ylabel('Y'),zlabel('f') |
评分
-
1
查看全部评分
-
|