C#中的函数编程

目录

通过函数表示数据

集合

空集

全集

单元素集

其他集合

二元操作

并集

交集

笛卡尔积

差集

对等差分

其他操作

对于那些想要更进一步的人

欧几里得平面

绘制磁盘

绘制水平和垂直半平面

函数

对于那些想要更进一步的人

分形

复数和绘图

牛顿分形

对于那些想要更进一步的人

延迟简介


  • 下载源码 - 121.6 KB

C#中的函数编程_第1张图片

 

通过函数表示数据

S是任何元素abc...(例如,桌子上的书,或者欧几里得平面的点)的集合,并让S'是这些元素的任意子集(例如,桌子上的绿皮书,或者半径为1的圆上以欧几里得平面原点为中心的点)。

集合S'特征函数S'(x)是一个函数,它将truefalseS的每个元素x相关联。

S'(x) = true if x is in S'
S'(x) = false if x is not in S'

S成为桌子上的一套书,让S'成为桌上的绿皮书。让ab是两个绿色的书,让cd是在表中的两个红色的本。然后:

S'(a) = S'(b) = true
S'(c) = S'(d) = false

S是欧几里德平面中的点的集合,并且让S'在半径为1的园中以欧几里得平面(0,0)原点为中心的点的集合(单位圆)。让ab在单位圆的两点,并让c并且d是半径为2的圆上以欧几里得平面原点为中心的点。然后:

S'(a) = S'(b) = true
S'(c) = S'(d) = false

因此,任何集合S'总是可以由其特征函数表示。一个函数,它将一个元素作为参数,并返回true如果该元素在S'中,否则返回false。换句话说,可以通过C#中的谓词来表示集合(抽象数据类型)。

Predicate set;

在接下来的部分中,我们将看到如何通过C#以函数方式表示集合代数中的一些基本集合,然后我们将在集合上定义泛型二进制运算。然后,我们将在欧几里德平面的子集上对数字应用这些操作。集合是抽象数据结构,数字的子集和欧几里得平面的子集是抽象数据结构的表示,最后二元操作是适用于抽象数据结构的任何表示的通用逻辑。

集合

本节通过C#介绍集合代数中一些基本集的表示。

空集

C#中的函数编程_第2张图片

 

E是空集和Empty是它的特征函数。在集合的代数中,E是没有元素的唯一集合。因此,Empty可以定义如下:

Empty(x) = false if x is in E
Empty(x) = false if x is not in E

因此,EC#中的表示可以定义如下:

public static Predicate Empty()
{
    return x => false;
}

在集合的代数中,Empty表示如下:

 

因此,运行以下代码:

Console.WriteLine("\nEmpty set:");
Console.WriteLine("Is 7 in {{}}? {0}", Empty()(7));

结果如下:

 

全集

 

C#中的函数编程_第3张图片

S是一集合和S'是包含所有要素的S的子集,All是其特色函数。在集合的代数中,S'是包含所有元素的完整集合。因此,All可以这样定义:

All(x) = true if x is in S

因此,S'C#中的表示可以定义如下:

public static Predicate All()
{
    return x => true;
}

在集合的代数中,All表示如下:

 

因此,运行以下代码:

Console.WriteLine("Is 7 in the integers set? {0}", All()(7));

结果如下:

 

单元素集

E是单元素集并且Singleton是它的特征函数。在集合的代数中,E也称为单元集合,或者1元组是具有恰好一个元素e的集合。因此,Singleton可以定义如下:

Singleton(x) = true if x is e
Singleton(x) = false if x is not e

因此,EC#中的表示可以定义如下:

public static Predicate Singleton(T e)
{
    return x => e.Equals(x);
}

因此,运行以下代码:

Console.WriteLine("Is 7 in the singleton {{0}}? {0}", Singleton(0)(7));
Console.WriteLine("Is 7 in the singleton {{7}}? {0}", Singleton(7)(7));

结果如下:

其他集合

本节介绍整数集的子集。

偶数

E是一个偶数集合,并且Even是它的特征函数。在数学中,偶数是一个2的倍数。因此,Even可以定义如下:

Even(x) = true if x is a multiple of 2
Even(x) = false if x is not a multiple of 2

因此,EC#中的表示可以定义如下:

