网站大量收购独家精品文档,联系QQ:2885784924

一元弱酸pH值的计算机数值求解.PPT

  1. 1、本文档共16页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 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)

hongbing2021 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档