第7章matlab解方程與函數(shù)極值71線性方程組求解72非線-免費(fèi)閱讀
【正文】 (2) 調(diào)用 fmin函數(shù)求最小值點(diǎn)。,x0) 這兩個(gè)函數(shù)的調(diào)用格式相似。 y139。tf=10。 tspan形式為 [t0,tf],表示求解區(qū)間。,39。 q(1)=*sin(x)*cos(y)。最優(yōu)化工具箱提供了 20多個(gè)選項(xiàng)，用戶可以使用 optimset命令將它們顯示出來(lái)。 function fx=funx(x) fx=x10.^x+2。0]) 非線性方程數(shù)值求解 單變量非線性方程求解 在 MATLAB中提供了一個(gè) fzero函數(shù)，可以用來(lái)求單變量非線性方程的根。2,2,1]。 在命令中調(diào)用函數(shù)文件 ，命令如下： A=[10,1,0。 f=(DL)\b。 [x,n]=jacobi(A,b,[0,0,0]39。 y=B*x0+f。 elseif nargin3 error return end D=diag(diag(A))。1,6,1,4]。當(dāng) X為對(duì)稱正定的，則 p=0， R與上述格式得到的結(jié)果相同；否則 p為一個(gè)正整數(shù)。 x=R\(Q\b) 或采用 QR分解的第 2種格式，命令如下： [Q,R,E]=qr(A)。 實(shí)現(xiàn) QR分解后，線性方程組 Ax=b的解 x=R\(Q\b)或x=E(R\(Q\b))。1,6,1,4]。注意，這里的矩陣 X必須是方陣。0,2,1,1。1,6,1,4]。 [L,U,P]=lu(X)：產(chǎn)生一個(gè)上三角陣 U和一個(gè)下三角陣 L以及一個(gè)置換矩陣 P，使之滿足 PX=LU。 b=[13,9,6,0]39。 例 73 用 QR分解求解例 71中的線性方程組。 x=E*(R\(Q\b)) (3) Cholesky分解 如果矩陣 X是對(duì)稱正定的，則 Cholesky分解將矩陣 X分解成一個(gè)下三角矩陣和上三角矩陣的乘積。如果 X為滿秩矩陣，則 R為一個(gè)階數(shù)為q=p1的上三角陣，且滿足 R39。 b=[13,9,6,0]39。 %求 A的對(duì)角矩陣 L=tril(A,1)。 n=n+1。,) 2． GaussSerdel迭代法 在 Jacobi迭代過(guò)程中，計(jì)算時(shí)，已經(jīng)得到，不必再用，即原來(lái)的迭代公式 Dx(k+1)=(L+U)x(k)+b可以改進(jìn)為Dx(k+1)=Lx(k+1)+Ux(k)+b，于是得到： x(k+1)=(DL)1Ux(k)+(DL)1b 該式即為 GaussSerdel迭代公式。 y=G*x0+f。1,10,2。 b=[9。該函數(shù)的調(diào)用格式為： z=fzero(39。 (2) 調(diào)用 fzero函數(shù)求根。如果想改變其中某個(gè)選項(xiàng)，則可以調(diào)用
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1