Microsoft MVP성태의 닷넷 이야기
Math: 50. C# - MathNet.Numerics의 Matrix(행렬) 연산 [링크 복사], [링크+제목 복사],
조회: 28384
글쓴 사람
정성태 (techsharer at outlook.com)
홈페이지
첨부 파일
(연관된 글이 4개 있습니다.)

C# - MathNet.Numerics의 Matrix(행렬) 연산

ML.NET의 경우 기본 라이브러리에는 행렬 타입이 없습니다. 그나마 기존에 Infer.NET으로 알려진 Microsoft.ML.Probabilistic에 Matrix 타입이 제공되므로 다음과 같이 추가하면 됩니다.

Install-Package Microsoft.ML.Probabilistic

<PackageReference Include="Microsoft.ML.Probabilistic" Version="0.3.1810.501" />

연산자 재정의가 되어 있으므로 행렬 곱도 자연스럽게 연산할 수 있고,

double[] data1 = Range.Create(1.0, 10.0, 1.0).ToArray();
Matrix mat1 = new Matrix(3, 3, data1);

double[] data2 = Range.Create(1.0, 10.0, 1.0).ToArray();
Matrix mat2 = new Matrix(3, 3, data2);

Matrix mat3 = mat1 * mat2;

Console.WriteLine(mat3);
/*
30  36  42
66  81  96
102 126 150
*/

요소 간 곱셈도 별도의 SetToElementwiseProduct 메서드를 이용하면 됩니다.

mat3.SetToElementwiseProduct(mat1, mat2);
Console.WriteLine(mat3);
/*
1  4  9
16 25 36
49 64 81
*/

하지만, 기능이 많지는 않습니다. 예를 들어 역행렬을 구하는 것조차 아직 미구현이니,

// 예외 발생: System.NotImplementedException: 'The method or operation is not implemented.'
Matrix invMat = mat3.Inverse();

의사 역행렬을 구하는 기능은 껍데기 메서드조차 없는 상태입니다.




이런 빈약한 행렬 처리를 보완하려면 MathNet.Numerics 라이브러리를 함께 연동하는 것이 좋습니다.

Install-Package MathNet.Numerics

<PackageReference Include="MathNet.Numerics" Version="4.7.0" />

어차피 2차원 배열로 행렬이 이뤄지기 때문에 ML.NET의 Matrix와 MathNet의 Matrix를 다음과 같이 서로 변환할 수 있습니다.

double[,] matData = ...;
Matrix mlMatrix = new Matrix(matData);

// ML.NET의 Matrix를 MathNet Matrix로 변환
Matrix<double> mathMatrix = CreateMatrix.DenseOfArray(mlMatrix.ToArray());

// MathNet Matrix를 ML.NET Matrix로 변환
mlMatrix = new Matrix(mathMatrix.ToArray());

따라서 이제부터는 MathNet의 Matrix를 이용해 의사 역행렬, 행렬 곱, 요소 간 행렬 곱 등의 연산을 할 수 있습니다.

Matrix<double> mat1 = CreateMatrix.DenseOfColumnMajor(3, 3, data1).Transpose();
Matrix<double> mat2 = CreateMatrix.DenseOfColumnMajor(3, 3, data2).Transpose();

Matrix<double> mat3 = mat1.Multiply(mat2); // 행렬 곱
               mat3 = mat1 * mat2; // 연산자도 지원
Matrix<double> mat4 = mat1.PointwiseMultiply(mat2); // 요소 간 곱
Matrix<double> invMat = mat4.Inverse(); // 역행렬

Matrix<double> mat5 = CreateMatrix.DenseOfArray(data.ToArray2D());
Matrix<double> pseudoInvMat = mat5.PseudoInverse(); // 의사 역행렬

(첨부 파일은 이 글의 예제 코드를 포함합니다.)




[이 글에 대해서 여러분들과 의견을 공유하고 싶습니다. 틀리거나 미흡한 부분 또는 의문 사항이 있으시면 언제든 댓글 남겨주십시오.]

[연관 글]






[최초 등록일: ]
[최종 수정일: 8/29/2021]

Creative Commons License
이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.
by SeongTae Jeong, mailto:techsharer at outlook.com

비밀번호

댓글 작성자
 



2021-07-09 12시12분
asc-community/AngouriMath
 - Open-source cross-platform symbolic algebra library for C# and F#. One of the most powerful in .NET. Can be used for both production and research purposes.
; https://github.com/asc-community/AngouriMath

