数学实验6选编.docxVIP

实验6-非线性方程求解 实验目的 1.掌握用MATLAB 软件求解非线性方程和方程组的基本用法，并对结果作初步分析。 2.练习用非线性方程和方程组建立实际问题的模型并进行求解。 实验内容 Q6 给定4种物质对应的参数ai,bi,ci和交互作用矩阵Q如下：? a1=18.607?a2=15.841?a3=20.443?a4=19.293?b1=2643.31?b2=2755.64?b3=4628.96?b4=4117.07?c1=239.73? c2=219.16? c3=252.64? c4=227.44? Q?=? ????1.000????0.192????2.169????1.611? ????0.316????1.000????0.477????0.524? ????0.377????0.360????1.000????0.296? ????0.524????0.282????2.065????1.000? 在压强p=760mmHg下，为了形成均相共沸混合物，温度和组分分别是多少？请尽量找出所有可能的解。 【问题分析】 设组分i所占比例为xi，则有: j=1nxj=1 xj≥0 组分i饱和气相压强为pi，活度系数为γi有，据饱和气相压强与活度系数的计算式(设qij为组分i对组分j的交互作用系数)： lnpi=ai-biT+ci lnγi=1-ln(j=1nxjqij)-j=1nxjqijk=1nxkqjk i=1,2,?,n P=γipi i=i,2,?,n 将上述式子整合为： biT+ci+lnj=1nxjqij+j=in(xjqijk=1nxkqij)-1-ai+lnp=0 i=1,2,?,n 将式子变形为： xibiT+ci+lnj=1nxjqij+j=inxjqijk=1nxkqij-1-ai+lnp=0 【运算程序】 首先编写如下的函数M文件： function f=azeofun(XT,n,P,a,b,c,Q) x(n)=1; for i=1:n-1 x(i)=XT(i); x(n)=x(n)-x(i); end T=XT(n); p=log(P); for i=1:n d(i)=x*Q(i,1:n); dd(i)=x(i)/d(i); end for i=1:n f(i)=x(i)*(b(i)/(T+c(i))+log(x*Q(i,1:n))+dd*Q(1:n,i)-a(i)-1+p); end 然后用所给数据编程，作如下计算： n=4; P=760; a=[18.607,15.841,20.443,19.293]; b=[2643.31,2755.64,4628.96,4117.07]; c=[239.73,219.16,252.64,227.44]; Q=[1.0 0.192 2.169 1.611 0.316 1.0 0.477 0.524 0.377 0.360 1.0 0.296 0.524 0.282 2.065 1.0]; XT0=[0.25,0.25,0.25,50]; [XT,Y]=fsolve(@azeofun,XT0,[],n,P,a,b,c,Q) 得到 XT =[0.0000 0.5858 0.4142 71.9657] Y = 1.0e-06 * [-0.0009 -0.0423 0.4428 -0.4700] 即4种物质组成均相共沸混合物时的比例分别为0,58.58%，41.42%，0，温度为71.9657°C。 改变温度进行多次尝试，可以得到如下结果 初值解XT0x1x2x3x4T[0.25,0.25,0.25,50]0.0000 0.5858 0.41420.000071.9657[0.25,0.4,0.1,70]0.00000.78030.00000.219776.9613[0.5,0,0.5,100]0.00000.00000.00001.000097.7712【结果分析】 本题与课本例题非常相似，原理及运算程序均比较简单。在做本题的过程中遇到最大的bug是设立不同的初值有可能得到相同的结果，所以在找不同的结果设立初值时花费了很长的时间。 Q8 假设商品在t时刻的市场价格为p（t），需求函数为 D(p(t))=c-dp(t)（c，d≥0）。而生产方的期望价格为q（t），供应函数为S（q（t））。当供销平衡时S(q(t))=D（q(t)）。 若期望价格与市场价格不符，商品市场不均衡，生产方t+1时期的期望价格将会调整，方式为q(t+1)-q(t)=