新·第2部分 MATLAB在信号处理中的应用-精.pptVIP

MATLAB在通信中的应用 一、 MATLAB基础知识 二、 MATLAB与信号处理 三、 MATLAB与通信仿真 四、 Simulink应用 第二部分 信号处理MATLAB仿真 信号及其运算表示 一、信号及其运算表示 连续信号 严格来说，MATLAB并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。 在MATLAB可视化绘图中 对于以t为自变量的连续信号，在绘图时统一用plot函数； 而对n为自变量的离散序列，在绘图时统一用stem函数。 单位冲激函数、单位冲激序列 t = -5:0.01:5;y = (t==0);subplot(1,2,1);plot(t, y, r);n = -5:5;x = (n==0);subplot(1,2,2);stem(n, x); 单位阶跃函数、单位阶跃序列 t =-5:0.01:5; y=(t=0); subplot(1,2,1); plot(t, y, r) n = -5:5;x = (n=0);subplot(1,2,2);stem(n, x); 生成上述三种信号 t =-5:0.01:5; subplot(2,2,1); a=2 y1=a.^t plot(t, y1, r) subplot(2,2,2); a=2; theat=pi/3; y2=sin(2*pi*t+theat) plot(t, y2) subplot(2,1,2); w=4; y3=exp((a+j*w)*t); plot(t, y3, ‘y) clc;clf; N=24;M=7;m=3; n1=0:8; x1=(0.8).^n1; %生成指数序列 x2=x1; n2=n1+m; subplot(4,1,1);stem(n1,x1);ylabel(x(n));axis([0 N-1 0 1]); subplot(4,1,2);stem(n2,x2);ylabel(x(n-3));axis([0 N-1 0 1]); n=0:N-1; if Mlength(x1) for i=1:length(x1)-M x1(i)=x1(i)+x1(length(x1)-i) end end xc=x1(mod(n,M)+1); %产生x(n)周期延拓x((n))8 xcm=x1(mod(n-m,M)+1); %产生移位序列x(n-3)周期延拓x((n-3))8 subplot(4,1,3);stem(n,xc);ylabel(x((n))_8);axis([0 N-1 0 max(xc)]); subplot(4,1,4);stem(n,xcm);ylabel(x((n-3))_8);axis([0 N-1 0 max(xcm)]); 已知某线性移不变系统单位抽样相应为 f1(t)=e-t[ε(t)- ε(t-2)]、 f2(t) = [ε(t)- ε(t-3)]，编制一个m文件，绘出f1(t)、f2(t)和f(t)=f1(t)* f2(t)的波形。 p=0.001; k1=0:p:2 f1=exp(-k1) k2=0:p:3; f2=rectpuls(k2-1.5,3); f=conv(f1,f2); %计算序列f1,f2的卷积和f f=f*p; k0=k1(1)+k2(1); %计算序列f非零样值的起点位置 k3=length(f1)+length(f2)-2; %计算卷积和f的非零样值的宽度 k=k0:p:k3*p; %确定卷积和f非零样值的时间向量 subplot(2,2,1); plot(k1,f1); title(f1(t)); xlabel(t); ylabel(f1(t)); axis([0 5 0 1.2]); subplot(2,2,2); plot(k2,f2); title(f2(t)); xlabel(t); ylabel(f2(t)); axis([0 5 0 1.2]); subplot(2,1,2); plot(k,f); %画出卷积f(t)的波形 title(f(t)=f1(t)*f2(t))); xlabel(t); ylabel(f(t)); axis([0 5 0 1.2]); 二、线性时不变系统时域分析 系统描述 系统表示方法转换 时域响应 连续系统 离散系统 3、线性时不变系统时域响应 连续系统 离散系统 冲激响应、阶跃响应 零状态响应、零输入响应、全响应 连续系统的时域响应 连续系统的时域