3楼: Originally posted by 杨幸遇 at 2014-05-28 08:54:54
万分感谢！金币送上！...

function Lag_interp
clear all;clc
Ti=[12 33 66];
T=[ 0          20          40          60          80         100];
p=[999.9        998.2        992.2        983.2        971.8        958.4];

for i=1:length(Ti);
    p_answer(i)=polyinterp(T,p,Ti(i));
end
Result=[Ti' p_answer']
   


function v=polyinterp(x,y,u)
n=length(x);
v=zeros(size(u));
for k=1:n;
    w=ones(size(u));
    for j=[1:k-1 k+1:n];
        w=(u-x(j))./(x(k)-x(j)).*w;
    end
v=v+w*y(k);
end

2楼: Originally posted by 月只蓝 at 2014-05-27 18:33:19
function Lag_interp
clear all;clc
Ti=;
T=;
p=;

for i=1:length(Ti);
    p_answer(i)=polyinterp(T,p,Ti(i));
end
Result=
   


function v=polyinterp(x,y,u)
n=length(x);
v=zeros(size(u)) ...