MATLAB實(shí)驗(yàn)三代數(shù)方程的近似解.ppt

《MATLAB實(shí)驗(yàn)三代數(shù)方程的近似解.ppt》由會(huì)員分享，可在線閱讀，更多相關(guān)《MATLAB實(shí)驗(yàn)三代數(shù)方程的近似解.ppt（24頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、實(shí) 驗(yàn) 三求 代 數(shù) 方 程 的 近 似 根 (解 )數(shù) 學(xué) 實(shí) 驗(yàn) q 問(wèn) 題 背 景 和 實(shí) 驗(yàn) 目 的實(shí) 驗(yàn) 三 、 近 似 求 解 代 數(shù) 方 程u 解 方 程 （ 代 數(shù) 方 程 ） 是 最 常 見(jiàn) 的 數(shù) 學(xué) 問(wèn) 題 之 一 ， 也 是眾 多 應(yīng) 用 領(lǐng) 域 中 不 可 避 免 的 問(wèn) 題 之 一 。u 目 前 還 沒(méi) 有 一 般 的 解 析 方 法 來(lái) 求 解 非 線 性 方 程 ， 但 如果 在 任 意 給 定 的 精 度 下 ， 能 夠 解 出 方 程 的 近 似 解 ， 則 可以 認(rèn) 為 求 解 問(wèn) 題 已 基 本 解 決 ， 至 少 可 以 滿 足 實(shí) 際 需 要 。u

2、本 實(shí) 驗(yàn) 主 要 介 紹 一 些 有 效 的 求 解 方 程 的 數(shù) 值 方 法 ： 對(duì)分 法 ， 迭 代 法 和 牛 頓 法 。 同 時(shí) 要 求 大 家 學(xué) 會(huì) 如 何 利 用Matlab 來(lái) 求 方 程 的 近 似 解 。 相 關(guān) 概 念0( )f x u 如 果 f(x) 是 一 次 多 項(xiàng) 式 ， 稱 上 面 的 方 程 為 線 性 方程 ； 否 則 稱 之 為 非 線 性 方 程 。q 線 性 方 程 與 非 線 性 方 程本 實(shí) 驗(yàn) 主 要 討 論 非 線 性 方 程 的 數(shù) 值 求 解 q 基 本 思 想 對(duì) 分 法將 有 根 區(qū) 間 進(jìn) 行 對(duì) 分 ， 判 斷 出 解 在 某

3、 個(gè) 分 段 內(nèi) ， 然 后再 對(duì) 該 段 對(duì) 分 ， 依 次 類 推 ， 直 到 滿 足 給 定 的 精 度 為 止 。q 適 用 范 圍求 有 根 區(qū) 間 內(nèi) 的 單 根 或 奇 重 實(shí) 根 。q 數(shù) 學(xué) 原 理 ： 介 值 定 理設(shè) f(x) 在 a, b 上 連 續(xù) ， 且 f(a) f(b)0， 則 由 介 值 定理 可 得 ， 在 (a, b) 內(nèi) 至 少 存 在 一 點(diǎn) 使 得 f()=0。 q 具 體 步 驟 對(duì) 分 法設(shè) 方 程 在 區(qū) 間 a,b 內(nèi) 連 續(xù) ， 且 f(a)f(b)0， 給 定精 度 要 求 ， 若 有 |f(x)| ， 則 x 就 是 我 們 所 需 要

