Question

matlab 程序问题

已知向量x=[1 2 3 4 5],y=[7 11 13 15 19].试分别用线性插值法和三次多项式插值法求x=1.3，2.4，3.7，4.8时所应对的y值为多少。

Answer 1

不知道你的程序是要什么样的结果，单纯从代码来看，就是全角和半角的问题。我已经修改如下：

建立M文件eq1．m

     function dy=eq1(x,y)

                  dy=zeros(2,1);

                  dy(1)=y(2);

                  dy(2)=1/5*sqrt(1+y(1)^2)/(1-x);

取x0=0，xf=0．9999，建立主程序如下： 

    x0=0;

       xf=0.9999;

      [x,y]=ode15s('eq1',[x0 xf],[0 0]);

      plot(x,y(:,1),'b.') 

      hold on

      y=0:0.01:2;

       plot(1,y,'b*')

得到的图像是

Answer 2

x=[1 2 3 4 5];
y=[7 11 13 15 19];
x1=[1.3,2.4,3.7,4.8];
y1=interp1(x,y,x1)%线性
y2=interp1(x,y,x1,'cubic')%三次多项式
结果：
y1 =

8.2000 11.8000 14.4000 18.2000

y2 =

8.4310 11.8960 14.3020 18.0293

Answer 3

设z=[1.3 2.4 3.7 4.8]
则执行以下命令：
z=[1.3 2.4 3.7 4.8];
x=[1 2 3 4 5];
y=[7 11 13 15 19];
yy=interp1(x,y,z,'linear') '线性插值
yyy=interp1(x,y,z,'cubic') '三次多项式插值

Answer 4