艾滋病疗法的评价及疗效的预测模型

1、艾滋病疗法的评价及疗效的预测模型摘要本文讨论了对艾滋病疗法的评价机制和疗效的预测方法，目的是预测治疗艾滋病的疗效以及终止治疗的最佳时间，从而降低成本并及时寻找到更好的疗法， 以提高人体免疫能力。对于问题一：通过对艾滋病机理分析，我们建立了CD4浓度、HIV浓度、受感染的CD4细胞浓度与时间关系的微分方程组，并利用 matlab求得方程的解 析解，再用非线性拟合法求得系数，最终确定函数。另外我们定义了一个评判指 标Y :人体内CD4浓度和HIV浓度的比值。由于CD4的浓度越高、HIV的浓度 越低反映治疗的效果越好，即丫越大治疗效果越好。根据对 丫的趋势分析，最终 得出结论：针对初始CD4浓度不同

2、的病人，最佳治疗终止时间在 2427周内。对于问题二：采取模式识别中的最大隶属原则，通过引入隶属函数来确定四 个疗法的最终隶属类型，最终结果为：优为第 4疗法；中为第3疗法；差为第1 疗法和第2疗法。并对最优疗法展开进一步的讨论，引入年龄因素影响产生的折 减系数.，采用问题一的方法，最终得出如下结论：最佳治疗终止时间为40周。对于问题三：建立优化模型，利用层次分析法计算CD4生成速率和疗法费用的相对重要性的权数利用这权重来作出判断，可将病人分为两类：当疗法费 用权重小于0.457时，不影响问题二中的评价和预测；当疗法费用权重大于0.457 时，问题二中的最优疗法改为第一种疗法。通过对模型中各个

3、方程的系数进行灵敏度分析，很好的验证了模型的稳定 性。说明我们的模型具有很强的可靠性。关键词：机理分析疗效评判指标最大隶属原则层次分析法灵敏度分析1#问题重述1.1. 基本内容艾滋病是当前人类社会最严重的瘟疫之一， 从 1981年发现以来的 20多年间， 它已经吞噬了近 3000 万人的生命。 HIV 病毒破坏人的免疫系统，使人体丧失抵 抗各种疾病的能力，从而严重危害人的生命。人类免疫系统的 CD4 细胞在抵御 HIV 的入侵中起着重要作用， 当 CD4 被 HIV 感染而裂解时， 其数量会急剧减少， HIV 将迅速增加，导致 AIDS 发作。艾滋病治疗的目的，是尽量减少人体内 HIV 的数量

4、，同时产生更多的 CD4，至少要有效地降低CD4减少的速度，以提高人体免疫能力。迄今为止人 类还没有找到能根治 AIDS 的疗法，目前的一些 AIDS 疗法不仅对人体有副作用， 而且成本也很高。许多国家和医疗组织都在积极试验、寻找更好的 AIDS 疗法。现在题目给出两组数据。ACTG320 （见附件1）是同时服用zidovudine （齐多夫定），lamivud ine （拉美夫定）和in di navir （茚地那韦）3种药物的300多名 病人每隔几周测试的CD4和HIV的浓度（每毫升血液里的数量）。193A （见附 件 2）是将 1300多名病人随机地分为 4组，每组按下述 4种疗法中的一

5、种服药， 大约每隔8周测试的CD4浓度（这组数据缺HIV浓度，它的测试成本很高）。4 种疗法的日用药分别为：600mg zidovudine或400mg didanosine （去羟基苷），这 两种药按月轮换使用；600 mg zidovudine加2.25 mg zalcitabine（扎西他滨）；600 mg zidovudine 力口 400 mg didanosing 600 mg zidovudine 力口 400 mg didanosine 再 加 400 mg nevirapine （奈韦拉平）。1.2. 拟解决的问题（1）利用附件 1 的数据，预测继续治疗的效果，或者确定最佳治

