11.2 Matlab 插值拟合小作业

目录

主要内容

定义的拉格朗日函数

执行的内容

原始数据

---

主要内容

10组dat文件，任选一组数据，运用拉格朗日插值法，拟合函数曲线。

定义的拉格朗日函数

function y0 = Lagrange(x,y,x0)
n = length(x); m = length(x0);
for i = 1:m
  z = x(i); 
  s = 0.0;
    for k = 1:n
        p = 1.0;
        for j= 1:n
            if j ~= k  
                p = p * (z - x(j)) / (x(k)-x(j));
            end
        end 
        s = p * y(k) + s;  
    end
    y0(i) = s
end
end

执行的内容

x = A(:,1); %A是将dat文件导入数据后，从新定义的矩阵名称
y = A(:,2);
xp=-80:2:640; %以1度为间隔，没写出来，所以写成了2为间隔，不知道有没有大神教一下
yp=Lagrange(x,y,xp);

plot(x,y,'*',xp,yp,'b','linewidth',1.5);
xlabel('转角度数（°）');
ylabel('缸压(Mpa）');
legend('原始离散数据','拉格朗日插值');

原始数据

欢迎小可爱们一起来练习

链接：https://pan.baidu.com/s/1UDWSi_v_W4xT__qxWvNZQQ 
提取码：fr3p