examp03.m

来自「ADSP TOOLBOX: Version 2.0 and gui m-file」· M 代码 · 共 9 行

clc,echo on 
%EXAMPLE 03
ts=0.05;t=0:ts:6;l=length(t);          % Choose time axis with ts=0.05
x=10*exp(-t).*sin(2*t);                % function x(t)
z=diff(x)/ts;                          % numerical derivative. NOTE: ts=0.05
ze=10*exp(-t).*(2*cos(2*t)-sin(2*t));  % Exact derivative
plot(t,x,t(2:l),z,t,ze),grid           % Plot numerical and exact result
echo off %end of example