6、疗终止时 间（继续治疗指在测试终止后继续服药， 如果认为继续服药效果不好， 则可选择 提前终止治疗）。（2）利用附件2的数据，评价4种疗法的优劣（仅以CD4为标准），并对 较优的疗法预测继续治疗的效果，或者确定最佳治疗终止时间。（3） 艾滋病药品的主要供给商对不发达国家提供的药品价格如下：600mg zidovudine 1.60美元， 400mg didanosine 0.85美元， 2.25 mg zalcitabine 1.85美元， 400 mg nevirapine 1.20 美元。如果病人需要考虑 4 种疗法的费用，对（ 2）中的 评价和预测（或者提前终止）有什么改变。基本假设1)

7、2).不发生新的机会性感染，原有的机会性感染不再复发。.病人服用药物按照医生的嘱咐，不存在漏服药物的情况 .人体内的资源所产生的健康 CD4 细胞的速率为常数。三、符号说明Ti第i阶段健康CD4细胞的浓度Si第i阶段2030之间CD4增长速率VID具体一个人的CD4增长的平均速率Vi第i种疗法整个组的CD4增长的平均速率第i种疗法的总人数A第i种评价向量Wi第i种特征向量四、模型的建立与求解4.1. 模型一 4. 1.1 问题一分析我们首先对艾滋病进行机理分析，梳理疗效与两种细胞之间的关系，根据它 建立初步模型。我们首先根据所查的资料：CD4细胞计数是衡量机体免疫功能的重要指标。美国已将CD4

8、<200个/ mm3作为艾滋病的一项指标。CD4+细胞的数量和功能的 降低是引起人体免疫功能缺陷的主要原因。CD4+细胞的数量是确定诊断、疾病 分期、制定抗病毒治疗和预防机会性感染方案的标准试验指标。CD4+细胞计数配合病毒载量测定还是预测疾病进程的可靠指标。及CD4计数和HIV相关并发症。通过散点图以及对数据的聚类分析，于是我们再根据题目把CD4分为三个评价阶段：040 个/ mm3 ； 40100个/ mm3 ； 100 个/mm3 以上。还有通过治疗尽可能早地将病毒抑制到检测水平以下 （<50拷贝/mL），以保 存免疫功能，延长患者的生命和改善生活质量。所以我们将 HIV的评

9、价是以50 拷贝/mL为标准。4.1.2 机理分析根据题意以及与医学实际更加符合， 需要考虑病毒产生的过程：HIV病毒感 染CD4细胞，在病毒整合酶催化下与细胞宿主染色体 DNA整合，会复制产生出 新的HIV病毒。为了用数学方法描述这一过程，我们引人健康的 CD4细胞的浓3度变量T,受感染的CD4细胞浓度变量I。假设没有HIV病毒时，T服从下面的 微分方程1 ：5这里s是由人体内的资源所产生的健康 CD4细胞的速率，是健康CD4 细胞的自然死亡率。当HIV病毒活动及有药物治疗时，我们除了考虑到上述的条件之外，还有治疗的效果和健康T细胞和自由病毒在结合时两者都在减少，因此第三个方程出现I项及第一

10、个方程和第三个方程各出现了 -kT项。设细胞的感染率服从双线性函数，其中:-为治疗系数，感染细胞的死亡率为， ki是健康T细胞和自由病毒的结合率,感染细胞的死亡率,k；表示病毒和T细胞结合后转化成感染T细胞的速率，卜I是%是溶解性感染T细胞的死亡率，N是每个溶解性感染T细胞产生病毒的数量于是有初步模型里的三个方程11:巴=s _ ST -k；VTdtd - k；VT - fldt dV = N.bI k；VT ： l vV。dt4.1.3 模型确定我们首先利用Matlab中的dsolve函数来对上面初步模型里的三个函数进行解析式求解，再通过对系数的整理、归纳成非线性指数函数T = c 1 Qe

11、31，再利用Matlab对上面的初步模型进行系数的确定，具体是将所给出的数据按照之 前的三个评价阶段及初步模型的函数分别拟合，最终确定系数。于是我们就用附录一里的数据代人方程中，按照CD4分为三组进行模拟，下图就是我们模拟出来图表： 第一组（CD4初始浓度：040个/mm3）：350口恍I1IIILJI.D610152D253D35ID466D图表1第一组CD4浓度与时间的关系第二组（CD4初始浓度：40100个/mm3）图表2第二组CD4浓度与时间的关系第三组（CD4初始浓度：100个/mm3以上）上面三个图形包含了艾滋病病人整个治疗过程，每个过程都基本符合我们通 过机理分析所得到的函数。再

