2012-04-19 62 views
4

假设我有一个列矩阵pols包含[theta,rho,z]的向量。这意味着,如果我有9个这样的向量,它将是一个9x3的矩阵。 这是很方便的把他们布置成这样,因为我可以只喂它们中的任何一个功能,如pol2cartMatlab:我可以通过唯一名称引用数组索引吗?

cart3 = pol2cart(pols(3,:)); 

和一定矢量,我可以通过索引找到其组件1,2 ,3:

rho5 = pols(5,2); 

但有时矩阵实际上是另一个更广泛的基质中,并且可以在中间,而不是开始,使得上述可能变成:

rho5 = pols(5,6); 

为了让代码更易读,以防其他人需要维护它,是否有通过唯一名称引用索引?像

rho5 = pols(5).rho; 

它可以早些时候.rho映射定义为具有RHO值的列​​。

我已经冒险将矩阵转换为单元格,然后使用mat2cellcell2struct来排列,但它似乎并不实用。或者,我可以做结构的数组,但后来我失去了做pol2cart(pols)的能力,而是必须做

pol2cart(pols.theta, pols.rho, pols.z); 

如此重复这个问题:我可以映射指数以独特的名字呢?

+1

你的第一行代码是行不通的。 [POL2CART](http://www.mathworks.com/help/techdoc/ref/pol2cart.html)需要2(或3)个输入参数,而不是单个向量参数。您正在使用哪个版本的MATLAB? – gnovice 2012-04-19 16:05:52

+0

其实我只在Octave上测试过它,因为我没有Matlab,但我必须帮助我的朋友使用它。所以如果pol2cart不接受这样的参数,那么第一个参数不成立。尽管如此,能够通过名称来处理索引也是很方便的。 – syockit 2012-04-20 10:39:55

回答

3

对于默认的MATLAB数据类型,不,你不能这样做。但是,您可以使用create your own new data type (i.e. class object)来存储您的数据。在类定义中,您会重载subsref方法来定义下标引用(即使用(),{}.)表现为您的新对象。这在处理对象数组方面可能会相当棘手,但这是可能的。

请注意,您还必须为您希望在新数据类型上使用的所有现有函数创建重载方法。具体而言,您必须为您的对象创建一个pol2cart方法,该方法可以在内部调用内置的pol2cart函数,并将来自对象的适当数据作为参数传递。

...这让我对目前的情况更简单的解决方案。您可以创建一个结构数组(或数组的标量结构)来存储数据,而不是创建一个全新的类对象,只需创建一个专门用于struct data types的新重载函数,这将简化调用语法。

我在另外两个答案herehere中讨论了内置数据类型的重载函数的更多细节。简而言之,您将创建一个名为@struct的文件夹,并将其放置在您的MATLAB path的文件夹中。在此@struct文件夹,你会然后把这个重载函数:

function varargout = pol2cart(polarCoordinates) 
    [varargout{1:nargout}] = pol2cart(polarCoordinates.theta, ... 
            polarCoordinates.rho, ... 
            polarCoordinates.z); 
end 

请注意,这是功能的流线型的版本,没有错误检查输入,等等。现在,让我们做一些样本数据:

pols = rand(9, 3); %# A 2-D array of data 
polStruct = struct('theta', pols(:, 1), ... %# Convert it to a scalar 
        'rho', pols(:, 2), ... %# structure of arrays 
        'z', pols(:, 3)); 

而且如下,你可以访问第五行的rho值:

rho5 = pols(5,2); 
rho5 = polStruct.rho(5); 

如果你想从极性转换为笛卡尔N个坐标,这里是你会怎么做每个数据类型:

[X,Y,Z] = pol2cart(pols(:,1), pols(:,2), pols(:,3)); %# Calls the built-in one 
[X2,Y2,Z2] = pol2cart(polStruct); %# Calls the overloaded one 

你可以检查它们各自产生相同的结果如下:

>> isequal([X Y Z],[X2 Y2 Z2]) 

ans = 

    1 %# True! 
+0

感谢您的详细解答!我认为类是要走的路......尽管恐怕我不能在代码中使用它,因为害怕让代码的所有者惊诧不已,他们不熟悉编程和面向对象的范例。 – syockit 2012-04-20 10:43:15

0

不,你不能这样做。就如此容易。

1

好的,正式答案可能是“否”,如上面的木片给出的。但是,如果你真的想做这样的事情,你可能会使用半黑客。具体而言,您可以定义一个类并重载运算符以实现(几乎)您想要的。 不幸的是,我看到Matlab不允许超载.,所以你必须使用一些其他的操作符。(请参阅下面的编辑)

只是为了给您提供一个想法,下面是一个返回矩阵的第i行M的类。

classdef Test 
    properties 
    M; 
    end 

    methods 
    function this = Test(M) 
     this.M = M; 
    end 

    function res = mpower(this, i) 
     res = this.M(i, :); 
    end 
    end 
end 

它可以这样运行:您自担风险

>> tmp = Test([1 2; 3 4]); 
>> tmp^1 

ans = 

    1  2 

>> tmp^2 

ans = 

    3  4 

使用! :)

编辑:

我错了上面。正如gnovice的答案中所提到的,您实际上可以使用方法subsref为自定义类定义.运算符。

+0

哇,这是生病!很高兴知道我能做到这一点:)虽然也许不适合与其他人一起工作(我认为这就是为什么有些人避免运营商超载?) – syockit 2012-04-20 10:38:32

相关问题