实验1 离散时间信号的matlab实现.pdfVIP

试验1 离散时间信号的MATLAB实现 第一部分 序列的生成 主要 第二部分 离散序列的基本运算与波形变换的MATLAB实现 内容 第三部分 线性时不变系统求解的MATLAB实现 第四部分 序列的傅里叶变换与信号重构MATLAB实现 试验 Matlab 平台 平台 试验 以第四部分序列的傅里叶变换与信号重构MATLAB实现为内容写一份试验 报告 报告 第一部分 序列的生成 1．正弦序列 离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出 p 序列的波形。下面就是正弦序列sin k 的MATLAB源程序。程序运行结果如图1.1所示。 6 例1-1 %正弦序列实现程序 k=0:39; fk=sin(pi/6*k); stem(k,fk) 图1.1 正弦序列波形 2．指数序列 离散指数序列的一般形式为cak ，可用MATLAB中的数组幂运算（即点幂运算）c*a.^k 来实 现。下面为用MATLAB编写绘 离散时间实指数序列波形的函数。 实指数序列 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 3 5 3 利用上述函数，实现实指数波形MATLAB程序如下（其中a 值分别为 ）。 4 4 4 4 例1-2 %离散时间实指数序列实现程序 subplot 221; dszsu(1,5/4,0,20); xlabel(k); title(f1[k]); subplot 222 dszsu(1,3/4,0,20); xlabel(k); title(f2[k]); subplot 223; dszsu(1,-5/4,0,20); xlabel(k); title(f3[k]); subplot 224; dszsu(1,-3/4,0,20); xlabel(k); title(f4[k]); 程序运行结果如图2所示。如图可知，对于离散时间实指数序列cak ，当a 的绝对值大于1 时，序列为随时间发散的序列，当a 的绝对值小于1时，序列为随时间收敛的序列。同时可 见，当a 的值小于零时，其波形在增长或衰减的同时，还交替地改变序列值的符号。 图1.2 不同底数的实指数序列 虚指数序列 对于离散时间虚指数序列，可用通过调用下列绘制虚指数序列时域波形的MATLAB函数。 function[]=dxzsu(n1,n2,w) %n1:绘 波形的虚指数序列的起始时间序号 %n2：绘 波形的虚指数序列的终止时间序号 %w：虚指数序列的角频率 k=n1:n2; f=exp(i*w*k); Xr=real(f) Xi=imag(f) Xa=abs(f) Xn=angle(f) subplot(2,2,1), stem(k,Xr,filled),title(实部); subplot(2,2,3), stem(k,Xi,filled),title(虚部); subplot(2,2,2), stem(k,Xa,filled),title(模); su