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

... 151  152  153  154  [155]  156  157  158  159  160  161  162  163  164  165  ...
NoWriterDateCnt.TitleFile(s)
1268정성태3/27/201237018오류 유형: 148. WCF svc 호출 시 HTTP Error 404.17 - Not Found [1]
1267정성태3/16/201235676.NET Framework: 314. C++의 inline asm 사용을 .NET으로 포팅하는 방법 [1]파일 다운로드1
1266정성태3/14/201238388개발 환경 구성: 149. RAID 1 구성 시 하드 디스크 장애 발생 해결에 대한 경험담
1265정성태3/13/201228713VC++: 61. 아이태니엄(IA64: Itanium) 에서 겪은 C++ 포인터 연산 문제 [2]
1264정성태3/10/201248513.NET Framework: 313. WELL512 난수 발생 알고리즘 - C# [5]파일 다운로드1
1263정성태3/9/201226809개발 환경 구성: 148. tinyget 사용법
1262정성태3/8/201247554개발 환경 구성: 147. .keystore 파일에 저장된 개인키 추출 방법과 인증기관으로부터 온 공개키를 합친 pfx 파일 만드는 방법 [1]
1261정성태3/7/201227845Phone: 2. 개발자용 윈도우 폰 7 기기 등록하는 방법
1260정성태3/6/201228424Phone: 1. 윈도폰 7 개발자 (회사) 등록하는 방법 [3]
1259정성태3/4/201239889Windows: 57. 새로 추가된 네트워크 커널 디버깅 및 PowerShell 3.0 [1]
1258정성태3/3/201241377개발 환경 구성: 146. SQL Server 2012에 포함된 LocalDB 기능 소개 [3]파일 다운로드1
1257정성태3/3/201229965.NET Framework: 312. Native 스레드와 Managed 스레드 개체의 상태 관계 [1]파일 다운로드1
1256정성태3/3/201233023Math: 5. Euler's totient function - C#파일 다운로드1
1255정성태3/3/201235267Math: 4. 소수 판정 및 소인수 분해 소스 코드 - C# [1]파일 다운로드1
1254정성태3/1/201230837Windows: 56. Windows 8 Consumer Preview를 사용해 보고... [1]
1253정성태3/1/201232063VS.NET IDE: 71. Visual Studio 11 Ultimate 베타 설치 [3]
1252정성태3/1/201229342Windows: 55. 윈도우 8 베타 설치 과정 [1]
1251정성태2/27/201233633VC++: 60. C/C++ Native 스레드 콜 스택 덤프를 얻는 공개 라이브러리 [2]파일 다운로드1
1250정성태2/27/201235631VC++: 59. C/C++ 프로젝트 빌드 속도 개선 - UnityBuild를 아세요? [3]
1249정성태2/26/201235876.NET Framework: 311. .NET 스레드 콜 스택 덤프 (5) - ICorDebug 인터페이스 사용법 [2]파일 다운로드3
1248정성태2/25/201246576.NET Framework: 310. C#의 Shift 비트 연산 정리파일 다운로드1
1247정성태2/25/201229346.NET Framework: 309. .NET 응용 프로그램에 기본 생성되는 스레드들에 대한 탐구 [1]파일 다운로드1
1246정성태2/25/201228687개발 환경 구성: 145. 한영 변환은 되지만, 정작 한글 입력이 안되는 경우
1245정성태2/25/201239983개발 환경 구성: 144. 윈도우에서도 유닉스처럼 명령행으로 원격 접속하는 방법
1244정성태2/24/201237010.NET Framework: 308. .NET System.Threading.Thread 개체에서 Native Thread Id를 구할 수 있을까? [1]파일 다운로드1
1243정성태2/23/201236998개발 환경 구성: 143. Visual Studio 2010 - .NET Framework 소스 코드 디버깅 - 두 번째 이야기 [1]
... 151  152  153  154  [155]  156  157  158  159  160  161  162  163  164  165  ...