lsqnonneg

음이 아닌 선형 최소제곱 문제 풀기

구문

x = lsqnonneg(C,d)

x = lsqnonneg(C,d,options)

x = lsqnonneg(problem)

[x,resnorm,residual] = lsqnonneg(___)

[x,resnorm,residual,exitflag,output] = lsqnonneg(___)

[x,resnorm,residual,exitflag,output,lambda] = lsqnonneg(___)

설명

다음 형식의 음이 아닌 최소제곱 곡선 피팅 문제 풀기

$$\underset{x}{\min }{\Vert C\cdot x-d\Vert }_2^2,\text{ where }x\ge 0.$$

x = lsqnonneg(C,d)는 x ≥ 0 조건하에 norm(C*x-d)를 최소화하는 벡터 x를 반환합니다. 인수 C와 d는 실수여야 합니다.

x = lsqnonneg(C,d,options)는 구조체 options에 지정된 최적화 옵션을 사용하여 최소화합니다. 이 옵션을 설정하려면 optimset (Optimization Toolbox)을 사용하십시오.

x = lsqnonneg(problem)은 problem의 최솟값을 구합니다. 여기서 problem은 구조체입니다.

[x,resnorm,residual] = lsqnonneg(___)는 이전의 모든 구문에 대해 잔차에 대한 2-노름(Norm)의 제곱 값인 norm(C*x-d)^2을 추가로 반환하고, 잔차 d-C*x를 반환합니다.

[x,resnorm,residual,exitflag,output] = lsqnonneg(___)는 lsqnonneg의 종료 조건을 설명하는 값 exitflag와 최적화 과정에 대한 정보를 포함하는 구조체 output을 추가로 반환합니다.

[x,resnorm,residual,exitflag,output,lambda] = lsqnonneg(___)는 라그랑주 승수 벡터 lambda를 추가로 반환합니다.

예제

음이 아닌 선형 최소제곱

선형 최소제곱 문제에 대해 음이 아닌 해를 계산하고, 그 결과를 비제약 문제의 해와 비교합니다.

문제 $$\min ||Cx-d||$$에 대해 C 행렬과 d 벡터를 준비합니다.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];
 
d = [0.8587
     0.1781
     0.0747
     0.8405];

제약 조건이 있는 해와 제약 조건이 없는 해를 계산합니다.

x = lsqnonneg(C,d)

x = 2×1

         0
    0.6929

xunc = C\d

xunc = 2×1

   -2.5627
    3.1108

x의 모든 요소는 음수가 아니지만, xunc의 일부 요소는 음수입니다.

두 해에 대한 잔차의 노름을 계산합니다.

constrained_norm = norm(C*x - d)

constrained_norm = 0.9118

unconstrained_norm = norm(C*xunc - d)

unconstrained_norm = 0.6674

제약 조건은 잔차 노름을 오직 증가시키만 하므로 제약 조건이 없는 해의 잔차 노름이 더 작습니다.

디폴트가 아닌 옵션을 사용한 음이 아닌 최소제곱

lsqnonneg 완료 시 출력값을 보려면 Display 옵션을 'final'로 설정하십시오.

옵션을 만듭니다.

options = optimset('Display','final');

문제 $$\min ||Cx-d||$$에 대해 C 행렬과 d 벡터를 준비합니다.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];

d = [0.8587
     0.1781
     0.0747
     0.8405];

options 구조체로 lsqnonneg를 호출합니다.

x = lsqnonneg(C,d,options);

Optimization terminated.

음이 아닌 최소제곱에서 잔차 구하기

출력값과 함께 lsqnonneg를 호출하여 해, 잔차 노름, 잔차 벡터를 구합니다.

문제 $$\min ||Cx-d||$$에 대해 C 행렬과 d 벡터를 준비합니다.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];

d = [0.8587
     0.1781
     0.0747
     0.8405];

해와 잔차에 대한 정보를 얻습니다.

 [x,resnorm,residual] = lsqnonneg(C,d)

x = 2×1

         0
    0.6929

resnorm = 0.8315

residual = 4×1

    0.6599
   -0.3119
   -0.3580
    0.4130

반환되는 잔차 노름이 반환되는 잔차 벡터 노름의 제곱인지 확인합니다.

 norm(residual)^2

ans = 0.8315

음이 아닌 최소제곱의 결과 검사하기

lsqnonneg 완료 후 모든 출력 인수를 요청하여 해와 풀이 과정을 검토합니다.

문제 $$\min ||Cx-d||$$에 대해 C 행렬과 d 벡터를 준비합니다.

C = [0.0372    0.2869
     0.6861    0.7071
     0.6233    0.6245
     0.6344    0.6170];

d = [0.8587
     0.1781
     0.0747
     0.8405];

문제를 풀고 모든 출력 인수를 요청합니다.

[x,resnorm,residual,exitflag,output,lambda] = lsqnonneg(C,d)

x = 2×1

         0
    0.6929

resnorm = 0.8315

residual = 4×1

    0.6599
   -0.3119
   -0.3580
    0.4130

exitflag = 1

output = struct with fields:
    iterations: 1
     algorithm: 'active-set'
       message: 'Optimization terminated.'

lambda = 2×1

   -0.1506
   -0.0000

exitflag는 1이며, 이는 올바른 해임을 나타냅니다.

x(1) = 0이고 그에 대응하는 lambda(1) $$\ne$$ 0이며, 이는 올바른 쌍대성을 표시합니다. 마찬가지로, x(2) > 0이고 그에 대응하는 lambda(2) = 0입니다.

