Microsoft MVP성태의 닷넷 이야기
Math: 17. C# - 복소수 타입의 승수를 지원하는 Power 메서드 [링크 복사], [링크+제목 복사],
조회: 27249
글쓴 사람
정성태 (techsharer at outlook.com)
홈페이지
첨부 파일
(연관된 글이 1개 있습니다.)

C# - 복소수 타입의 승수를 지원하는 Power 메서드

Math.PI를 오일러 공식에 대입하면 다음의 식을 얻게 됩니다.

eπi + 1 = 0;

그런데, 저 값을 실제로 코드로 구현하려면 어떻게 해야 할까요? 직관적인 코드로는 다음과 같이 작성해야 하는데,

Complex i = new Complex(0, 1);

double negOne = Math.Pow(Math.E, (Math.PI * i)); // 컴파일 에러!

애석하게도 닷넷 BCL에 기본 제공되는 Math.Pow 메서드는 복소수 타입의 제곱을 지원하지 않기 때문에 컴파일 오류가 발생합니다.

다행히 검색해 보면, Math.NET 라이브러리에 구현된 Complex32 타입은 복소수 타입의 제곱 함수를 제공해 줍니다.

Math.NET Numerics 3.11.1 
; https://www.nuget.org/packages/MathNet.Numerics/

Math.NET은 Nuget으로도 제공하기 때문에 다음과 같이 참조 추가를 해주고,

Install-Package MathNet.Numerics 

Complex32 타입에서 제공하는 Power 메서드를 사용해 복소수 승수를 사용해 주면 됩니다. 이렇게!

Complex32 i = new Complex32(0, 1);
Complex32 e = new Complex32((float)Math.E, 0);
Complex32 pi = new Complex32((float)Math.PI, 0);

Complex32 negOne = e.Power(pi * i);

Console.WriteLine(negOne); // (-1, 1.509958E-07)




Complex32 타입의 소스 코드를 참조하면 당연히 System.Nuemrics의 Complex 타입에도 복소수를 받는 Pow 메서드를 구현하는 것이 가능합니다. 이를 위해 다음과 같은 도우미 함수를 Math.NET으로부터 복사해 가져오고,

static class Helper
{
    public static readonly Complex Zero = new Complex(0, 0);
    public static readonly Complex One = new Complex(1, 0);
    public static readonly Complex NaN = new Complex(double.NaN, double.NaN);

    public static bool IsReal(this Complex complex)
    {
        return (complex.Imaginary == 0.0);
    }

    public static bool IsZero(this Complex complex)
    {
        return ((complex.Real == 0.0) && (complex.Imaginary == 0.0));
    }

    public static bool IsRealNonNegative(this Complex complex)
    {
        return ((complex.Imaginary == 0.0) && (complex.Real >= 0.0));
    }

    public static Complex Power(this Complex src, Complex exponent)
    {
        if (src.IsZero())
        {
            if (exponent.IsZero())
            {
                return One;
            }

            if (exponent.Real > 0f)
            {
                return Zero;
            }

            if (exponent.Real >= 0f)
            {
                return NaN;
            }

            if (exponent.Imaginary != 0f)
            {
                return new Complex(double.PositiveInfinity, double.PositiveInfinity);
            }

            return new Complex(double.PositiveInfinity, 0f);
        }

        Complex complex = exponent * src.NaturalLogarithm();
        return complex.Exponential();
    }

    public static Complex NaturalLogarithm(this Complex src)
    {
        if (src.IsRealNonNegative())
        {
            return new Complex(Math.Log(src.Real), 0);
        }

        return new Complex(0.5 * (Math.Log(src.MagnitudeSquared())), src.Phase());
    }

    public static double MagnitudeSquared(this Complex src)
    {
        return ((src.Real * src.Real) + (src.Imaginary * src.Imaginary));
    }

    public static double Phase(this Complex src)
    {
        if ((src.Imaginary == 0f) && (src.Real < 0f))
        {
            return Math.PI;
        }

        return Math.Atan2(src.Imaginary, src.Real);
    }

    public static Complex Exponential(this Complex src)
    {
        double real = Math.Exp(src.Real);

        if (src.IsReal())
        {
            return new Complex(real, 0f);
        }

        return new Complex(real * (Math.Cos(src.Imaginary)),
                            real * (Math.Sin(src.Imaginary)));
    }
}

System.Numerics.Complex 타입을 이용해 Math.Pow 기능을 다음과 같이 구현할 수 있습니다.

Complex i = new Complex(0, 1);
Complex e = new Complex(Math.E, 0);

Complex negOne = e.Power(Math.PI * i);

Console.WriteLine(negOne); // (-1, 1.22460635382238E-16)

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




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

