曲线是圆形的开始和结束两个小数点,而不是更准确的值

14日视图(30天)
这是我的兆瓦:
clc;
关闭所有;
%有输入常数
cP = 0.05;
cN = 0.05;
没有发生= 0.005;
logb0 = 2.88;
logb1 = 1.17;
logb2 = 4.85;
logb3 = 2.58;
logKs1 = -12.6;
pKa1 = -2.143;
pKa2 = -7.205;
pKa3 = -12.34;
pKb1NH3 = -4.753355133;
%有一系列的轴(0 - 14)
pH = 0:0.01:14;
%有计算cH2PO4 cHPO4 cPO4,他们必须为黑色曲线方程计算(他们是最终的变量曲线方程和依赖pH值)
pKa1 Ka1 =权力(10日);
pKa2 Ka2 =权力(10日);
pKa3卡=权力(10日);
cH2PO4 = cP。* Ka1。*(权力(ph)) ^(2) /((电力(ph)。^ (3)) + Ka1。*(电力(ph)。^ (2)) + Ka1。* Ka2。*(电力(ph)。^ (1) + Ka1。* Ka2。*卡);
cHPO4 = cP。* * Ka2 Ka1。*电源(10,ph)。/((电力(ph)。^ (3)) + Ka1。*(电力(ph)。^ (2)) + Ka1。* Ka2。*(电力(ph)。^ (1) + Ka1。* Ka2。*卡);
cPO4 = cP。* * Ka2 Ka1。* Ka3. /((电力(ph)。^ (3)) + Ka1。*(电力(ph)。^ (2)) + Ka1。* Ka2。*(电力(ph)。^ (1) + Ka1。* Ka2。*卡);
%计算cNH4,它有黑色曲线方程计算(这是最终的曲线方程的变量也取决于pH值)
pOH = 14-pH;
pKb1NH3 Kb1 =权力(10日);
cNH4 = cN。* Kb1. /(权力(poh) + Kb1);
%有一些偏黑色曲线的计算方程
logb0 b0 =权力(10日);
logb1 b1 =权力(10日);
logb2 b2 =权力(10日);
logb3 b3 =权力(10日);
logKs1 Ks1 =权力(10日);
pH-14寇=权力(10日);
%这是黑色的曲线方程
logcNP = Ks1. / (cPO4 * cNH4)。* (1 + b1。* cH2PO4 + b0。* cHPO4 + b2。* cPO4 + b3 *正);
%有一个黑色的横截面曲线计算红色虚线(康联值为0.005的代码)
Index_NP =找到(diff(签署(logcNP -康联)));
i = 1:元素个数(Index_NP);
IndexRange_NP = max (1, Index_NP (i) 1): min(元素个数(pH) Index_NP (i) + 1);
Crossection_NP (i) = interp1 (logcNP (IndexRange_NP), pH值(IndexRange_NP),没有发生);
结束
%以下边界consitions红色曲线(不是冲)
pH_NP = (Crossection_NP (1) < = pH&pH < = Crossection_NP (2));
%这是最后的情节曲线(黑),crossections(青色点),红色虚线和洋红色的曲线
logcNP semilogy (pH值,“- - -”,“颜色”,“k”,“MarkerFaceColor”,' w ',“MarkerEdgeColor”,“米”);
盒子;
持有;
semilogy (pH值(pH_NP) logcNP (pH_NP),“- - -”,“颜色”,“米”,“MarkerFaceColor”,' w ',“MarkerEdgeColor”,“米”);
盒子;
持有;
情节([0 Crossection_NP(1)),(没有发生发生),“——”,“颜色”,“米”,“线宽”,0.5);
盒子;
持有;
情节(Crossection_NP(1),没有发生,“度”,“MarkerFaceColor”,' w ',“MarkerSize”5);
盒子;
持有;
情节(Crossection_NP(2),都没有发生,“度”,“MarkerFaceColor”,' w ',“MarkerSize”5);
盒子;
持有;
轴([0 14 0.000000000000000001 - 1000000000000000000]);
包含(“pH值(-)”);
ylabel (“c{(摩尔/ dm ^ {3})}”);
广场;
网格;
;
盒子;
持有;
兆瓦导致全尺寸图:
如果你放大为更详细的视图的红色曲线可以看出:开始和结束
似乎开始和结束点的红色曲线在两个最近的小数点,而不是圆形的更精确的数据存储在变量Crossection_NP(1)和Crossection_NP (2)。为什么开始和结束的红色曲线的两个最近的小数在策划和如何固定?
我不得不提醒,红色曲线应该开始和结束在青色点。
谢谢你的回答。

答案(0)

类别

找到更多的在计算几何帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!