实验,线性代数方程组的数值解法非线性方程组求解..docxVIP

数学实验郭明钊 2012011880 化21实验5种群繁殖问题分析：这道题和课本上的例题比较像，不同点在于多了收获量还有限制条件，即各个年龄的种群数量保持不变，(K=1,2…,n)。这道题的思路就是像例题那样列出线性代数方程组然后解之。模型的建立：若已知，且给定收获量。根据和可以得到以下结果： …… 稍作变换可得：……于是可以得到模型如下： L=当年各个年龄段的数量（向量）为：再令 于是就有： 对于第二问，运用matlab实现：clcL=[-1 0 5 3 0;0.4 -1 0 0 0;0 0.6 -1 0 0;0 0 0.6 -1 0;0 0 0 0.4 -1];%构造L矩阵H=[0; 500; 400; 200; 100];x=L\H;x1=vpa(x,5) %求解当前时段各年龄组的数量，保留5位有效数字。 得到结果如下：x1 = 8481.0 2892.4 1335.4 601.27 140.51 于是可以知道分别是8481.0 2892.4 1335.4 601.3 140.5对于第三问，只和第二问有一点差别，就是均是500.用matlab实现时和第二问的方法一样。clcL=[-1 0 5 3 0;0.4 -1 0 0 0;0 0.6 -1 0 0;0 0 0.6 -1 0;0 0 0 0.4 -1];%构造L矩阵H=[0; 500; 500; 500; 500];x=L\H;x1=vpa(x,5) %求解当前时段各年龄组的数量，保留5位有效数字。得到的结果如下：x1 = 10981.0 3892.4 1835.4 601.27 -259.49于是分别是10981.0 3892.4 1835.4 601.3 -259.5 时可以实现各个年龄段的收获量都为500. 但是从结果看出，x5为-259.49，但种群数量不可能为负数，在本题所给条件下，无法使h1~h5均为500。问题总结：这道题用matlab实现时十分的容易，就是解基本的线性代数方程组，这道题的关键就是如何正确地建立模型，只要找准数量之间的关系就能很好的建立模型。实验6均相共沸混合物的组分1、问题分析：所谓共沸混合物，是指有两种或两种以上物质组成的液体混合物，在某种压力下被蒸馏或局部汽化时，在气体状态下和在液体状态下保持相同的组分。课本中已经建立了模型如下：组分xi（i=1，…, n）应满足：?（1）在压强P不大的情况下，稳定条件可以表示，其中是组分i的饱和汽相压强，与温度T有关，由确定，其中为常数。是组分i的液相活度系数，由?确定，其中qij表示组分i与组分j的交互作用参数，qij构成交互作用矩阵Q，Q不一定是对称矩阵。从而经过对数变换，稳定条件可以改写成?（2）则给定组成均相共沸混合物的n种物质，参数ai, bi, ci和交互作用矩阵Q是可以通过实验得到的，可以作为已知系数。在一定的压强P下，模型（1），(2) 描述了确定均相共沸混合物的组分xi的条件。且（1），（2）是含有个未知数非线性方程组，可以利用MATLAB优化工具包的fsolve求解，注意到（1）式是一个简单的线性等式，可以从中消去1个未知数，这通常会使解的效果更好。例如，我们可以从（1）中解出将它带入（2）式，得到含有n个未知数XT=（）的非线性方程组。2、matlab实现，根据题目中所给的数据可以matlab编写如下：首先编写如下的函数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=[1/4,1/4,1/4,70];[XT,Y]=fsolve(@azeofun,XT0,[],n,P,a,b,c,Q)得到XT =[0.0000 0.5858 0.4142 71.9657]Y =