exm040622_1.m

MATLAB 6.5 从基础讲解其功能

y=inline('sin(t)^2*exp(-a*t)-b*abs(t)','t','a','b');				%<1>  

a=0.1;b=0.5;t=-10:0.01:10;	%对自变量采样，采样步长不宜太大。
y_char=vectorize(y);		%为避免循环，把y改写成适合数组运算形式。	 <3>
Y=feval(y_char,t,a,b);		%在采样点上计算函数值。
clf,plot(t,Y,'r');hold on,plot(t,zeros(size(t)),'k');	%画坐标横轴
xlabel('t');ylabel('y(t)'),hold off  
 

zoom on						%在MATLAB指令窗中运行，获局部放大图
[tt,yy]=ginput(5);zoom off	%在MATLAB指令窗中运行，用鼠标获5个零点猜测值。
tt			%显示所得零点初始猜测值（该指令可在Notebook中运行）。  

[t4,y4,exitflag]=fzero(y,tt(4),[],a,b)								%<11>  

[t3,y3,exitflag]=fzero(y,tt(3),[],a,b)  

op=optimset('fzero')		%提取fzero所采用的options缺省设置  

op=optimset('tolx',0.01);	%把终止计算的相对误差阈值设置得较大
op.TolX						%观察新设置值。注意TolX字母的大小写。  
[t4n,y4n,exitflag]=fzero(y,tt(4),op,a,b)	%采用新的op设置参数。