Microsoft MVP성태의 닷넷 이야기
Math: 50. C# - MathNet.Numerics의 Matrix(행렬) 연산 [링크 복사], [링크+제목 복사],
조회: 30334
글쓴 사람
정성태 (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
정성태

... 106  107  [108]  109  110  111  112  113  114  115  116  117  118  119  120  ...
NoWriterDateCnt.TitleFile(s)
11314정성태9/28/201724892디버깅 기술: 98. windbg - 덤프 파일로부터 닷넷 버전 확인하는 방법
11313정성태9/25/201723360디버깅 기술: 97. windbg - 메모리 덤프로부터 DateTime 형식의 값을 알아내는 방법파일 다운로드1
11312정성태9/25/201726792.NET Framework: 685. C# - 구조체(값 형식)의 필드를 리플렉션을 이용해 값을 바꾸는 방법파일 다운로드1
11311정성태9/20/201718899.NET Framework: 684. System.Diagnostics.Process 객체의 명시적인 해제 권장
11310정성태9/19/201724675.NET Framework: 683. WPF의 Window 객체를 생성했는데 GC 수집 대상이 안 되는 이유 [3]
11309정성태9/13/201721359개발 환경 구성: 335. Octave의 명령 창에서 실행한 결과를 복사하는 방법
11308정성태9/13/201723532VS.NET IDE: 121. 비주얼 스튜디오에서 일부 텍스트 파일을 무조건 메모장으로만 여는 문제파일 다운로드1
11307정성태9/13/201726293오류 유형: 421. System.Runtime.InteropServices.SEHException - 0x80004005
11306정성태9/12/201724392.NET Framework: 682. 아웃룩 사용자를 위한 중국어 스팸 필터 Add-in
11305정성태9/12/201725641개발 환경 구성: 334. 기존 프로젝트를 Visual Studio를 이용해 Github의 신규 생성된 repo에 올리는 방법 [1]
11304정성태9/11/201722269개발 환경 구성: 333. 3ds Max를 Hyper-V VM에서 실행하는 방법
11303정성태9/11/201726373개발 환경 구성: 332. Inno Setup 파일의 관리자 권한을 제거하는 방법
11302정성태9/11/201721906개발 환경 구성: 331. SQL Server Express를 위한 방화벽 설정
11301정성태9/11/201719932오류 유형: 420. SQL Server Express 연결 오류 - A network-related or instance-specific error occurred while establishing a connection to SQL Server.
11300정성태9/10/201725644.NET Framework: 681. dotnet.exe - run, exec, build, restore, publish 차이점 [3]
11299정성태9/9/201724022개발 환경 구성: 330. Hyper-V VM의 Internal Network를 Private 유형으로 만드는 방법
11298정성태9/8/201727607VC++: 119. EnumProcesses / EnumProcessModules API 사용 시 주의점 [1]
11297정성태9/8/201724153디버깅 기술: 96. windbg - 풀 덤프에 포함된 모든 닷넷 모듈을 파일로 저장하는 방법
11296정성태9/8/201725809웹: 36. Edge - "이 웹 사이트는 이전 기술에서 실행되며 Internet Explorer에서만 작동합니다." 끄는 방법
11295정성태9/7/201724641디버깅 기술: 95. Windbg - .foreach 사용법
11294정성태9/4/201723738개발 환경 구성: 329. 마이크로소프트의 CoreCLR 프로파일러 예제 빌드 방법 [1]
11293정성태9/4/201725141개발 환경 구성: 328. Visual Studio(devenv.exe)를 배치 파일(.bat)을 통해 실행하는 방법
11292정성태9/4/201723391오류 유형: 419. Cannot connect to WMI provider - Invalid class [0x80041010]
11291정성태9/3/201723461개발 환경 구성: 327. 아파치 서버 2.4를 위한 mod_aspdotnet 마이그레이션
11290정성태9/3/201728008개발 환경 구성: 326. 아파치 서버에서 ASP.NET을 실행하는 mod_aspdotnet 모듈 [2]
11289정성태9/3/201725728개발 환경 구성: 325. GAC에 어셈블리 등록을 위해 gacutil.exe을 사용하는 경우 주의 사항
... 106  107  [108]  109  110  111  112  113  114  115  116  117  118  119  120  ...