nuget - AngouriMath
; https://www.nuget.org/packages/AngouriMath

AngouriMath 1.3 update
; https://habr.com/en/post/565996/

What's new in AngouriMath 1.2?
; https://habr.com/en/post/545436/

------------------------------

Numpy.NET - C#/F# bindings for NumPy - a fundamental library for scientific computing, machine learning and AI
; https://github.com/SciSharp/Numpy.NET
정성태

1  2  3  4  [5]  6  7  8  9  10  11  12  13  14  15  ...
NoWriterDateCnt.TitleFile(s)
13862정성태1/13/20255092Linux: 112. Linux - 데몬을 위한 SELinux 보안 정책 설정
13861정성태1/11/20255304Windows: 276. 명령행에서 원격 서비스를 동기/비동기로 시작/중지
13860정성태1/10/20254783디버깅 기술: 216. WinDbg - 2가지 유형의 식 평가 방법(MASM, C++)
13859정성태1/9/20255699디버깅 기술: 215. Windbg - syscall 이후 실행되는 KiSystemCall64 함수 및 SSDT 디버깅
13858정성태1/8/20255848개발 환경 구성: 738. PowerShell - 원격 호출 시 "powershell.exe"가 아닌 "pwsh.exe" 환경으로 명령어를 실행하는 방법
13857정성태1/7/20256082C/C++: 187. Golang - 콘솔 응용 프로그램을 Linux 데몬 서비스를 지원하도록 변경파일 다운로드1
13856정성태1/6/20254729디버깅 기술: 214. Windbg - syscall 단계까지의 Win32 API 호출 (예: Sleep)
13855정성태12/28/20246561오류 유형: 941. Golang - os.StartProcess() 사용 시 오류 정리
13854정성태12/27/20246467C/C++: 186. Golang - 콘솔 응용 프로그램을 NT 서비스를 지원하도록 변경파일 다운로드1
13853정성태12/26/20245185디버깅 기술: 213. Windbg - swapgs 명령어와 (Ring 0 커널 모드의) FS, GS Segment 레지스터
13852정성태12/25/20246495디버깅 기술: 212. Windbg - (Ring 3 사용자 모드의) FS, GS Segment 레지스터파일 다운로드1
13851정성태12/23/20245494디버깅 기술: 211. Windbg - 커널 모드 디버깅 상태에서 사용자 프로그램을 디버깅하는 방법
13850정성태12/23/20246891오류 유형: 940. "Application Information" 서비스를 중지한 경우, "This file does not have an app associated with it for performing this action."
13849정성태12/20/20246881디버깅 기술: 210. Windbg - 논리(가상) 주소를 Segmentation을 거쳐 선형 주소로 변경
13848정성태12/18/20246294디버깅 기술: 209. Windbg로 알아보는 Prototype PTE파일 다운로드2
13847정성태12/18/20246188오류 유형: 939. golang - 빌드 시 "unknown directive: toolchain" 오류 빌드 시 이런 오류가 발생한다면?
13846정성태12/17/20246901디버깅 기술: 208. Windbg로 알아보는 Trans/Soft PTE와 2가지 Page Fault 유형파일 다운로드1
13845정성태12/16/20245497디버깅 기술: 207. Windbg로 알아보는 PTE (_MMPTE)
13844정성태12/14/20247500디버깅 기술: 206. Windbg로 알아보는 PFN (_MMPFN)파일 다운로드1
13843정성태12/13/20245583오류 유형: 938. Docker container 내에서 빌드 시 error MSB3021: Unable to copy file "..." to "...". Access to the path '...' is denied.
13842정성태12/12/20245773디버깅 기술: 205. Windbg - KPCR, KPRCB
13841정성태12/11/20246502오류 유형: 937. error MSB4044: The "ValidateValidArchitecture" task was not given a value for the required parameter "RemoteTarget"
13840정성태12/11/20245696오류 유형: 936. msbuild - Your project file doesn't list 'win' as a "RuntimeIdentifier"
13839정성태12/11/20246808오류 유형: 936. msbuild - error CS1617: Invalid option '12.0' for /langversion. Use '/langversion:?' to list supported values.
13838정성태12/4/20246412오류 유형: 935. Windbg - Breakpoint 0's offset expression evaluation failed.
13837정성태12/3/20247226디버깅 기술: 204. Windbg - 윈도우 핸들 테이블 (3) - Windows 10 이상인 경우
1  2  3  4  [5]  6  7  8  9  10  11  12  13  14  15  ...