实验2利用MATLAB分析信号频谱及系统的频率特性

1广东技术师范学院实验报告学院：自动化学院专业：班级：成绩：姓名：学号：计算机编号：实验地点：微机6室实验日期：指导教师签名：实验（二）项目名称：利用MATLAB分析信号频谱及系统的频率特性一．实验目的1．深入理解信号频谱的概念，掌握典型信号的频谱以及Fourier变换的主要性质及其matlab实现；2．学习和掌握连续时间系统的频率特性及其幅度特性、相位特性的物理意义及其matlab实现；3．掌握抽样定理。二．实验原理1.对于非周期信号)(tf，其傅立叶变换及其反变换式定义如下：dejFtfdtetfjFtjtj)(21)()()(式中，)(jF是原函数)(tf的傅立叶变换，称为频谱函数，它是一个复函数，可以写成)()()(jejFjF。它的模量|)(|jF是频率的函数，代表信号中各频率分量的相对大小；相角)(也是频率的函数，代表有关频率分量的相位。为了与周期信号的频谱相一致，人们习惯上把|)(|jF~与)(~曲线分别称为非周期信号的幅度频谱与相位频谱。容易看出，它们在形状上与相应的周期信号频谱包络线相同。通过典型信号频谱以及Fourier变换性质的研究，可以初步掌握Fourier分析方法的应用，同时验证一些典型信号的频谱以及傅立叶变换的主要性质，使实验者能够直观地了解信号的时域、频域波形对照，加深对信号频谱的理解。Matlab提供了能直接求解傅立叶变换和反变换的函数fourier（）、ifourier（）。调用格式分别为：F=fourier(f)f=ifourier(F)2．频域分析法与时域分析法的不同之处主要在于信号分解的单元函数不同。在频域分析法中，信号分解成一系列不同幅度、不同频率的等幅正弦函数，通过求取对每一单元激励产生的响应，并将响应叠加，再转换到时域以得到系统的总响应。所以说，频域分析法是一种变域分析法。它把时域中求解响应的问题通过Fourier级数或Fourier变换转换成频域中的问题；在频域中求解后再转换回时域从而得到最终结果。预习情况良好操作情况良好考勤情况全勤数据处理情况良好2所谓频率特性，也称频率响应特性，是指系统在正弦信号激励下稳态响应随频率变化的情况，包括幅度随频率的响应和相位随频率的响应两个方面。利用系统函数也可以确定系统频率特性，公式如下：)(|)(|HjjsejHsHjH(3–1)幅度响应用jH表示，相位响应用)(H表示。Matlab提供了专门对连续时间系统频率响应H(jω)进行分析的函数freqs()。该函数可以求出系统频率响应的数值解，并可绘出系统的幅频和相频响应曲线。一般调用格式：[h,w]=freqs(b,a,n)其中h为返回w所定义的频率点w上系统频率响应的幅值；b为系统频率响应分子多项式系数，a为系统频率响应分母多项式系数，n为输出频率点个数。3．一个频谱受限的信号f(t)，如果频谱只占据-ωm~ωm的范围，则信号f(t)可以用等间隔的抽样值唯一的表示。而抽样间隔必须不大于1/(2fm)。也就是说：对于带限信号，当ωS≥2ωm时，频谱不发生混叠，可用理想低通滤波器将原信号从抽样信号中无失真地恢复；否则，频谱将会混叠。三．实验内容1.请用MATLAB绘制函数t/tsin)t(f、傅立叶变换)]t(f[F及其)]t(f[FF1的波形。（注：请调用函数()ifourier(),fourier）symst;x=sin(t)/t;F=fourier(x);subplot(2,1,1);ezplot(F,[-1.5,1.5]);title('傅里叶变换F[f(t)]的波形');legend('pi*heaviside(1-w)-pi*heaviside(-w-1)');xlabel('t轴');ylabel('F轴');f=ifourier(F);subplot(2,1,2);ezplot(f,[-1515]);title('傅里叶反变换的波形');legend('((pi*(cos(x)*i+sin(x)))/x-(pi*(cos(x)*i-sin(x)))/x)/(2*pi)');xlabel('t轴');ylabel('f轴');32.一RLC二阶高通滤波器如下图所示。已知2,0.4,0.05RLHCF，请用MATLAB求其频率响应并绘制幅频响应和相频响应曲线。(注：请先计算出相应的参数)+CLRu2(t)u1(t)--1.5-1-0.500.511.50123t轴傅里叶变换F[f(t)]的波形F轴-15-10-505101500.51t轴傅里叶反变换的波形f轴pi*heaviside(1-w)-pi*heaviside(-w-1)((pi*(cos(x)*i+sin(x)))/x-(pi*(cos(x)*i-sin(x)))/x)/(2*pi)（第1题）4解：CRjwCLjwLRjwLRjwLRjwCRCLjwLRRjwLjwLRjwCRjwLjwLRRjwLjwCRjwLjwH)()()(11111111)(22代入数据2,0.4,0.05RLHCF，得：40)(8)(8.0)(8.005.02)(05.04.0)(24.0)(24.0)(2222jwjwjwjwjwjwjwHB=[0.800];A=[0.8840];[h,w]=freqs(B,A,100);h1=abs(h);h2=angle(h);subplot(2,1,1);plot(w,h1,'-r','LineWidth',1.5);gridon;xlabel('角频率(W)');ylabel('幅度');title('H(jw)的幅频特性');subplot(2,1,2);plot(w,h2*180/pi,'c','LineWidth',1.5);gridon;xlabel('角频率(W)');ylabel('相位(度)');title('H(jw)的相频特性');53．设有限频带信号f(t)=5+2cos(2πt)+cos(4πt)。（1）计算该信号的奈奎斯特频率；（2）以不同的采样频率对该信号进行采样，画出采样前后信号的频谱，对比和分析信号临界采样、过采样和欠采样情况下，信号频谱有何变化010203040506070809010000.51角频率(W)幅度H(jw)的幅频特性0102030405060708090100050100150200角频率(W)相位(度)H(jw)的相频特性（第2题）（1）解：该信号的临界抽样角频率：82s临界抽样频率：HZfss42临界抽样周期：sTs25.0display('奈奎斯特周期1秒，Ts1，过采样，Ts1，欠采样');display('请输入采样周期：');Ts=input('Ts=');t=0:0.01:40;y=5+2*cos(2*pi*t)+cos(4*pi*t);subplot(221);plot(t,y,'r');axis([03310]);6012346810t单位：sf(t)-50500.511.52x104F(j)单位:pi012346810t单位：sfs(t)-202050100150200单位:piFs(j)7、012346810t单位：sf(t)-50500.511.52x104F(j)单位:pi012346810t单位：sfs(t)-202050100150200单位:piFs(j)012346810t单位：sf(t)-50500.511.52x104F(j)单位:pi012346810t单位：sfs(t)-2020200400600800单位:piFs(j)（过抽样）8