X=[120:130];
temp=load('~mchung/public/html/teaching/data/unemployment.data');
unemployment=temp(:,3);
>> Y=unemployment(X)
Y =
6.0000
6.3000
6.3000
6.3000
6.9000
7.5000
7.6000
7.8000
7.7000
7.5000
7.5000
XX=[120:130];
>> S=csaps(X,Y,1,XX)'
S =
6.0000
6.3000
6.3000
6.3000
6.9000
7.5000
7.6000
7.8000
7.7000
7.5000
7.5000
XX=[120:0.1:130];
S=csaps(X,Y,1,XX)';
plot(X,Y,'o')
hold on
plot(XX,S,'-')
S=csaps(X,Y,0,XX)';
plot(X,Y,'o')
hold on
plot(XX,S,'-')
S=csaps(X,Y,1/2,XX)';
plot(X,Y,'o')
hold on
plot(XX,S,'-')
X=[1:252];
XX=[1:0.2:252];
plot(X,unemployment,'o')
hold on
S=csaps(X,unemployment,0,XX);
plot(XX,S,'-')
S=csaps(X,unemployment,1/256,XX)';
plot(XX,S,'-')
S=csaps(X,unemployment,1/2,XX)';
plot(XX,S,'-')
S=csaps(X,unemployment,1,XX)';
plot(XX,S,'-')