数学建模微分方程实验作业.docxVIP

实验1实验报告2012326601042 曹骏森 12信科1班一、完成教材（2013高教版）P129的 1,2,3,4实验；1、dsolve(Dy+3*y=8,y(0)=2,x)ans =8/3-2/3*exp(-3*x)2、dsolve((1+x^2)*D2y=2*x*Dy,y(0)=1,Dy(0)=3,x)ans =1+x^3+3*x3、dsolve(D4y-2*D3y+D2y=0,x)ans =C1*exp(x)+C2*exp(x)*x+C3+C4*x4[x,y]=dsolve(2*Dx+4*x+Dy-y=exp(t),Dx+3*x+y=0,x(0)=3/2,y(0)=0,t)x =-3/2*exp(7^(1/2)*t)*(1/3*7^(1/2)+1/3)-3/2*exp(-7^(1/2)*t)*(-1/3*7^(1/2)+1/3)+1/6*exp(t)+1/2*7^(1/2)*exp(7^(1/2)*t)*(1/3*7^(1/2)+1/3)-1/2*7^(1/2)*exp(-7^(1/2)*t)*(-1/3*7^(1/2)+1/3)y =exp(7^(1/2)*t)*(1/3*7^(1/2)+1/3)+exp(-7^(1/2)*t)*(-1/3*7^(1/2)+1/3)-2/3*exp(t)二、完成教材（2013高教版）P205的 1,2,3,4实验；1、function exampleclearode45(@fun,[0,1],2)function f=fun(x,y)f=8-3*y;2、function exampleclearode45(@fun,[0,1],[1,3])function f=fun(x,t)f=[t(1);(1+x^2)/(2*x*t(1))];3、function exampleclearode45(@fun,[0,10],[1,2,3])function f=fun(t,x)f=[x(2);x(3);2*x(3)-x(2)];4、function exampleclearode45(@fun,[0,5],[1,2])function f=fun(t,x)f=[exp(t)+3*x(2)+2*x(1);-3*x(1)-x(2)]3.从制冰厂购买了一块立方体的冰块，在运输途中发现，第一小时大约融化了1/4,（1） 求冰块全部融化要多长时间（设气温不变）（2） 如运输时间要2.5小时，问：运输途中冰块大约融化掉多少？解：设冰块永远保持立方体，冰块的融化速度与其表面积成比例， 在t时刻，体积为v(t),面积为S（t）；1)即 dv/dt =- kS = -kv^(2/3)dv/v^(2/3) = kdt两块积分得： v(t)=k(c-t)^3v(1)/v(0) = (c-1)^3/(c)^3 = 1/4 c=2^(2/3)== 3.85 小时； 所以3.85小时后融化完.2)v(2.5)/v(0) = (c-2.5）^3/c^3 = 0.043;即只剩下4.3%4、正常人身上也有癌细胞，一个癌细胞直径约为10μm，重约0.001μg，（1）当患者被查出患有癌症时，通常直径已有1cm以上（即已增大1000倍），由此容易算出癌细胞转入活动期已有300多天，故如何在早期发现癌症是攻克癌症的关键之一。（2）手术治疗常不能割去所有癌细胞，故有时需进行放射疗法。射线强度太小无法杀死癌细胞，太强病人身体又吃不消且会使病人免疫功能下降。一次照射不可能杀死全部癌细胞。请设计一个可行的治疗方案（医生认为当体内癌细胞数小于100000个时即可凭借体内免疫系统杀灭)。解：假设在不考虑每次放射性射线杀死的免疫细胞，从而影响人体免疫系统功能的前提下，有一位癌症病人的癌细胞的直径为1cm，此癌症病人的体内含有的癌细胞的含量为106个，现在医院要对病人进行治疗，打算用17/2σ天使癌症病人的体内的癌细胞数目从原来的106个减少到体内免疫细胞可杀死的水平，每两次放射性治疗的时间间隔为1/2σ。决定用9次级放射性治疗来使病人的癌细胞得到控制，为达到此目的每次治疗要用射线治疗的射线强度为杀死癌细胞的x所需的射线量。?根据假设条件，编写matlab程序：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个癌细胞的射线量。