python - f 在 MATLAB 代码中似乎总是正数，但在我的 Python 代码中，f 也得到负值

Question

如何转换following MATLAB代码转为Python？ Here is my solution ，但它并不完全产生相同的结果。例如，f 在 MATLAB 代码中似乎始终为正值，但在我的 Python 代码中，f 也得到负值。

有什么想法可以修复该程序吗？

大多数情况下，我担心这些:
MATLAB:

    for k = 1 : nx  
            j = k+2;   

Python:

for k in range(1,nx+1):  
    j = k+2  

MATLAB:

    [V,D] = eig(A, B);  
    DD = diag(D);  
    keep_idxs = find( ~isinf(DD) );  
    D = diag( DD(keep_idxs) );  
    V = V(:, keep_idxs);  
    [lambda, idx] = min(diag(D));  
    f = V(:,idx);  

Python:

w,vr = scipy.linalg.decomp.eig(A,B)  
w = w.real  
vr = vr.real  
w = w[2:-1-2]  
lambda_ = w.min()  
idx = w.argmin()  
f = vr[:,idx]  

MATLAB:

    f = f(3:end-2);  
    [nf, nf_idx] = max(abs(f)); % L_infty norm  
    n2 = f(nf_idx); % normalize sign away, too  
    f = f ./ n2;  

Python:

f = f[2:-1-1]  
nf = max(np.absolute(f))  
nf_idx = np.absolute(f).argmax()  
nf_idx = np.ma.argmax(f)  
n2 = f[nf_idx]  
f = f/n2  

MATLAB:

    xx = -kappa:h:kappa;  

Python:

xx = np.arange(-kappa, kappa+h, h)  

它们彼此等效吗？如果是，那么为什么它们不产生完全相同的结果？

Answer 1

我不了解 matlab，但对于 python 代码

for k in range(1,nx+1):  
  j = k+2

与

相同

j = nx+2

这不是您的主要问题，但令人担忧。