입력 인수

C — 선형 승수

실수 행렬

선형 승수로, 실수 행렬로 지정됩니다. 다음 문제에서 변수 C를 나타냅니다.

$$\underset{x}{\min }{\Vert C\cdot x-d\Vert }_2^2,\text{ where }x\ge 0.$$

호환성을 위해, C의 행 개수는 d의 길이와 같아야 합니다.

예: C = [1,2;3,-1;-4,4]

데이터형: double

d — 부가 항

실수형 벡터

부가 항으로, 실수형 벡터로 지정됩니다. 다음 문제에서 변수 d를 나타냅니다.

$$\underset{x}{\min }{\Vert C\cdot x-d\Vert }_2^2,\text{ where }x\ge 0.$$

호환성을 위해, d의 길이는 C의 행 개수와 같아야 합니다.

예: d = [1;-6;5]

데이터형: double

options — 최적화 옵션

optimset 등이 반환하는 구조체

최적화 옵션으로, optimset 등이 반환하는 구조체로 지정됩니다. optimset을 사용하여 options 구조체에서 다음 필드의 값을 설정하거나 변경할 수 있습니다. 자세한 내용은 최적화 옵션 설정하기 항목을 참조하십시오.

Display	표시 수준: 'notify'(디폴트 값)는 함수가 수렴하지 않는 경우에만 출력값을 표시합니다. 'off' 또는 'none'은 출력값을 표시하지 않습니다. 'final'은 최종 출력값만 표시합니다.
TolX	x에 대한 종료 허용오차로, 양의 스칼라입니다. 디폴트 값은 10*eps*norm(C,1)*length(C)입니다. 허용오차와 중지 조건 항목을 참조하십시오.

예: options = optimset('Display','final')

데이터형: struct

problem — 문제 구조체

구조체

문제 구조체로, 다음 필드를 가진 구조체로 지정됩니다.

필드 이름	항목
C	실수 행렬
d	실수형 벡터
solver	'lsqnonneg'
options	optimset에서 반환되는 options 구조체

데이터형: struct

출력 인수

x — 해(Solution)

실수형 벡터

해로, 실수형 벡터로 반환됩니다. x의 길이는 d의 길이와 같습니다.

resnorm — 잔차 노름의 제곱

음이 아닌 스칼라

잔차 노름의 제곱으로, 음이 아닌 스칼라로 반환됩니다. norm(C*x-d)^2과 같습니다.

residual — 잔차

실수형 벡터

잔차로, 실수형 벡터로 반환됩니다. 잔차는 d - C*x입니다.

exitflag — lsqnonneg가 중지된 이유

정수

lsqnonneg가 중지된 이유로, 정수로 반환됩니다.

1	함수가 해 x로 수렴되었습니다.
0	반복 횟수가 options.MaxIter를 초과했습니다.

output — 최적화 과정에 대한 정보

구조체

최적화 과정에 대한 정보로, 다음 필드를 가진 구조체로 반환됩니다.

iterations	수행된 반복 횟수
algorithm	'active-set'
message	종료 메시지

lambda — 라그랑주 승수

실수형 벡터

라그랑주 승수로, 실수형 벡터로 반환됩니다. 라그랑주 승수의 요소는 상보성 조건 x'*lambda = 0을 충족합니다. 즉, x(i)가 0에 가까우면 lambda(i) < 0이고, x(i) > 0인 경우 lambda(i)가 0에 가까워진다는 의미입니다.

알고리즘

lsqnonneg에 사용되는 알고리즘은 [1]에 설명되어 있습니다. 이 알고리즘은 가능한 일련의 기저 벡터로 시작하며 관련 쌍대 벡터(Dual Vector) lambda를 계산합니다. 그런 다음 다른 가능한 후보와 기저를 교환하기 위해 lambda의 최댓값에 해당하는 기저 벡터를 선택합니다. 이 작업은 lambda ≤ 0이 될 때까지 계속됩니다.

대체 기능

앱

최적화 라이브 편집기 작업은 lsqnonneg에 대한 시각적 인터페이스를 제공합니다.

참고 문헌

[1] Lawson, C. L. and R. J. Hanson. Solving Least-Squares Problems. Upper Saddle River, NJ: Prentice Hall. 1974. Chapter 23, p. 161.

확장 기능

C/C++ 코드 생성

MATLAB® Coder™를 사용하여 C 코드나 C++ 코드를 생성할 수 있습니다.

C/C++ 코드 생성의 경우:

• 가변 크기 배열에 대한 지원을 활성화해야 합니다.

• 출력 구조체의 종료 메시지는 변환되지 않습니다.

• 코드 생성 시 이 함수에 대해 희소 행렬 입력값은 지원되지 않습니다.

스레드 기반 환경

MATLAB®의 backgroundPool을 사용해 백그라운드에서 코드를 실행하거나 Parallel Computing Toolbox™의 ThreadPool을 사용해 코드 실행 속도를 높일 수 있습니다.

이 함수는 스레드 기반 환경을 완전히 지원합니다. 자세한 내용은 스레드 기반 환경에서 MATLAB 함수 실행하기 항목을 참조하십시오.

버전 내역

R2006a 이전에 개발됨

참고 항목

mldivide | optimset | 최적화

×

예제 열기

이 예제의 수정된 버전이 있습니다. 사용자가 편집한 내용을 반영하여 이 예제를 여시겠습니까?

아니요, 수정된 버전을 덮어씁니다. 예