基于Mathemtica绘制一元酸滴定碱的滴定曲线（含V--PH曲线方程的推导）

基于Mathemtica绘制一元酸碱滴定曲线（含V–PH曲线方程的推导）

本文主要是关于Mathemtica绘制一元酸碱滴定曲线的介绍，并且用多种方法来绘制滴定曲线

一、酸滴定碱的V-PH方程的推导

酸滴定碱一共有四种情况，并且这四种情况已经被我同意起来了，如果是强酸或者强碱那么a或b是1（另外一种角度来解读的话就是k≥10^5），这是因为如此，我们就进行了统一。

二、基于MMA绘制滴定曲线的方法

1、 ContourPlot函数绘制（推荐）

（1）强酸滴定强碱模型
0.1 MHCl滴定0 0.1 M，20mlNaOH

首先，我们输入如下代码：

Ka = 10^5; Kb = 10^5; Kw = 10^(-14); C1 = 0.1; C2 = 0.2 ; V2 = 20;
a = Ka/(10^(-x) + Ka); b = Kb/(Kw/10^(-x) + Kb); c = 
 Kw/10^(-x) - 10^(-x);
ContourPlot[V == (b*C1*V2 - c*V2)/(a*C1 + c), {V, 0, 40}, {x, 0, 14}]

绘制图片如下：

（2）强酸滴定弱碱模型：

0.1MHCl滴定0.1M，20ml氨水

输入如下代码：

Ka = 10^5; Kb = 10^(-4.75); Kw = 
 10^(-14); C1 = 0.1; C2 = 0.2 ; V2 = 20;
a = Ka/(10^(-PH) + Ka); b = Kb/(Kw/10^(-PH) + Kb); c = 
 Kw/10^(-PH) - 10^(-PH);
ContourPlot[V == (b*C1*V2 - c*V2)/(a*C1 + c), {V, 0, 40}, {PH, 0, 14}]

绘制如下图：

（3）弱酸滴定强碱模型：

0.1M醋酸滴定0.1M，20ml氢氧化钠

输入如下代码：

Ka = 10^(-4.75); Kb = 10^5; Kw = 
 10^(-14); C1 = 0.1; C2 = 0.2 ; V2 = 20;
a = Ka/(10^(-x) + Ka); b = Kb/(Kw/10^(-x) + Kb); c = 
 Kw/10^(-x) - 10^(-x);
ContourPlot[V == (b*C1*V2 - c*V2)/(a*C1 + c), {V, 0, 40}, {x, 0, 14}]

绘制如图：

（4）弱酸滴定弱碱：

0.1M醋酸滴定0.1M，20ml氨水

输入如下代码：

Ka = 10^(-4.75); Kb = 10^(-4.75); Kw = 
 10^(-14); C1 = 0.1; C2 = 0.2 ; V2 = 20;
a = Ka/(10^(-x) + Ka); b = Kb/(Kw/10^(-x) + Kb); c = 
 Kw/10^(-x) - 10^(-x);
ContourPlot[V == (b*C1*V2 - c*V2)/(a*C1 + c), {V, 0, 40}, {x, 0, 14}]

绘制图像如下：