matlab課程設(shè)計(jì)word版(存儲(chǔ)版)
【正文】 ? ? ? ? 43/32 22 ??? yx 算法說(shuō)明： 此題相當(dāng)于求兩一個(gè)二元二次方程組的解，故為便于清楚地顯示出兩橢圓的相對(duì)位置，用 ezplot 函數(shù)把兩個(gè)橢圓畫在同一個(gè)坐標(biāo)系中，然后利用 fsolve 函數(shù)解方程組得到兩橢圓的交點(diǎn)即方程組的解。,[1,5, 8,8])。 [x,y]=solve(f1,f2,39。 %合并 x,y兩個(gè)矩陣 intersection_x_y=double(middle) %將符號(hào)解轉(zhuǎn)換成數(shù)值解 程序運(yùn)行結(jié)果： 開始 依照所給方程畫出兩個(gè)橢圓 用符號(hào)解法求兩橢圓的交點(diǎn) 將符號(hào)解轉(zhuǎn)換成數(shù)值解 結(jié)束 intersection_x_y = + + + + + 。2*(x3)^2+(y/3)^2=439。 %畫第一個(gè)橢圓 hold on ezplot(39。主應(yīng)力平面方向角（角度）： 39。or39。)。 text(Smin,0,39。.r39。 %畫出 x坐標(biāo)軸 hold on plot(Sa,0,39。 Tmin=min(Tau)。A39。)。.r39。Tau_xy(MPa)=39。clc。RungKutta法解 39。b39。 %設(shè)定函數(shù)初值 PointNum=5。 %真值解 fun=inline(39。 %得公式 k4 y(k+1,:)=y(k,:)+(f1+2*(f2+f3)+f4)/6。 %得 x向量值 輸入求解的自變量范圍 求出待求簡(jiǎn)單微分方程的真值解 用 MATLAB 自帶函數(shù) ode23求解待求微分方程 結(jié)束 用自編函數(shù)四階龍格－庫(kù)塔（ RK）方法求解待求微分方程 開始 for k=1:(PointNum) %迭代計(jì)算 f1=h*feval(fun,x(k),y(k,:),varargin{:})。 確定求解范圍內(nèi)的步長(zhǎng) k = 取點(diǎn)數(shù)？ 否 求解：???????????????????) ,()21 ,21()21 ,21(),(3423121hkyhxfkhkyhxfkhkyhxfkyxfknnnnnnnn 求解并輸出：)22(6 3211 kkkhyy nn ????? 是 結(jié)束
畢業(yè)設(shè)計(jì)相關(guān)推薦
鄂ICP備17016276號(hào)-1