我正在通过移植我在matlab / octave中编写的一些代码来尝试numpy。在matlab中,我可以定义python切片的等价物,然后根据需要增加它。例如,在我的matlab代码中我有
HXx_range = 1:NHXx;
HXy_range = 1:NHXy;
blah blah blah
Hx(HXx_range, HXy_range) = Da(Hx_media(HXx_range, HXy_range)).*Hx(HXx_range, HXy_range) + Db(Hx_media(HXx_range, HXy_range)).*(Ez(HXx_range,HXy_range) - Ez(HXx_range,**HXy_range+1**));
Hy(HYx_range, HYy_range) = Da(Hy_media(HYx_range, HYy_range)).*Hy(HYx_range, HYy_range) + Db(Hy_media(HYx_range, HYy_range)).*(Ez(**HYx_range+1**,HYy_range) - Ez(HYx_range,HYy_range));
Ez(EZx_range, EZy_range) = Ca(Ez_media(EZx_range, EZy_range)).*Ez(EZx_range, EZy_range) + Cb(Ez_media(EZx_range, EZy_range)).*(Hy(EZx_range,EZy_range) - Hy(**EZx_range-1**,EZy_range) + Hx(EZx_range,**EZy_range-1**) - Hx(EZx_range,EZy_range));
'**'中的术语(如'HXy_range + 1')是关键部分; HXy_range + 1等于2:(NHXy + 1)。在python中,我可以用类似的方式定义切片:
HXx_range = slice(0, NHXx)
然而,HXx_range + 1给了我一个错误。当然,我可以为此做一个新的切片,但它不是那么干净。有办法解决这个问题吗?
感谢。
答案 0 :(得分:3)
如果将HXy_range
定义为numpy数组,则可以根据需要增加它。如果您希望将其用作切片,可以形成slice(*HXy_range)
:
In [26]: HXy_range = np.array([1,10])
In [27]: HXy_range+1
Out[27]: array([ 2, 11])
In [28]: slice(*(HXy_range+1))
Out[30]: slice(2, 11, None)
答案 1 :(得分:1)
不,Python slice
个实例是immutable。要使用标准slice
实例,您必须每次都创建一个新实例。不幸的是,you can't subclass slice
either。