打靶法(含matlab程序)(參考版)
2025-07-10 13:58本頁(yè)面
【正文】 b39。r39。 n=m.*exp(1/2*m)。r39。 x=ans(:,1)。 h= dbf(f,x0l,x0u,y0l,y0u,h,1e6)。 alfa=0。 x0l=0。+(h/6)*(L1+2*L2+2*L3+L4)。+h*L3)。+(h/2)*L2)。+(h/2)*L1)。%迭代初值for i=1:m1 L1=f(t(i),x(:,i))。%區(qū)間長(zhǎng)度t(1)=t0。 function x=rk4(f,t0,x0,h,a,b)%rungkuta法求每個(gè)點(diǎn)的近似值（參考大作業(yè)一）t=a:h:b。,yvalue39。 yvalue=y1(1,:)。 y1=y2。 s1=s2。 y2=rk4(ff,a,x0,h,a,b)。 end end if flag~=1 while abs(y1(1,n)beta)eps s2=s1(y1(1,n)beta)*(s1s0)/(y1(1,n)y0(1,n))。 y1=rk4(ff,a,x0,h,a,b)。 else s1=s0+1。 if abs(y0(1,n)beta)=eps flag=1。%迭代初值flag=0。%x取值范圍n=length(xvalue) s0=。打靶法 源代碼：創(chuàng)建M 文件：function ys=dbf(f,a,b,alfa,beta,h,eps) ff=(x,y)[y(2),f(y(1),y(2),x)]。,y,wucha39。),p=[k39。)title(39。),legend(39。 ylabel(39。軸\it x39。mp39。 wucha=wucha(1:k,:)。 y=y(1:n+1,1)。 k=k+1。beta]。 y=[alpha。Y=A\B1。A=A1(1:n1,1:n1)+A2(1:n1,1:n1)+A3(1:n1,1:n1)。endB(1,1)=(h.^2)*k3(1)(1+h*k1(1)/2)*alpha。 k3(k)=feval(q3,X(k))。A2(k,k)=
高考資料相關(guān)推薦
鄂ICP備17016276號(hào)-1