我正在尝试重现Dirac Delta功能:
我的代码:
x = -30:1:30;
y = zeros(1,numel(x)); %sets all values initially to zero
y(x==0)= inf; % the point corresponding to x=0 is set to inf
plot(x,y,'d')
axis([-40 40 0 inf])
我的代码产生:
答案 0 :(得分:10)
您可以使用stem
执行此操作,并将其'Marker'
指定为向上箭头...
% Open figure
figure;
% Blue stem plot at x=0, to y=75. Marker style is up arrow
stem(0, 75,'color','b','linewidth',2,'marker','^')
% Add infinity label at x=0, y = 82 = 75 + fontsize/2, where we plotted up to 75
text(0,82,'∞','FontSize',14)
% Set axes limits
xlim([-40,40])
ylim([0,90])
您可以看到output plot here,但请参阅下面的修改以获得改进版本。
注意,当然您应该选择相对于绘图上任何其他数据较大的y值。在这个例子中,我选择了75来粗略匹配你想要的示例图。 MATLAB不能在inf
处绘制一个值,因为无穷大位于y轴的哪个位置?
编辑:您可以在评论中指出由Marco建议的其他'≈'
字符打破y轴。将xlim
和ylim
组合到一个axis
调用中,并更改y轴标记以帮助指示轴中断,我们得到以下结果:
stem(0, 80,'color','b','linewidth',2,'marker','^')
text([-42,0,38], [80,87,80], {'≈','∞','≈'}, 'Fontsize', 14)
axis([-40, 40, 0, 100])
yticks(0:20:60)
答案 1 :(得分:3)
要显示无穷大,您不应将y
设置为无穷大。为此,您可以将y
设置为与轴值成比例的较大值。例如,如果轴类似于[min_x max_x min_y max_y]
,则可以设置y(x==0) = max_y*10
。
在您的情况下,您将拥有:
x = -30:1:30; min_x = min(x) - 10; max_x = max(x) + 10;
y = zeros(1,numel(x));
% compute values of y here
% ...
min_y = min(y) - 10; max_y = max(y) + 10;
y(x==0)= 10 * max_y;
plot(x,y,'d');
axis([min_x max_x min_y max_y]);
答案 2 :(得分:-2)