信号与系统仿真试验试验指示书要点-read.pptVIP

《信号与系统》仿真实验 实验指示书要点 实验三 信号的频谱分析 系统仿真的要点 关于快速傅立叶变换的fft函数的调用 频谱分析举例 一、fft函数的调用 fft(x)—DFT的快速算法 （1） 数据长度为N； 若N为2的整数次幂，则是基-2的 DFT，运算速度较快； 若N不为2的整数次幂，则运算速 度较慢。 （2）若数据实际长度小于N，则自动补 一些零，使之长度等于N. 复指数座标 正变换 反变换 其中 fft(x) 或fft(x,N)的调用编程 利用fft计算功率谱的函数演示——fftdemo (1)产生信号 运行软件Matlab, 进入Matlab Comand Windows 运行fftdemo命令，该程序是按照一步一步的断点执行的，它能一边显示注解，一边给出每一步的结果，所以我们根据这个演示，能学会fft函数的调 用及其有关应用技巧。 要求产生一个时间从0到250毫秒的，含有噪声的，频率为50到120Hzde 时域信号, 噪声的标准差为2。离散信号的时间间隔为1毫秒。 程序如下： t = 0:.001:.25; x = sin(2*pi*50*t) + sin(2*pi*120*t); y = x + 2*randn(size(t)); plot(y(1:50)), title(Noisy time domain signal), （2）函数调用计算频率特性计算功率谱以及显示结果 Y = fft(y,256); Pyy = Y.*conj(Y)/256; f = 1000/256*(0:127); plot(f,Pyy(1:128)), title(Power spectral density), ... ZOOM plot(f(1:50),Pyy(1:50)), title(Power spectral density), ... echo off End 以上程序是Matlab 提供的fftdemo演示举例，同学们跟着它执行一遍，有助于理解和掌握。 * t=(0:1/255:1); x=exp(-0.9)*sin(2*pi*60*t); y=fft(x); ry=real(y); a=abs(y); plot(t,ry,r-,t,a,b--); grid 这里是归一化频率 （1）产生信号 （2）调用函数计算频率特性(fft) （3）计算功率谱(Pyy) （4）显示结果 To get started, type one of these commands: helpwin, helpdesk, or demo. For information on all of the MathWorks products, type tour. ?fftdemo clc % This example shows the use of the FFT function for spectral % analysis. A common use of FFTs is to find the frequency % components of a signal buried in a noisy time domain signal. % % First we need to create some data. Consider data sampled at % 1000 Hz. We start by forming a time axis for our data, running % from t=0 until t=.25 in steps of 1 millisecond: t = 0:.001:.25; pause % Strike any key to continue. pause % Strike any key to continue. % Next we can form a signal containing 50 Hz and 120 Hz: x = sin(2*pi*50*t) + sin(2*pi*120*t); pause % Strike any key to continue. % and add some random noise with a standard deviation of 2 to % produce a noisy