scatter的图和plot的图分开显示

jerrynigel

借鉴了论坛大佬的code加上自己的修改添加，想要请问，最后怎么把scatter和plot的图分开显示在两个不同的图里面？

1. x=[0 6.8 18.6 21.3 3 12.5 13 20.7 17.7 20.8 1.1 6.3 13.2 3.4 3.5 18.9 11.3 13.1];
   
2. y=[1.7 1.1 0.8 1.2 3.1 3 3.1 4.2 4.1 5.8 6.4 5.8 7.4 5.3 5.9 9.7 7.6 10.7];
   
3. z=[0.59 0.58 0.53 1.25 1.24 0.2 0.18 0.14 0.64 0. 1.06 0.73 0.15 0.15 0.1 0.08 0.17 0.8];
   
4. scatter3(x,y,z)
   
5. x=x';
   
6. y=y';
   
7. z=z';
   
8. hold on
   
9. Z=[ones(length(x),1),x,y,x.^2,x.*y,y.^2,x.^3,x.^2.*y,x.*y.^2,y.^3];
   
10. A=Z\z;
    
11. X=min(x)-1:0.2:max(x)+1;
    
12. Y=min(y)-1:(max(y)-min(y)+2)/(length(X)+1):max(y)+1;
    
13. [x y]=meshgrid(X,Y);
    
14. z=A(1)+A(2)*x+A(3)*y+A(4)*x.^2+A(5)*x.*y+A(6)*y.^2+A(7)*x.^3+A(8)*x.^2.*y+A(9)*x.*y.^2+A(10)*y.^3;
    
15. mesh(x,y,z)
    
16. max_where=find(z==max(z));
    
17. h = max(max(z)) %h deflection
    
18. P = 10; % P Pressure, Unit: KPa
    
19. a = 10; % a Bulge window size: nm
    
20. v = 0.25; % v Possion's ratio
    
21. t = 500; % t Thickness of the sample
    
22. Sigma_xx = (P*R)/t % Stress
    
23. E = ((P*4*(a^4))*(1-(v^2)))/(3*(h^3)*t) % E Young's modulus
    
24. Epsilon = ((P*R)/t)*(1-(v^2))/E %Strain
    
25. plot(Sigma_xx, Epsilon,'*')

复制代码

20141303

删去hold on在scatter和plot前加一句figure试试

jerrynigel

20141303 发表于 2019-12-3 08:48
删去hold on在scatter和plot前加一句figure试试

没看懂，能否写出来吗？

jerrynigel

20141303 发表于 2019-12-3 08:48
删去hold on在scatter和plot前加一句figure试试

OK了，太感谢了！

jerrynigel

20141303 发表于 2019-12-3 08:48
删去hold on在scatter和plot前加一句figure试试

我还想请教一个问题，如果把xyz的散点和在一起形成一个x的矩阵，然后用surf（x）的语法，那怎么更改后面的code得到这个surface的拟合呢？

20141303

好像可以使用interp函数和griddata函数，查查帮助文档吧