- clc;
- clear;
- % syms x y
- % y=1
- % x=solve(y+x==2,x)
- % disp(x)
- % global h5;
- global h6;
- h1=figure
- h2=axes('parent',h1,...
- 'unit','normalized',...
- 'position',[0.2,0.3,0.65,0.6],...
- 'box','on')
- % ezsurf('x^2*exp(-x^2-y^2)')
- x=[-5, -2.5, -2.5, -2.5, -0.5, -2.5, -0.5, -2.5, -0.5, -2.5, -0.5, -2.5, -5.5, 1.5, 1.5, -2.5, 1.5, 2.5, -0.5, -2.5, -3, -3, -3, -3, -3, -3, -3, -3, -3, -3, 2.1, 2.1, 2.1, 2.1, 2.1, 2.1, 2.1, 2.1, -2, 5.5, 25, 45, -2, 7, 31, -1.5, -1.5, -1.5, 3.5, 27.52, -1.5, -1.5, -1.5, 3.5, 24.25, -2.5, -2.5, -2.5, 8, 28, -2.5, -1.5, 9, 28, -3, -3, 10, -3.89, 7, 18, -5.24, 7, 17, 25, -6.59, 7, 17, 27, 37, -2.92, -2.92, 11, -4.49, 8, 20, -6.29, 7, 7, 18, 29, -3, 10, 20, 30, -3, 8, 17, 27, -3, -3, -3, -3, 10, 10, 10, 10, 20, 20, 20, 20, 30, 30, 30, 30, -3, -3, -3, -3, 8, 8, 8, 8, 17, 17, 17, 17, 27, 27, 27, 27]
- y=[21, 21, 41, 41, 61, 61, 76, 76, 92, 92, 102, 102, 119, 119, 136, 136, 153, 153, 167, 167, 28, 35, 49, 56, 85, 113, 157, 162, 173, 179, 31, 51, 66, 89.5, 101, 108, 127.5, 144.5, 84.25, 84.25, 84.25, 84.25, 60, 60, 60, 84.25, 84.25, 84.25, 84.25, 84.25, 60, 60, 60, 60, 60, 84.25, 84.25, 84.25, 84.25, 84.25, 60, 60, 60, 60, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 84.25, 84.25, 84.25, 84.25, 60, 60, 60, 60, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 84.25, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60]
- z=[1151, 1141.75, 1151, 1134, 1151, 1125.1, 1151, 1124, 1151, 1122.25, 1151, 1124, 1151, 1130.5, 1151, 1130.5, 1151, 1130.5, 1151, 1138.75, 1128.13, 1120.48, 1111.2, 1102.8, 1093.5, 1106.75, 1112.5, 1112.5, 1132.5, 1144.5, 1124.05, 1108, 1098.4, 1092.7, 1096.2, 1100.53, 1109.2, 1109.2, 1088.5, 1092, 1092, 1088.5, 1094.2, 1094.2, 1094.2, 1131, 1122, 1113, 1103, 1106, 1138, 1127, 1116, 1107, 1111, 1120, 1120, 1120, 1096.3, 1096.3, 1122, 1122, 1098.5, 1098.5, 1128, 1128, 1128, 1120, 1120, 1120, 1111, 1111, 1111, 1111, 1102, 1102, 1102, 1102, 1102, 1140, 1128, 1128, 1116, 1116, 1116, 1104, 1104, 1104, 1104, 1104, 1102.2, 1102.2, 1102.2, 1102.2, 1105.2, 1105.2, 1105.2, 1105.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1102.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2, 1105.2]
- type=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 5, 5, 5, 5, 5, 5, 5, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7]
- % x=[0,0,1,1]
- % y=[0,0,0,0]
- % z=[0,1,0,1]
- subplot(1,2,1)
- leg=linspace(1,1,8)
- for ii = 1:length(x)
- hold on
- if type(ii)==0
- res1=scatter3(x(ii),y(ii),z(ii),'r', 'fill')
-
- end
- if type(ii)==5
- res8=scatter3(x(ii),y(ii),z(ii),'g', 'fill')
- end
- if type(ii)==1
- res2=scatter3(x(ii),y(ii),z(ii),'b', 'fill')
- end
- if type(ii)==2
- res3=scatter3(x(ii),y(ii),z(ii),'c', 'fill')
- end
- if type(ii)==3
- res4=scatter3(x(ii),y(ii),z(ii),'m', 'fill')
- end
- if type(ii)==4
- res5=scatter3(x(ii),y(ii),z(ii),'y', 'fill')
- end
- if type(ii)==6
- res6=scatter3(x(ii),y(ii),z(ii),'k', 'fill')
- end
- if type(ii)==7
-
- res7=scatter3(x(ii),y(ii),z(ii),36,[0.5,0.1,0.23], 'fill')
- end
-
- end
复制代码
|