4、的 f(x) 在 區(qū) 間 (a,b) 內(nèi) 的 近 似 根 。；， 計(jì) 算令 )( 2/)( )1( 00 xfbax ；輸 出 結(jié) 果停 止 計(jì) 算 ， 的 近 似 根 ，就 是 我 們 所 要， 則若 000 |)(| )2( xxxxf ；否 則 令， 令若 bbxaxbaaxfaf 1010110 , ;, 0)()( )3( ；輸 出 結(jié) 果， 則 停 止 計(jì) 算 ，若令 11111 |)(|,2/)( )4( xxxfbax ；否 則 令， 令若 1212121211 , ;, 0)()( bbxaxbaaxfaf . .Matlab程 序 見(jiàn) fulu1.m q 收 斂 性 分 析

5、 對(duì) 分 法 收 斂 性 = =1 1 11 1 1 1| | ( ) ( ) ( )2 2 2 2k k k k k kx b a b a b a 設(shè) 方 程 的 根 為 x* (ak , bk ) ， 又 ， 所 以2k kk a bx 0(k )對(duì) 分 法 總 是 收 斂 的u 但 對(duì) 分 法 的 收 斂 速 度 較 慢u 通 常 用 來(lái) 試 探 實(shí) 根 的 分 布 區(qū) 間 ， 或 給 出 根 的 一 個(gè) 較 為 粗 糙 的 近 似 。根 據(jù) 上 面 的 算 法 ， 我 們 可 以 得 到 一 個(gè) 每 次 縮 小 一 半 的區(qū) 間 序 列 ak , bk ， 在 (ak , bk ) 中

6、 含 有 方 程 的 根 。 迭 代 法q 基 本 思 想u 構(gòu) 造 f (x) = 0 的 一 個(gè) 等 價(jià) 方 程 ： ( )x xu 從 某 個(gè) 近 似 根 x0 出 發(fā) ， 計(jì) 算得 到 一 個(gè) 迭 代 序 列 0k kx 1 ( )k kx x k = 0, 1, 2, . . (x) 的 不 動(dòng) 點(diǎn)f (x) = 0 x = (x)等 價(jià) 變 換f (x) 的 零 點(diǎn) u 若 收 斂 ， 即 ， 假 設(shè) (x) 連 續(xù) ， 則q 收 斂 性 分 析迭 代 法 的 收 斂 性 1lim lim ( ) limk k kk k kx x x lim *kk x x *x ( *)x kx

7、* ( *)x x ( *) 0f x 即 注 ： 若 得 到 的 點(diǎn) 列 發(fā) 散 ， 則 迭 代 法 失 效 ！ q 定 義 ： 迭 代 法 收 斂 性 判 斷q 定 理 2： 如 果 定 理 1 的 條 件 成 立 ， 則 有 如 下 估 計(jì) 1 0| *| | |1 kk qx x x xq 11| *| | |1k k kx x x xq 如 果 存 在 x* 的 某 個(gè) 鄰域 =(x*- , x* + ), 使得 對(duì) x0 開(kāi) 始 的 迭 代 xk+1 = (xk) 都 收 斂 , 則 稱 該 迭 代 法 在 x* 附 近 局 部 收 斂。q 定 理 1： 設(shè) x* =(x*)， 的

8、 某 個(gè) 鄰 域 內(nèi) 連 續(xù) ， 且 對(duì) x 都 有 |(x)|q 1, 則 對(duì) x0 ， 由 迭代 xk+1 = (xk) 得 到 的 點(diǎn) 列 都 收 斂 。 迭 代 法 收 斂 性 判 斷 1 0| *| | |1 kk qx x x xq q 定 理 3： 已 知 方 程 x =(x)， 且(1) 對(duì) xa, b，有 (x)a, b；(2) 對(duì) xa, b，有 |(x)|q syms x f=sin(x)+3*x2; g=diff(f,x) g=diff(sin(x)+3*x2,x) f=inline(函 數(shù) 表 達(dá) 式,變 量 1,變 量 2,.)y=f(數(shù)值列表)代 入 的 數(shù) 值

9、列 表 順 序 應(yīng) 與 定 義 時(shí) 的 變 量 名 順 序 一 致例 ： 附 錄 ： inlinel inline 命 令 可 以 用 來(lái) 定 義 一 個(gè) 內(nèi) 聯(lián) 函 數(shù)l 調(diào) 用 方 式 ： u 這 種 函 數(shù) 定 義 方 式 是 將 f 作 為 一 個(gè) 內(nèi) 部 函 數(shù) 調(diào) 用 。 其 特點(diǎn) 是 ： 調(diào) 用 方 式 最 接 近 于 我 們 平 時(shí) 對(duì) 函 數(shù) 的 定 義 ， 使 程 序更 具 可 讀 性 。 同 時(shí) 由 于 它 是 基 于 Matlab 的 數(shù) 值 計(jì) 算 內(nèi) 核的 ， 所 以 它 的 運(yùn) 算 速 度 較 快 ， 程 序 更 有 效 率 。u 這 種 定 義 方 式 的 缺 點(diǎn) ：l 定 義 一 個(gè) 內(nèi) 聯(lián) 函 數(shù) 用 去 的 內(nèi) 存 空 間 比 相 同 條 件 下 其 他的 方 法 要 大 得 多 。l 該 方 法 只 能 對(duì) 數(shù) 值 進(jìn) 行 代 入 ， 不 支 持 符 號(hào) 代 入 ， 并 且對(duì) 于 定 義 后 的 函 數(shù) 不 能 進(jìn) 行 求 導(dǎo) 等 符 號(hào) 運(yùn) 算 。自 定 義 函 數(shù) 方 式 （ 三 ） 自 定 義 函 數(shù) 方 式 （ 三 ） 教 材 ： P87, 4q 作 業(yè) （ 要 求 寫 實(shí) 驗(yàn) 報(bào) 告 ）上 機(jī) 作 業(yè)