数学实验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，则有:组分i饱和气相压强为pi，活度系数为有，据饱和气相压强与活度系数的计算式(设qij为组分i对组分j的交互作用系数)：将上述式子整合为：将式子变形为：【运算程序】首先编写如下的函数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);endT=XT(n);p=log(P);for i=1:n d(i)=x*Q(i,1:n);dd(i)=x(i)/d(i);endfor 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。改变温度进行多次尝试，可以得到如下结果初值解XT0T[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)=r[p(t)-q(t)](0r1)，以p(t)=[c-D(p(t))]/d=[c-S(q(t))]/d带入,得到关于q(t)的递推方程. 设S(x)=arctan(μx),μ=4.8,d=0.25,r=0.3,以c为可变参数,讨论期望价格q(t)的变化规律,是否有混沌现象出现,并找出前几个分岔点,观察分岔点的极限趋势是否符合Feigenbaum常数揭示的规律。【问题分析】依题意可得：q(t+1)=(1-r)q(t)+r(c-arctan(μq(t)))/d;代入数据得：q(t+1)=0.7q(t)+1.2(c-arctan(4.8q(t)));由此可构建函数并寻找规律。【运算程序】首先编写如下运算程序：function chaos(iter_fun,x0,r,n)kr=0;for