graf1.asv

该代码为原创

x=[5.226278,5.218761,5.211189,5.205401,5.190243,5.185260,5.179887,5.159825,5.150973,5.137968,5.125573,5.086965,5.070354,5.061931,5.033913,4.989862,4.978391,4.945518,4.873425,4.848972,4.827147,4.772400,4.731393,4.707590,4.622537,4.613050,4.563600,4.509898,4.466851,4.454617,4.415051,4.384782,4.391827,4.340787,4.301199,4.342952,4.275477,4.034977,4.166478,4.201476,4.033207,3.728374,3.159405,3.168083]
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,4.800002]
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,'.')