我用Mathematica画的碱滴定酸的分布曲线（长期更新）

本文主要记录Mathemtica画酸碱滴定曲线的图像。
一、 一元强碱滴定一元酸的V-PH图像
方法一:数据法
首先我们看推导过程:

我们得到的V-H+的关系式，接下来我们Mathemtica里面的表格功能拿到数据，在利用Origin画图。

1. 一元强碱滴定一元强酸
在这里我以氢氧化钠滴定滴定盐酸为例子，并且酸碱浓度均为0.1M,酸起始体积为20ml（下面f-PH图条件也是如此，不在进行讲述）
我们先上代码：

c = 0.1; C2 = 0.1; V2 = 20; Kw = 10^(-14);
b = Kw/10^(-PH) - 10^(-PH);
Table[{PH, (C2*V2 + V2*b)/(0.1 - b)}, {PH, 0, 12.5, 0.1}] // TableForm
TableForm

接下来我们运行，拿到数据：
记住我们在Mathemtica里面复制数据时，选择文本格式复制：

我们将拿到的数据，放在Origin当中去，下面这个是数据：

接下来，我们画图即可：

这就是我们画出来的图像（当然这个过程对图像进行处理，比如曲线加粗）

** 2. 一元强碱滴定一元弱酸**
在这里我以氢氧化钠滴定醋酸为例子,酸碱浓度均为0.1，酸体积为20ml
我们先上代码：

c = 0.1; C2 = 0.1; V2 = 20; Ka = 10^(-4.75);
Kw = 10^(-14);
a = Ka/(10^(-PH) + Ka);
b = Kw/10^(-PH) - 10^(-PH);
Table[{PH, (a*C2*V2 + V2*b)/(0.1 - b)}, {PH, 0, 12.5, 
   0.1}] // TableForm
TableForm

用Origi做图与上面操作是一样的，我们看图片：

这就是利用Mathemtica画一元强碱滴定一元酸的V-PH图的方法了。

3 总结

根据我们上面的推断，强酸和弱酸的区别在于a，而我们知道Ka=10^5，我们就可以认为是强酸了，所以最后我们可以统一一个模型，该模型的代码如下：

c = 0.1; C2 = 0.1; V2 = 20; Ka = 10^(-4.75);
Kw = 10^(-14);
a = Ka/(10^(-PH) + Ka);
b = Kw/10^(-PH) - 10^(-PH);
Table[{PH, (a*C2*V2 + V2*b)/(0.1 - b)}, {PH, 0, 12.5, 
   0.1}] // TableForm
TableForm

PS：代码的第一排，我们需要自己输入，如果是强酸Ka=10^5,。

方法二、解方程得到反函数法

方法三、数组法画图

二、一元强碱滴定一元酸的f-PH图

首先，我们看推导过程：

方法一：数据法

1、一元强碱滴定一元强酸

代码如下：

c = 0.1; C2 = 0.1; V2 = 20; Kw = 10^(-14);
b = Kw/10^(-PH) - 10^(-PH);
Table[{PH, (c*(C2*V2 + V2*b))/((c - b) C2*V2)}, {PH, 0, 14, 
   0.1}] // TableForm
TableForm

拿到数据之后到Origin里面画图，图像如下：

2、一元强碱滴定一元弱酸

c = 0.1; C2 = 0.1; V2 = 20; Ka = 10^(-4.75);
Kw = 10^(-14);
a = Ka/(10^(-PH) + Ka);
b = Kw/10^(-PH) - 10^(-PH);
Table[{PH, (c*(a*C2*V2 + V2*b))/((c - b) C2*V2)}, {PH, 0, 14, 
   0.1}] // TableForm
TableForm

拿到数据之后，我们去Origin画图，图像如下：