matlab實(shí)驗(yàn)二次插值法

《matlab實(shí)驗(yàn)二次插值法》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)名稱：二次插值法院（系）：機(jī)電學(xué)院專業(yè)班級：機(jī)械制造及其自動(dòng)化姓名：學(xué)號：2013年5月13日實(shí)驗(yàn)一：二次插值法實(shí)驗(yàn)日期：2013年5月13日一、實(shí)驗(yàn)?zāi)康牧私釳ATLAB的基本運(yùn)用了解MATLB在優(yōu)化中的使用二、實(shí)驗(yàn)原理利用目標(biāo)函數(shù)飛f（x）上極小點(diǎn)附近若干點(diǎn)的函數(shù)值構(gòu)造一個(gè)與原函數(shù)相近的插函數(shù)平p（x），以插值函數(shù)的最優(yōu)解作為目標(biāo)函數(shù)的最優(yōu)解的近似解。以插值函數(shù)平p（x）的極小點(diǎn)作為新的中間插入點(diǎn)，多次迭代，逐步縮小區(qū)間。三、實(shí)驗(yàn)內(nèi)容二次插值法程序：function[xmin,fmin]=erf(f,a0,b0,epsilon)%二分法一維

2、割搜索計(jì)算%a0，b0初始區(qū)間%fx2測試函數(shù)a=a0;%初始化b=b0;x1=a;f1=f(x1);x3=b;f3=f(x3);x2=0.5*(x1+x3);f2=f(x2);c1=(f3-f1)/(x3-x1);c2=((f2-f1)/(x2-x1)-c1)/(x2-x3);xp=0.5*(x1+x3-c1/c2);fp=f(xp);while(abs(xp-x2)>=epsilon)%迭代求解ifx2fpf1=f2;x1=x2;x2=xp;f2=fp;elsef3=fp;x3=xp;endelseiff2>fpf3=f2;x3=x2

3、;f2=fp;x2=xp;elsef1=fp;x2=xp;endendc1=(f3-f1)/(x3-x1);c2=((f2-f1)/(x2-x1)-c1)/(x2-x3);xp=0.5*(x1+x3-c1/c2);fp=f(xp);endiff2>fp%輸出xmin=xp;fmin=f(xp);elsexmin=x2;fmin=f(x2);endend函數(shù)程序：function[zhi]=fx3(x)%32aê?oˉêyzhi=exp(x)-2*x;end調(diào)用執(zhí)行程序：[xmin3,fmin3]=erf(@fx3,0,2,0.1)執(zhí)行結(jié)果：xmin3=0

4、.6621fmin3=0.6147一、實(shí)驗(yàn)小結(jié)通過本實(shí)驗(yàn)了解了了matlab的基本操作方法，了解二次插值法的原理與基本運(yùn)用