exm045_1.m

来自「北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压」· M 代码 · 共 26 行

%exm045_1.m
%（1）使用字符串表示被处理函数
P1=0.1;P2=0.5;		%按泛函指令要求，这里参数必须用P1,P2,...表示
y_C='sin(x).^2.*exp(-P1*x)-P2*abs(x)';	%这里自变量必须用x表示	<2>  

%（2）作图法观察函数零点分布
x=-10:0.01:10;	    %对自变量采样，采样步长不宜太大。
Y=eval(y_C);		%在采样点上计算函数值。
clf,shg
plot(x,Y,'r');hold on,plot(x,zeros(size(x)),'k');	%画坐标横轴
xlabel('t');ylabel('y(t)'),hold off  

disp('按任意键继续！')
disp('下面将通过鼠标操作从图上获取 5 个零点的近似值。') 
pause

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

disp('按任意键继续！')
disp('下面求靠近tt(4) 的精确零点。') 
pause

[t4,y4,exitflag]=fzero(y_C,tt(4),[],P1,P2)