Predicate even = i => i % 2 == 0;

因此,运行以下代码:

Console.WriteLine("Is {0} even? {1}", 99, even(99));
Console.WriteLine("Is {0} even? {1}", 998, even(998));

结果如下:

奇数

E是一个奇数集合并且Odd是它的特征函数。在数学中,奇数是一个不是2的倍数的数字。因此,Odd可以定义如下:

Odd(x) = true if x is not a multiple of 2
Odd(x) = false if x is a multiple of 2

因此,EC#中的表示可以定义如下:

Predicate odd = i => i % 2 == 1;

因此,运行以下代码:

Console.WriteLine("Is {0} odd? {1}", 99, odd(99));
Console.WriteLine("Is {0} odd? {1}", 998, odd(998));

结果如下:

3的倍数

E是一个3的倍数的集合兵器MultipleOfThree是它的特征函数。在数学中,3的倍数是可被3整除的数。因此,MultipleOfThree可以定义如下:

MultipleOfThree(x) = true if x is divisible by 3
MultipleOfThree(x) = false if x is not divisible by 3

因此,EC中的表示可以定义如下:

Predicate multipleOfThree = i => i % 3 == 0;

因此,运行以下代码:

Console.WriteLine("Is {0} a multiple of 3? {1}", 99, multipleOfThree(99));
Console.WriteLine("Is {0} a multiple of 3? {1}", 998, multipleOfThree(998));

结果如下:

5的倍数

E5的倍数的集合兵器MultipleOfFive是它的特征函数。在数学中,5的倍数是可被5整除的数。因此,MultipleOfFive可以定义如下:

MultipleOfFive(x) = true if x is divisible by 5
MultipleOfFive(x) = false if x is not divisible by 5

因此,EC中的表示可以定义如下:

Predicate multipleOfFive = i => i % 5 == 0;

因此,运行以下代码:

Console.WriteLine("Is {0} a multiple of 5? {1}", 15, multipleOfFive(15));
Console.WriteLine("Is {0} a multiple of 5? {1}", 998, multipleOfFive(998));

结果如下:

质数

很久以前,当我接触Project Euler问题时,我不得不解决以下问题:

By listing the first six prime numbers: 2, 3, 5, 7, 11, and 13, 
we can see that the 6th prime is 13.
What is the 10 001st prime number?

为了解决这个问题,我首先必须编写一个快速算法来检查给定的数字是否为素数。编写算法后,我编写了一个迭代算法,遍历素数直到找到第10 001个素数。

E是素数的集合并且Prime是它的特征函数。在数学中,素数是大于1的自然数,除了1和自身之外没有正除数。因此,Prime可以定义如下:

Prime(x) = true if x is prime
Prime(x) = false if x is not prime

因此,EC中的表示可以定义如下:

Predicate prime = IsPrime;

IsPrime是检查给定数字是否为素数的方法。

static bool IsPrime(int i)
{
    if (i == 1) return false; // 1 is not prime
    if (i < 4) return true;   // 2 and 3 are primes
    if ((i >> 1) * 2 == i) return false; // multiples of 2 are not prime
    if (i < 9) return true; // 5 and 7 are primes
    if (i % 3 == 0) return false; // multiples of 3 are not primes

    // If a divisor less than or equal to sqrt(i) is found
    // then i is not prime
    int sqrt = (int)Math.Sqrt(i);
    for (int d = 5; d <= sqrt; d += 6)
    {
        if (i % d == 0) return false;
        if (i % (d + 2) == 0) return false;
    }

    // Otherwise i is prime
    return true;
}

因此,运行下面的代码来解决我们的问题:

int p = Primes(prime).Skip(10000).First();
Console.WriteLine("The 10 001st prime number is {0}", p);

其中Primes定义如下:

static IEnumerable  Primes(Predicate prime)
{
    yield return 2;

    int p = 3;
    while (true)
    {
        if (prime(p)) yield return p;
        p += 2;
    }
}

结果如下:

二元操作

本节介绍了从给定集合和操作集合构造新集合的几个基本操作。下面是集合代数中的维恩图

C#中的函数编程_第4张图片

并集

C#中的函数编程_第5张图片

EF两个集合。EF的并集,用E u F表示的是EF的所有元素的集合。

