matlab实验二次插值法

《matlab实验二次插值法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验报告实验名称：二次插值法院（系）：机电学院专业班级：机械制造及其自动化姓名：学号：2013年5月13日实验一：二次插值法实验日期：2013年5月13日一、实验目的了解MATLAB的基本运用了解MATLB在优化中的使用二、实验原理利用目标函数飞f（x）上极小点附近若干点的函数值构造一个与原函数相近的插函数平p（x），以插值函数的最优解作为目标函数的最优解的近似解。以插值函数平p（x）的极小点作为新的中间插入点，多次迭代，逐步缩小区间。三、实验内容二次插值法程序：function[xmin,fmin]=erf(f,a0,b0,epsilon)%二分法一维

2、割搜索计算%a0，b0初始区间%fx2测试函数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函数程序：function[zhi]=fx3(x)%3²âÊÔº¯Êýzhi=exp(x)-2*x;end调用执行程序：[xmin3,fmin3]=erf(@fx3,0,2,0.1)执行结果：xmin3=0

4、.6621fmin3=0.6147一、实验小结通过本实验了解了了matlab的基本操作方法，了解二次插值法的原理与基本运用