Bingung dengan implementasi ridge MATLAB

Saya memiliki dua implementasi berbeda ridgedi MATLAB. Hanya satu

1. $\mathbf x = (\mathbf{A}'\mathbf{A}+\mathbf{I}\lambda)^{-1}\mathbf{A}'\mathbf b$

   (seperti yang terlihat di halaman regresi ridge Wikipedia ), dengan$\mathbf{I}$ menjadi matriks identitas kolom ukuran ($\mathbf{A}$) $\times$ kolom ($\mathbf{A}$), dan

2. Saya hanya memanggil "punggungan" Matlab dengan

   x = ridge(A, b, lambda)

Masalah saya adalah keduanya menghasilkan hasil yang berbeda. (1) mengembalikan hasil yang saya inginkan (saya tahu ini dengan membandingkan hasil dengan orang lain) tetapi mengapa (2) tidak mengembalikan hasil yang sama?

Matriks saya $\mathbf A$jarang, diisi dengan 1% 1 dan 99% 0. Beberapa kolom mengandung hampir no 1. Perbedaan terbesar tampaknya adalah bahwa koefisien untuk kolom tersebut dengan sangat sedikit 1 sangat dekat dengan 0 in (1), tetapi bisa sangat jauh dari 0 in (2)

Adakah yang tahu mengapa ini berbeda dan bagaimana saya bisa memodifikasi panggilan di (2) untuk memberikan hasil yang sama seperti (1)?

Ini adalah program matlab untuk memvalidasi apa yang dikatakan kardinal, sebenarnya karena pemusatan dan penskalaan

% Create A(10 by 3 matrix) and b(10 by 1 matrix)
A=rand(10,3);
b=rand(10,1);
lambda=0.01
% centering and scaling A 
s=std(A,0,1);
s=repmat(s,10,1);
A=(A-repmat(mean(A),10,1))./s;

%check the result
X1=inv(A'*A+eye(3)*lambda)*A'*b;
X2=ridge(b,A,lambda,1);

x1 maka sama dengan x2

Anda harus mengkhususkan skala pada ridge sebagai 0, sehingga terlihat seperti x = ridge (A, b, lambda, 0). Dalam hal ini, baris pertama x adalah konstan dan sisanya adalah koefisien. Dengan kata lain, x (2: end, :) harus hasil yang sama seperti yang Anda dapatkan dengan menggunakan (1). Ini jelas dinyatakan dalam dokumentasi mathlab. Semoga ini bisa membantu.