Union成为并操作。因此,可以在C#中实现如下Union操作:

public static Predicate Union(this Predicate e, Predicate f)
{
    return x => e(x) || f(x);
}

如您所见,Union是一组特征函数的扩展函数。所有操作将被定义为集合的特征函数上的扩展函数。从而,运行以下代码:

Console.WriteLine("Is 7 in the union of Even and Odd Integers Set? {0}", Even.Union(Odd)(7));

结果如下:

交集

C#中的函数编程_第6张图片

EF两个集合。EF的交集,由E n F表示的是所有元素的集合,他们是EF的成员。

Intersection成为交叉操作。因此,可以在C#中实现如下Intersection操作:

public static Predicate Intersection(this Predicate e, Predicate f)
{
    return x => e(x) && f(x);
}

如您所见,Intersection是一组特征函数的扩展函数。从而,运行以下代码:

Predicate multiplesOfThreeAndFive = multipleOfThree.Intersection(multipleOfFive);
Console.WriteLine("Is 15 a multiple of 3 and 5? {0}", multiplesOfThreeAndFive(15));
Console.WriteLine("Is 10 a multiple of 3 and 5? {0}", multiplesOfThreeAndFive(10));

结果如下:

笛卡尔积

C#中的函数编程_第7张图片

EF两个集合。EF笛卡儿积,用E × F表示的是所有有序对(e, f) 集合,因此eE的成员,fF的成员。

CartesianProduct成为笛卡尔积操作。因此,CartesianProduct可以在C#中实现如下操作:

public static Func CartesianProduct(this Predicate e, Predicate f)
{
    return (x, y) => e(x) && f(y);
}

如您所见,CartesianProduct是一组特征函数的扩展函数。从而,运行以下代码:

Func cartesianProduct = multipleOfThree.CartesianProduct(multipleOfFive);
Console.WriteLine("Is (9, 15) in MultipleOfThree x MultipleOfFive? {0}", cartesianProduct(9, 15));

结果如下:

差集

C#中的函数编程_第8张图片

EF两个集合。FE差集,用E \ F表示的是所有元素的集合,其中的元素是E的成员但不是F的成员。

Complement成为差集操作。因此,Complement可以在C#中实现如下操作:

public static Predicate Complement(this Predicate e, Predicate f)
{
    return x => e(x) && !f(x);
}

如您所见,Complement是一个集合的特征函数的扩展方法。从而,运行以下代码:

Console.WriteLine("Is 15 in MultipleOfThree \\ MultipleOfFive set? {0}", 
          multipleOfThree.Complement(multipleOfFive)(15));
Console.WriteLine("Is 9 in MultipleOfThree \\ MultipleOfFive set? {0}", 
          multipleOfThree.Complement(multipleOfFive)(9));

结果如下:

对等差分

C#中的函数编程_第9张图片

EF两个集合。EF对等差分,用▲ F表示的是所有元素的集合,其中的元素是EF的成员,但不是EF的交集成员。

SymmetricDifference成为对等差分操作。因此,SymmetricDifference可以在C#中以两种方式实现操作。一个简单的方法是使用并和差操作,如下:

public static Predicate SymmetricDifferenceWithoutXor(this Predicate e, Predicate f)
{
    return Union(e.Complement(f), f.Complement(e));
}

另一种方法是使用XOR二进制操作如下:

public static Predicate SymmetricDifferenceWithXor(this Predicate e, Predicate f)
{
    return x => e(x) ^ f(x);
}

如您所见,SymmetricDifferenceWithoutXorSymmetricDifferenceWithXor是集合的特征函数的扩展方法。从而,运行以下代码:

// SymmetricDifference without XOR
Console.WriteLine("\nSymmetricDifference without XOR:");
Predicate sdWithoutXor = prime.SymmetricDifferenceWithoutXor(even);
Console.WriteLine
("Is 2 in the symetric difference of prime and even Sets? {0}", sdWithoutXor(2));
Console.WriteLine
("Is 4 in the symetric difference of prime and even Sets? {0}", sdWithoutXor(4));
Console.WriteLine
("Is 7 in the symetric difference of prime and even Sets? {0}", sdWithoutXor(7));

