王 珊 尹云旺2
(1.山東科技大學測繪科學與工程學院,山東 青島266590;2.西安大地測繪股份有限公司,陜西 西安710000)
文章開始之前,需要特殊說明一點,在其他編程語言中,也有元胞數(shù)組的概念,但是本文中所闡述的元胞數(shù)組,是基于MATLAB軟件介紹的,它為 MATLAB中的一種數(shù)據(jù)類型[1],相關(guān)特性可能與其他軟件有所不同。
軟件開發(fā)過程中,我們有時需要一次傳遞多維數(shù)據(jù)到一個循環(huán)中,普通數(shù)據(jù)傳遞一次只能傳遞一個數(shù)據(jù),但是使用元胞,我們可以一次傳遞多維數(shù)據(jù)。與結(jié)構(gòu)體通過域名訪問元素值不同,元胞數(shù)組通過索引方式訪問包含在其中的元素值,正是由于元胞數(shù)組的這種特性,使得其與循環(huán)聯(lián)合使用時,可以解決很多軟件開發(fā)過程中遇到的問題。
元胞數(shù)組由元胞組成,是一個能夠存放任何數(shù)據(jù)類型、任意大小數(shù)組的數(shù)據(jù)容器。同一個元胞數(shù)組中各元胞的內(nèi)容可以不同,可以是MATLAB或用戶自定義的數(shù)據(jù)類型[2]。每個數(shù)據(jù)類型可以有任意維數(shù)。能夠混合存儲任意類型和大小的數(shù)組是元胞數(shù)組一個重要功能,另外一個經(jīng)常使用元胞的方式是存儲不同長度的字符串。像MATLAB中其他數(shù)組一樣,元胞數(shù)組必須為矩形形狀,也就是說元胞數(shù)組的任何一維的每行每列都必須相等。了解結(jié)構(gòu)體的人可能知道,結(jié)構(gòu)體和元胞數(shù)組設(shè)計和目的都很相似,都提供了一種存儲不同大小和不同類型數(shù)據(jù)的方式,兩者最大的區(qū)別在于結(jié)構(gòu)體使用域名來訪問存儲在其中的數(shù)據(jù),而元胞數(shù)組則使用數(shù)字索引,所以我們可以利用索引的方法,一次傳遞多維數(shù)據(jù)。
為方便理解本文中的相關(guān)例子,我們先來講解一下元胞數(shù)組的一些基本操作。元胞數(shù)組的操作包括元胞數(shù)組創(chuàng)建、元胞數(shù)組連接、元胞數(shù)組索引、元胞數(shù)組賦值和從元胞數(shù)組返回值等內(nèi)容[3],這里僅對元胞數(shù)組創(chuàng)建、賦值和索引等內(nèi)容做簡要介紹。
2.2.1 創(chuàng)建元胞數(shù)組
在MATLAB中創(chuàng)建元胞數(shù)組和創(chuàng)建其他類型數(shù)組類似,主要的區(qū)別是創(chuàng)建元胞數(shù)組時需要使用花括號將元胞數(shù)組的內(nèi)容或索引括起來,而創(chuàng)建其他類型數(shù)組時使用的為中括號[4]。元胞數(shù)組中的每行元素用逗號或空格符做間隔,各列之間用分號做間隔。在賦值運算符“=”右邊用花括號“{}”來為一個元胞數(shù)組賦值。例如使用命令A(yù)={B C D E}就可以創(chuàng)建一個包含B C D E不同數(shù)據(jù)類型的元胞數(shù)組,可以在一個元胞數(shù)組中嵌套另一個元胞數(shù)組。也可以使用函數(shù)cell來創(chuàng)建元胞數(shù)組,運行cell(m,n)命令就創(chuàng)建了一個大小為m×n的元胞數(shù)組。
2.2.2 元胞數(shù)組索引
當處理元胞數(shù)組時,你可以選擇元胞數(shù)組中整個元胞來進行操作,也可以選擇元胞中的內(nèi)容進行操作,第一種方式是元胞索引,第二種方式是內(nèi)容索引[5]。x=C(s)命令返回元胞數(shù)組C中下標為s的元胞,x=C{s}命令返回元胞數(shù)組C中下標為s的內(nèi)容,假設(shè)一個元胞數(shù)組包含多個元胞,每個元胞又為包含一個或多個元素的數(shù)組,則x=C{s}(t)命令返回元胞數(shù)組中由下標s指定的元胞中的下標為t的數(shù)組元素。
在這里列舉一個簡單的例子來說明上述操作,首先在MATLAB軟件中輸入下面兩條語句:
data(1,:)= {'張 三'49 {58 98.3 [103 72]}};
data(2,:)= {'李 四'25 {60 98.6 [105 75]}};
這兩個語句創(chuàng)建了一個2×3的元胞數(shù)組data,其中每一行中又嵌套了一個元胞數(shù)組,嵌套的元胞數(shù)組中最后一個元素為一個向量。使用a=data(2,1)命令時,命令窗口返回a=’李四’;使用b=data{2,1}命令時,命令窗口返回b=李四;使用c=data{2,3}(1)命令時,命令窗口返回c=[60],使用d=data{2,3}{1}命令時,命令窗口返回d=60,各運算結(jié)果如圖1所示。
使用whos命令來查看各個變量的類型,如圖2所示,變量a、c和data的類型為元胞,變量b的類型為字符型,變量d的類型為數(shù)值型。
在使用MATLAB編程時,我們可能遇到這樣的情況,需要將一個以上的數(shù)據(jù)同時傳遞到一個循環(huán)中。例如,在編程繪制礦區(qū)工作面示意圖時,假設(shè)一個礦區(qū)有三個矩形工作面。
如果直接使用Plot命令繪圖這三個工作面,那么三個工作面首尾則由一條直線連接在一起,如下圖圖3所示:
為了解決上述問題帶來的不便,一種較理想的方法是使用一個循環(huán)每次繪制一個工作面。但是,這就需要一次傳遞給循環(huán)一個工作面的數(shù)據(jù)。使用其他類型數(shù)據(jù)傳遞給循環(huán)時,一次只能傳遞一個數(shù)據(jù)。由于元胞值可以為一個數(shù)組,每個元胞又可以整體看做一個數(shù)據(jù),同時可以使用索引尋訪每個元素值,所以我們可以在循環(huán)中一次傳遞多維數(shù)據(jù),從而實現(xiàn)上述需求。編寫MATLAB程序如下:
[filename,pathname]=uigetfile(...
{'*.txt';'*.xls';'*.xlsx';'*.xlsb';'*.xlsm';...
'*.*'},...
'坐標數(shù)據(jù)導入');
if isequal([filename,pathname],[0,0])
return
else
File=fullfile(pathname,filename);
end %讀取工作面坐標文件
gzmzb=importdata(File);%導入工作面坐標數(shù)據(jù)到變量gzmzb中
x_gzm=gzmzb(:,1);%將變量gzmzb第一列數(shù)據(jù)賦值給x_gzm
y_gzm=gzmzb(:,2);%將變量gzmzb第二列數(shù)據(jù)賦值給y_gzm
k=length(x_gzm)/4;%根據(jù)坐標求工作面數(shù)目
gzm=cell(k,1);%創(chuàng)建和工作面數(shù)目相等的元胞數(shù)組
gzm(:)={zeros(5,2)};%對元胞數(shù)組中每個元胞初始化
jj=1;
for mm=1:length(x_gzm)
if mod(mm,4)==0
gzm{jj}(1,1)=x_gzm(mm-3);
gzm{jj}(2,1)=x_gzm(mm-2);
gzm{jj}(3,1)=x_gzm(mm-1);
gzm{jj}(4,1)=x_gzm(mm);
gzm{jj}(5,1)=x_gzm(mm-3);
gzm{jj}(1,2)=y_gzm(mm-3);
gzm{jj}(2,2)=y_gzm(mm-2);
gzm{jj}(3,2)=y_gzm(mm-1);
gzm{jj}(4,2)=y_gzm(mm);
gzm{jj}(5,2)=y_gzm(mm-3);
jj=jj+1;
end
end %將4個坐標組成一個工作面完整坐標,為了使工作面閉合,%添加第五個坐標,其值和第一個坐標相等
for nn=1:k
plot(gzm{nn}(:,1),gzm{nn}(:,2),'-r','linewidth',2.5);
hold on
end %繪制工作面示意圖
程序中使用cell(k,1)命令創(chuàng)建了一個大小為k×1的元胞數(shù)組,然后使用gzm(:)={zeros(5,2)}命令,將元胞數(shù)組中每個元胞的值初始化為5×2的全零數(shù)組,這樣做可以提高程序運行速度,并防止發(fā)生內(nèi)存溢出錯誤,最后使用循環(huán)將工作面坐標賦值給每個元胞并繪圖。運行該程序如圖4所示:
MATLAB不僅功能強,而且使用方便,易于掌握[6]。使用其他方式也可以實現(xiàn)同樣效果,但是使用元胞數(shù)組傳遞多維數(shù)據(jù)更方便、高效、直觀。元胞數(shù)組設(shè)計的目的就是為了存儲不同類型、不同大小的數(shù)據(jù),并通過索引來尋訪其中的值,所以元胞數(shù)組和循環(huán)的聯(lián)合使用,必定成為傳遞多維數(shù)據(jù)的一種理想方案。