不知道这个怎么解决,可以详细讲一下吗,本人不太懂
function y=lglrcz(x0,y0,x)
n=length(x0);
m=length(x);
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-x0(j))/(x0(k)-x0(j));
end
end
s=p*y0(k)+s;
end
y(i)=s;
end
format long
clc;clf;
x=[0,4.74,9.05,19,38,76,95,114,133,152,171,190];
y=[0,5.23,8.1,11.97,16.15,16.34,14.63,12.16,6.69,7.03,3.99,0];
x0=0:0.1:190;
lglr=lglrcz(x,y,x0);
laglr=lglrcz(x,y,12)
fdxx=interp1(x,y,x0);
fendxx=interp1(x,y,12)
scyt=interp1(x,y,x0,'spline');
sancyt=interp1(x,y,12,'spline')
plot(x,y,'*',x0,lglr,'r',x0,fdxx,'g'x0,scyt,'b')
legend('观测点','拉格朗日插值曲线','分段线性插值曲线','三次样条插值曲线')
plot(x,y,'*',x0,lglr,'r',x0,fdxx,'g',x0,scyt,'b')
少了个逗号