// SymmetricDifference with XOR
Console.WriteLine("\nSymmetricDifference with XOR:");
Predicate sdWithXor = prime.SymmetricDifferenceWithXor(even);
Console.WriteLine("Is 2 in the symetric difference of prime and even Sets? {0}", sdWithXor(2));
Console.WriteLine("Is 4 in the symetric difference of prime and even Sets? {0}", sdWithXor(4));
Console.WriteLine("Is 7 in the symetric difference of prime and even Sets? {0}", sdWithXor(7));

结果如下:

C#中的函数编程_第10张图片

其他操作

本节介绍集合上其他有用的二进制操作。

包含

Contains是检查元素是否在集合中的操作。此操作是一个集合的特征函数的扩展函数,它将元素作为参数,如果元素在集合中则返回true,否则返回false

因此,此操作在C#中定义如下:

public static bool Contains(this Predicate e, T x)
{
    return e(x);
}

因此,运行以下代码:

Console.WriteLine("Is 7 in the singleton {{0}}? {0}", Singleton(0).Contains(7));
Console.WriteLine("Is 7 in the singleton {{7}}? {0}", Singleton(7).Contains(7));

结果如下:

Add是将一个元素添加到集合中的操作。此操作是一个集合的特征函数的扩展函数,它将元素作为参数并将其添加到集合中。

因此,此操作在C#中定义如下:

public static Predicate Add(this Predicate s, T e)
{
    return x => x.Equals(e) || s(x);
}

因此,运行以下代码:

Console.WriteLine("Is 7 in {{0, 7}}? {0}", Singleton(0).Add(7)(7));
Console.WriteLine("Is 0 in {{1, 0}}? {0}", Singleton(1).Add(0)(0));
Console.WriteLine("Is 7 in {{19, 0}}? {0}", Singleton(19).Add(0)(7));

结果如下:

C#中的函数编程_第11张图片

删除

Remove是从集合中删除元素的操作。此操作是对集合的特征函数的扩展函数,该函数将元素作为参数并将其从集合中移除。

因此,此操作在C#中定义如下: 

public static Predicate Remove(this Predicate s, T e)
{
    return x => !x.Equals(e) && s(x);
}

因此,运行以下代码:

Console.WriteLine("Is 7 in {{}}? {0}", Singleton(0).Remove(0)(7));
Console.WriteLine("Is 0 in {{}}? {0}", Singleton(7).Remove(7)(0));

结果如下:

对于那些想要更进一步的人

您可以通过函数式编程看到我们在C#中使用集合代数是多么容易。在前面的部分中显示了最基本的定义。但是,如果你想进一步,你可以考虑:

  • 关系集
  • 抽象代数,如单倍体,群,场,环,K-矢量空间等
  • 包含排除原则
  • 罗素的悖论
  • 康托尔的悖论
  • 双向量空间
  • 定理和推论

欧几里得平面

在上一节中,集合的基本概念是在C#中实现的。在本节中,我们将练习在平面点集(欧几里德平面)上实现的概念。

绘制磁盘

C#中的函数编程_第12张图片

磁盘是由圆圈限定的平面的子集。有两种类型的磁盘。封闭的磁盘是包含构成其边界的圆的点的磁盘,而打开的磁盘是不包含构成其边界的圆的点的磁盘。

在本节中,我们将设置特征函数盘,并在WPF应用程序绘制。

要设置特征函数,首先需要一个计算平面中两点之间欧氏距离的函数。该函数实现如下:

public static double EuclidianDistance(Point point1, Point point2)
{
    return Math.Sqrt(Math.Pow(point1.X - point2.X, 2) + Math.Pow(point1.Y - point2.Y, 2));
}

其中Point是在System.Windows命名空间中定义的struct。这个公式是基于毕达哥拉斯定理。

C#中的函数编程_第13张图片

其中c欧几里德距离(point1.X - point2.X)²(point1.Y - point2.Y)²

Disk特征函数的闭盘。在集合的代数中,实数集中的闭盘的定义如下:

其中ab是中心和R半径的坐标。

因此,DiskC#中的实现如下:

public static Predicate Disk(Point center, double radius)
{
    return p => EuclidianDistance(center, p) <= radius;
}

