MATLAB在信号与系统中的应用(精编).pptVIP

MATLAB在信号与系统中的应用 7.1连续信号和系统 7.1.1连续信号的表示和运算 1.向量表示法： t=-10:1.5:10 F=sin(t)./t plot(t,f) 程序运行结果： 2.符号运算表示法 如果信号可以用一个符号表达式来表示它，则可用ezplot命令绘制出信号的波形。 例： f=sin(pi/4*t); ezplot(f,[-16,16]) 用MATLAB实现连续信号的可视化： （1）单位阶跃信号 单位阶跃函数的数学表达式 MATLAB程序为： t0=0;tf=5;dt=0.05;t1=1;%定义时间的起始值，步长 t=[t0:dt:tf];st=length(t); n1=floor((t1-t0)/dt);%求t1对应的样本序号 x1=[zeros(1,n1),ones(1,st-n1)]; %产生阶跃信号 stairs(t,x1),grid on%绘图 axis([0,5,0,1.1])%为使方波顶部避开图框，改变图框坐 标 阶跃信号 (2)单位冲激函数：可以看作是宽度为△（程序中用dt表示)，幅度为1/ △的矩形脉冲，即 axis([0,5,0,22]) %使脉冲顶部避开图框，改变图框坐标 程序运行结果如下： 7.1.2连续系统的冲激响应、阶跃响应及MATLAB实现 描述连续系统的微分方程为： impulse()调用格式： impulse(b,a) impulse(b,a,t) impulse(b,a,t1:p:t2):t1-t2时间范围内，以时间 间隔p为均匀采样冲激响应的时域波形。 y=impulse(b,a,t1:p:t2)：该调用格式并不绘制出 系统的阶跃冲激响应的波形，而是求出其数值解。 2.step()函数 step()函数将绘出由向量a和b表示的连续系统的阶 跃响应g(t)在制定时间范围内的波形图，并能求出 其数值解。 其调用格式和impulse()函数一样。 例：已知描述某连续系统的微分方程为 7.2傅里叶分析 傅里叶变换的定义 傅里叶反变换的定义 注意：在调用这两个函数之前，要用sym命令对所有用到的变量进行说明。对于fourier(）中的函数F，也要用符号定义符sym将它们说明为符号表达式。 例：设 ，试画出f(t)及其幅频图 MATLAB程序： syms t x; x=1/2*exp(-2*t)*sym(heaviside(t)); F=fourier(x); subplot(2,1,1); ezplot(x); subplot(2,1,2); ezplot(abs(F)) 7.3 离散信号和系统 一、离散时间信号（stem()) 1.单位序列（单位脉冲序列） 二、离散系统的单位响应 当输入为单位序列是离散系统产生的零状态响应称 为系统的单位响应，MATLAB中求离散系统的单位响 应，并绘制其时域波形的函数impz()。 设描述离散系统的差分方程为 例：已知描述某离散系统的差分方程如下： 2y(k)-2y(k-1)+ (k-2)=f(k)+ 3f(k-1)+ 2f(k-2) 用MATLAB绘制出系统在0～100时间范围内单位响 应波形。 MATLAB程序： a=[2 -2 1]; b=[1 3 2]; impz(b,a) 7.4利用MATLAB求连续系统的响应 一、线性定常连续系统可用如下所示的线性常系数 为微分方程来描述 解： MATLAB程序： a=[1 -0.25 0.5]; b=[1 1]; t=0:20; x=(1/2).^t; y=filter(b,a,x); subplot(2,1,1) stem(t,x) title(输入序列) subplot(2,1,2) stem(t,y) title(输出序列) 7.5用MATLAB实现拉普拉斯逆变换和逆Z变换 7.5.1拉普拉斯逆变换 用MATLAB函数来实现连续时间信号的拉普拉斯逆 变换。 K：为包含F（s）部分分式展开的多项式项的系数cj的行 向量,若M≤N,则K为空阵。 例：已知连续信号的拉普拉斯变换为 试用MATLAB求其拉普拉斯逆变换f(t)。 MATLAB程序： a=[1 0 4 0]; b=[2 4]; [r p k]=residue(b,a) 程序运行结果： r = -0.5000 - 0.5000i -0.5000 + 0.5000i 1.0000 p = 0 + 2.0000i 0 - 2.0000i 0 k = [ ] abs(r) ans = 0.7071 0.7071 1.0000 angle