graf1.m

该代码为原创

x=[2.043254,1.989856,1.954187,1.930405,1.887477,1.865163,1.845817,1.829529,1.809862,1.794058,1.784073,1.759142,1.743078,1.723457,1.714270,1.699557,1.685074,1.680042,1.667581,1.661690,1.653694,1.639933,1.633348,1.617274,1.599312,1.595121,1.592699,1.570442,1.544275,1.530192,1.529304,1.393409,1.417193,1.380351,1.254166,1.150996,1.164059,1.152916,1.084194,1.161304,0.818581,1.096061,0.739569];
y=[0.301030,0.477121,0.602060,0.698970,0.903090,1.000000,1.113943,1.204120,1.301030,1.397940,1.491362,1.602060,1.698970,1.799341,1.903090,2.000000,2.096910,2.198657,2.301030,2.399674,2.499687,2.599883,2.699838,2.800029,2.899821,3.000000,3.099681,3.200029,3.299943,3.400020,3.499962,3.599992,3.700011,3.799961,3.899985,4.000000,4.099991,4.200002,4.300008,4.400002,4.500003,4.600003,4.700002];
function yy=pline(x,y)
x=x(1:44);
y=y(1:44);
x=reshape(x,44,1);
y=reshape(y,44,1);
A=[x ones(44,1)]
b=y;
b=A'*A;
b=A'*b;
yy=B\b;
yy=yy';
p=pline(M,S)
%S1=polyval(p,M)
%plot(M,S1)
%hold on
%plot(M,S,'.')