例6-1(2).m

来自「matlab基础与提高」· M 代码 · 共 12 行

dt=0.01;
t=[0:dt:1.5];
w=4*sqrt(3);
y=e3b(t,w);
subplot(2,1,1),plot(t,y),grid;
s=length(t)-1
Dy=diff(y)/dt;
subplot(2,1,2),plot(t(1:s),Dy),grid;
%标注的程序略去，另外建立一个子函数ex1bs.m，其内容为：
function xvalues=ex51bs(tvalues,w)
xvalues= sqrt(3)/3*exp(-2* tvalues).*sin(w* tvalues +pi/6);