西北工业大学信号与系统实验五.docxVIP

西北工业大学 《信号与系统》实验报告 学 院: 学 号: 姓 名: 专 业: 实验时间: 实验地点: 指导教师: 软件与微电子学院 毅字楼 陈勇 西北工业大学 200年 月 、实验目的 掌握离散时间傅里叶变换的综合公式和分析公式，以及离散时间傅里叶变换 的性质，并学会使用fft和ifft的数值计算实现分析公式和综合公式。 二、实验要求 仔细阅读实验指导书，练习对fft、ifft、residue等函数的使用，理解对DTFT 的数值实现的两个主要限制，分析实际应用中的一些问题，理解全通系统，用 DTFT来设计和辨识LTI系统； 三、实验设备(环境) 操作系统：Win10 系统内存：8GB MATLAB 版本：R2014a 四、实验内容与步骤 5.1计算dtft的样本 基本题： (1)问题描述与分析 题目中要求使用解析方法计算矩形脉冲的 DTFT ,通过计算可以得出 X=(1-exp(-1*j*w*10)).*1./(1-exp(-1*j*w)) ;转换后，画出频谱图的幅值图和相位图； 用fftshift对DTFT进行重新排列，再次画出 X对w的幅值和相位 (2)代码实现 % (b) N=100; k=[0:N-1]; w=2*pi*k/N; X=(1-exp(-10*j*w))./(1-exp(-j*w)); figure; subplot(2,1,1) stem(w,abs(X)) ylabel( abs(X) ) subplot(2,1,2) stem(w,angle(X)) ylabel( angle(X) ) % (c) N=100; k=[0:N-1]; w=2*pi*k/N; w1=w-pi; x=ones(1,10); X=fft(x,N); X1=fftshift(X); figure; subplot(2,1,1) stem(w1,abs(X1)) ylabel( abs(XI) ) subplot(2,1,2) stem(w1,angle(X1)) ylabel( angle(XI) ) (3) Result(图形) MVASDaL 1 1. L LN L刑kMHkj血匕疝怖;面fffrh MVASDa L 1 1. L L N L刑kMHkj血匕疝怖;面fffrh缶加M kl L - \)x(^Fgna 4 2 o-2 0 1 2 3 4 5 6 7 4Y^qx^sDa2 4 Y^qx^sDa 0 -2 -4 -3 -2 -10 12 3 4-4（4）实验结论 -4 -3 -2 -1 0 1 2 3 4 -4 由图（b）和图（c）两个图像中可以看出，图（c）是由图（b）后半部 分移到前面所得到的图像，根据理论，理论上图（ c）前半部分的频率与图 （b）后半部分的频率相差一个2冗，而函数的周期即为2冗，所以理论上该 部分图形相同，实验图像结果与理论设想相符，所以该代码是正确的。 5.2电话拨号音 基本题： 程序1如何使用计时函数tic、toc ; clc; clear; REPS = 1000; tic; %启动计时器 total=0; for i=1:REPS % 做累加 total = total + i; end Time =double(toc);% 关闭计时器并转换浮点数，64位 REPS1 = 10000; tic;%重新启动计时器 total1=0; for i=1:REPS1 % 做累加 total1 = total1 + i; end Timel =double(toc);% 关闭计时器并转换浮点数，64位 prompt={运算时间比}; titla=显示窗口 ’; a1=num2str(Time1/Time);%num2str把非整数数组转换为字符串 def={a1}; lines=1; sita=inputdlg(prompt,titla,lines,def); 程序2 function dtmf %DTMF generation n=[0:999]; d0=sin(0.7217*n)+sin(1.0247*n); d1=sin(0.5346*n)+sin(0.9273*n); d2=sin(0.5346*n)+sin(1.0247*n); d3=sin(0.5346*n)+sin(1.1328*n); d4=sin(0.5906*n)+sin(0.9273*n); d5=sin(0.5906*n)+sin(1.0247*n); d6=sin(0.5906*n)+sin(1.1328*n); d7=sin(0.6535*n)+sin(0.9273*n); d8=sin(0.6535*n)+sin(1.0247*n); d9=sin(0.6535*n)+sin(1.1328*n); ko=