信号与系统的频域分析专题研讨

PAGE5信号与系统的频域分析专题研讨【目的】(1)加深对信号与系统频域分析基本原理和方法的理解。(2)学会利用信号抽样的基本原理对信号抽样过程中出现的一些现象的进行分析。(3)通过实验初步了解频谱近似计算过程中产生误差的原因。(4)学会用调制解调的基本原理对系统进行频域分析。【研讨题目】1．信号的抽样 频率为f0Hz的正弦信号可表示为按抽样频率fsam=1/T对x(t)抽样可得离散正弦序列x[k]在下面的实验中，取抽样频率fsam=8kHz。(1)对频率为2kHz,2.2kHz,2.4kHz和2.6kHz正弦信号抽样1秒钟，利用MATLAB函数sound(x,fsam)播放这四个不同频率的正弦信号。(2)对频率为5.4kHz,5.6kHz,5.8kHz和6.0kHz正弦信号抽样1秒钟，利用MATLAB函数sound(x,fsam)播放这四个不同频率的正弦信号。(3)比较(1)和(2)的实验结果，解释所出现的现象。【题目分析】【信号抽样过程中频谱变化的规律】【比较研究】 连续的播放两段音频信号，比较函数sound和wavplay的异同。【仿真结果】【结果的理论分析和解释】【自主学习内容】【阅读文献】【发现问题】(专题研讨或相关知识点学习中发现的问题)：【问题探究】【仿真程序】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，分别用两个子程序近似计算信号的频谱，比较两种方法的计算时间和误差；(3)若将信号的时域有效宽度定义为其中表示信号在时域的最大值。试分析时域有效宽度对近似计算的影响。给出一个由信号时域有效宽度估计近似计算中所需信号长度的经验公式。(4)定义信号频域有效宽度为其中表示信号在频域的最大值。给出一个由信号频域有效宽度估计近似计算中所需抽样频率的经验公式。(5)用计算机录分别一段男生和女生的语音信号，计算其频谱并比较其特点。（4）

先要找各找一段男女生的音频，把它以WAV的格式保存在你所编辑程序的文件夹中，再用wavread()函数把信号载入程序中来运行。对于女生：源程序：fs=8000;[x,fs,hits]=wavread('C:\DocumentsandSettings\n\桌面\新建文件夹\女声.wav');wavplay(x,fs);N=length(x);df=fs/N;n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;figure(4);subplot(2,2,1);plot(x);title('原始女生信号的时域波形');ylabel('幅值/A');grid;y0=fft(x);figure(4);subplot(2,2,2);plot(f,abs(y0(n1)));title('原始女生信号的频谱图');xlabel('频率w/pi');ylabel('幅值');grid;[y,fs,bits]=wavread('C:\DocumentsandSettings\n\桌面\新建文件夹\男声.wav');wavplay(y,fs);N=length(y);df=fs/N;n1=1:N/2;f=[(n1-1)*(2*pi/N)]/pi;subplot(2,2,3);plot(y);title('原始男生信号的时域波形');ylabel('幅值/A');grid;y0=fft(y);subplot(2,2,4);plot(f,abs(y0(n1)));title('原始男生信号的频谱图');xlabel('频率w/pi');ylabel('幅值');grid;仿真结果：【结果分析】由以上分析可以得到男生的频谱高于女生的频谱【自主学习内容】fft以及wavread函数的使用。3．幅度调制和连续信号的Fourier变换 本题研究莫尔斯码的幅度调制与解调。本题中信号的形式为 LISTNUMchf其中信号x(t)由文件ctftmod.mat定义(该文件在硬盘上)，可用命令Loadctftmod将文件ctftmod.mat定义的变量装入系统内存。运行命令Loadctftmod后，装入系统的变量有afbfdashdotf1f2tx其中bfaf：定义了一个连续系统H(s)的分子多项式和分母多项式。可利用freqs(bf,af,w)求出该系统的频率响应，也可用sys=tf(bf,af)得到系统的模型，从而用lsim求出信号通过该系统的响应。dashdot：给出了莫尔斯码中的基本信号dash和dot的波形f1f2：载波频率t:信号x(t)的抽样点x:信号x(t)的在抽样点上的值 信号x(t)含有一段简单的消息。Agend007的最后一句话是Thefutureoftechnologyliesin···还未说出最后一个字，Agend007就昏倒了。你(Agend008)目前的任务就是要破解Agend007的最后一个字。该字的信息包含在信号x(t)中。信号x(t)具有式REF_Ref45974728\r\h(1)的形式。式中的调制频率分别由变量f1和f2给出，信号m1(t)，m2(t)和m3(t)对应于字母表中的单个字母，这个字母表已用国际莫尔斯码进行编码，如下表所示：A·H····OV···B···I··P··W·C··J·Q·X··D··K·R··Y·E·L···S···Z··F···MTG·N·U··字母B可用莫尔斯码表示为b=[dashdotdotdot]，画出字母B莫尔斯码波形；(2)用freqs(bf,af,w)画出由bf和af定义的系统的幅度响应；(3)利用lsim求出信号dash通过由sys=tf(bf,af)定义的系统响应，解释你所获得的结果；(4)用解析法推导出下列信号的Fourier变换(5)利用(4)中的结果，设计一个从x(t)中提取信号m1(t)的方案，画出m1(t)