數(shù)值計算實驗word版-資料下載頁
2025-01-07 16:08本頁面
【正文】 function values y % max1 is the maximum number of iterations %Output p1 is the secant method approximation to a zero % err is the error estimate for p1 % k is the number of iterations % y is the function value f(p1) for k=1:max1 p2=p1feval(f,p1)*(p1p0)/( feval(f,p1)feval(f,p0))。 err=abs(p2p1)。 relerr=2*err/(abs(p2)+delta)。 p0=p1。 p1=p2。 y=feval(f,p1)。 if(errdelta)|(relerrdelta)|(abs(y)epsilon),break,end end 二、數(shù)值試驗內(nèi)容 1）用 Newton 迭代或加速 Newton 迭代法求下列方程的根 1． 100?p ，求 91 的近似值（精確到小數(shù)點后 10 位）． 2． 20 ??p ，求 3 7? 的近似值（精確到小數(shù)點后 10 位）． 3． )ln()1()( xxxf ?? ， 1,2 * ?? xm 初始值 20?p ． 4．設(shè) xxexf ??)( ，分別取初始近似值 ?p 和 ?p ．討論其結(jié)果． 2）用割線迭代算法求下列方程的根 設(shè) xxexf ??)( ，初始近似值 , ?? pp ． 數(shù)值實驗 6 求矩陣的特征對 一、 方法與程序（或 Matlab 命令） Hessenberg 陣 H 的 QR 算法 （ 1）使用 Givens 變換對 H 陣進行 QR 分解； （ 2）令 1HH? 則 kkk RUH ? ％ kH 的 QR 分解 kkk URH ??1 對稱矩陣 A 的 Lanczos 算法 （１）利用 Lanczos 迭代將 A 三對角化 Tnnn QTQA? （ 2）求 0)det( ?? nTI? 的根 i? （ 3）求 A 的特征對 ),( )(ii z? 有關(guān)命令 （ 1）［ Q,R］ =qr(A) % 輸出正交陣 Q 和上三角陣 R （ 2） [P,R]=eig(A) % 輸出 A 的特征值對角陣 R 和特征向量 P 二、實驗內(nèi)容 利用 Hatlab 中的命令，編寫 QR 算法程序 求矩陣 5H 的所有特征對 ?????????????????01000982008763076544654325H 利用 Matlab 語言編寫 Lanczos 算法和特征方程求根的程序，計算 5A 的所有特征對． ?????????????????42100242101242101242001245A 利用“ eig”命令求 5H 和 5A 的特征對，并與 QR 算法與 Lanczos 算法得到的特征對進行比較．
公司管理相關(guān)推薦
鄂ICP備17016276號-1