2006年数学建模B题题目即答案
B题目: 艾滋病疗法的评价及疗效的预测
艾滋病是当前人类社会最严重的瘟疫之一,从1981年发现以来的20多年间,它已经吞噬了近3000万人的生命。
艾滋病的医学全名为“获得性免疫缺损综合症”,英文简称AIDS,它是由艾滋病毒(医学全名为“人体免疫缺损病毒”, 英文简称HIV)引起的。这种病毒破坏人的免疫系统,使人体丧失抵抗各种疾病的能力,从而严重危害人的生命。人类免疫系统的CD4细胞在抵御HIV的入侵中起着重要作用,当CD4被HIV感染而裂解时,其数量会急剧减少,HIV将迅速增加,导致AIDS发作。
艾滋病治疗的目的,是尽量减少人体内HIV的数量,同时产生更多的CD4,至少要有效地降低CD4减少的速度,以提高人体免疫能力。
迄今为止人类还没有找到能根治AIDS的疗法,目前的一些AIDS疗法不仅对人体有副作用,而且成本也很高。许多国家和医疗组织都在积极试验、寻找更好的AIDS疗法。
现在得到了美国艾滋病医疗试验机构ACTG公布的两组数据。 ACTG320(见附件1)是同时服用zidovudine(齐多夫定),lamivudine(拉美夫定)和indinavir(茚地那韦)3种药物的300多名病人每隔几周测试的CD4和HIV的浓度(每毫升血液里的数量)。193A(见附件2)是将1300多名病人随机地分为4组,每组按下述4种疗法中的一种服药,大约每隔8周测试的CD4浓度(这组数据缺HIV浓度,它的测试成本很高)。4种疗法的日用药分别为:600mg zidovudine或400mg didanosine(去羟基苷),这两种药按月轮换使用;600 mg zidovudine加2.25 mg zalcitabine(扎西他滨);600 mg zidovudine
加400 mg didanosine;600 mg zidovudine加400 mg didanosine,再加400 mg nevirapine
(奈韦拉平)。
请你完成以下问题:
(1)利用附件1的数据,预测继续治疗的效果,或者确定最佳治疗终止时间(继续治疗指在测试终止后继续服药,如果认为继续服药效果不好,则可选择提前终止治疗)。
(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)中的评价和预测(或者提前终止)有什么改变。
艾滋病疗法的评价及疗效的预测
摘要:人类最宝贵的生命时刻受着疾病的威胁,新的治疗
的发现给人们的健康增加了一份保障。对一种治疗方案的疗效评价与预测更是对这种疗法的发展前景起着关键的决定作用。疗效不好甚至有严重负面作用的药物就应当及时制止应用于临床治疗,以免给人类造成重大损失。
本文以艾滋病的几种疗法为研究对象,对其做出疗法的评价及疗效的预测。首先,我们对检测患病人群所得数据进行整理,按实际情况抽取其中的有用数据,然后根据这些数据模拟血液中CD4、HIV的含量关于时间t的函数关系。并根据所得函数进行疗效预测,得到治疗方案的优劣以及最佳终止治疗时间的判定。
然后,对较优的艾滋病疗法采用灰色系统理论,建立GM(1,1)预测模型,来确定较优方案继续治疗的疗效预测。考虑到实际中人体健康状况随年龄不同的差异,我们将患病人群按年龄的大小分为若干组来进行预测,并将预测结果拟合成光滑曲线,从图象观察走势。
最后,将价格因素纳入到考虑范围,综合考虑价格与疗效的结合。类比于普通商品的性价比,我们可以将某种治疗方案疗效的量化函数与治疗费用的比值作为指标函数,可看作该治疗
的“性价比”,通过这一指标的高低来判定相应治疗方案的优劣。
关键字: 曲线拟合 灰色系统 灰色预测 疗效评价与预测 指标函数
1
一、 问题重述
艾滋病是当前人类社会最严重的瘟疫之一,从1981年发现以来的20多年间,它已经吞噬了近3000万人的生命。
艾滋病的医学全名为“获得性免疫缺损综合症”,英文简称AIDS,它是由艾滋病毒(医学全名为“人体免疫缺损病毒”, 英文简称HIV)引起的。这种病毒破坏人的免疫系统,使人体丧失抵抗各种疾病的能力,从而严重危害人的生命。人类免疫系统的CD4细胞在抵御HIV的入侵中起着重要作用,当CD4被HIV感染而裂解时,其数量会急剧减少,HIV将迅速增加,导致AIDS发作。
艾滋病治疗的目的,是尽量减少人体内HIV的数量,同时产生更多的CD4,至少要有效地降低CD4减少的速度,以提高人体免疫能力。
迄今为止人类还没有找到能根治AIDS的疗法,目前的一些AIDS疗法不仅对人体有副作用,而且成本也很高。许多国家和医疗组织都在积极试验、寻找更好的AIDS疗法。
现在得到了美国艾滋病医疗试验机构ACTG公布的两组数据。 ACTG320(见附件1)是同时服用zidovudine(齐多夫定),lamivudine(拉美夫定)和indinavir(茚地那韦)3种药物的300多名病人每隔几周测试的CD4和HIV的浓度(每毫升血液里的数量)。193A(见附件2)是将1300多名病人随机地分为4组,每组按下述4种疗法中的一种服药,大约每隔8周测试的CD4浓度(这组数据缺HIV浓度,它的测试成本很高)。4种疗法的日用药分别为:600mg zidovudine或400mg didanosine(去羟基苷),这两种药按月轮换使用;600 mg zidovudine加2.25 mg zalcitabine(扎西他滨);600 mg zidovudine
加400 mg didanosine;600 mg zidovudine加400 mg didanosine,再加400 mg nevirapine
(奈韦拉平)。
请你完成以下问题:
预测继续治疗的效果,或者确定最佳治疗终止时间(继续(1)利用附件1的数据,
治疗指在测试终止后继续服药,如果认为继续服药效果不好,则可选择提前终止治疗)。
(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(题目中提供了患病者在不同时间点的治疗效果,即CD4和HVI的数量。仔细分析本题所给数据就会发现每个患者测试的次数不相同,测试的时间点的间距也没有规律,而且由于测试成本昂贵,所提供的时间点不超过6个。所以对于测试次数较少的患者没有什么参考价值,对于此类数据应剔除。有个别变化明显区别大多数变化规律的也予以剔除。
2(对于附件1在测试时间初始处,由于个体差异和健康状况不同,每个人的CD4的数量差距比较大但HVI波动较小,所以根据CD4的初始状态把患者分成若干段,分别
2
进行分析和预测。同理由于年龄的波动也较大,于是根据年龄将附件2中的患者分成若干段,采用灰色系统预测CD4的走势。
三、模型假设
1 CD4、HIV数量随时间的变化是连续的,对数据点做曲线拟合。
2 CD4与HIV数量的变化只与治疗的药物有关,外界状况对患者体内CD4和HIV数量的 影响忽略。
3 最佳治疗终止时间后患者患病状况不予考虑。
4 患者在治疗过程中没有间断停药和漏服。
四、模型的建立与求解
一 问题一的模型:
(一)模型一的建立:
基于前面叙述的原则一,首先对附件1中的不利数据即考察时间点较少及明显不符合大致走向的数据从数据
中剔除,从而获得较有规律的数据表,然后在结合CD4和HIV两个因素参量来观察和预测。
由于个人在接受治疗初期时CD4的含量因人而异同,而HIV的数量变化较小,从而依据CD4的初始状态将测试人群划分为七段,分别是[0,30)、[30,60)、[60,90)、[90,120)、[120,150)、[150,180)、[180,230)七个组。然后分别各分段的测试人群CD4和HIV数量的变化,利用数学软件matlab分别对其做四次多项式拟合。
下图是各个组测试人群CD4数量与时间的变化关系:
图 1
3
下图是各个组测试人群HIV数量与时间t的变化关系:
图 2
下图是CD4数量在总体下关于t 的拟和效果
图 3
4
下图是HIV数量在总体下关于时间t的拟合效果:
图 4
从以上CD4与时间t变化关系图中可以看出:在患者接受治疗初期,CD4有一段大约十周时间的明显增长,这是由于药物对体内的HIV病毒产生了作用,从而使体内的CD4的数量增加。但此后的增长放缓,这是由于HIV的产生了抗药作用,抑制了CD4的增长。 到了20周左右又达到了第二次峰值,出现了治疗时期的反弹现象,直到35周左右,HIV产生了足够强的抗药性,使得患者体内的CD4急剧的减少,不再产生反弹现象。
从以上HIV与时间t变化关系图中看出,在治疗初期,HIV有大约10周时间的减少,受到了药物的作用。但HIV产生了抗药性,产生大约10周的缓慢增长。接着又被药物抑制增长,使数量有了明显的减少。大约在35周附近,HIV的抗药性变强,药物对其不再产生作用,从而HIV的数量急剧增加。
从图3中可看出CD4与HIV呈现此消彼长的现象,而且二者的峰值基本在同一时间点附近达到对立状态,从而此图为在确定何时为最佳治疗终止时间提供了有利的说明和预测。
综上所述,在此种治疗方案下,患者体内的CD4的数量随时间的波动出现了两次波峰值(此时正好对应了HIV数量波动两次波谷值)并且第二次的峰值较第一次振动强烈。在以后的时段不会在再有反弹增加。
因而我们预测的结果为接受治疗后患者将在10周附近出现第一次的增长的峰值(在此几乎同时出现HIV减少的第一次峰值),之后出现较平缓的减少趋势,在20周附近又出现增长并且增长速度不断增加,至35周附近将达到第二次增长的峰值(在此也同时出现了HIV减少的第二次峰值),并且比前一次峰值的出现要剧烈。之后由于HIV的抗药性很强使得CD4的数量比较快的减少,并不再有反弹增加的现象。
5
(二)模型结果的分析与验证
根据处理过的数据表中每个测试患者的CD4的数量求出在测试时间的每个时间段内的CD4的变化率。然后在总体上利用数学软件matlab对CD4的变化率分别拟合曲线。
下图是CD4变化率与时间t的变化关系:
图 5
根据参考文献[5]中提到的:如果以CD4细胞的增长率30%以上为治疗有效,既无间断停药也无漏服的治疗有效率达到65.71%。
所谓的最佳治疗终止时间即要找到一个时间使得在这个时间附近治疗效率尽量大。从图 4中我们能够得到CD4的变化率大于30%的时间段大致为30周至37周,而且在35周附近增长率达到峰值; 而对于前面模型预测的最佳治疗时间25周附近落在了这个区间内,而且增长率正好达到了峰值。因此根据参考文献和模型的提前假设知前面预测的最佳治疗终止时间的精确度是很高的。
二 问题二的模型
(一)评价模型
1、模型的建立
本题关键在于确定评价的指标,通常采用模糊评价模型。而题目中要求仅以CD4为目标,则转化为单指标函数。通过单指标函数来评价4中疗法的的优劣。
从统计学角度出发,1300多名病人随机的分为4组,因此可以认为每组的年龄分布,CD4数量的初始状态分布是大致相同的。因此从总体出发来寻找指标函数。
我们利用数学软件matlab从总体上拟合每种疗法中CD4的数量关于时间t变化的曲线,将函数记为:
f(x) (i=1,2,3,4) i
将单指标函数定义为:
f(x) F=max() (i=1,2,3,4) i
6
我们取从第0周到第40周作为观察时段,四种疗法的CD4数量关于时间t的变化如下图所示:
图 6
从图中我们可以清楚的看到,在第0周到第40周这段时间内,峰值最大的是疗法四,即用疗法四是指标函数取到了最大值。
因而在仅以CD4为标准,我们评价疗法四最优。
2、模型的分析与验证:
从总体上来拟合CD4和时间t的变化规律会产生误差,因为患者在接受治疗后体内的CD4的数量有可能直接增大,也有可能直接减小,尽管总体拟合效果来看CD4的数量较高,但具体到每个患者波动向哪个方向(即增大还是减小)未能体现出来。
在此论证角度下,我们来考察四种疗法的稳定性。
所谓稳定性,就是在患者吃药后的短期内体内CD4向增大方向发展的趋势。我们将剔除不利结点后的数据表作统计,记:
m(i):表示疗法i患者接受药物治疗初期的短时间内CD4出现负增长的
人数
n(i):表示疗法i接受治疗的总人数
h(i): 表示疗法i的稳定性
则 h(i)=1-m(i)/n(i)
通过对已知数据表的分析得到了以下数据:
h(1)=0.6172
h(2)=0.5
h(3)=0.6121
h(4)=0.7540
以上数据表明,疗法四的稳定性最高。此结果既弥补了上述模型的误差性,同时有
7
验证了上述模型的正确性。
(二)预测模型
“灰色预测”主要是以少数数据建模并加以预测的方式,透过对于少数且离散的数据加以攫取、灰生成、灰建模、还原之后,最后以其灰预测模型的解求得灰数列预测。灰色预测的主要目的得乃针对少数据、离散、不确定性的数列掌握其趋势并建立出预测模型来。
现在预测艾滋病患者采用最优疗法(即第四种疗法)治疗测试结束后一年内血液中CD4的数量。为简化数据处理的过程,避免反常数据对预测结果的负面影响,根据原则二现将人群以年龄为界限分成21~30、31~40、41~50、51~60四组,将时间分为第0周,第8周,第16周,第24周,第32周,第40周六个时间点。处理数据时,我们采用灰色系统模型。本文预测患者血液中CD4的含量考虑的因素为血液中HIV的含量,但是实际情况是该因素中影响参数既有已知参数,又有未知参数和不确定参数,为一灰色系统。运用GM(1,1)模型对各年龄段患者血液中CD4的含量 进行预测分析,数据处理效果是合理有效的。
模型假设: A
1(以每个年龄段患者血液中CD4含量的平均水平为数据处理对象;
2(将每个时间点附近的点归为这个时间点进行处理;
3(用题目中所给的数据模拟出的患者血液中CD4含量的曲线能很好的反应实际情
况;
4(未来采取与以前同等的治疗方案和治疗力度;
参数说明:
(0)表示有i列原始数列 (其中i=1,2,3,4,5,6); xi
(0)(n)表示第i列原始数列中的每一个元素 xi
(其中n=1,2,3,4,5,6);
B 模型的建立:
采用灰色预测模型GM(1,1),对测试后一年中CD4数量的变化进行预测 ,结合离散模型建立一个按时间做逐段分析的模型。其中原始数据是一个序列
(t)={ (1),(2),„,(n)} (t=1,2,...,n) (1) XXXX0000
从而对(t)建立微分方程形式,即GM(1,1)模型(其中a,u为待辩史系数)。 X1
(1)dx(1) () (2) ,ax,udt
根据各年龄段患者血液中CD4含量的数据统计如下表:
[20,30] [30,40] [40,50] [50,60]
0 2.550361 2.949066 3.117051 2.932679
8 2.900796 3.330598 3.543438 3.769477
16 3.274816 3.029413 3.423345 3.446007
24 2.695796 2.934388 3.317191 3.410415
32 2.653804 2.791325 3.268072 3.192738
40 2.566657 2.584073 3.072224 2.950325
根据上表把21~30、31~40、41~50、51~60这四个年龄段患者血液中CD4含量的原始数据分别列写成数据形式:
(0)=(2.550361,2.900796,3.274816,2.695796,2.653804,2.566657) X1
8
(0)=(2.949066,3.330598,3.029413,2.934388,2.791325,2.584073) X2(0)=(3.117051,3.543438,3.423345,3.317191,3.268072,3.072224) X3
(0)=(2.932679,3.769477,3.446007,3.410415,3.192738,2.950325) X4
建立四个时间响应函数方程:
(1)dx(1)1? ,ax,u111dt(1)dx(1)2 ? ,ax,u222dt
(1)dx(1)3 ? ,ax,u333dt(1)dx(1)4 ? ,ax,u444dt
T设,其中; ,,i,1,2,3,4,5,6A,a,uiiiT,,0,,0,,0,,0; ,,,,,,,,,,Y,x2,x3,x4,x5iiiii
按最小二乘法的超定方程组解法可得
,1TT ,,A,BBBYiiiii11,,xx,1/2(1,2)1ii,,,,,,,,,,11xx,1/22,31,,,,,,ii,,,,,,,,其中 B= i、、、、,,
,,,,11,,,,,,,,xx,1/25,61,,ii
各类时间响应函数方程的解为:
u,(a*k),a*(k,1)(0)(0)111 X(k,1),(X(1),)*(e,e)11a1
u,(a*k),a*(k,1)(0)(0)222 X(k,1),(X(1),)*(e,e)22a2
u,(a*k),a*(k,1)(0)(0)333 X(k,1),(X(1),)*(e,e)33a3
u,(a*k),a*(k,1)(0)(0)444 X(k,1),(X(1),)*(e,e)44a4
利用各类白化形式的微分方程来对第七个时间点的各类人血液中CD4的平均含量
进行求解,也就是让依次取7,8,9,10,11,12,代入时间响应函数方程求出在未k
来一年里不同年龄人所占的比重,根据最后所得到的数据就能对未来患者血液中CD4浓
度的发展趋势作出预测分析所有的结果如下表所示:
20-30 30-40 40-50 50-60
48 2.456786 2.447104 3.008729 2.822889
56 2.348538 2.306382 2.911369 2.668309
64 2.245059 2.173752 2.817113 2.522194
72 2.14614 2.048749 2.725909 2.384081
80 2.051579 1.930934 2.637657 2.25353
88 1.961184 1.819895 2.552262 2.130128
9
根据原始数据和预测一年以后得到的数据拟合成曲线,针对不同的年龄段拟合成4条曲线,如下图所示:
图 7
从图象上可以看出,CD4的走向是缓慢减少的而且不同年龄层上的患者到达峰值是不同的,峰值所对应的时间也不相同。
年龄在20到30之间时,到达峰值的时间在16周附近,此也即为最佳治疗终止时间。
年龄在30到40之间时,到达峰值的时间在10周附近,此也即为最佳治疗终止时间。
年龄在40到50之间时,到达峰值的时间在15周附近,此也即为最佳治疗终止时间。
年龄在50到60之间时,到达峰值的时间在15周附近,此也即为最佳治疗终止时间。 C 模型的检验和分析
回代检验
有实测值与模拟制得比较结果得出,用GM(1,1)模型所预测的CD4含量与实际值
最大误差为-8.0%,平均误差0.16%,吻合较好,精度较高。
关联度检验
模拟值与实测值之间的关联度(R)下式确定
,,,(t),minmin|x(t),x(t)|,,maxmax|x(t),x(t)|/[x(t),x(t)|,,maxmax|x(t),x(t)|]01010101
式中R为关联度; 为实测值于摸拟值在(t)的关联系数; 为分辨系数,0-1之,(t),
间,一般取0.5, 为模拟和实测值的残差值的绝对最小minmin|x(t),x(t)|01
值;maxmax|x(t),x(t)|为模拟和实测值的残差值的绝对最大值。经计算得关联度01
R=0.62>0.5,则模型有足够的精度。
经上述检验,可以看出,该模型预测精度较高,故不需作残差分析修定,此模型可以作为该测百分比的预测模型,预测结果见上表。
三 问题三的模型
问题三是在问题二的基础上又考虑经济对治疗艾滋病的影响,即在仅以CD4为标准
10
的基础上再增加对治疗费用的考虑,我们建立一个新的指标函数模型来求解。
类比于普通商品的性价比,我们可以将某种治疗方案疗效的量化函数与治疗费用的比值作为指标函数,可看作该治疗方法的“性价比”,通过这一指标的高低来判定相应治疗方案的优劣
为了简单起见,将患者在治疗过程中CD4的含量看作是治疗效果,它是关于一个关于时间t的函数。而治疗费用则是根据题目中每种疗法的日用药的费用与治疗时间的乘积,于是建立起指标函数为:
指标函数=治疗效果/治疗费用;
在这里治疗效果我们采用问题二中四种疗法拟合的曲线来代替,则可以作出相应的指标函数:
f1=
(-0.00000047190507*x^4+0.00006592932538*x^3-0.00273175795027*x^2+0.02200991
333399*x+2.98485748047478)/7*1.225*x
f2=
-0.00000205450268*x^4+0.00017137897680*x^3-0.00462182638939*x^2+0.02793036(
537014*x+3.10257653254036)/7*3.45*x
f3=
(-0.00000361036121*x^4+0.00038497021508*x^3-0.01370990968307*x^2+0.16995320
464564*x+3.02677739250685)/7*2.45*x
f4=
(-0.00000185266104*x^4+0.00021706946957*x^3-0.00828791234678*x^2+0.10222444
504744*x+2.95333264000659)/7*4.3*x
利用matlab 将f1,f2,f3,f4关于时间t的变化规律反映到下图:
图 8
11
从图象观察,指标函数在相当长的一段时间里是呈上升趋势,而且图中上升最快的曲线是疗法四的指标函数。
结合问题二中的图6我们可以发现在疗法4指标函数增长的过程中,CD4的数量达到了峰值,而且其峰值超过其它疗法。即疗法4的效果和“性价比”(即指标函数)相比较其他疗法较优。
所以我们得到结论:疗法四最优。
五、模型的评价与推广
模型的评价:
1(本题的关键在于处理数据,在建立模型处理数据时规定了两个原则,使得能够
充分的利用有效数据,避免了特殊结点对模型的影响。
2(针对本题所给数据离散的特点,在做出相应的预测后,用曲线拟合的方法将离
散数据点用光滑的曲线来反映。
3(定下合理的假设,使对模型结果的分析和验证作出更好的解释,使人更容易理
解。
4(本题充分挖掘数据中的关系,对数据的特点从图象的角度来观看,更容易得出
合理的结论。
模型的改进:
1(可建立更多的因素变量来反映数据的变化,例如患者的停药或漏服现象对数据
的影响。
2(在处理数据时,尽量科学化,在避免特殊结点影响模型的同时,也尽量不要将
能反映规律的结点删掉。
参考文献
[1] 姜启源,数学模型,北京:高等教育出版社,1987
[2] 郝红伟,MATLAB6实例
,北京:中国电力出版社,2001.9 [3] 张志涌等,精通Matlab6.5版教程,北京:北京航空航天大学出版社,2003.3 [4] 李庆扬,关治,白峰杉,数值计算原理,北京:清华大学出版社,2000.9 [5] 于兰,豆智慧,曲淑霞等,HAART药物依从性对CD4细胞增长率的影响研究,中国
艾滋病性病,11(4):255-257,2005
12
附录:
1(问题二的灰色预测的程序代码实现:
x0=[2.550361 2.900796 3.274816 2.695796 2.653804 2.566657];
for i=1:6
s=0;
for j=1:i
s=s+x0(j);
end
x1(i)=s;
end
z1(1)=0;
for i=2:6
z1(i)=0.5*x1(i)+0.5*x1(i-1);
end
B=[-z1(2) 1;-z1(3) 1;-z1(4) 1;-z1(5) 1;-z1(6) 1]; Y=[x0(2),x0(3),x0(4),x0(5),x0(6)]';
U=(B'*B)\B'*Y;
for i=7:12
x(i-6)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2)));
x0(i)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2)));
end
x
t=[0 8 16 24 32 40 48 56 64 72 80 88]; a=polyfit(t,x0,4)
ezplot('-0.00000024249800*x^4+0.00004891300192*x^3-0.00328255524646*x^2+0.0
6786915594827*x+2.57788513753247',[1,90]) hold on
x0=[2.949066 3.330598 3.029413 2.934388 2.791325 2.584073];
for i=1:6
s=0;
for j=1:i
s=s+x0(j);
end
x1(i)=s;
end
z1(1)=0;
for i=2:6
z1(i)=0.5*x1(i)+0.5*x1(i-1);
end
B=[-z1(2) 1;-z1(3) 1;-z1(4) 1;-z1(5) 1;-z1(6) 1]; Y=[x0(2),x0(3),x0(4),x0(5),x0(6)]';
U=(B'*B)\B'*Y;
13
for i=7:12
x(i-6)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2))); x0(i)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2))); end
x
t=[0 8 16 24 32 40 48 56 64 72 80 88];
a=polyfit(t,x0,4)
ezplot('-0.00000014783411*x^4+0.00003024751139*x^3-0.00203852222453*x^2+0.0
3204116804071*x+3.01572219979800',[1,90])
hold on
x0=[3.117051 3.543438 3.423345 3.317191 3.268072 3.072224]; for i=1:6
s=0;
for j=1:i
s=s+x0(j);
end
x1(i)=s;
end
z1(1)=0;
for i=2:6
z1(i)=0.5*x1(i)+0.5*x1(i-1);
end
B=[-z1(2) 1;-z1(3) 1;-z1(4) 1;-z1(5) 1;-z1(6) 1]; Y=[x0(2),x0(3),x0(4),x0(5),x0(6)]';
U=(B'*B)\B'*Y;
for i=7:12
x(i-6)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2)));
x0(i)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2))); end
x
t=[0 8 16 24 32 40 48 56 64 72 80 88];
a=polyfit(t,x0,4)
ezplot('-0.00000014608628*x^4+0.00003076679028*x^3-0.00216698705891*x^2+0.0
4515032663197*x+3.17374575047441',[1,90])
hold on
x0=[2.932679 3.769477 3.446007 3.410415 3.192738 2.950325]; for i=1:6
s=0;
for j=1:i
s=s+x0(j);
end
x1(i)=s;
14
end
z1(1)=0;
for i=2:6
z1(i)=0.5*x1(i)+0.5*x1(i-1);
end
B=[-z1(2) 1;-z1(3) 1;-z1(4) 1;-z1(5) 1;-z1(6) 1];
Y=[x0(2),x0(3),x0(4),x0(5),x0(6)]';
U=(B'*B)\B'*Y;
for i=7:12
x(i-6)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2)));
x0(i)=(x0(1)-U(2)/U(1))*(exp(-U(1)*(i-1))-exp(-U(1)*(i-2))); end
x
t=[0 8 16 24 32 40 48 56 64 72 80 88];
a=polyfit(t,x0,4)
ezplot('-0.00000029417459*x^4+0.00006022911938*x^3-0.00411803817585*x^2+0.0
8555977073217*x+3.04205656177807',[1,90])
hold on
grid on
2(问题三指标函数画图的代码:
ezplot('(-0.00000047190507*x^4+0.00006592932538*x^3-0.00273175795027*x^2+0.
02200991333399*x+2.98485748047478)/7*1.225*x',[1,70]) hold on
ezplot('(-0.00000205450268*x^4+0.00017137897680*x^3-0.00462182638939*x^2+0.
02793036537014*x+3.10257653254036)/7*3.45*x',[1,70]) hold on
ezplot('(-0.00000361036121*x^4+0.00038497021508*x^3-0.01370990968307*x^2+0.
16995320464564*x+3.02677739250685)/7*2.45*x',[1,70]) hold on
ezplot('(-0.00000185266104*x^4+0.00021706946957*x^3-0.00828791234678
*x^2+0.10222444504744*x+2.95333264000659)/7*4.3*x',[1,70])
15