- 1、本文档共16页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
一元弱酸pH值的计算机数值求解;一、问题的提出:;代入[HA]的型体公式:
得到一元三次方程
[H+]3+Ka[H+]2-(Kac+Kw)[H+]-KaKw=0
二元弱酸,精确式为四次方程。
三元弱酸,精确式为五次方程。;二、问题的分析:;二分法:;;牛顿法:;;#includeiostream.h
#includeiomanip.h
#includemath.h
doubleFx(double,double,double);
constdoubleKw=1.0E-14;
intmain()
{
doublek,c,r,x0,A,B,PH;
doublepka;
coutpKa=;
cinpka;
;coutKa=endl;
cink;
coutc=;
cinc;
cout精度:;
cinr;
k=pow(10,-pka);
A=1.0E-14.0;
B=1;//确定扫描范围
x0=0.5(A+B);
while(abs(B-A)r)
{
if(f(x0)*f(A))0
A=x0;
else
B=x0;
x0=0.5(A+B);
}
;PH=-log10(B);
cout该一元酸的pH=setprecision(2)
setiosflags(ios::fixed|ios::showpoint)Pendl;
return0;
}
//functiondefinition
doubleFx(doubles,doubleK,doubleC)
{
return(((s+K)*s)*s)-(K*C+Kw)*s-K*Kw);
};牛顿法;coutpKa=;
cinpka;
coutc=;
cinc;
cout精度:;
cinr;
k=pow(10,-pka);
x0=sqrt(k*c);
x1=x0-Fx(x0,k,c)/DFx(x0,k,c);
while(fabs((x1-x0)/x1)r)
{
x0=x1;
x1=x0-Fx(x0,k,c)/DFx(x0,k,c);
}
PH=-log10(x1);
cout该一元酸的pH=setprecision(2)
setiosflags(ios::fixed|ios::showpoint)PHendl;
return0;
}
;//functiondefinition
floatFx(floats,floatK,floatC)
{
return(((s+K)*s)*s-(K*C+Kw)*s-K*Kw);
}
floatDFx(floatt,floatK,floatC)
{
return(3*t+2*K)*t-(K*C+Kw);
};三、讨论:;谢谢
文档评论(0)