dogsurve.m

用matlab绘制三次样条曲线。系数由c++程序计算得出。

%first

x=1:0.05:2;
y=3+(x-1)-0.34681*(x-1).^2+0.04681*(x-1).^3;
plot(x,y);
hold on;
x=2:0.05:5;
y=3.7+0.44681*(x-2)-0.20638*(x-2).^2+0.0265552*(x-2).^3;
plot(x,y);
x=5:0.05:6;
y=3.9-0.0744797*(x-5)+0.0326168*(x-5).^2+0.341863*(x-5).^3;
plot(x,y);
x=6:0.05:7;
y=4.2+1.01634*(x-6)+1.05821*(x-6).^2-0.574548*(x-6).^3;
plot(x,y);
x=7:0.05:8;
y=5.7+1.40911*(x-7)-0.665439*(x-7).^2+0.156329*(x-7).^3;
plot(x,y);
x=8:0.05:10;
y=6.6+0.54722*(x-8)-0.19645*(x-8).^2+0.0239201*(x-8).^3;
plot(x,y);
x=10:0.05:13;
y=7.1+0.0484602*(x-10)-0.0529297*(x-10).^2-0.00255607*(x-10).^3;
plot(x,y);
x=13:0.05:17;
y=6.7-0.338131*(x-13)-0.0759343*(x-13).^2+0.00574178*(x-13).^3;
plot(x,y);


%second

x=17:0.05:20;
y=4.5+3*(x-17)-1.10071*(x-17).^2+0.126162*(x-17).^3;
plot(x,y);
x=20:0.05:23;
y=7-0.197875*(x-20)+0.0347502*(x-20).^2-0.0229307*(x-20).^3;
plot(x,y);
x=23:0.05:24;
y=6.1-0.608501*(x-23)-0.171626*(x-23).^2+0.280127*(x-23).^3;
plot(x,y);
x=24:0.05:25;
y=5.6-0.111371*(x-24)+0.668756*(x-24).^2-0.357385*(x-24).^3;
plot(x,y);
x=25:0.05:27;
y=5.8+0.153987*(x-25)-0.403398*(x-25).^2+0.0882022*(x-25).^3;
plot(x,y);
x=27:0.05:27.7;
y=5.2-0.401178*(x-27)+0.125815*(x-27).^2-2.568*(x-27).^3;
plot(x,y);


%third

x=27.7:0.05:28;
y=4.1+0.33(x-27.7)+2.26205*(x-27.7).^2-3.79941*(x-27.7).^3;
plot(x,y);
x=28:0.05:29;
y=4.3+0.661386*(x-28)-1.15743*(x-28).^2+0.29604*(x-28).^3;
plot(x,y);
x=29:0.05:30;
y=4.1-0.765347*(x-29)-0.269307*(x-9).^2-0.0653465*(x-29).^3;
plot(x,y);