Matlab---基本

连续信号

b = 5;
a = 0.8;
t = 0:0.001:5;
x = b*exp(-a*t).*sin(pi*t);
plot(t,x);

离散信号

c = 2;
d = 0.8;
k = -5:5;
y = c*d.^k; 
stem(k,y);

求零状态响应

系统的微分方程为
y”(t) + 2y’(t) + 77y(t) = f(t)
在t≥0接入激励f(t)=10sin(2πt)，求零状态响应。

sys = tf([1],[1 2 77]);
t = 0:0.01:5;
f = 10*sin(2*pi*t);
y = lsim(sys, f, t);
plot(t,y);
xlabel('Time(sec)')
ylabel('y(t)')

系统的冲激响应和阶跃响应

求以下系统的冲激响应和阶跃响应。
7y”(t) + 4y’(t) + 6y(t) = f’(t)+ f(t)

a=[7 4 6]; 
b=[1 1]; 
subplot(2,1,1)
impulse(b,a); 
subplot(2,1,2)
step(b,a); 

卷积

t11=0;t12=1;t21=0; t22=2;
t1 = t11:0.001:t12;
ft1= 2*rectpuls(t1-0.5,1);    %幅值为1，宽度为1，t1=0.5点对称的方波信号
t2 = t21:0.001:t22;
ft2 = t2;
t3 = t11+t21:0.001:t12+t22;
ft3 = conv(ft1,ft2);
ft3 = ft3*0.001;
plot(t3,ft3)
title('ft1(t)*ft2(t)')

离散系统的零状态响应

R = 51;                      %输入信号长度
d = rand(1,R)-0.5;           %产生离散随机数
k = 0:R-1;
s = 2*k.*(0.9.^k);
f = s+d;
figure(1);stem(k,f,'.');      %显示加噪信号
M=5;b=ones(M,1)/M;            %f前系数均为1/M
a=1;
y=filter(b,a,f);            %求零状态响应
figure(2);
stem(k,y,'.');              %显示平滑滤波结果