12、利用Matlab把各个阶段的数据代入非线性函数中进行模拟，确定函数的 未知参数，求得结果如下（具体程序见附录）：第一组：T. = 29.5932e0.0405t第二组：T2 =167.7509 - 99e 加342机第三组：T3 =165.3167 -0.0096e。现在对继续服药的效果进行预测：定义衡量疗效的评判指标丫 y ：人体内CD4浓度和HIV浓度的比值。由于 CD4的浓度越高、HIV的浓度越低反映治疗的效果越好，即丫越大治疗效果越好。 并对疗效指标丫与时间t进行分析：通过T与V的比值，得出Y =at2 bt c ；我们首先利用Matlab中的dsolve函数将按照同上的分组的数据代人

13、上述方 程，得到结果如下：第一组：£ =-0.1462t26.77386.1691 ；模拟得到下图：图表4第一组CD4浓度与HIV浓度的关系图上显示当匕=25时，函数值为最大，表明此时刻治疗效果最佳。另外从图 上可以看出，在25周以后继续进行治疗，疗效反而变差，因此最佳治疗终止时 间为25周。第二组：2Y2 二-0.1054t4.4408t 25.9585模拟得到下图：图表4第二组CD4浓度与HIV浓度的关系显示当t2=27时，函数值为最大。表明此时刻治疗效果最佳。另外从图上可 以看出，在27周以后继续进行治疗，疗效反而变差，因此最佳治疗终止时间为25周。第三组:Y3 =-0.086

14、9t23.4752t 55.8551 ;模拟得到下图：图表4第三组CD4浓度与HIV浓度的关系显示当t3=24时，函数值为最大。表明此时刻治疗效果最佳。另外从图上可以看出，在24周以后继续进行治疗，疗效反而变差，因此最佳治疗终止时间为 24周。纵上所述，针对初始CD4浓度不同的病人，最佳治疗终止时间在2427周 内。42 模型二4.2.1问题二分析对四种疗法我们考虑采用模糊数学中的模式识别的最大隶属原则进行识别， 看其属于哪个评价程度。我们把疗法的评价程度分为优、中、劣三个等级，其中 以每种疗法治疗后CD4的增长速度作为标准，即为积分高低的标准。接下来对 疗法的优劣采用最大隶属度进行综合评价，

15、得到评价矩阵来判断。 再利用模型一中的函数进行预测和评价。4.2.2综合评价四个疗法我们把所给的数据按照各个疗法分组，将一些只进行一次测试的病人分别剔 除掉，再把剩余的按照个人进行统计分析。一个疗法的积分计算方法为整个组 CD4增长速度的平均速度，其中CD4增 长速度是一个人一次测试得到的 CD4减去之前测试的，再除以这次测试的时间 与上次时间的差即二.VidVi- i =1,2,3,4niVID二CD4k 1 -cD4k （k各个人实际测试次数） ty tk得到的具体数据如下表：疗 法亍 1234积只-0.0-0.0-0.00.0分5549417325252273由上表可知最终的总积分都分布