为了查看集合,我决定实现一个在欧几里得平面中绘制集合的函数Draw。我选择了WPF,因此使用System.Windows.Controls.Image作为画布和Bitmap作为上下文。

因此,我通过Draw方法建立了下面说明的欧几里德平面

C#中的函数编程_第14张图片

以下是该方法的实现:

public static void Draw(this Predicate set, Image plan)
{
    Drawing.Bitmap bitmap = new Drawing.Bitmap((int)plan.Width, (int)plan.Height);

    //
    // Graph drawing
    //
    double semiWidth = plan.Width / 2;
    double semiHeight = plan.Height / 2;

    double xMin = -semiWidth;
    double xMax = +semiWidth;
    double yMin = -semiHeight;
    double yMax = +semiHeight;

    for (int x = 0; x < bitmap.Height; x++)
    {
        double xp = xMin + x * (xMax - xMin) / plan.Width;

        for (int y = 0; y < bitmap.Width; y++)
        {
            double yp = yMax - y * (yMax - yMin) / plan.Height;

            if (set(new Point(xp, yp)))
            {
                bitmap.SetPixel(x, y, Drawing.Color.Black);
            }
        }
    }

    plan.Source = Imaging.CreateBitmapSourceFromHBitmap(
                            bitmap.GetHbitmap(),
                            IntPtr.Zero,
                            System.Windows.Int32Rect.Empty,
                            BitmapSizeOptions.FromWidthAndHeight(bitmap.Width, bitmap.Height));

}

Draw方法中,创建具有与欧几里德平面容器相同的宽度和相同高度的bitmap 然后,如果bitmap每个像素点(x,y)都属于set,则将其替换为一个黑点。xMinxMaxyMinyMax是在上方图中所示的欧几里得平面的边界值。

如您所见,Draw是一组点的特征函数的扩展函数。因此,运行以下代码:

Plan.Disk(new Point(0, 0), 20).Draw(plan);

结果如下:

C#中的函数编程_第15张图片

绘制水平和垂直半平面

水平垂直半平面或者是平面分割欧几里得空间的两个子集之一。水平半平面是两个子集的任意一个子集,其中一个平面通过与Y轴垂直的线将欧几里得空间分割成两个子集。垂直半平面是两个子集的任意一个子集,其中一个平面通过与x轴垂直的线将欧几里得空间分割成两个子集。

在本节中,我们将设置水平垂直半平面的特征函数,在WPF应用程序中绘制它们,看看如果我们将它们与磁盘子集组合,我们可以做些什么。

HorizontalHalfPlane水平半平面的特征函数HorizontalHalfPlanEC中的实现如下:

public static Predicate HorizontalHalfPlane(double y, bool lowerThan)
{
    return p => lowerThan ? p.Y <= y : p.Y >= y;
}

因此,运行以下代码:

Plan.HorizontalHalfPlane(0, true).Draw(plan);

结果如下:

C#中的函数编程_第16张图片

VerticalHalfPlane垂直半平面的特征函数VerticalHalfPlanEC中的实现如下:

public static Predicate VerticalHalfPlane(double x, bool lowerThan)
{
    return p => lowerThan ? p.X <= x : p.X >= x;
}

因此,运行以下代码:

Plan.VerticalHalfPlane(0, false).Draw(plan);

结果如下:

C#中的函数编程_第17张图片

在本文的第一部分中,我们在集合上设置了基本的二元操作。因此,通过组合diskhalf-plane的交集,我们可以绘制半磁盘子集。

因此,运行以下示例:

Plan.VerticalHalfPlane(0, false).Intersection(Plan.Disk(new Point(0, 0), 20)).Draw(plan);

结果如下:

C#中的函数编程_第18张图片

函数

本节介绍欧几里德平面上的集合的函数。

转变

C#中的函数编程_第19张图片

Translate是转变在平面上的点的函数。在欧几里德几何中,Translate 是一个将给定点在指定方向上移动恒定距离的函数。因此,C#中的实现如下:

public static Func Translate(double deltax, double deltay)
{
    return p => new Point(p.X + deltax, p.Y + deltay);
}

其中(deltax, deltay)是转变的常量向量。

TranslateSet是转化平面中集合的函数。此函数在C#中简单实现如下:

