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,'-')