信号与系统与MATLAB实现.docxVIP

《信号与系统》实验指导书 本题也可先建立一个画单位阶跃序列 (k k ) 的 M 函数文件，画图时调用。 0 function jyxulie(k1,k2,k0) k=k1:-k0-1; kk=-k0:k2; n=length(k); nn=length(kk) u=zeros(1,n); uu=ones(1,nn); stem(kk,uu,filled) hold on stem(k,u,filled) hold off title( 单位阶跃序列 ) axis([k1 k2 0 1.5]) 保存文件名 jyxulie.m 画图时在命令窗口调用，例： jyxulie(-3,8,0) 例 2-4 实指数序列信号 k f (k) ca ，c、 a 是实数。 建立一个画实指数序列的 M 函数文件，画图时调用。 function dszsu(c,a,k1,k2) %c：指数序列的幅度 %a：指数序列的底数 %k1：绘制序列的起始序号 %k2：绘制序列的终止序号 k=k1:k2; x=c*(a.^k); stem(k,x,filled) hold on plot([k1,k2],[0,0]) hold off 5 k 调用该函数画信号： ) ( ) f1 (k) ( k 4 dszsu(1,5/4,0,40) 3 k ， f 2(k) ( ) (k) 4 波形。 dszsu(1,-3/4,0,40) 2 .离散信号的卷积和 两个有限长序列 f1，f2 卷积可调用 MATLAB 函数 conv，调用格式是 f=conv(f1,f2), f 是卷积结果，但不显示时间序号，可自编一个函数 dconv 给出 f 和 k，并画图。 function [f,k]=dconv(f1,f2,k1,k2) %The function of compute f=f1*f2 % f: 卷积和序列 f(k) 对应的非零样值向量 8 《信号与系统》实验指导书 % k： 序列 f(k) 的对应序号向量 % f1: 序列 f1(k) 非零样值向量 % f2: 序列 f2(k) 的非零样值向量 % k1: 序列 f1(k) 的对应序号向量 % k2: 序列 f2(k) 的对应序号向量 f=conv(f1,f2) %计算序列 f1 与 f2 的卷积和 f k0=k1(1)+k2(1); %计算序列 f 非零样值的起点位置 k3=length(f1)+length(f2)-2; %计算卷积和 f 的非零样值的宽度 k=k0:k0+k3 %确定卷积和 f 非零样值的序号向量 subplot(2,2,1) stem(k1,f1) %在子图 1 绘序列 f1(k) 时域波形图 title(f1(k)) xlabel(k) ylabel(f1(k)) subplot(2,2,2) stem(k2,f2) %在图 2 绘序列 f2(k) 时波形图 title(f1(k)) xlabel(k) ylabel(f2(k)) subplot(2,2,3) stem(k,f); %在子图 3 绘序列 f(k) 的波形图 title(f(k)f1(k) 与 f2(k) 的卷积和 f(k)) xlabel(k) ylabel(f(k)) h=get(gca,position); h(3)=2.5*h(3); set(gca,position,h) %将第三个子图的横坐标范围扩为原来的 2.5 倍 例 2-5 求卷积和， f1 (k ) (k 1) 2 (k) (k 1) f2 (k ) (k 2) (k 1) (k) (k 1) (k 2) f1=[1 2 1]; k1=[-1 0 1]; f2=ones(1,5); k2=-2:2; [f, k]=dconv(f1,f2,k1,k2) 由运行结果知， f 的长度等于 f1 和 f2 长度之和减一， f 的起点是 f1 和 f2 的起点之和， f 的终点是 f1 和 f2 的终点之和。 3． 离散系统的单位响应 MATLAB 提供画系统单位响应函数 impz，调用格式是 impz(b,a) 式中 b 和 a 是表示离散系统的行向量； impz(b,a,n) 式中 b 和 a 是表示离散系统的行向量，时间范围是 0~n； impz(b,a,n1,n2) 时间范围是 n1~n2 ；y=impz(b,a,n1,n2) 由 y 给出数值序列； 9 《信号与系统》实验指导书 % k： 序列 f(k) 的对应序号向量 % f1: 序列 f1(k) 非零样值向量 % f2: 序列 f2(k) 的非零样值向量 % k1: 序列 f1(k) 的对应序号向量 % k2: 序列 f2(k) 的对应序号向量 f=conv(f1,f2) %计算序列 f1 与 f2