打靶法(含matlab程序)(編輯修改稿)
2024-08-03 13:58 本頁面
【文章內(nèi)容簡介】 (k,:)39。+k1*h/2。k2=feval(dydx1,x2,y2)。 k3=feval(dydx1,x2,Y1(k,:)39。+k2*h/2)。 k4=feval(dydx1, X(k)+h,Y1(k,:)39。+k3*h)。 Y1(k+1,:)=Y1(k,:)+h*(k139。+2*k239。+2*k339。+k439。)/6,k=k+1。endu=Y1(:,1)for k=1:nk1=feval(dydx2,X(k),Y2(k,:)) x2=X(k)+h/2。y2=Y2(k,:)39。+k1*h/2。k2=feval(dydx2,x2,y2)。 k3=feval(dydx2,x2,Y2(k,:)39。+k2*h/2)。 k4=feval(dydx2, X(k)+h,Y2(k,:)39。+k3*h)。 Y2(k+1,:)=Y2(k,:)+h*(k139。+2*k239。+2*k339。+k439。)/6,k=k+1。endv=Y2(:,1)Y=u+(betau(n+1))*v/v(n+1)for k=2:n+1wucha(k)=norm(Y(k)Y(k1))。 k=k+1。endX=X(1:n+1)。Y=Y(1:n+1,:)。k=1:n+1。wucha=wucha(1:k,:)。P=[k39。,X39。,Y,wucha39。]。plot(X,Y(:,1),39。ro39。,X,Y1(:,1),39。g*39。,X,Y2(:,1),39。mp39。)xlabel(39。軸\it x39。)。 ylabel(39。軸\it y39。)legend(39。是邊值問題的數(shù)值解y(x)的曲線39。,39。是初值問題1的數(shù)值解u(x)的曲線39。, 39。是初值問題2的數(shù)值解v(x)的曲線39。)title(39。用線性打靶法求線性邊值問題的數(shù)值解的圖形39。)%有限差分法function [k,A,B1,X,Y,y,wucha,p]=yxcf(q1,q2,q3,a,b,alpha,beta,h)n=fix((ba)/h)。 X=zeros(n+1,1)。 Y=zeros(n+1,1)。 A1=zeros(n,n)。
高考資料相關推薦
鄂ICP備17016276號-1