157阶吧 关注:38,858贴子:4,080,917
  • 10回复贴,共1

急!!!!!求计算方法作业,明天就要交了,完全不会啊,题目如下

只看楼主收藏回复

插值方法的应用



1楼2011-10-27 22:29回复
    求求各位大神了,给跪了


    2楼2011-10-27 22:30
    回复
      2025-08-26 23:10:48
      广告
      不感兴趣
      开通SVIP免广告
      看不懂,跪了


      3楼2011-10-27 22:36
      回复
        用埃尔米特插值多项式的方法做的吧、、、、


        4楼2011-10-27 22:36
        回复
          用matlb实现计算
          方法步骤
          1:建立m脚本,其代码如下:
          function [f,f0] = Hermite(x,y,y_1,x0)
          %输入x,y,y_1,x0,4个参数
          %求已知数据点的埃尔米特插值多项式
          %已知数据点的x坐标向量:x
          %已知数据点的x坐标向量:y
          %已知数据点的x坐标向量:y_1
          %插值点的x坐标:x0
          %求得的埃尔米特插值多项式:f
          %x0处的插值:f0
          syms t;%定义符号变量t
          f = 0.0;
          if(length(x) == length(y))
          if(length(y) == length(y_1))
          n = length(x);
          else
          disp('y和y的导数的维数不相等!');
          return; %提起终止程序运行
          end
          else
          disp('x和y的维数不相等!');
          return;
          end %检错
          for i=1:n %赋初值
          h = 1.0;
          a = 0.0;
          for j=1:n
          if( j ~= i)
          h = h*(t-x(j))^2/((x(i)-x(j))^2);
          a = a + 1/(x(i)-x(j));
          end
          end
          f = f + h*((x(i)-t)*(2*a*y(i)-y_1(i))+y(i));%计算插值多项式
          end
          f0=subs(f,'t',x0); %计算插值点的函数值2:输入对应数据求解可得答案。
          方法步骤
          1:建立m脚本,其代码如下:
          function [f,f0] = Hermite(x,y,y_1,x0)
          %输入x,y,y_1,x0,4个参数
          %求已知数据点的埃尔米特插值多项式
          %已知数据点的x坐标向量:x
          %已知数据点的x坐标向量:y
          %已知数据点的x坐标向量:y_1
          %插值点的x坐标:x0
          %求得的埃尔米特插值多项式:f
          %x0处的插值:f0
          syms t;%定义符号变量t
          f = 0.0;
          if(length(x) == length(y))
          if(length(y) == length(y_1))
          n = length(x);
          else
          disp('y和y的导数的维数不相等!');
          return; %提起终止程序运行
          end
          else
          disp('x和y的维数不相等!');
          return;
          end %检错
          for i=1:n %赋初值
          h = 1.0;
          a = 0.0;
          for j=1:n
          if( j ~= i)
          h = h*(t-x(j))^2/((x(i)-x(j))^2);
          a = a + 1/(x(i)-x(j));
          end
          end
          f = f + h*((x(i)-t)*(2*a*y(i)-y_1(i))+y(i));%计算插值多项式
          end
          f0=subs(f,'t',x0); %计算插值点的函数值
          2:输入对应数据求解可得答案。


          5楼2011-10-27 22:37
          回复
            复制了两遍 不好意思


            6楼2011-10-27 22:38
            回复
              完全不懂。。。。还是谢谢你


              7楼2011-10-27 22:40
              回复
                额。。这个 我建议你还是去图书馆借本关于matlb使用的参考书吧。。。那里面会详细介绍各种插值计算的方法实现。。。


                8楼2011-10-27 22:49
                回复
                  2025-08-26 23:04:48
                  广告
                  不感兴趣
                  开通SVIP免广告
                  留扣扣吧 给你个我们以前的作业 熄灯之前~~~~


                  9楼2011-10-27 22:52
                  回复
                    472045480


                    10楼2011-10-27 22:55
                    回复
                      回复9楼:
                      1458409707这位大哥给我也整一份,真心感谢了


                      来自掌上百度11楼2011-10-27 23:57
                      回复