微分和差分方程作业.docVIP

微分方程和差分方程作业题 一．Apollo卫星的运动轨迹的绘制 解： M文件的内容： function dx=appollo(t,x) u=1/82.45; u1=1-u; r1=sqrt((x(1)+u)^2+x(3)^2); r2=sqrt((x(1)-u1)^2+x(3)^2); dx=[x(2) 2*x(4)+x(1)-u1*(x(1)+u)/r1^3-u*(x(1)-u1)/r2^3 x(4) -2*x(2)+x(3)-u1*x(3)/r1^3-u*x(3)/r2^3]; 命令窗口内容： x0=[1.2;0;0;-1;%x0(i)对应与xi的初值 [t,y]=ode45(@appollo,[0,20],x0,options); plot(y(:,1),y(:,3))%绘制图形 title(Appollo卫星运动轨迹) xlabel(X) ylabel(Y) 运行结果： 二．每户人家都希望有一套属于自己的住房，但又没有足够的资金一次买下。这就产生了贷款买房问题。某新婚夫妇急需一套属于自己的住房。他们看到一则理想的房产广告：“名流花园之高尚住宅公寓，供工薪阶层选择。一次性付款优惠价40.2万元。若不能一次性付款也没关系，只付首期款为15万元，其余每月1977.04元等额偿还，15年还清(公积金贷款月利息为3.675‰）。表示第k个月还款后的欠款 模型： 代入n=180、 r=0.003675、 B=1977.04 结果：A=260000（元）则原价为41万元 三．2周后成熟，平均产卵100个，四周龄的成虫平均产卵150个。假设每个卵发育成2周龄成虫的概率为0.09，（称为成活率），2周龄成虫发育成4周龄成虫的概率为0.2。 （1）假设开始时， 0~2，2~4，4~6周龄的昆虫数目相同，计算2周、4周、6周后各种周龄的昆虫数目； （2）讨论这种昆虫各种周龄的昆虫数目的演变趋势：各周龄的昆虫比例是否有一个稳定值？昆虫是无限地增长还是趋于灭亡？ （3）假设使用了除虫剂，已知使用了除虫剂后各周龄的成活率减半，问这种除虫剂是否有效？,2-4周虫的数量为,4-6周虫数量为 又开始时，0-2，2-4，4-6周龄的昆虫数目相同,将其数目设为100.即，依题意可得模型： 则可解出两周后0-2周龄昆虫为25000，2-4周龄为9,4-6周龄为20；四周后昆虫为3900,2-4周龄为2250,4-6周龄为2；六周后0-2周龄昆虫为225270,2-4周龄为351,4-6周龄为450 （3）由各周龄的成活率减半可建立模型： 由MATLAB作图可得： 则由图可分析得：当年份足够长时，整个昆虫的数量相比未使用除虫剂时有了一定的减少，但还在无限增长，故这种除虫剂的效果不明显。 四．早期肿瘤的体积增长满足Malthus模型（，其中λ为常数），（1）求肿瘤的增倍时间σ。根据统计资料，一般有σ∈(7,465)（单位为天），肺部恶性肿瘤的增倍时间大多大于70天而小于465天（发展太快与太慢一般都不是恶性肿瘤），故σ是确定肿瘤性质的重要参数之一（2）为方便起见，医生通常用肿瘤直径来表示肿瘤的大小，试推出医生用来预测病人肿瘤直径增大速度的公式. (3) 正常人身上也有癌细胞，一个癌细胞直径约为10μm，重约0.001μg.，当患者被查出患有癌症时，通常直径已有1cm以上（即已增大1000倍），由此容易算出癌细胞转入活动期已有30σ天，故如何在早期发现癌症是攻克癌症的关键之一。手术治疗常不能割去所有癌细胞，故有时需进行放射疗法。射线强度太小无法杀死癌细胞，太强病人身体又吃不消且会使病人免疫功能下降。一次照射不可能杀死全部癌细胞，请设计一个可行的治疗方案（医生认为当体内癌细胞数小于105个时即可凭借体内免疫系统杀灭。 解： （1），其中λ为常数 （其中c为常数） 当 （2） 因为 t=0时， 所以 （3）假设，不考虑杀死免疫细胞从而影响人类免疫功能，一位病人癌细胞直径为1cm，癌细胞含量为 ，治疗打算用17/2 天使癌细胞减少到可杀死水平，每两次治疗的间隔为1/2 ,用9次治疗使癌细胞得到控制。 程序： clear; clc; for x=1:900000 k=1000000-x; for i=1:9 k=sqrt(2)*k-x; end if k=100000 break end end x 解得，x=301006.每隔1/2 进行一次放射性治疗，强度为杀死301006个癌细胞。