概述
1、时域响应
当一个系统是线性时不变系统时,则输出可由输入与单位冲激响应卷积来表示,连续LTI系统响应为:
离散LTI系统响应为:
然后MATLAB中提供了卷积函数conv。
clear
clc
dt=input('输入时间间隔 dt=');
x=2*ones(1,fix(10/dt));%fix函数为朝0方向最近值,取整数
h=exp(-0.5*[0:fix(10/dt)]*dt);
y=conv(x,h);
t=dt*([1:length(y)]-1);
plot(t,y)
grid
clear
clc
x=2*ones(1,10);
n=[0:14];
h=2.^n;
y=conv(x,h);
stem(y);
xlabel('n');
ylabel('y[n]');
2、连续LTI系统的零输入响应函数initial
格式为:[y,t,x]=initial(a,b,c,d,x0)
计算出连续时间LTI系统由于初始状态x0所引起的零输入响应y,x为状态记录,t为仿真所用的时间间隔。
clear
clc
a=[-0.6 -0.8;0.9 0];
b=[1;0];
c=[2 7];
d=[0];
x0=[1;0];
t0=0:0.1:20;
initial(a,b,c,d,x0,t0);
grid
title('LTI系统零输入响应曲线');
3、离散系统的零输入响应函数dinitial
格式为:[y,x,n]=dinitial(a,b,c,d,x0)
计算出离散时间LTI系统由于初始状态x0所引起的零输入响应y,x为状态响应,n为仿真所用的点数。
clear
clc
a=[-0.6 -0.8;0.9 0];
b=[1;0];
c=[2 7];
d=[0];
x0=[1;0];
dinitial(a,b,c,d,x0);
title('离散系统零输入响应曲线');
4、滤波函数filter
从频域角度上来说,无论是连续或离散LTI系统,系统对于输入信号的响应,本质上对输入信号的频谱进行不同选择处理过程,这个过程称为滤波。
格式:y=filter(B,A,x)
对向量x中的数据进行滤波处理,即求解差分方程,产生输出序列向量y。
clear
clc
B=1;
A=[1,-0.5];
n=0:31;
x=0.8.^n;
y=filter(B,A,x);
subplot(121);stem(x);
subplot(122);stem(y);
最后
以上就是可靠砖头为你收集整理的离散信号及离散系统的matlab编程实现_【信号处理工具箱】— 线性时不变系统...的全部内容,希望文章能够帮你解决离散信号及离散系统的matlab编程实现_【信号处理工具箱】— 线性时不变系统...所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
发表评论 取消回复