函數(shù)零點(diǎn)_第1頁(yè)
已閱讀1頁(yè),還剩11頁(yè)未讀 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

1、,,建模實(shí)驗(yàn):購(gòu)房貸款的利率和最佳訂貨量,不難算出,你向銀行總共借了25.2萬(wàn),30年內(nèi)共要還51.696萬(wàn),這個(gè)案例中貸款年利率是多少呢?,《新民晚報(bào)》2000年3月30日上的一則房產(chǎn)廣告:,設(shè)xk為第k個(gè)月的欠款數(shù), a為月還款數(shù), r為月利率,模型假設(shè),模型建立,根據(jù) a=0.1436, x0=25.2, x360=0得到 25.2(1+r)360 – 0.1436[(1+r)360-1]/r=0,模型求解,常識(shí)上,r

2、應(yīng)比當(dāng)時(shí)活期存款月利率略高一些。我們用活期存款月利率0.0198/12 作為迭代初值,用fzero求解>>clear; fun=inline('25.2*(1+r)^360-((1+r)^360-1)/r*0.1436' ,'r')>>r=fzero(fun,0.0198/12);>>R=12*r得年利率為5.53%,一、預(yù)備知識(shí):零點(diǎn)、極值和最小二乘法,一元

3、非線性方程的一般形式為 f (x) = 0若對(duì)于數(shù)?有f (?) = 0, 則稱?為方程的解或根,也稱為函數(shù)f (x)的零點(diǎn),若f (?) = 0, f ’(?)?0 則?稱為單根。若有k >1, f (?) = f ’(?) = …= f (k-1)(?) = 0但f (k)(?)?0 , 稱為k重根,非線性方程(組)求解通常用數(shù)值方法,,零點(diǎn),求近似解,常見的有二分法、牛頓法等,極值,最小二乘法,假設(shè)已知經(jīng)驗(yàn)公式y(tǒng)

4、=f(c,x)(c和x均可為向量), 要求根據(jù)一批有誤差的數(shù)據(jù)(xi,yi), i=0,1,…,n, 確定參數(shù)c.這樣的問題稱為數(shù)據(jù)擬合。,當(dāng)f關(guān)于c是線性函數(shù),問題轉(zhuǎn)化為一個(gè)線性方程組求解,且其解存在唯一。 如果f關(guān)于c是非線性函數(shù),問題轉(zhuǎn)化為一個(gè)函數(shù)極值問題,1 多項(xiàng)式  y=polyval(p,x) 求得多項(xiàng)式p在x處的值y, x可以是一個(gè)或多個(gè)點(diǎn) p3=conv(p1,p2) 返回多項(xiàng)式p1和p2的乘

5、積,[p3,r]=deconv(p1,p2) p3返回多項(xiàng)式p1除以p2的商,r返回余項(xiàng),x=roots(p) 求得多項(xiàng)式p的所有復(fù)根.p=polyfit(x,y,k) 用k次多項(xiàng)式擬合向量 數(shù)據(jù)(x, y),返回多項(xiàng) 式的降冪系數(shù),二、函數(shù)零點(diǎn)、極值和擬合的MATLAB命令,MATLAB中一個(gè)多項(xiàng)式用系數(shù)降冪排列向量來(lái)表示,例2.用2

6、次多項(xiàng)式擬合下列數(shù)據(jù). x 0.1 0.2 0.15 0 -0.2 0.3 y 0.95 0.84 0.86 1.06 1.50 0.72,例1.求多項(xiàng)式x3 + 2 x2 - 5的根,x=[0.1 0.2 0.15 0 -0.2 0.3];>> y=[0.95 0.84 0.86 1.06 1.50 0.72];>> p=polyfi

7、t(x,y,2),p = 1.7432 -1.6959 1.0850,» roots([1 2 -5]),2. 定義非線性函數(shù),fun=inline(‘funstr’,’var’) 定義一個(gè)Inline函數(shù),其中funstr是函數(shù)的表達(dá)式, var是變量名,例、f=inline(‘x.^2-3’); %建立一元函數(shù),默認(rèn)x為自變量; u=inline(‘x.

8、^y-3’, ‘x’,‘y’); %建立二元函數(shù),x為第一自變量, y為第二自變量;,fun=@Mfun 定義一個(gè)函數(shù)句柄,這里Mfun是函數(shù)的M文件表達(dá)方式,x0為向量[a, b]時(shí),返回在[a,b]中的零點(diǎn),3.函數(shù)零點(diǎn),x=fzero(Fun, x0) 返回一元函數(shù)Fun的一 個(gè)零點(diǎn),其中Fun為函數(shù)

9、句柄、內(nèi)嵌函數(shù)或字符串表達(dá)方式。,x0為標(biāo)量時(shí), 返回函數(shù)在x0附近的零點(diǎn);,[x,f,h]=fsolve(Fun, x0) x返回一元或多元函數(shù)Fun在x0附近的一個(gè)零點(diǎn),其中x0為迭代初值; f返回Fun在x的函數(shù)值, 應(yīng)該接近0; h返回值如果大于0, 說(shuō)明計(jì)算結(jié)果可靠,否則計(jì)算結(jié)果不可靠。,例3 求函數(shù)y=xsin(x2-x-1)在(-2, -0.1)內(nèi)的 零點(diǎn),解(解法一) 首

10、先定義函數(shù):f=inline('x*sin(x^2-x-1)','x')f = Inline function: f(x) = x*sin(x^2-x-1)fzero(f,[-2,-0.1])??? Error using ==> fzeroThe function values at the interval endpoints must differ in sig

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 眾賞文庫(kù)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論