數(shù)值分析課程課程設計-閱讀頁
2025-02-02 15:57本頁面
【正文】 三次多項式中的四個系數(shù)；（2）用所求出的三次多項式函數(shù)繪制出飛機降落曲線。a2=[1 ,x2,x2^2,x2^3]。a4=[0,1,2*x4,3*x4^2]。 a2。a4]。 y2。y4]。x=12000:250:0。d39。x39。y39。因而發(fā)表論文，提出了大大有名的摩爾定律（Moore’s Law），并預測未來這種增長仍會延續(xù)下去。這些數(shù)據是推出摩爾定律的依據：年代19591962196319641965增加倍數(shù)13456試從表中數(shù)據出發(fā)，推導線性擬合的函數(shù)表達式。y=[1,3,4,5,6]。39。r*39。x39。y39。運行結果：p1 = +003 * y1 = 線性擬合的函數(shù)表達式： Y=+003解法二：運行結果：xs = +003 * x=[1959 1962 1963 1964 1965]。3。5。for i=1:length(x) for j=1:2 A(i,j)=x(i)^(j1)。*A)。*y))從而年代y與增加倍數(shù)x之間的關系為：y=+5．1 用幾種不同的方法求積分的值。解：syms xi1=int(4/(1+x^2),x,0,1)運行結果：i1 =pii2 = 3i3 = i4 = a=0。h=ba。h=(ba)/M。for k=1:(M1) x=a+h*k。endi4=h*(4/(1+a^2)+4/(1+b^2))/2+h*i4結果分析：牛頓萊布尼茲公式得到精確結果，采用梯形公式得到的結果比采用Simpson公式的精確度要低，采用復化梯形公式在步長取得越來越小的狀態(tài)下可以提高精度。x=diff(cos(t))。z=diff(2cos(t)sin(t))。 %符號函數(shù)積分digits(14)。t=t0+[0:N]39。x(1,:)=x039。 f=f39。end以39。保存function f=Euler_fun(t,x)f=*x./(1+2*t)。39。Euler_fun39。fh=dsolve(39。,39。)。 fx(k)=subs(fh,ft(k))。39。=f(x,y)%a,b為左右端點%N為迭代步長%h為步長%ya為初值h=(ba)/N。Y=zeros(1,N+1)。Y(1)=ya。 k2=h*feval(f,T(j)+h/2,Y(j)+k1/2)。 k4=h*feval(f,Y(j)+h,Y(j)+k3)。endans=[T39。] 以39。保存function z=f(x,y)z=*y./(1+2*x)。39。f39。其中，是初值問題 的解。Dv*log(x)=2*x39。v(0)=039。x39。f=(1log(v))*v f =2*(1log(2*Ei(1,2*log(x))))*Ei(1,2*log(x)) ezplot(39。)輸出結果：f =2*(1log(2*Ei(1,2*log(x))))*Ei(1,2*log(x))的系數(shù)為[1 a1 a2 …a9 a10]=[1 –55 1320 –181 50 157 773 –902 055 341 693 0 840 950 0 127 535 76 106 286 40 632 880 0]試用多項式的求根指令roots求出該10次方程的10個根，然后修改9次項的系數(shù)55為56，得新的10次方程，求解新的方程，觀察根的變化是否很顯著。 roots(p)ans = + + + + + p=[1 56 1320 18150 157773 902055 3416930 8409500 12753576 10628640 6328800]。解析：ode**函數(shù)無能為力Matlab中，提供了 bvp解算器。close。odefun=inline(39。y(1)1]39。t39。y39。bcfun=inline(39。yb(1)]39。ya39。yb39。sol=bvp4c(odefun,bcfun,sinit)t=linspace(0,1,101)。plot(t,y(1,:),(1,:),39。)legend(39。,39。,2)xlabel(39。),ylabel(39。)運行結果：sol = x: [0 1] y: [2x5 double] yp: [2x5 double] solver: 39。8．3 通過平面上的三個點P1（x1,y1）,P2(x2,y2),P3(x3,y3)可以確定一個圓的方程。解：（1）將圓上三個點的坐標數(shù)據分別代入可得：（2）整理上述方程組得：其中，； 運行結果為：x = 2 4 20r = 5O = 1 2（3）（4）MATLAB代碼：A=[1 3 1。6 2 1]。50。x=A\br=sqrt((x(1)^2+x(2)^2)/4x(3))x0=(x(1)/2)。O=[x0,y0]結果：半徑 為5，圓心坐標（1，2）
環(huán)評公示相關推薦
鄂ICP備17016276號-1