S_S频域分析研究性学习报告.doc

5信号与系统的频域分析专题研讨【目的】(1) 加深对信号与系统频域分析基本原理和方法的理解。(2) 学会利用信号抽样的基本原理对信号抽样过程中出现的一些现象的进行分析。(3) 通过实验初步了解频谱近似计算过程中产生误差的原因。(4)学会用调制解调的基本原理对系统进行频域分析。【研讨题目】1信号的抽样频率为f0 hz的正弦信号可表示为按抽样频率fsam=1/t对x(t)抽样可得离散正弦序列xk在下面的实验中，取抽样频率fsam=8khz。(1)对频率为2khz, 2.2 khz, 2.4 khz和 2.6 khz正弦信号抽样1 秒钟，利用matlab函数 sound(x, fsam)播放这四个不同频率的正弦信号。(2)对频率为5.4 khz, 5.6khz, 5.8 khz和 6.0khz正弦信号抽样1 秒钟，利用matlab函数 sound(x, fsam)播放这四个不同频率的正弦信号。(3)比较(1)和(2)的实验结果，解释所出现的现象。【题目分析】利用计算机读出正弦信号，通过声音的尖锐程度和响声大小分别判断音频大小和响度大小。【信号抽样过程中频谱变化的规律】【比较研究】利用系统的help，阅读函数sound和wavplay的使用方法。 连续播放两段音频信号，比较函数sound和wavplay的异同。【仿真结果】利用计算机读出声音信号，发现：（1） 以（2）中的各频率的抽样后信号比以（1）中各频率的抽样后信号音频要高，并且还伴有明显的高低频的混叠。（2） 以（1）中的各频率的抽样后的信号随着抽样频率的增大音频升高；以（2）中的各频率的抽样后信号随着抽样频率的增大音频有所降低。【结果的理论分析和解释】由信号抽样定理知，当抽样频率大于等于信号最高频率的两倍时，才没有混叠现象。也就是说，当取样频率从一个大于信号最高频率的两倍的数值逐渐降低时，得到的频谱在刚开始的一段范围内（）都没有混叠；而抽样频率继续降低时（时），开始出现混叠，并且随着抽样频率的降低，混叠越来越严重。【自主学习内容】【阅读文献】【发现问题】 (专题研讨或相关知识点学习中发现的问题)：【问题探究】【仿真程序】fsam=8000; % f0=2000f0=2000;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=2200f0=2200;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=2400f0=2400;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=2600f0=2600;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=5400f0=5400;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=5600f0=5600;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=5800f0=5800;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam); % f0=6000f0=6000;t=0:1/fsam:1;x=sin(2*pi*f0*t);sound(x,fsam);2 连续时间信号fourier变换的数值近似计算计算连续信号频谱是对信号和系统进行频域分析的基础，由于实际信号大多无简单的解析表达式，所以要用数值方法进行近似计算。本题要求对频谱近似计算中误差的原因进行初步的分析，希望能在计算实际信号频谱的近似计算中起一定的指导作用。若信号x(t)的非零值在区间，则可用下面提供的函数ctft1或ctft2近似计算其频谱。函数ctft的调用形式为x,f=ctft1(x,fsam,n)x,f=ctft2(x,fsam,n)其中调用变量x 存放信号x(t)的抽样值，fsam表示对连续信号x(t)的抽样频率(hz)，n表示用dft进行近似计算时dft 的点数，为了能高效的进行计算，n最好取2的整数次幂，如512, 1024等。返回变量x是计算出的信号频谱的抽样值，f(单位hz)表示对应的频率抽样点。返回变量x一般是复数，可用函数abs(x)计算出幅度谱，函数angle(x) 计算出相位谱。(1)阅读程序ctft2，叙述该程序的基本原理。该程序中有一处需要产生一个大的2维矩阵，指出该行程序，并评价该方法的优缺点。(2)取抽样频率fsam=100hz, 信号抽样长度n=1024， 分别用两个子程序近似计算信号的频谱，比较两种方法的计算时间和误差；(可用tic, toc 计算程序运行时间)(3)若将信号的时域有效宽度定义为其中表示信号在时域的最大值。试分析时域有效宽度对近似计算的影响。给出一个由信号时域有效宽度估计近似计算中所需信号长度的经验公式。(4)定义信号频域有效宽度为其中表示信号在频域的最大值。给出一个由信号频域有效宽度估计近似计算中所需抽样频率的经验公式。(5)用计算机录分别一段男生和女生的语音信号，计算其频谱并比较其特点。(6)讨论：计算误差产生的主要原因？如果不知信号的解析表达式，如何分析计算误差？%近似计算连续信号频谱的函数function x,f=ctft1(x,fs,n)x=fftshift(fft(x,n)/fs;f=-fs/2+(0:n-1)*fs/n;function x,f=ctft2(x,fs,n)tk=(0:n-1)/fs;% 时域抽样点df=fs/n; %频域抽样间隔fm=(0:n/2)*df; %频域抽样点x=x*exp(-j*2*pi*tk*fm)/fs;% 近似计算信号频谱f=-fliplr(fm(2:end) fm; % 增添负频率点x=-conj(fliplr(x(2:end) x;% 增添频率点对应的频谱【题目分析】【仿真结果】【结果分析】分析误差产生的原因及解决问题的方法。【自主学习内容】【阅读文献】【发现问题】 (专题研讨或相关知识点学习中发现的问题)：【问题探究】【仿真程序】3幅度调制和连续信号的fourier变换本题研究莫尔斯码的幅度调制与解调。本题中信号的形式为其中信号x(t)由文件ctftmod.mat定义(该文件在硬盘上)，可用命令load ctftmod 将文件ctftmod.mat定义的变量装入系统内存。运行命令load ctftmod后，装入系统的变量有af bf dash dot f1 f2 t x其中bf af： 定义了一个连续系统h(s)的分子多项式和分母多项式。可利用freqs(bf,af,w)求出该系统的频率响应，也可用sys=tf(bf,af)得到系统的模型，从而用lsim求出信号通过该系统的响应。dash dot： 给出了莫尔斯码中的基本信号dash和dot的波形f1 f2： 载波频率t: 信号x(t)的抽样点x: 信号x(t)的在抽样点上的值信号x(t)含有一段简单的消息。agend 007的最后一句话是the future of technology lies in 还未说出最后一个字，agend 007就昏倒了。你(agend 008)目前的任务就是要破解agend 007的最后一个字。该字的信息包含在信号x(t)中。信号x(t)具有式0的形式。式中的调制频率分别由变量f1和f2给出，信号m1(t)，m2(t)和m3(t)对应于字母表中的单个字母，这个字母表已用国际莫尔斯码进行编码，如下表所示：a - h o - - -v -b - i p -w -c -j - - -q - -x -d -k -r -y -e l -s z -f -m - -t -g - -n -u -(1) 字母b可用莫尔斯码表示为b=dash dot dot dot，画出字母b莫尔斯码波形；(2) 用freqs(bf,af,w)画出由bf和af定义的系统的幅度响应；(3) 利用lsim求出信号dash通过由sys=tf(bf,af)定义的系统响应，解释你所获得的结果；(4)用解析法推导出下列信号的fourier变换(5)利用(4)中的结果，设计一个从x(t)中