实验四数字信号源的设计(精选).docVIP

实验四、数字信号源的设计 实验目的 学习和掌握MATLAB程序设计基本方法。 学习并掌握数字信号处理中常用时域离散序列信号的设计。 了解MATLAB信号处理工具箱提供的常用信号使用方法。 实验内容及步骤 （1）通过编程实现以下信号序列： 单位采样序列； 参考代码： function y impseq n1,n2,n0 if n0 n1 | n0 n2 | n1 n2 error ‘Arguments must satisfy n1 n0 n2’ end n n1:n2; y [ n-n0 0]; 单位阶跃序列； 参考代码： function y stepseq n1,n2,n0 if n0 n1 | n0 n2 | n1 n2 error ‘Arguments must satisfy n1 n0 n2’ end n n1:n2; y [ n-n0 0]; 矩形序列； 参考代码 function y srcjy n1,n2,n3,n4 n n1:n2; y [ n-n3 0 n-n4 0]; 实指数序列； 参考代码： function y srcexp a, n y a.^n （2）通过MATLAB帮助系统学习以下函数的用法： 方波——square 语法格式： x square t x square t,duty 参考代码： t 0:0.1*pi:6*pi; y square t ； axis [0 7*pi -1.5 1.5] plot t，y ； 正弦波——sin 锯齿波——sawtooth 语法格式： x sawtooth t x sawtooth t,width 参考代码： t 0:0.001:1.5; % 抽取长度1.5s，抽样频率1000Hz x sawtooth t ; %产生信号频率为50Hz plot t,x axis [0 0.2 -1 1] sinc信号 语法格式： y sinc x （3）序列的基本操作 信号加 function [y,n] sigadd x1,n1,x2,n2 % Implements y n x1 n +x2 n % % [y,n] sigadd x1,n1,x2,n2 % y sum sequence over n, which includes n1 and n2 % x1 first sequence over n1 % x2 second sequence over n2 n2 can be different from n1 % n min min n1 ,min n2 :max max n1 ,max n2 ; y1 zeros 1,length n ; y2 y1; y1 find n min n1 n max n1 1 x1; y2 find n min n2 n max n2 1 x2; y y1+y2; 信号乘 function [y,n] sigmult x1,n1,x2,n2 % Implement y n x1 n *x2 n % % [y,n] sigmult x1,n1,x2,n2 % y product sequence over nm which includes n1 and n2 % x1 first sequence over n1 % x2 second sequence over n2 n2 can be different from n1 % n min min n1 ,min n2 :max max n1 ,max n2 ; y1 zeros 1,length n ; y2 y1; y1 find n min n1 n max n1 1 x1; y2 find n min n2 n max n2 1 x2; y y1.*y2; 改变比例 移位 折叠——用fliplr函数，格式：fliplr x 取样和——用sum函数，格式：y sum x n1:n2 取样积——用prod函数，格式：y prod x n1:n2 信号能量——用sum函数，格式：Ex sum x.*conj x 或Ex sum abs x .^2 信号功率——用sum函数 卷积——用conv函数，格式：y conv a,b （4）综合应用 用MATLAB设计一个程序，调用（1）和（2）中的函数，并以一个Figure同时显示4个函数的方式进行显示。 Examples Perform ideal bandlimited interpolat