exa10_9_1.m

来自「有关matlab r2007的实用教程」· M 代码 · 共 19 行

%设定自变量数组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('绘图示例'),ylabel('y(t)')  
%求导数
Dy=diff(y)/0.01;
%求导后数组长度减去1
s=length(t)-1;
subplot(2,1,2),plot(t(1:s),Dy),grid 
%加标注
xlabel('t /s')
ylabel('Dy(t)')