连续时间信号的表示LTI系统的时域分析
例一:
代码如下:
首先构建阶跃函数:
function f=u(t)
f=(t>=0);
end
试题代码:
t=-10:0.01:10;
y1=2*exp(-1.5.*t).*u(t);
y2=2*sin(2*pi*t+pi/4);
y3=2*(u(t)-u(t-1));
subplot(2,2,2),plot(t,y1,'r'),axis([-2 4 0 3]),xlabel('t'),ylabel('y=2exp(-1.5t)*u(t)'),grid on
subplot(2,2,3),plot(t,y2),axis([-2*pi 2*pi -2 2]),xlabel('t'),ylabel('y=2sin(2pi*t+pi/4'),grid on
subplot(2,2,4),plot(t,y3),axis([-1 2 -1 3]),xlabel('t'),title('分段函数'),grid on
例二:
代码如下:
阶跃函数:
function f=funct(t)
f=u(t+2)-u(t)+(-t+1).*(u(t)-u(t-1));
end
原函数:
function f=funct(t)
f=u(t+2)-u(t)+(-t+1).*(u(t)-u(t-1));
end
试题代码:
t=-10:0.001:10;
subplot(2,3,1),plot(t,funct(t)),xlabel('t'),ylabel('y1'),axis([-4 4 0 1]),title('原函数'),grid on
subplot(2,3,2),plot(t,funct(t-2)),xlabel('t1'),ylabel('y1'),axis([-2 4 0 1]),title('y=f(t-2)'),grid on
subplot(2,3,3),plot(t,funct(3*t)),xlabel('t2'),ylabel('y1'),axis([-2 2 0 1]),title('y=f(3t)'),grid on
subplot(2,3,4),plot(t,funct(-t)),xlabel('t3'),ylabel('y1'),axis([-2 3 0 1]),title('y=f(-t)'),grid on
subplot(2,3,5),plot(t,funct(-3*t-2)),xlabel('t4'),ylabel('y1'),axis([-2 1 0 1]),title('y=f(-3t-2)'),grid on
例题三:
代码如下:
t=[0:0.1:5];
a=[1 4 4];
b=[1 3];
x=exp(-t).*stepfun(t,0);
y1=impulse(b,a,t);
y2=lsim(b,a,x,t);
subplot(2,2,1),plot(t,y1,'-r*'),xlabel('t'),ylabel('h(t)'),title('冲激响应'),grid on,legend('impulse h(t)')
subplot(2,2,2),plot(t,y2),xlabel('t'),ylabel('yf(t)'),title('零状态响应'),grid on,legend('lism yf(t))')
例题四:
代码如下:
dt=0.01; t=-3:dt:3;
f1=2*(heaviside(t+1)-heaviside(t-1));
f2=heaviside(t+2)-heaviside(t-2);
f=conv(f1,f2)*dt;
n=length(f);
tt=(0:n-1)*dt-2;
subplot(221),plot(t,f1),grid on;
axis([-3,3,-0.2,1.2]);
title('f1(t)');
xlabel('t')
subplot(222);
plot(t,f2);
grid on;
axis([-3,3,-0.2,1.2]);
title('f2(t)');
xlabel('t')
subplot(212);
plot(tt,f);
grid on;
title('f(t)=f1(t)*f2(t)');
xlabel('t')