19次方系數(shù)對(duì)方程根的影響課程設(shè)計(jì)報(bào)告書(shū)(編輯修改稿)
2024-09-25 05:35 本頁(yè)面
【文章內(nèi)容簡(jiǎn)介】 (i)=f(x(i))。 end a=polyfit(x,y,3) syms x f1=a(1)*x^3+a(2)*x^2+a(3)*x+a(4) x=[1 0 1]。 for i=1:5 E(i)=(f(x(i))(a(1)*x(i)^3+a(2)*x(i)^2+a(3)*x(i)+a(4)))^2。 end sum(E) 輸出結(jié)果為 a = 8 f1 = 4878849915647781/1298074214633706907132624082305024*x^31600/2639*x^25348847520430703/649037107316853453566312041152512*x+1514/2639 （擬合的多項(xiàng)式） ans = （均方誤差） 九個(gè)點(diǎn)的時(shí)候， matlab 編碼為： x=[1 0 1]。 for i=1:9 y(i)=f(x(i))。 end a=polyfit(x,y,3) syms x f2=a(1)*x^3+a(2)*x^2+a(3)*x+a(4) x=[1 0 1]。 for i=1:5 E1(i)=(f(x(i))(a(1)*x(i)^3+a(2)*x(i)^2+a(3)*x(i)+a(4)))^2。 end sum(E1) 輸出結(jié)果為： a = f2 = 728732707776039/2535301202056458802993406410752*x^31263030908712921/2251799813685248*x^2+4915246442354361/20282409603651670423947251286016*x+1093229300764671/2251799813685248 （最小二乘擬合多項(xiàng)式） ans = （均方誤差） 9 用復(fù)合梯形公式求積分 120 41 dxx??的值。 function i=biraprl(f,a,b,eps) %復(fù)化梯形公式求函數(shù) f 在區(qū)間 [a， b]上的定積分 %函數(shù)名： f %積分下限： a %積分上 限： b %積分精度： eps %積分值： i %積分劃分的子區(qū)間個(gè)數(shù)： step n=1。 h=(ba)/2。 i1=0。 i2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h。 while abs(i2i1)eps n=n+1。 h=(ba)/n。 i1=i2。 i2=0。 for i=0:n1 x=a+h*i。 x1=x+h。 i2=i2+(h/2)*(subs(sym(f),findsym(sym(f)),x)+... subs(sym(f),findsym(sym(f)),x1))。 end end 四階龍格 庫(kù)塔法分別求解下列初值問(wèn)題； 2(1 ) 39。 , ( 0) 1 。 [ 0 , 1 ]12( 2) 39。 , ( 0) 2 。 [ 0 , 1 ]1yy y xxxyy y xx? ? ??? ? ? ?? function[x,y]=runge_kutta(fun,x0,xt,y0,pointnum,varargin) if nargin3 y0=0。 end y(1,:)=y0(:)39。 h=(xtx0)/(pointnum1)。 x=x0+[0:pointnum]39。*h。 for k=1:pointnum f1=h*feval(fun,x(k),y(k,:),varargin{:})。 f1=f1(:)39。 10 f2=h*feval(fun,x(k)+h/2,y(k,:)+f1/2,varargin{:})。 f2=f2(:)39。 f3=h*feval(fun,x(k)+h/2,y(k,:)+f2/2,varargin{:})。 f3=f3(:)39。 f4=h*feval(fun,x(k)+h,y(k,:)+f3,varargin{:})。 f4=f4(:)39。 y(k+1,:)=y
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1