exm040104.m

来自「matlab电子书籍」· M 代码 · 共 20 行

clear
d=pi/8;					
t=0:d:pi/2;				
y=0.2+sin(t);				
s=sum(y);					
s_sa=d*s;					
s_ta=d*trapz(y);			
disp(['sum求得积分',blanks(3),'trapz求得积分'])
disp([s_sa, s_ta])
t2=[t,t(end)+d];					
y2=[y,nan];						
stairs(t2,y2,':k')				
hold on
plot(t,y,'r','LineWidth',3)	
h=stem(t,y,'LineWidth',2);		
set(h(1),'MarkerSize',10)
axis([0,pi/2+d,0,1.5])			
hold off
shg