matlab實驗二次插值法

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

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

2、割搜索計算%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í)行結果：xmin3=0

4、.6621fmin3=0.6147一、實驗小結通過本實驗了解了了matlab的基本操作方法，了解二次插值法的原理與基本運用