최소 제곱 회귀의 범위 제한

Question

A(k1,k2) 데이터 집합에 대한 회귀를 원했지만 회귀를 -K1<k1<K1 및 -K2<k<K2 범위 내로 제한하고 싶습니다. A은 60x60 행렬로 구성된 이미지 단계이며, 크기는 MxN입니다. A의 정규화 된 주파수 영역의 중심으로부터 0.1N/2만큼 최소 제곱 근사를 수행합니다.

A=rand(60); 
[m, n]=size(A); 
[M,N] = meshgrid(1:m,1:n); 
X = [M(:), N(:)]; 
B=regress(A(:), X); %regression will be done on all the values of A, part where adjustment needed 
hat=reshape(X*B,m,n); 

Answer 1

당신은 먼저 회귀 수행 할에 행렬의 부분 집합을 선택할 수 있습니다 : 여기

코드의 일부이다 호기심에서

% generate the full image 
A_full=rand(60); 
[m, n]=size(A_full); 

% select the part you want, 
% it is not very clear to me if this is really the part you want, 
% but I think you will be able to change it to your needs 
A=A_full(floor(m/2-0.1*m/2):ceil(m/2+0.1*m/2), floor(n/2-0.1*n/2):ceil(n/2+0.1*n/2)); 

% perform the regression on the selected part of A (like you did it) 
[m, n]=size(A); 
[M,N] = meshgrid(1:m,1:n); 
X = [M(:), N(:)]; 
B=regress(A(:), X); %regression will be done on all the values of A, part where adjustment needed 
hat=reshape(X*B,m,n);