【MATLAB】一維搜尋之牛頓法
clear;clc;formatcompact
syms x
f =3*x^4-16*x^3+30*x^2-24*x+8;
% f =x^4-4*x^3-6*x^2-16*x+4; %書上P99例題
x0 = 3; %設定初始點
i = 1;
while(i)
x0 = x0-subs(diff(f,x),x0)/subs(diff(diff(f,x),x),x0);%牛頓法迭代
if abs(subs(diff(f,x),x0))<1e-4 %迭代停止的條件
i = 0;
end
end
disp(sprintf('用牛頓法求得的近似解為%9.15f',double(x0))) %列印結果
程式輸出結果：
用牛頓法求得的近似解為2.000000002785312
相關推薦
【MATLAB】一維搜尋之牛頓法
clear;clc;formatcompact syms x f =3*x^4-16*x^3+30*x^2-24*x+8; % f =x^4-4*x^3-6*x^2-16*x+4; %書上P99例題 x0 = 3; %設定初始點 i = 1; while(i) x
【MATLAB】一維搜尋之平分法
clc;clear; syms x a = -10; b = 10; f =3*x^4-16*x^3+30*x^2-24*x+8; dfdx =diff(f,x); while(1) x0 = 1/2*(a+b); dfdx0 = subs(dfdx,x,x
【工程優化】一維搜尋方法
一維搜尋方法的分類如下： 這篇文章主要講解黃金分割法、二分法、牛頓法這三種一維搜尋方法。黃金分割法只用到原函式，二分法用到函式的一階導，牛頓法用到函式的二階導。由於本文主要對研一上學期的課程中的部分演算法進行程式實現，理論部分大多參考上課的課件。 黃金分
【Java】一維數組
div bsp exceptio 異常 read exc linear void ava /* 1.演示初始化一維數組的三種方法 2.演示數組的 length 屬性，用 length 獲取數組長度 */ public class linear_array { p
【 MATLAB 】訊號處理工具箱之 fft 案例分析
上篇博文：【 MATLAB 】訊號處理工具箱之fft簡介及案例分析介紹了MATLAB訊號處理工具箱中的訊號變換 fft 並分析了一個案例，就是被噪聲汙染了的訊號的頻譜分析。 這篇博文繼續分析幾個小案例： Gaussian Pulse 這個案例是將高斯脈衝從時域變換到頻域，高斯脈衝的資
【 MATLAB 】訊號處理工具箱之fft簡介及案例分析
目錄 Syntax Description Y = fft(X) Y = fft(X,n) Y = fft（X，n，dim） Examples Noisy Signal Syntax Y = fft(X) Y = fft(X
【 MATLAB 】訊號處理工具箱之 ifft 簡介及案例分析
這篇博文和上篇博文對應：【 MATLAB 】訊號處理工具箱之fft簡介及案例分析 目錄 ifft Syntax Description 案例分析 Inverse Transform of Vector Padded Inverse Transform of Matri
【 MATLAB 】nextpow2 函式用法之 Optimize FFT with Padding
您可以使用nextpow2來填充傳遞給fft的訊號。 這樣做可以在訊號長度不是2的精確冪次時加速FFT的計算。 Optimize FFT with Padding 下面這個例子展示了 使用填充優化FFT的案例，通過使用函式nextpow2完成： clc c
【 MATLAB 】訊號處理工具箱之 dct 簡介及案例分析
dct Discrete cosine transform Syntax y = dct(x) y = dct(x,n) y = dct(x,n,dim) y = dct(___,'Type',dcttype) Description y = dct（x）
【 MATLAB 】訊號處理工具箱之波形產生函式 pulstran
前兩篇博文和這篇博文有些許聯絡： MATLAB幫助文件稱pulstran函式為：脈衝串（Pulse train）產生函式。 pulstran從連續函式或採樣原型脈衝生成脈衝序列。 語法格式：
一維搜尋方法/黃金分割法（附matlab程式碼）
一維搜尋方法中的黃金分割法（附matlab程式碼） 一維搜尋方法：一維搜尋，又稱一維優化，是指求解一維目標函式 f(X) 最優解的過程，分為試探法和插值法。 黃金分割法：屬於一維搜尋方法中的試探法，適用於[a，b]區間上的任何單谷函式求極小值問題。 證明r=0.
【matlab】MarkDown Letex 編碼 之 隨機過程及應用（三）
**Provement of Gaussian Distribution：** 設正態分佈概率密度函式是 $$f(x)=\frac{1}{\sqrt{2π}\sigma}*e^{\
【原創】運維基礎之Ansible（1）簡介、安裝和使用
ets 安裝 yum ant gem get 結構 ges describe 官方：https://www.ansible.com/ 一 簡介 Ansible is a radically simple IT automation engine that automate
【原創】運維基礎之Nginx（1）簡介、安裝、使用
官方：http://nginx.org nginx [engine x] is an HTTP and reverse proxy server, a mail proxy server, and a generic TCP/UDP proxy server, originally written by
【原創】運維基礎之Redis（1）簡介、安裝、使用
lists 腳本 分享 ngs 參考 ports eos 運維基礎 lru redis 5.0.3 官方：https://redis.io/ 一 簡介 Redis is an open source (BSD licensed), in-memory data str
【MATLAB】二維矩陣視覺化 MATLAB繪圖
基本繪圖函式 1、plot(y) 如果y是向量，則採用向量的索引值作為橫座標值，以向量元素的值作為縱座標值。 如果y是實數矩陣，則相當於對y的每一列進行二維繪圖。 如果y是複數組成的向量，則相當於plot(real(y),imag(y))。 2、plot(x,y)
【原創】運維基礎之Nginx（3）location
大小 uri 規則 ati 第一個 基礎 匹配規則 最大 cati nginx location =：精確匹配（必須全部相等） ~：大小寫敏感，正則匹配 ~*：忽略大小寫，正則匹配 ^~：只需匹配uri部分，精確匹配 @：內部服務跳轉，精確匹配 規則
【 MATLAB 】filter 函式介紹（一維數字濾波器）
filter 1-D digital filter Syntax y = filter(b,a,x) y = filter(b,a,x,zi) y = filter(b,a,x,zi,dim) [y,zf] = filter(___) Description
第十四屆華中科技大學程序設計競賽 J Various Tree【數值型一維BFS/最小步數】
hust there ini amp rep targe ati 數值 set 鏈接：https://www.nowcoder.com/acm/contest/106/J 來源：牛客網 題目描述 It’s universally acknowledged that th
【51nod-1521】一維戰艦
不能 編號 格子 每次 .html ray 長方形 http main 愛麗絲和鮑博喜歡玩一維戰艦的遊戲。他們在一行有n個方格的紙上玩這個遊戲（也就是1×n的表格）。 在遊戲開始的時候，愛麗絲放k個戰艦在這個表格中，並不把具體位置