f9_6.m

matlab的一些例子,对入门的新手有帮助

% f9_6   
% Copyright S. Nakamura, 2001
set(gcf, 'NumberTitle','off','Name', 'Figure 9.6')

clear; clf; hold off
f = [  1 2 3 5 1 2 4 2  ];  m = length(f);
s=1:m;
plot([-2 14], [0 6], '.'); hold on
xlabel('s','fontsize',18); ylabel('f','fontsize',18);  
plot(s,f,'o','linewidth',2)
for k=1:m
  z=int2str(k+2); sk = s(k); fk=f(k); text(sk+0.2,fk-0.2,z,'fontsize',18)
end
t = 0:0.1:1; t2=t.^2; t3=t.^3;
for i=2:m-2
   fb = 1/6*((1-t).^3*f(i-1)+(3*t3-6*t2+ 4)*f(i) + ...
        (-3*t3+3*t2 + 3*t + 1)*f(i+1) + t3*f(i+2) );
   plot(s(i)+t,fb,'linewidth',2)
end
%text (-1,6.3, ' With 3 repeating points at edges', ...
% 'FontSize',[24])
%title (' With no repeating points at the edges')
%print ../fig9d50.ps