[연관 글]






[최초 등록일: ]
[최종 수정일: 5/21/2016]

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

비밀번호

댓글 작성자
 




... 46  47  48  49  50  51  52  53  54  55  56  57  [58]  59  60  ...
NoWriterDateCnt.TitleFile(s)
12580정성태3/28/202120625오류 유형: 708. SQL Server Management Studio - Execution Timeout Expired.
12579정성태3/28/202121393오류 유형: 707. 중첩 가상화(Nested Virtualization) - The virtual machine could not be started because this platform does not support nested virtualization.
12578정성태3/27/202122056개발 환경 구성: 560. Docker Desktop for Windows 기반의 Kubernetes 구성 (2) - WSL 2 인스턴스에 kind가 구성한 k8s 서비스 위치
12577정성태3/26/202122758개발 환경 구성: 559. Docker Desktop for Windows 기반의 Kubernetes 구성 - WSL 2 인스턴스에 kind 도구로 k8s 클러스터 구성
12576정성태3/25/202121490개발 환경 구성: 558. Docker Desktop for Windows에서 DockerDesktopVM 기반의 Kubernetes 구성 (2) - k8s 서비스 위치
12575정성태3/24/202120508개발 환경 구성: 557. Docker Desktop for Windows에서 DockerDesktopVM 기반의 Kubernetes 구성 [1]
12574정성태3/23/202125887.NET Framework: 1030. C# Socket의 Close/Shutdown 동작 (동기 모드)
12573정성태3/22/202123075개발 환경 구성: 556. WSL 인스턴스 초기 설정 명령어 [1]
12572정성태3/22/202122240.NET Framework: 1029. C# - GC 호출로 인한 메모리 압축(Compaction)을 확인하는 방법파일 다운로드1
12571정성태3/21/202119389오류 유형: 706. WSL 2 기반으로 "Enable Kubernetes" 활성화 시 초기화 실패 [1]
12570정성태3/19/202126176개발 환경 구성: 555. openssl - CA로부터 인증받은 새로운 인증서를 생성하는 방법
12569정성태3/18/202125754개발 환경 구성: 554. WSL 인스턴스 export/import 방법 및 단축 아이콘 설정 방법
12568정성태3/18/202118283오류 유형: 705. C# 빌드 - Couldn't process file ... due to its being in the Internet or Restricted zone or having the mark of the web on the file.
12567정성태3/17/202121546개발 환경 구성: 553. Docker Desktop for Windows를 위한 k8s 대시보드 활성화 [1]
12566정성태3/17/202121160개발 환경 구성: 552. Kubernetes - kube-apiserver와 REST API 통신하는 방법 (Docker Desktop for Windows 환경)
12565정성태3/17/202117796오류 유형: 704. curl.exe 실행 시 dll not found 오류
12564정성태3/16/202118475VS.NET IDE: 160. 새 프로젝트 창에 C++/CLI 프로젝트 템플릿이 없는 경우
12563정성태3/16/202121411개발 환경 구성: 551. C# - JIRA REST API 사용 정리 (3) jira-oauth-cli 도구를 이용한 키 관리
12562정성태3/15/202121542개발 환경 구성: 550. C# - JIRA REST API 사용 정리 (2) JIRA OAuth 토큰으로 API 사용하는 방법파일 다운로드1
12561정성태3/12/202120663VS.NET IDE: 159. Visual Studio에서 개행(\n, \r) 등의 제어 문자를 치환하는 방법 - 정규 표현식 사용
12560정성태3/11/202122154개발 환경 구성: 549. ssh-keygen으로 생성한 PKCS#1 개인키/공개키 파일을 각각 PKCS8/PEM 형식으로 변환하는 방법
12559정성태3/11/202122792.NET Framework: 1028. 닷넷 5 환경의 Web API에 OpenAPI 적용을 위한 NSwag 또는 Swashbuckle 패키지 사용 [2]파일 다운로드1
12558정성태3/10/202121681Windows: 192. Power Automate Desktop (Preview) 소개 - Bitvise SSH Client 제어 [1]
12557정성태3/10/202119751Windows: 191. 탐색기의 보안 탭에 있는 "Object name" 경로에 LEFT-TO-RIGHT EMBEDDING 제어 문자가 포함되는 문제
12556정성태3/9/202116621오류 유형: 703. PowerShell ISE의 Debug / Toggle Breakpoint 메뉴가 비활성 상태인 경우
12555정성태3/8/202121347Windows: 190. C# - 레지스트리에 등록된 DigitalProductId로부터 라이선스 키(Product Key)를 알아내는 방법파일 다운로드2
... 46  47  48  49  50  51  52  53  54  55  56  57  [58]  59  60  ...