資源描述：
《matlab數(shù)學(xué)建模實(shí)例》由會(huì)員上傳分享，免費(fèi)在線(xiàn)閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、第四周3.functiony=mj()forx0=0:0.01:8x1=x0^3-11.1*x0^2+38.79*x0-41.769;if(abs(x1)<1.0e-8)x0endend4.分別用簡(jiǎn)單迭代法、埃特金法、牛頓法求解方程，并比較收斂性與收斂速度（e分別取10-3、10-5、10-8）。簡(jiǎn)單迭代法：functiony=jddd(x0)x1=(20+10*x0-2*x0^2-x0^3)/20;k=1;while(abs(x1-x0)>=1.0e-3)x0=x1;x1=(20+10*x0-2*x0^2-x0^3)/20;k=k+
2、1;endx1k埃特金法：functiony=etj(x0)x1=(20-2*x0^2-x0^3)/10;x2=(20-2*x1^2-x1^3)/10;x3=x2-(x2-x1)^2/(x2-2*x1+x0);k=1;while(abs(x3-x0)>=1.0e-3)x0=x3;x1=(20-2*x0^2-x0^3)/10;x2=(20-2*x1^2-x1^3)/10;x3=x2-(x2-x1)^2/(x2-2*x1+x0);k=k+1;endx3k牛頓法：functiony=newton(x0)x1=x0-fc(x0)/df(x0)
3、;k=1;while(abs(x1-x0)>=1.0e-3)x0=x1;x1=x0-fc(x0)/df(x0);k=k+1;endx1kfunctiony=fc(x)y=x^3+2*x^2+10*x-20;functiony=df(x)y=3*x^2+4*x+10;第六周1.解例6-4（p77)的方程組，分別采用消去法（矩陣分解）、Jacobi迭代法、Seidel迭代法、松弛法求解，并比較收斂速度。消去法：x=ad或[L,U]=lu(a);x=inv(U)inv(L)dJacobi迭代法：functions=jacobi(a,d,x
4、0)D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1);C=inv(D);B=C*(L+U);G=C*d;s=B*x0+G;n=1;whilenorm(s-x0)>=1.0e-8x0=s;s=B*x0+G;n=n+1;endnSeidel迭代法：functions=seidel(a,d,x0)D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1);C=inv(D-L);B=C*U;G=C*d;s=B*x0+G;n=1;whilenorm(s-x0)>=1.0e-5x0=s
5、;s=B*x0+G;n=n+1;endn松弛法：functions=loose(a,d,x0,w)D=diag(diag(a));U=-triu(a,1);L=-tril(a,-1);C=inv(D-w*L);B=C*((1-w)*D+w*U);G=w*C*d;s=B*x0+G;n=1;whilenorm(s-x0)>=1.0e-8x0=s;s=B*x0+G;n=n+1;endn2.練習(xí)MATLAB的常用矩陣語(yǔ)句，就龍格現(xiàn)象函數(shù)（p88）練習(xí)插值語(yǔ)句interp,spline，并比較。3.測(cè)得血液中某藥物濃度隨時(shí)間的變化值為:t(h)
6、0.250.51.01.52.03.04.06.08.010.0C(mg/L)19.3018.1515.3614.1012.899.327.555.243.862.88求t=0.45,1.75,5.0,6.0時(shí)的濃度C.分別用n=4,5,9的拉格朗日插值計(jì)算；并用樣條函數(shù)插值計(jì)算,并比較結(jié)果。拉格朗日插值:functions=lagr(n)x=[0.250.51.01.52.03.04.06.08.010.0];y=[19.3018.1515.3614.1012.899.327.555.243.862.88];x0=[0.451.75
7、5.06.0];m=length(x0);fori=1:mD=abs(x-x0(i));I=1;whileI<=n+1fora=1:length(x)ifD(a)==min(D)c(I)=a;D(a)=max(D)+1;breakendendI=I+1;endb=sort(c);z=x0(i);t=0.0;fork=1:length(b)u=1.0;forj=1:length(b)ifj~=ku=u*(z-x(b(j)))/(x(b(k))-x(b(j)));endendt=t+u*y(b(k));ends(i)=t;end樣條函數(shù)差
8、值：Interp1(x,y,x0,’spline’)Spline(x,y,x0)第八周1.給定某藥物濃度隨時(shí)間的變化值(作業(yè)3)，1）分別采用樣條函數(shù)和三點(diǎn)公式（設(shè)h=0.1)求結(jié)點(diǎn)處的導(dǎo)數(shù)值，并比較結(jié)果。2）求該時(shí)間段
