matlab與科學(xué)工程計(jì)算-資料下載頁
2025-02-13 14:22本頁面
【正文】 2 0],6,... [0 0 0],[],... ’cons’... )。 function v = fun(x) v =3*x(1)x(2)^2+x(3)。 function [c,ceq] = cons(x) c(1)= x’*x16。 c(2)= x(2)*(x(1)+x(3))12。 ceq = []。 8 s . t . 321 ??? xxx0, 321 ?xxx62 21 ?? xx16232221 ??? xxx12)( 312 ?? xxx最優(yōu)解 [6 0 0]’ 74 微分方程工具箱 (ode45) [t,y] = ode45(odefun,tspan,y0,options,p1,p2...) 求解函數(shù) 工具箱 函數(shù) 算法參數(shù) help odeset 求解區(qū)間 初值 解的向量 可用 plot(x,y)畫出解 目標(biāo)函數(shù) 其他參數(shù) 75 微分方程工具箱 (ode45) 39。 1 ,(0 ) 1 , 0 2y y xyx? ? ?? ? ? ( ) e xy x x??真 解function dydx =f(x,y) dydx = y – x + 1。 [x,y] = ... ode45(’f’,[0,2],1)。 plot(x,y,’r’)。 76 多項(xiàng)式計(jì)算 ? 多項(xiàng)式表示：向量 22 ?xp = [1 0 2] 表示 ? 多項(xiàng)式求值 polyval(p,1:5) ? 多項(xiàng)式求根 r = roots(p) ? 用根生成多項(xiàng)式 p = poly(r) ?多項(xiàng)式相乘 p = [1 0 2]。 q = [2 1]。 s = conv(p,q)。 ?多項(xiàng)式相除 [q,r] = deconv(s,p)。 ?多項(xiàng)式求導(dǎo) d = polyder(p) 77 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 0 . 200 . 20 . 40 . 60 . 811 . 2多項(xiàng)式計(jì)算 ? 多項(xiàng)式擬合 (polyfit) x = linspace(0,pi,8) x = 0 y = sin(x) y = 0 p = polyfit(x,y,2) p = xx = linspace(0,pi,200)。 yy = polyval(p,xx)。 plot(xx,yy,39。r39。,x,y,39。b39。,39。linewidth39。,2) 78 插值 (interp1) yi=interp1(x， y， xi， 39。method39。) 插值方法 被插值點(diǎn) 插值節(jié)點(diǎn) 插值結(jié)果 ‘ nearest’: 最鄰近插值‘linear’: 線性插值 ‘ spline’: 三次樣條插值 ‘ cubic’: 立方插值 缺省時(shí) : 分段線性插值 79 插值 (interp2) z=interp2(x0,y0,z0,x,y,’method’) 被插值點(diǎn) 插值方法 插值節(jié)點(diǎn) ‘ nearest’ 最鄰近插值 ‘ linear’ 雙線性插值 ‘ cubic’ 雙三次插值 缺省時(shí) , 雙線性插值 被插值點(diǎn)的函數(shù)值 80 統(tǒng)計(jì)工具箱 ?參數(shù)估計(jì) ?betafit, binofit, mle, … ?概率密度函數(shù) ?betapdf, binopdf, hygepdf, … ?各種分布的均值、方差 ?betastat, binostat, chi2stat, … ?線性、非線性的各種回歸模型 ?anova1, anova2, polyfit, regress, … ?lsqnonneg, nlinfit, nlintool, … 81 matlab上機(jī)操作 ?遵守機(jī)房各項(xiàng)規(guī)定 ?保持機(jī)房衛(wèi)生整潔 ?matlab系統(tǒng)若有問題，請(qǐng)重新開機(jī)，在選擇進(jìn)入系統(tǒng) 前 按 CtrlR 82 matlab練習(xí)題 ? 求下面表達(dá)式的值 其中 76242/2412321xxxxxxxxxxxy???????????????????????????????????????????? ? 6 1 7 , , , , , , , ?? ?x83 matlab練習(xí)題 ? 已知矩陣 A,B如下 , 計(jì)算并比較 A*B, B*A, A.*B, B.*A, A/B, A\B, A./B, A.\B ?????????????????????????????11032132101 ,13221026 BA84 matlab練習(xí)題 ? 已知矩陣 A,B如下 , 給出這兩個(gè)矩陣的 (1)逆 (2)行列式 (3)跡 (4)特征值和特征向量 (5)最簡(jiǎn)階梯形矩陣 ?????????????????????????????11032132101 ,13221026 BA85 matlab練習(xí)題 ? 已知矩陣 A及多項(xiàng)式如下 , 給出多項(xiàng)式在矩陣上的值； 比較結(jié)果矩陣和原矩陣的特征值和特征向量。 1)( ,13221024??????????????????xxxfA86 matlab綜合訓(xùn)練 ?矩陣 : 寫一個(gè)程序，對(duì)某一輸入的對(duì)稱矩陣 A進(jìn)行分解 , A=LDL’, 其中 L為單位下三角矩陣， D為對(duì)角矩陣。 參考 《 數(shù)值計(jì)算 》 中的 Cholesky分解。 ?畫圖 : 給出某一曲面方程及一垂直于 XOY平面的平面方程，利用平面把曲面割成兩片，畫出這兩片，中間的縫隙要逐漸變大。 ?算法 : 輸入平面上的至少四個(gè)點(diǎn)，挑選當(dāng)中的部分點(diǎn)，并把它們連成一個(gè)多邊形，要求面積盡量大，給出其面積之值。 87 matlab命令速查 : .* ./ .^ \ abs all and any asin asinh bar3 ceil cell char chol et3 conj contour conv deconv det diag disp doc else elseif end eig exp eval eye false feval fix floor fmincon fminu for fprintf function fzero help hilb i if imag inf input interp1 interp2 inv j linprog linspace log log10 lookfor magic mesh meshgrid NaN not num2str ode45 ones or pascal pi pie3 plot plot3 polar poly polyder polyfit polyval prod rand rank real roots round sin sinh sqrt str2num struct surf title trace tril triu true view what which while who whos xlabel xor ylabel zeros 88 作業(yè) : 三維畫圖 在某海域測(cè)得一些點(diǎn) (x,y)處的水深 z由下表給出，船的吃水深度為 5英尺，在矩形區(qū)域（ 75， 200） *（ 50， 150）里的哪些地方船要避免進(jìn)入。 x y z 4 8 6 8 6 8 8 x y z 9 9 8 8 9 4 9 89 作業(yè) : 三維畫圖 90 THANK YOU
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1