ones(33,1)表示矩阵33行1列的1矩阵，即每个元素均为1
GNP0'表示矩阵GNP0的转置
X=[ones(33,1),GNP0']表示把两个矩阵拼成一个矩阵，并赋值给X
参考代码：
GNP0=1:3;X=[ones(3,1), GNP0'];disp(X);
输出
1 1
1 2
1 3
