文档介绍：线性系统时域分析的Matlab实现
1
信号的表示及可视化
信号的时域运算及变换
线性系统的时域分析
连续系统的冲激响应、阶跃响应
离散系统的单位响应
连续系统的响应
2
向量表示法
>> t = -10::10;
>> f = sin(t)./t;
>> plot(t,f);
>> title('f(t)=Sa(t)')
>> xlabel('t')
>> axis([-10,10,-,])
3
符号向量表示法
>> f = sym('sin(pi/4*t)');
>> ezplot(f,[-16,16])
4
单位阶跃信号
t = -1::3;
f = heaviside(t);
plot(t,f)
axis([-1 3 - ])
5
f = sym('heaviside(t+3)-2*heaviside(t)');
ezplot(f,[-5,4]);
hold on,plot([0,0],[-1,1]);
axis([-5,4,- ]);hold off
6
t = 0::3;
f = exp(-*t);
plot(t,f);
axis([0,3,0,])
7
单位冲击信号
function chongji(t1,t2,t0)
dt = ;
t = t1:dt:t2;
n = length(t);
x = zeros(1,n);
x(1,(-t0-t1)/dt+1) = 1/dt;
stairs(t,x);
axis([t1,t2,0,])
chongji(-1,5,0)
8
离散时间信号
>> k=[-3,-2,1,0,1,2,3];
>> f=[1,2,-1,3,2,4,-1];
>> stem(k,f,'filled'),axis([-4,4,-,])
9
单位采样序列
function [x,n]=mm601_1(n0,n1,n2)
%x(n) = delta(n-n0); n1<= n <=n2
%-------------------------------------------
%[x,n] = mm601_1(n0,n1,n2)
n = [n1:n2];
x = [(n-n0) == 0 ];
stem(n,x)
xlabel('n')
ylabel('x(n)')
title('Delta Sequence');
grid
[x,n]=mm601_1(0,-10,10);
10