게시판 즐겨찾기
편집
드래그 앤 드롭으로
즐겨찾기 아이콘 위치 수정이 가능합니다.
[잡담성글]답없는 Cholesky 쓰끼 by.Matlab
게시물ID : programmer_14231짧은주소 복사하기
작성자 : 니콜캣우먼
추천 : 0
조회수 : 716회
댓글수 : 3개
등록시간 : 2015/11/02 21:29:29


촐 스키를 처음 매트랩으로 구현해 봤습니다.

u_ii 성분과 u_ij 성분으로 변환된 윗삼각행렬 하나와 그 전치행렬이 
반대로 행렬곱 되면 원본놈이 나오는 것으로 알고 있는데요


매트랩에서는 세상에, 너무도 간단하게



>>B=rand(3)
B =

    0.3290    0.1587    0.4877
    0.0324    0.9230    0.0216
    0.8427    0.7151    0.9285

아무렇게나 랜덤한 3바이3 행렬을 하나 만들었구요

운이 좋게도 

>>C=chol(B)
C =

    0.5736    0.2766    0.8502
         0    0.9201   -0.2321
         0         0    0.3896


변환이 바로 가능한 행렬이었네요. C에 촘스키를 넣었습니다. 
>>CC=C.'
CC =

    0.5736         0         0
    0.2766    0.9201         0
    0.8502   -0.2321    0.3896

이렇게 전치행렬로 혹은, 
>>CC=transpose(C)
해도 결과는 같습니다.

이제 이걸
CC*C 하면 B가 나오겠지요

>>CC*C 
ans =

    0.3290    0.1587    0.4877
    0.1587    0.9230    0.0216
    0.4877    0.0216    0.9285


아, 뭔가 나왔네요
확인 해 볼까요?

>>B
B =

    0.3290    0.1587    0.4877
    0.0324    0.9230    0.0216
    0.8427    0.7151    0.9285


.... ...

... .....?????



>>>>CC*C 
ans =

    0.3290    0.1587    0.4877
    0.1587    0.9230    0.0216
    0.4877    0.0216    0.9285
>>>>B
B =

    0.3290    0.1587    0.4877
    0.0324    0.9230    0.0216
    0.8427    0.7151    0.9285

.... ....???

전체 추천리스트 보기
새로운 댓글이 없습니다.
새로운 댓글 확인하기
글쓰기
◀뒤로가기
PC버전
맨위로▲
공지 운영 자료창고 청소년보호