Question

我试图将一些Matlab代码转换为Python。我遇到切片问题。

Matlab代码：

demod_1_b=-1*mod_noisy*2.*sin(2*pi*Fc*t+phi);
y = filter(Hd,demod_1_b);
y2=conv(y,raised)/ConvFac;
%% till this line the length in python and Matlab are same
y2=y2(sa/2:end-(sa/2));
%%%% when i write this line in Python it gives me wrong answer it should come out as   26     but in python it gives me 33 i think i havnt converted it in a rigth way 
demod_3_b=y2(sa/2:sa:end);

Python代码：

demod_1_b=-1*mod_noisy*2*sin((2*pi*Fc*t)+phi)

N=10
Fc=40
Fs=1600
d=firwin(numtaps=N,cutoff=40,nyq=Fs/2)
print(len(d))
Hd=lfilter( d, 1.0, demod_1_b)
y2=(convolve(Hd,raised))/Convfac
print(len(y2))
y2=y2[(sa/2)-1:-sa/2]
print(len(y2))
 # problem starts here
demod_3_b=y2[(sa/2)-1:sa:,]
print(len(demod_3_a))

我只是想问一下，Matlab中的demod_3_b=y2(sa/2:sa:end);和Python中的demod_3_v=y2[(sa/2)-1:sa:,]是一样的吗？

Answer 1

是的，你的索引是错误的。在NumPy中，以下indexing适用：

基本切片语法是i：j：k其中i是起始索引，j是停止索引，k是步（k≠0）。

因此，您在Python中寻找的是：

y2[(sa/2)-1::sa]

与Matlab不同，stepsize是最后一个输入。如果您想处理数组的整个长度，请不要在两个:之间放置任何内容。

Matlab与Python：切片

1 个答案: