exm04111_1.m

matlab mfiles实例分析

t=linspace(0,5,100);y=1-cos(3*t).*exp(-t);

plot(t,y,'b');grid;hold on,plot(t,0.95*ones(size(t)),'r');hold off
t_ginput=ginput(1)			%它的第一个元素就是穿越时刻。  

it=min(find(y>0.95));		%确定原始数据中第一次使y>0.95的元素“下标”
T=(it-3):(it+3);		%在第一穿越点两侧，各取3点，作为插值“基准”		<2>
t_nearst=interp1(y(T),t(T),0.95,'nearst');						%	<3>
t_linear=interp1(y(T),t(T),0.95);								%	<4>
t_cubic=interp1(y(T),t(T),0.95,'cubic');						%	<5>
t_spline=interp1(y(T),t(T),0.95,'spline');						%	<6>
disp(['  t_nearst  t_linear   t_cubic   t_spline'])
disp([t_nearst  t_linear   t_cubic   t_spline])  

t_zero=fzero('1-cos(3*x).*exp(-x)-0.95',0.5)