public static Predicate TranslateSet(this Predicate set, 
double deltax, double deltay)
{
    return x => set(Translate(-deltax, -deltay)(x));
}

TranslateSet是集合上的扩展函数。它以deltax作为参数,deltax是第一个欧几里得维度中的delta距离,deltay是第二个欧几里得维度中的delta距离。如果点Pxy在集合S被转变,则其坐标将变为x'y'=x + delatxy + deltay。因此,点X ' - delatxY' - DELTAY将始终属于集合S。在集合代数中,TranslateSet称为同构,换句话说,所有转变的集合形成转变组T,其与空间本身同构。这解释了函数的主要逻辑。

因此,在我们的WPF应用程序中运行以下代码:

TranslateDiskAnimation();

其中TranslateDiskAnimation描述如下:

private const double Delta = 50;
private double _diskDeltay;
private readonly Predicate _disk = Plan.Disk(new Point(0, -170), 80);

private void TranslateDiskAnimation()
{
    DispatcherTimer diskTimer = new DispatcherTimer { Interval = new TimeSpan(0, 0, 0, 1, 0) };
    diskTimer.Tick += TranslateTimer_Tick;
    diskTimer.Start();
}

private void TranslateTimer_Tick(object sender, EventArgs e)
{
    _diskDeltay = _diskDeltay <= plan.Height ? _diskDeltay + Delta : Delta;
    Predicate translatedDisk = _diskDeltay <= plan.Height ? 
                                      _disk.TranslateSet(0, _diskDeltay) : _disk;
    translatedDisk.Draw(plan);
}

结果如下:

C#中的函数编程_第20张图片

相似扩大

C#中的函数编程_第21张图片

Scale是将任何点M发送到另一个点N的函数,这样段SNSM在同一条线上,但按系数lambda进行缩放。在集合的代数中,Scale表述如下:

因此,C#中的实现如下:

public static Func Scale
(double deltax, double deltay, double lambdax, double lambday)
{
    return p => new Point(lambdax * p.X + deltax, lambday * p.Y + deltay);
}

其中(deltax, deltay)是转变的常量向量,(lambdax, lambday)是?向量。

ScaleSet是对计划中的集合上应用相似扩大的函数。此函数在C#中简单实现如下:

public static Predicate ScaleSet(this Predicate set, 
       double deltax, double deltay, double lambdax, double lambday)
{
    return x => set(Scale(-deltax / lambdax, -deltay / lambday, 1 / lambdax, 1 / lambday)(x));
}

