數(shù)值計算大作業(yè)——劉

1、課 程 設 計課程名稱：數(shù)值分析設計題目：數(shù)值計算大作業(yè)學 號：S315070064姓 名：劉峰完成時間：2015年10月25日題目一、非線性方程求根1.題目 假設人口隨時間和當時人口數(shù)目成比例連續(xù)增長，在此假設下人口在短期內的增長建立數(shù)學模型。（1）如果令表示在時刻的人口數(shù)目，表示固定的人口出生率，則人口數(shù)目滿足微分方程，此方程的解為；（2）如果允許移民移入且速率為恒定的，則微分方程變成，此方程的解為；假設某地區(qū)初始有1000000人，在第一年有435000人移入，又假設在第一年年底該地區(qū)人口數(shù)量1564000人，試通過下面的方程確定人口出生率，精確到；且通過這個數(shù)值來預測第二年年末的人口數(shù)

2、，假設移民速度保持不變。2.數(shù)學原理采用牛頓迭代法，牛頓迭代法的數(shù)學原理是，對于方程，如果是線性函數(shù)，則它的求根是很容易的，牛頓迭代法實質上是一種線性化方法，其基本思想是將非線性方程逐步歸結為某種線性方程來求解。設已知方程有近似根(假定)，將函數(shù)在點進行泰勒展開，有于是方程可近似地表示為這是個線性方程，記其根為，則的計算公式為，這就是牛頓迭代法，簡稱牛頓法。3.程序設計作出函數(shù)的圖像，大概估計出根的位置fplot(1000*exp(x)+(435*x)*(exp(x)-1)-1564,0 3);grid大概估計出初始值x=0.5function p1,err,k,y=newton(f,df,p

3、0,delta,max1)% f是非線性系數(shù)% df是f的微商% p0是初始值% dalta是給定允許誤差% max1是迭代的最大次數(shù)% p1是牛頓法求得的方程近似解% err是p0誤差估計% k是迭代次數(shù)p0,feval(f,p0)for k=1:max1 p1=p0-feval(f,p0)/feval(df,p0); err=abs(p1-p0); p0=p1; p1,err,k,y=feval(f,p1) if(err a=-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,0 b=1,1,1,1,1,1,1,1,1,1 c=-0.5,-0.5,-0

4、.5,-0.5,-0.5,-0.5,-0.5,-0.5,-0.5,0 f=0.5,0,0,0,0,0,0,0,0,0得到此題中的a,b,c,f矩陣：a = -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 b = 1 1 1 1 1 1 1 1 1 1c = -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 -0.5000 0f = 0.5000 0 0 0 0 0 0 0 0 0然后在MATLAB中調用之前保存的迭代法函數(shù)

5、function，在命令窗口中輸入：chase(a,b,c,f)回車得到結果： x=chase(a,b,c,f)x =0.9000 0.8000 0.7000 0.6000 0.5000 0.4000 0.3000 0.2000 0.1000 0追趕法為一種特殊的LU分解法。追趕法是求解三對角矩陣的常用方法，但從整體編程角度分析，其程序編寫較迭代法復雜，但通用性較好。追趕法求解三對角矩陣不但節(jié)省存儲單元，而且可以減少計算量，是工程技術中比較常用的數(shù)學工具。三、數(shù)值積分1、題目衛(wèi)星軌道是一個橢圓，橢圓周長的計算公式是, 這里是橢圓的半長軸, 是地球中心與軌道中心(橢圓中心)的距離, 記為近地點距

6、離, 為遠地點距離, 公里為地球半徑,則, 某人造衛(wèi)星近地點距離公里,遠地點距離公里, 試用Romberg方法求衛(wèi)星軌道的周長，精確到。2.數(shù)學原理龍貝格方法是在梯形公式、辛普森公式和柯特斯公式之間的關系的基礎上，構造出一種加速計算積分的方法。 作為一種外推算法, 它在不增加計算量的前提下提高了誤差的精度。 龍貝格方法的主要過程是將粗糙的梯形公式逐步加工成精度較高的辛普森公式和科特斯公式的方法稱為龍貝格方法。復化梯形公式 在復化梯形公式中，每個內節(jié)點既是前一個小區(qū)間的終點，又是后一個小區(qū)間的起點，因此上式可以改寫為復化梯形公式余項 復化梯形公式的遞推公式為復化辛普森求積公式與復化梯形公式類似，

7、每個內節(jié)點需用兩次，因此有顯然復化辛普森公式在n趨于無窮大時，他的收斂速度比復化梯形公式更快。以表示二分k次后求得的梯形值，且以表示序列的m次加速度，理查森外推法的遞推公式可寫成龍貝格算法的計算過程如下：(1) 取求(2) 利用變步長梯形公式，其中k為區(qū)間的二分次數(shù)，即或(3) 依橫行次序求加速值，逐個求出的第k行其余各元素(4) 當相鄰對角元素之差的絕對值小于預先給定的精度時，終止計算。表3-1龍貝格算法遞推表kh0b-a12343.程序設計function R=romberg(f,a,b,n)format longR=zeros(n+1,n+1);R(0+1,0+1)=(b-a)/2*(f

8、eval(f,a)+feval(f,b);for i=1:n,h=(b-a)/2i; s=0; for k=1:2(i-1), s=s+feval(f,a+(2*k-1)*h); end R(i+1,0+1)=R(i-1+1,0+1)/2+h*s;endfor j=1:n,fac=1/(4j-1); for m=j:n, R(m+1,j+1)=R(m+1,j-1+1)+fac*(R(m+1,j-1+1)-R(m-1+1,j-1+1); endend4.結果分析與討論本題根據算法原理在matlab中編寫完龍貝格算法的自定義程序后，直接輸入符合格式的函數(shù)積分就可得到相應軌道周長。調用MATLAB龍貝格算法的函數(shù)后可算得 R = romberg(4*7800*sqrt(1-(973.5/7880)2*sin(x)2),0,pi/2,6)計算出來得出R=49136.836545由此可得精10-6確到的衛(wèi)星軌道