不知道怎样利用matlab进行分组运算，哪位大神给指导一下吧

adogpp

各位大神好，本人现在遇到了一个棘手的问题，不知如何使用matlab进行分组预算。
如图所示，本人已将数据进行分组完毕，运算部分的代码很简单，本人也已经完成，只是这个识别分组这里不太会，如果不能实现分组识别计算的话，几万组数据只保存在一个excel文件中，手动输入的话，确实有些繁琐，烦请各位老师指导一下。
9 V6 I# y; N7 Q7 L/ B4 k本人已完成的代码部分也一并提交，代码中的数据为其中一组数据，望各位老师予以添加修改并指导，多谢！

目前本人的代码只是将图片中的前两列数据手动复制进代码，本人想直接将图中的三列数据全部导进去，并根据第三列的数组进行分组计算（数字相同的为一组）

ps：本人使用的是Matlab r2018a，运行环境的macos

% g, ]. T1 v- D

/ y3 t- P4 U. e0 D6 m* C
x = [4.1441    0
) M1 N1 ?1 C6 D4.4684  0
5.7308  0.002
9.4417  0.0042
11.6335 0.0059
13.4882 0.0073
14.247  0.0078
14.6688 0.008
# Z7 ?0 _' g% [4 s# B13.7427 0.0068
11.8059 0.0044
0 B7 ~- Y0 Y9 S- K$ o! m8.6033  0.0015
5.2315  -0.0013
1.6184  -0.0036
$ q- @& {% y, k' c-0.8267 -0.0055
7 t2 m& S) Z* q6 g+ @" q-3.5249 -0.0076
, A/ \" n$ a; {, ]* e; J$ Q. b-5.0424 -0.0089
' D5 X/ r. G4 S; ]/ P+ P' h-5.97   -0.0096
$ H& v9 L6 P- u  I-5.9694 -0.0099
, ~0 W5 o1 N, j-5.5492 -0.0089
0 |; `- U1 K/ o- V-3.0212 -0.0062
-0.4088 -0.0035];
+ C* {3 ~$ X3 a, W# Z# [, h4 k  M %  p0=[1 1 1 1 1 1];
) I2 t; P( Y! o* |* U, \- O' H  p0=[0.005 0.005 0.005 0.005 0.005 0.005];
warning off
F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4)*x(:,1)+p(5)*x(:,2)+p(6);
% 拟合系数，最小二乘方法
  H# R9 c. y7 j$ wp=nlinfit(x,zeros(size(x,1),1),F,p0);
4 G7 i( Y+ T# N% Lp(1)
p(2)
- M1 S& j( P9 T( up(3)
p(4)
p(5)
p(6)
! N* t3 Z8 H; Y5 w1 J, W
0 J% h  ]" ?% q5 }* F+ P9 [
5 x9 h+ @5 X0 b4 W3 t$ ]A=p(1)/p(6);
4 V! P. o, |) G$ u; r$ x) bB=p(2)/p(6);
" }: h4 g( A! i/ `1 TC=p(3)/p(6);
D=p(4)/p(6);
& O- k9 J: `6 }" j) ]E=p(5)/p(6);
. h2 Z! g0 z5 n/ B9 b
  \" f  P; k/ h0 J
%%椭圆中心
6 f) o% v# u: y9 [; r+ K1 QX_center = (B*E-2*C*D)/(4*A*C - B^2);
Y_center = (B*D-2*A*E)/(4*A*C - B^2);
fprintf(' X_center=%g, Y_center=%g\n',X_center,Y_center);
- [4 `- d! ?8 E5 p%%长短轴
a= 2*sqrt((2*A*(X_center^2)+2*C*(Y_center^2)+2*B*X_center*Y_center-2)/(A+C+sqrt(((A-C)^2+B^2))));
b= 2*sqrt((2*A*(X_center^2)+2*C*(Y_center^2)+2*B*X_center*Y_center-2)/(A+C-sqrt(((A-C)^2+B^2))));


%%长轴倾角
* Z# q! n# n# I6 @& N/ bq=0.5 * atan(B/(A-C));
fprintf(' q=%g\n',q);
3 q" H9 O2 Z0 N# i
. G& `3 z( B8 k# r- k
. s# [  C% I! V7 J, `. l% Sfprintf(' a=%g, b=%g\n',a,b);
; ?) P/ E: f: D" ?% Zplot(x(:,1),x(:,2),'ro');

7 S# x; J4 A+ J
hold on;
xmin=min(x(:,1));
7 R7 C1 U& O9 C2 z# [0 }xmax=max(x(:,1));
ymin=min(x(:,2));
ymax=max(x(:,2));
xmid=(xmax+xmin)/2;
ymid=(ymax+ymin)/2;
! V; F  N6 p% lstri=(xmax-xmid)*(ymax-ymid)/2;
4 T  j2 |2 L' _3 DS=a*b*pi;
- R$ r  a2 m; c! a: [3 \* }: \z=S/4/pi/stri;
( L% j$ E" {2 M; s5 p' r% 将数据组集到data
data = [xmin',ymin',xmax',ymax',xmid',ymid',a', b', A', B',C',D',E',stri',S',z'];
% 将data切割成m*n的cell矩阵
[m, n] = size(data);
data_cell = mat2cell(data, ones(m,1), ones(n,1));   
% A" b5 h5 J5 i; x. Y9 n" Q/ @% 添加变量名称
title = {'xmin','ymin','xmax','ymax','xmid','ymid','a', 'b', 'A', 'B','C','D','E','stri','S','z'};
% 将变量名称和数值组集到result
# d$ C7 D$ i! ~7 t" tresult = [title; data_cell];
% 将result写入到wind.xls文件中
s = xlswrite('wind.xls', result);


% 作图
: \* L; c: O7 ^1 M$ P! l2 T%ezplot(@(x,y)F(p,[x,y]),[xmin,xmax,ymin,ymax]);
) e. f* c% q4 t8 G# `+ G2 h" Z* J; H6 L%title('曲线拟合');
5 |( B  _/ o! P1 f& k) q" d2 ?%legend('样本点','拟合曲线')


  t# l; D1 n2 F5 r4 Q, T


3 Z7 c8 t5 I" `1 K/ A5 Y
" }9 b5 Q" n' |3 V  e, c+ k

xixihahaheihei

帮你顶一下