exm02131_2.m

来自「matlab mfiles实例分析」· M 代码 · 共 9 行

t=-2*pi:pi/10:2*pi;	%该自变量数组中，存在0值。
y=sin(t)./t;			%在t=0处，按IEEE规则，计算将产生NaN
tt=t+(t==0)*eps;		%逻辑数组参与运算，使0元素被一个"机器零"小数代替。
yy=sin(tt)./tt;		%用数值可算的sin(eps)/eps近似替代sin(0)/0极限.
subplot(1,2,1),plot(t,y),axis([-7,7,-0.5,1.2]),
xlabel('t'),ylabel('y'),title('残缺图形')
subplot(1,2,2),plot(tt,yy),axis([-7,7,-0.5,1.2])
xlabel('t'),ylabel('yy'),title('正确图形')