查看: 1373|回复: 2

三点二次插值程序问题！

字体大小: 正常放大

wujiang1994

1 主题	12 听众	68 积分

升级 66.32%

TA的每日心情

	开心 2015-2-8 09:10

签到天数: 9 天

[LV.3]偶尔看看II

自我介绍: 学生

电梯直达

1^#

发表于 2015-1-9 16:12 |只看该作者 |正序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

做一个关于三点二次插值的程序设计，却总是出错，目标方程是f=t^3-3*t+2，初始点是t1=-1;t2=0;t3=2;
源代码
function sdeccz
syms t  ;
t1=-1;t2=0;t3=2;
f=t^3-3*t+2;
k=0;
u0=0.5*{f(t1)*(t2^2-t3^2)+f(t2)*(t3^2-t1^2)+f(t3)*(t1^2-t2^2)}/{(t2-t3)*f(t1)+(t3-t1)*f(t2)+(t1-t2)*f(t3)};
u=u0;
while(abs(u-t2)<0.05)
f1=subs(f,t,t1);
f2=subs(f,t,t2);
f3=subs(f,t,t3);
u=0.5*{f1*(t2^2-t3^2)+f2*(t3^2-t1^2)+f3*(t1^2-t2^2)}/{(t2-t3)*f1+(t3-t1)*f2+(t1-t2)*f3};
fu=subs(f,x,u);
      if u>t2
         if fu>=f2
               t2=u;
               f2=fu;
         else
               t1=t2;
               f1=f2;
               t2=u;
               f2=fu;
         end
      else
         if fu>=f2
               t1=u;
               f1=fu;
               t2=t1;
               f2=f1;
         else
               t2=u;
               f2=fu;
               t3=t2;
               f3=f2;
         end
      end
      k=k+1
      u;
      fu;
end
流程图：

liuchen.JPG (61.5 KB, 下载次数: 157)