例6-1(1).m

来自「对初学matlab的有很大帮助 认真学应该能学好」· M 代码 · 共 18 行

%设定自变量数组t
t=[0:0.01:1.5];   
%固定频率   					 
w=4*sqrt(3);
%使用数组运算式       					 
y= sqrt(3)/3*exp(-2*t).*sin(w*t+pi/6);  
%绘制曲线并且加上坐标网格	 
subplot(2,1,1),plot(t,y),grid	
%加标注				
title('绘图示例'),xlabel('时间t'),ylabel('y(t)')  
%求导数	 
Dy=diff(y)/0.01;     
%求导后数组长度减去1                    
s=length(t)-1;           		          	
subplot(2,1,2),plot(t(1:s),Dy),grid 
%加标注			  
ylabel('Dy(t)')