16、在-0.0560.023范围内。我们把疗法已分为优、中、差三类，每类都看作xm-0.10.11上的模糊集，分别记作A、B、C。其各自的隶属函数分别为'0X0.0 5 54 90-0.0 5 5 9-0.0 2 5 49A X22 x 0.45 4 9 -0.0 2 5 490.04 5 1、10.0 4 51x0.1'2(x + 0.5$x 兰0.055492(x+0.5259)2 0.05549<x 兰-0.02549"b(X ) = "21 - 2(x + 0.4549)- 0.02549 v x 兰 0.045100.0451 <x&quo

17、t;.11 2(x +0.5fx 兰-0.05549，1 -2(x +0.5259了-0.05549c-0.02549"c(X )= *0-0.02549c x 兰0.0451i00.0451 vx 兰0.1其中对于在集合范围内的任意x都有 xB XC X =1。现在我们对四种疗法的最终积分进行判断，并指出四种疗法分别属于三类中 的哪一类。由上表的数据我们可以知道：第一种疗法：max'% x1 ,b x1 ,七 x1 k'-max0,0.3951,0.6104-0.6049第二种疗法：max'% x1 , % x2 ,% x2 i； ；= max0,0.487

18、4,0.5126； = 0.5126第三种疗法：max'% X3 ,% X3 , "c X3 i；= max b.3688,0.6322,0； = 0.6322第四种疗法:maxiA x4，丄B x4 丄C x4 f = max1,0,0,1纵上，我们得出的疗法评价结果如下：优：第4疗法；中：第3疗法；差：第1疗法，第2疗法423对较优疗法的研究我们先挑选20、30、40、50四个年龄段的数据，由散点图（见附录）可见：四个图共同的特征是由斜率大小来分 CD4增加的速率。还有图中反映出处 于20到30岁之间的人CD4浓度增长是快的，减少是慢的，由此我们推出年龄 对CD4浓度的大

19、小是有一定影响的再加上资料所查得免疫能力随着年龄的变化 而变化的，其中青春期为发育阶段，青年为成熟阶段，中年期到老年期为衰老阶 段，呈现出一个有波峰的曲线。我们把曲线按青春期、青年、中年、更年期、老 年期五个阶段画出折线图如下：1Q O,0.00.60.40.2012345系列10.910.950.90.81:青春期2:青年3:中年4:更年期5:老年数系减折6图5五个阶段的折减系数图其中我们定义各个阶段的折减向量，即为二0.9,1,0.95,0.9,0.8/接下来求得第四种疗法的相应的指数方程：T二g - C2e X3t。附表二中对采用疗法四治疗的数据部分进行处理，结果如下:11#3j g J

20、GLIJIL a510IS20253035-10图6再根据方程及数据综合，拟合出方程的未知参数，由此得到具体方程：T =3.1869 -0.1424e"7576t。于是我们把上面的折减系数.考虑进来，加入方程中，总结得：T =3.1869 -0.1424e".7576t'。最后用此方程对四个年龄段的数据再次进行拟合，得到下图：下图充分表现出各个年龄段的CD4增长速率趋势基本相同，斜率相同3 22I11 1 -a ia-*i316-3L14-3112-31-3oe-3U05-nn*5101S3025303&40图7红代表50岁 绿代表20岁 黑代表10岁和40

21、岁 蓝代表30岁（由于数据十分接近导致部分数据点重叠）如图所示40周以后疗效趋于一个平稳值，加上理论分析中的免疫学指标， 得出：40周以后治疗，疗效不明显，建议停止用药。模型三4.3.1问题三分析根据题意，问题三目的是建立优化模型，将价格列入对影响疗法的因数之内, 通过层次分析法，确定它们的权重，从而作出评价。4.3.2模型三建立根据题意，问题三是在问题二的基础上加上费用因数作用后的评价及预测。 我们利用层次分析法先计算CD4生成速率和疗法费用的相对重要性的权数及分 析，得到的矩阵如下：|(1/213#通过Matlab对该矩阵进行处理，得到它的特征根为：-max =2，归一化的特征向量为W二0

22、6667 0.33331，再进行一致性指标的检验得 Cl = 0 , RI = 0 ,n 1一致性比率为CR =0 ： 0.1，则其通过一致性指标检验。所以此矩阵评估比较合理。我们再提出评价向量，它是把各个疗法费用与总的费用之比再经过归一化的。于是得到 CD4 生成速率的评价向量Al - '0.161243,0.201226,0.248998,0.38853 和 医疗价 格的评 价向量耳=-0.081731,-0.331731,-0.235577,-0.350962,而评价矩阵W和A的数据代人疗法的评价方程y = WA得到的结果如下表一:疗法序号1234费用1.6/0.853.452.

23、453.65评价结果0. 2407530. 0707220. 2624240. 426102（注：费用单位美元）表格数据明显显示虽然疗法4相对其它疗法稍微贵了点，但大多数病人还是 比较注重疗效的，尽量会去选择它的。于是四种疗法的费用对问题二中的评价和 预测没有什么改变。更进一步预测疗法1和疗法3由于在费用上占的优势，会有病人因此考虑并 使用选择它们。4.3.3进一步研究我们深入考虑到将病人分为普通人和穷人两种。4.3.1中的病人，我们认为为普通人，于是我们接下来对穷人进行讨论。对穷人：CD4生成速率和疗法费用的相对重要性的权数及分析得到的矩阵为Hr!B2 =3，其余步骤同上。归一化的特征向量为

24、 W2 = b1666 0.8334】，3 1 一CR =0 ： 0.1,则其通过一致性指标检验，因而较为合理。它的 CD4 生成速率的评价向量 A50.161243,0.201226,0.248998,0.3885341医疗费用的评价向量 甩 -_0.081731,_0.33173,_0.23557,_0.35092?，而 A 二 A5 !A6j根据疗法的评价方程结果如下表三:疗法序号1234费用1.883.452.453.65评价结果0. 2876410. 2271250. 2535560. 231678表三明显十分符合实际情况，穷人偏重于考虑物美价廉的疗法，经济方面制 约着他们的选择，于

25、是他们只能充分考虑自身情况进而选择疗法。 因而，问题二 中的最优疗法应为第一种疗法。当W3 =0.542 0.458，代入方程y二WA，可得X = 0.05 -0.0429 0.0271 0.04981.此时最优疗法为第一疗法，可见费用权重为0.458是最优疗法选择的阀值五、模型的灵敏度分析我们对模型中的方程进行灵敏度分析：以T29.5932e0.0405t 为例，进行灵敏度分析- 29,5932-0.0405，当变动1%时代人方程计算出，Ti变动：1%。当1变动1%时代人方程计算出，T1变动：1%。其余方程的步骤同上，结论如下：T2中、1、 各个变动一个百分点时代人方程计算出，T2分别变动：

26、1 %、0.5%和 0.5%。T3中、- > 各个变动一个百分点时代人方程计算出，T3分别变动：0.9%、0.1%和 0.01%。丫1中、1、各个变动一个百分点时代人方程计算出，£分别变动：0.12%、0.02%和 0.03%。丫2中、1、 各个变动一个百分点时代人方程计算出，丫2分别变动：0.9%、2%和 1.3%。丫3中、1、 各个变动一个百分点时代人方程计算出，丫3分别变动：0.5%、1.5%和 0.7%。说明我们所求得方程比较稳定。15六、模型的评价及改进1) .本模型在处理大量的问题时采用了统计方法， 对数据的处理结果都是以平均 值在处理，不会产生由于数据的偶然性而造

27、成较大偏差的问题，具有可靠的 稳定性。2) .在处理提问的问题时，我们对每个问题进行分析之后，将问题分类处理，保 证数据的反映问题的真实性。3) .通过灵敏度分析，模型比较稳定当然，我们的模型也存在改进之处：在模型二中，我们的方程只是对治疗在系数方面的体现，缺乏在函数上的考虑。 我们还可以将模型分组更细更优，还要考虑到方程的置信区间。七、参考文献1 马知恩等，传染病动力学的数学模型与研究，北京：科学出版社，2004。2 姜启源等，数学模型(第三版)，北京：高等教育出版社，2003.8。3 徐斌等，抗艾滋病治疗进展，!，2006.9.15。4 徐莲芝，艾滋病临床表现及诊断治疗， 2006.9.1

28、附录1) . T =-(-s*k*c-s*k*b+s*N*k1-a*b*c*C3*exp(-t*a-t*k)-k*b*c*C3*exp(-t*a-t*k)+C3*aA2* c*exp(-t*a-t*k)+2*C3*a*c*k*exp(-t*a-t*k)+C3*"2*c*exp(-t*a-t*k)-C3*aA2*b*exp (-t*a-t*k)-2*C3*a*b*k*exp(-t*a-t*k)-C3*kA2*b*exp(-t*a-t*k)+C3*aA3*exp(-t*a-t*k)+3*C3*aA2*k*exp(-t*a-t*k)+3*C3*kA2*a*exp(-t*a-t*k)+C3*

29、kA3*exp(-t*a-t*k)/(a+ k)/(k*c+k*b-N*k1)I =k1*(s*k*c+s*k*b-s*N*k1+exp(-b*t)*bA2*C2*a+exp(-b*t)*bA2*C2*k+exp(-b*t)*C2 *c*a*b+exp(-b*t)*C2*c*k*b+C3*aA2*b*exp(-t*a-t*k)+2*C3*a*b*k*exp(-t*a-t*k)+C3*kA2*b*exp(-t*a-t*k)+a*b*c*C3*exp(-t*a-t*k)+k*b*c*C3*exp(-t*a-t*k)/b/(a+k)/ (k*c+k*b-N*k1)-(-s*kA2*b*c-bA2*s

30、*kA2+2*s*k*N*k1*b+s*k*N*k1*c-s*NA2*k1A2-cA2*exp(c*t)*C1*k*a*b-cA2*exp(c*t)*C1*kA2*b-bA2*exp(c*t)*C1*k*a*c-bA2*exp(c*t)*C1*kA2*c+exp(c*t)*C1*N*k1*a*b*c+exp(c*t)*C1*N*k1*k*b*c+exp(-b*t)*N*k1*C2*a*b*c+exp(-b*t)*N*k1*C2*k*b*c-C3*a*bA2*k*exp(-t*a-t*k)*c-C3*kA2*bA2*exp(-t*a-t*k)*c+C3*N*k1*a*exp(-t*a-t*k)*

31、b*c+C3*N*k1*k*exp(-t*a-t*k)*b*c+C3*aA2*k*exp(- t*a-t*k)*b*c+2*C3*kA2*a*exp(-t*a-t*k)*b*c+C3*kA3*exp(-t*a-t*k)*b*c)/b/c/(a+k)/ (k*c+k*b-N*k1)T,I,V=dsolve( DT=s-a*T-k*T*V ','DI=k1*T*V+b*I ','DV=N*I-k*T*V+c*V','t')2) .a=0,4,8,25,40;>> b=88,138,142,164,210;>> beta

32、0=5,3,1;>> beta,r,J=nlinfit(a,b,'eq7',beta0)；>> aa,delta=nlpredci('eq7',a,beta,r,J);>> c=0,3,7,25,31,44;>> d=89,122,160,130,162,180;>> beta0=5,3,1;>> beta,r,J=nlinfit(c,d,'eq7',beta0)；>> cc,delta=nlpredci('eq7',c,beta,r,J);>

33、;> e=0,4,8,27,42;>> f=42,186,207,292,192;>> beta,r,J=nlinfit(e,f,'eq7',beta0)>> ee,delta=nlpredci('eq7',e,beta,r,J);>> g=0,4,9,25,43;>> h=61,93,115,226,269;>> beta,r,J=nlinfit(g,h,'eq7',beta0)>> gg,delta=nlpredci('eq7',g,be

34、ta,r,J);>> i=0,4,7,24,41;>> j=44,130,138,133,185;>> beta,r,J=nlinfit(i,j,'eq7',beta0)>> ii,delta=nlpredci('eq7',i,beta,r,J);>> k=0,4,8,24,40;>> l=35,69,108,168,307;>> beta,r,J=nlinfit(k,l,'eq7',beta0)>> kk,delta=nlpredci('eq

35、7',k,beta,r,J);>> m=0,4,8,25,42;>> n=47,103,114,139,192;>> beta,r,J=nlinfit(m,n,'eq7',beta0)；>> mm,delta=nlpredci('eq7',m,beta,r,J);>> o=0,5,9,25,40;>> p=68,114,139,209,209;>> beta,r,J=nlinfit(o,p,'eq7',beta0)>> oo,delta=nlp

36、redci('eq7',o,beta,r,J);>> plot(a,b,'b.',a,aa,'r',c,d,'b.',c,cc,'r',e,f,'b.',e,ee,'r',g,h,'b.',g,gg,'r',i,j,'b.',i,ii,'r',k,l,'b.',k,kk,'r', m,n,'b.',m,mm,'r',o,p,'b.'

37、;,o,oo,'black')；3) .a=0,4,8,25,40;b=178,228,126,171,99;beta0=5,3,1;beta,r,J=nlinfit(a,b,'eq7',beta0);aa,delta=nlpredci('eq7',a,beta,r,J);c=0,3,8,24,40;d=161,220,316,645,451;beta0=5,3,1;beta,r,J=nlinfit(c,d,'eq7',beta0);cc,delta=nlpredci('eq7',c,beta,r,J);e=0,4

38、,7,24,39;f=115,176,252,168,252;beta0=5,3,1;beta,r,J=nlinfit(e,f,'eq7',beta0);ee,delta=nlpredci('eq7',e,beta,r,J);g=0,3,8,24,40;h=161,220,316,645,451;beta0=5,3,1;beta,r,J=nlinfit(g,h,'eq7',beta0);gg,delta=nlpredci('eq7',g,beta,r,J);i=0,4,8,23,39;j=112,231,288,360,385;b

39、eta0=5,3,1;beta,r,J=nlinfit(i,j,'eq7',beta0)ii,delta=nlpredci('eq7',i,beta,r,J);k=0,4,8,24,39;l=123,155,189,226,203;beta0=5,3,1;beta,r,J=nlinfit(k,l,'eq7',beta0);kk,delta=nlpredci('eq7',k,beta,r,J);m=0,5,9,25,40;n=139,123,132,133,189;beta0=5,3,1;beta,r,J=nlinfit(m,n,&

40、#39;eq7',beta0);mm,delta=nlpredci('eq7',m,beta,r,J);o=0,5,9,24,40;p=247,281,351,356,382;beta0=5,3,1;beta,r,J=nlinfit(o,p,'eq7',beta0);oo,delta=nlpredci('eq7',o,beta,r,J);plot(a,b,'b.',a,aa,'r',c,d,'b.',c,cc,'r',e,f,'b.',e,ee,'bl

41、ack',g,h,'b.',g,gg,'r',i,j,'b.',i,ii,'r',k,l,'b.',k, kk,'r',m,n,'b.',m,mm,'r',o,p,'b.',o,oo,'r')4) .a=0,4,9,23,40;>> b=14,62,110,122,320;>> beta0=0,8,2;>> beta0=8,2;>> beta,r,J=nlinfit(a,b,

42、9;eq7',beta0);>> aa,delta=nlpredci('eq7',a,beta,r,J);>> plot(a,b,'b.',d,dd,'r')>> e=0,4,8,24,40;>> f=33,62,87,128,168;>> beta0=8,2;>> beta,r,J=nlinfit(e,f,'eq7',beta0);>> ee,delta=nlpredci('eq7',e,beta,r,J);>>

43、; d=0,4,7,24,40;>> c=10,11,30,137,225;>> beta0=8,2;>> beta,r,J=nlinfit(d,c,'eq7',beta0);>> dd,delta=nlpredci('eq7',d,beta,r,J);>> g=0,4,8,25,42;>> h=12,56,77,137,122;> beta0=8,2;>> beta,r,J=nlinfit(g,h,'eq7',beta0);>> gg,delt

44、a=nlpredci('eq7',g,beta,r,J);>> i=0,4,8,24,40;>> j=3,9,10,65,104;>> beta0=8,2;>> beta,r,J=nlinfit(i,j,'eq7',beta0);>> ii,delta=nlpredci('eq7',i,beta,r,J);>> k=0,5,8,24,46;>> l=21,86,149,117,127;>> beta0=8,2;>> beta,r,J=nli

45、nfit(k,l,'eq7',beta0);>> kk,delta=nlpredci('eq7',k,beta,r,J);>> m=0,3,8,24,40,48;>> n=8,34,67,45,88,99;>> beta0=8,2;>> beta,r,J=nlinfit(m,n,'eq7',beta0);>> mm,delta=nlpredci('eq7',m,beta,r,J);>> o=0,3,8,24,40;>> p=3,36,5

46、5,86,146;>> beta0=8,2;>> beta,r,J=nlinfit(o,p,'eq7',beta0);>> oo,delta=nlpredci('eq7',o,beta,r,J);>>plot(a,b,'b.',a,aa,'r',d,c,'b.',d,dd,'r',e,f,'b.',e,ee,'r',g,h,'b.',g,gg,'r',i,j,'b.',i,

47、ii,'r',k,l,'b.',k,kk,'r' ,m,n,'b.',m,mm,'r',o,p,'b.',o,oo,'black')5) .> u=0,4,8,25,40'>> v=1.86,2.76,45.90,72.94,30.29'>> x=o nes(5,1) u u.A2;>> b,bint,r,rint,stats=regress(v,x)6) .a=0,5,8,24,37;b=3.89,23.89,43.82,2

48、8.54,32.56;c=0,4,7,24;d=1.89,4.07,11.54,68.5;e=0,6,12,27,44;f=1.27,6.05,10,23.7,25.81;g=0,4,8,24,40;h=6.47,16.76,37.83,41.29,46.67;i=0,4,8,25,42;j=2.4,23.33,45.29,80.59,71.76; k=0,4,8,24,40;l=0.53,1.64,1.82,13.54,22.13;m=0,4,8,24,40;n=3.57,45.59,66.15,94.12,37.37;o=0,4,8,24,45;p=4.14,23.33,36.76,68.

49、48,37.21;q=0,4,8,24,39;r=7.02,10,12.86,22.17,45.13;s=0,3,8,24,40;t=0.61,17.14,23.91,40.95,85.88;u=0,4,8,25,40;v=1.86,2.76,45.90,72.94,30.29; w=0,4,8,24,40;x=7.12,35.29,35.29,54.12,88.26;y=0,3,8,23,39;z=5.56,60.59,73.08,76.92,167.64;plot(a,b,'r.',c,d,'r.',e,f,'r.',g,h,'r.&

50、#39;,i,j,'r.',k,l,'r.',m,n,'r.',o,p,'r.',q,r,'r.',s,t,'r.',u,v,'r.',w,x,'r.', y,z,'r.')7) .A=0,4,8,27,42'>> B=7.12,50.27,76.67,51.23,33.1'>> C=ones(5,1) A A.A2;>> b,bint,r,rint,stats=regress(B,C)>>

51、; Z=b(1)+b(2)*A+b (3) *A.A2;a=0,4,8,27,42;>> b=7.12,50.27,76.67,51.23,33.1;>> c=0,4,7,24,39;>> d=13.93,66.67,136.11,118.24,129.41;>> e=0,5,9,24,40;>> f=11.05,111.14,69.23,138.95,196.47;g=0,4,8,22,38;>> h=10.19,103.89,150,117.65,77.5;>> i=0,4,9,23,42;>>

52、 j=16.84,20.79,36.8,75.59,29.78;>> k=0,4,8,24,33;>> l=16.42,40.56,68.15,103.89,49.46;>> plot(A,Z,'black',a,b,'r.',c,d,'r.',e,f,'r.',g,h,'r.',i,j,'r.',k,l,'r.')8) .A=0,3,7,23,40'>> B=28.57,57.27,135,57.27,64.89'&g

53、t;> C=ones(5,1) A A.A2;>> b,bint,r,rint,stats=regress(B,C)>> Z=b(1)+b(2)*A+b(3)*A.A2;>> a=0,4,8,25,40;>> b=32.36,58.46,26.80,42.75,19.8;>> c=0,4,8,26,46;>> d=22.44,88.82,67.65,53.21,35.29;>> e=0,3,8,24,40;>> f=29.27,88,166.32,379.41,265.29;>>

54、g=0,4,8,24,39;>> h=26.23,67.21,45.43,48.57,38.43;>> i=0,3,7,23,40;>> j=28.57,57.27,135,57.27,64.89;>> k=0,4,8,24,39;>> l=21.8,40.33,74.62,79.41,108.82;>> m=0,4,8,24,39;>> n=34.17,81.58,111.18,132.94,106.84;>> plot(A,Z,'black',a,b,'r.',c

55、,d,'r.',e,f,'r.',g,h,'r.',i,j,'r.',k,l,'r.',m,n,'r.')>> g=0,8.2857,16,24,32,40;>> h=3.1869-0.1424*exp(-75.7576.*g);>> c=0,7.7143,15.7143,23.8571,27.8571,34.8571;>> d=3.1869-0.1424*exp(-71.96972.*c);>> a=0,8,16,24,31.8571;>> b=3.1869-0.1424*exp(-60.60608.*a);>> e=0,7.8571,15.8571,23.8571,31.8571;>> f=3.1869-0.1424*exp(-68.1804.*e);>> plot(a,b,'r.