ScaleSet是集合上的扩展函数。它以deltax作为参数,deltax是第一个欧几里得维度中的delta距离deltay是第二个欧几里得维度中的delta距离,以及(lambdax, lambday)是常数因子向量?如果点Pxy)在集合S中通过ScaleSet转换,则其坐标将更改为(x'y'=lambdax*x+delatxlambday*y+deltay)。因此,点((x-delatx/lambdax,(y-deltay/lambday)将始终属于集合S,如果?当然,与向量0不同。在集合代数中,ScaleSet称为同构,换句话说,所有同构的集合形成同构群H,同构于空间本身。这解释了函数的主要逻辑。

因此,在我们的WPF应用程序中运行以下代码:

ScaleDiskAnimation();

其中ScaleDiskAnimation描述如下:

private const double Delta = 50;
private double _lambdaFactor = 1;
private double _diskScaleDeltay;
private readonly Predicate _disk2 = Plan.Disk(new Point(0, -230), 20);

private void ScaleDiskAnimation()
{
    DispatcherTimer scaleTimer = new DispatcherTimer { Interval = new TimeSpan(0, 0, 0, 1, 0) };
    scaleTimer.Tick += ScaleTimer_Tick;
    scaleTimer.Start();
}

private void ScaleTimer_Tick(object sender, EventArgs e)
{
    _diskScaleDeltay = _diskScaleDeltay <= plan.Height ? _diskScaleDeltay + Delta : Delta;
    _lambdaFactor = _diskScaleDeltay <= plan.Height ? _lambdaFactor + 0.5 : 1;
    Predicate scaledDisk = _diskScaleDeltay <= plan.Height
                                        ? _disk2.ScaleSet(0, _diskScaleDeltay, _lambdaFactor, 1)
                                        : _disk2;
    scaledDisk.Draw(plan);
}

结果如下:

C#中的函数编程_第22张图片

旋转

C#中的函数编程_第23张图片

Rotation是以theta角度旋转一点的函数。在矩阵代数中,Rotation表述如下:

其中x'y'是旋转后点的坐标,x'y'的公式如下:

这个公式的演示非常简单。看看这个旋转。

C#中的函数编程_第24张图片

演示如下:

C#中的函数编程_第25张图片

因此,C#中的实现如下:

public static Func Rotate(double theta)
{
    return p => new Point(p.X * Math.Cos(theta) - p.Y * Math.Sin(theta), 
           p.X * Math.Cos(theta) + p.Y * Math.Sin(theta));
}

RotateSet是在平面的集合上应用的带角度旋转的函数?此函数在C#中简单实现如下。

public static Predicate RotateSet(this Predicate set, double theta)
{
    return p => set(Rotate(-theta)(p));
}

RotateSet是集合上的扩展函数。它以theta作为参数,即旋转角度。如果点Pxy)通过集合S中的RotateSet进行变换,则其坐标将变为(x', y') = (x * cos(?) - y * sin(?), x * cos(?) + y * sin(?))因此,点(x' * cos(?) + y' * sin(?), x' * cos(?) - y' * sin(?))将永远属于集合S。在集合的代数中,RotateSet称为同构,换句话说,所有旋转的集合形成旋转组R,其与空间本身同构。这解释了该函数的主要逻辑。

因此,在我们的WPF应用程序中运行以下代码:

RotateHalfPlaneAnimation();

其中RotateHalfPlaneAnimation描述如下:

private double _theta;
private const double TwoPi = 2 * Math.PI;
private const double HalfPi = Math.PI / 2;
private readonly Predicate _halfPlane = Plan.VerticalHalfPlane(220, false);

private void RotateHalfPlaneAnimation()
{
    DispatcherTimer rotateTimer = new DispatcherTimer 
              { Interval = new TimeSpan(0, 0, 0, 1, 0) };
    rotateTimer.Tick += RotateTimer_Tick;
    rotateTimer.Start();
}

private void RotateTimer_Tick(object sender, EventArgs e)
{
    _halfPlane.RotateSet(_theta).Draw(plan);
    _theta += HalfPi;
    _theta = _theta % TwoPi;
}

结果如下:

C#中的函数编程_第26张图片

对于那些想要更进一步的人

很简单,不是吗?对于那些想要更进一步的人,你可以探索这些:

  • 椭圆
  • 三维欧氏空间
  • 椭圆体
  • 抛物面
  • 双曲面
  • 球面谐波
  • 超椭球体
  • 妊神星
  • 同形体
  • 焦平面

分形

C#中的函数编程_第27张图片

分形是具有通常超过其拓扑维度并且可能落在整数之间的分形维数的集合。例如,Mandelbrot集是由一系列复数二次多项式定义的分形:

Pc(z) = z^2 + c

其中c是复数。Mandelbrot分形被定义为所有点的集合c,使得上述序列不逃逸到无穷大。在集合的代数中,这表达如下:

Mandelbrot集如上图所示。

分形(抽象数据类型)总是可以在C#中表示如下:

Func fractal;

复数和绘图

为了能够绘制分形,我需要操纵复数。因此,我使用了Meta.numerics库。我还需要一个实用程序来绘制Bitmap中的复数,因此我使用了类ColorMapClorTriplet类。

牛顿分形

我创建了一个牛顿分形(抽象数据类型表示)P(z) = z^3 - 2*z + 2,可在下面找到。

public static Func NewtonFractal()
{
    return z => z * z * z - 2 * z + 2;
}

为了能够绘制复数,我需要更新Draw函数。因此,我创建了一个使用ColorMapClorTriplet类的Draw函数的重载。下面是C#中的实现。

public static void Draw(this Func fractal, Image plan)
{
    var bitmap = new Bitmap((int) plan.Width, (int) plan.Height);

    const double reMin = -3.0;
    const double reMax = +3.0;
    const double imMin = -3.0;
    const double imMax = +3.0;

    for (int x = 0; x < plan.Width; x++)
    {
        double re = reMin + x*(reMax - reMin)/plan.Width;
        for (int y = 0; y < plan.Height; y++)
        {
            double im = imMax - y*(imMax - imMin)/plan.Height;

            var z = new Complex(re, im);
            Complex fz = fractal(z);

            if (Double.IsInfinity(fz.Re) || Double.IsNaN(fz.Re) || Double.IsInfinity(fz.Im) ||
                Double.IsNaN(fz.Im))
            {
                continue;
            }

            ColorTriplet hsv = ColorMap.ComplexToHsv(fz);

            ColorTriplet rgb = ColorMap.HsvToRgb(hsv);
            var r = (int) Math.Truncate(255.0*rgb.X);
            var g = (int) Math.Truncate(255.0*rgb.Y);
            var b = (int) Math.Truncate(255.0*rgb.Z);
            Color color = Color.FromArgb(r, g, b);

            bitmap.SetPixel(x, y, color);
        }
    }

    plan.Source = Imaging.CreateBitmapSourceFromHBitmap(
        bitmap.GetHbitmap(),
        IntPtr.Zero,
        Int32Rect.Empty,
        BitmapSizeOptions.FromWidthAndHeight(bitmap.Width, bitmap.Height));
}

因此,运行以下代码:

Plan.NewtonFractal().Draw(plan);

结果如下:

C#中的函数编程_第28张图片

对于那些想要更进一步的人

对于那些想要更进一步的人,你可以探索这些:

  • Mandelbrot分形
  • 朱莉娅分形
  • 其他牛顿分形
  • 其他分形

延迟简介

在本节中,我们将看到如何从.NET Framework 3.5版开始创建一个Lazy类型。

延迟评估是一种评估策略,它将表达式的评估延迟到需要它的值,并且还避免重复评估。与其他非严格的评估策略(如按名称调用)相比,共享可以通过指数因子减少某些函数的运行时间。下面列出了延迟评估的好处。

  • 通过避免不必要的计算以及评估复合表达式的错误条件来提高性能
  • 构造潜在无限数据结构的能力:我们可以轻松地创建一个无限的整数集,例如通过一个函数(参见集合部分中素数的例子)
  • 将控制流(结构)定义为抽象而不是基元的能力

我们来看看下面的代码:

public class MyLazy
{
    #region Fields

    private readonly Func _f;
    private bool _hasValue;
    private T _value;

    #endregion

    #region Constructors

    public MyLazy(Func f)
    {
        _f = f;
    }

    #endregion

    #region Operators

    //
    // Use objects of type MyLazy as objects of type T 
    // through implicit keyword
    //
    public static implicit operator T(MyLazy lazy)
    {
        if (!lazy._hasValue)
        {
            lazy._value = lazy._f();
            lazy._hasValue = true;
        }

        return lazy._value;
    }

    #endregion
}

MyLazy是一个包含以下字段的泛型类:

  • _f延迟评估的函数,返回T类型值
  • _valueT类型的值(冻结值)
  • _hasValue:一个布尔值,指示是否已计算该值

为了使用类型MyLazy的对象作为类型T的对象,使用implicit关键字。评估在类型铸造时完成,此操作称为解冻

因此,运行以下代码:

var myLazyRandom = new MyLazy(GetRandomNumber);
double myRandomX = myLazyRandom;
Console.WriteLine("\n Random with MyLazy: {0}", myRandomX);

其中GetRandomNumber返回随机double如下:

static double GetRandomNumber()
{
    Random r = new Random();
    return r.NextDouble();
}

给出以下输出:

.NET Framework 4引入了一个用于延迟评估的类System.Lazy。此类通过属性Value返回值。运行以下代码:

var lazyRandom = new Lazy(GetRandomNumber);
double randomX = lazyRandom;

给出编译错误,因为类型Lazy与类型double不同。

要使用类System.Lazy的值,必须按如下方式使用该属性Value

var lazyRandom = new Lazy(GetRandomNumber);
double randomX = lazyRandom.Value;
Console.WriteLine("\n Random with System.Lazy.Value: {0}", randomX);

输出如下:

.NET Framework 4还为延迟评估推出了ThreadLocalLazyInitializer

你可能感兴趣的:(CSharp.NET,函数编程)