二 连续信号的域分析.docVIP

实验二 连续信号的频域分析 学号： 姓名 成绩： 一、实验目的 （1）熟悉周期信号的傅立叶分解与合成。 （2）掌握傅立叶变换的数值计算。 二、实验原理 （1）周期信号的傅立叶级数 周期信号可展开为傅立叶级数，反映出周期信号可分解成直流分量、基波和各次谐波分量。周期为T，高度为E的方波，其傅立叶级数为 此处方波是奇函数，无直流分量，只包含基波和奇次谐波。由基波和各次谐波也可以合成方波，其误差取决于所取谐波的项数。所取谐波越多，其均方误差越小。 （2）连续信号的傅立叶变换 若f(t)是时限信号，上式中n的取值个数是有限的，设n是从n1~n2共取N个值，则 其中nτ用t表示；τ是t的增量，用dt表示；ω用w表示。利用MATLAB，则有 F=f*exp(-j*t′*w)*dt 其中， t′是N维列向量，w是M维行向量，t′*w是矩阵exp(-j*t′*w)也是N*M维矩阵。f同样傅立叶的逆变换变为： f= F*exp(j*t*w′)* dw/pi) 三、验证性实验 方波高度为2，周期为2π，分解其1、3、5、7、9次谐波。 程序： %方波的分解,基波角频率为1 t=0:0.01:2*pi; y=zeros(10,length(t)); % 预留10次谐波存储空间 for k=1:2:9 % 求1,3,5,7,9次谐波 y(k,:)=4 *sin(k*t)/k/pi; % k次谐波 end %figure(Position,[100,100,400,200]); plot(t,y(1:9,:)); grid; 问题：方波的角频率是多大？1、3、5、7、9次谐波的周期、幅度分别是多大？ 解：方波的角频率为1。一次谐波的周期2π，幅度为4/π二次谐波的周期2π/3，五次谐波的周期2π/5，幅度为4/5π。七次谐波的周期2π/7， 九次谐波的周期2π/9。 求单个矩形脉冲的傅立叶变换，画出其幅度频谱 程序： clear N=100; % 时间的分隔点数 Tf=10; % 信号的时间宽度 dt= Tf /N; t=[1:N]*dt; f=[ones(1,N/2),zeros(1,N/2)]; %矩形脉冲信号 Wf=10; % 需求的频谱宽度 M=100; % 需求的频谱点数 w1=linspace(0,Wf, M); F1=f*exp(-j*t*w1)*dt; % 求傅立叶变换 w=[-fliplr(w1),w1(2: M)]; % 补上负频率 F=[fliplr(F1),F1(2: M)]; % 补上负频率上的频谱 subplot(2,1,1),plot(t,f,linewidth,1.5),grid subplot(2,1,2),plot(w,abs(F),linewidth,1.5),grid 问题： 从幅度谱观察，矩形脉冲的带宽为多大？ 如果将矩形脉冲延迟10秒，得到的幅度谱有没有不同？带宽为多大？ 试求cos(5t)的幅度谱。 解：从幅度谱观察，矩形脉冲的带宽为2； 如果将矩形脉冲延迟10秒，得到的幅度谱如下图所示： cos(5t)的幅度谱为下图所示： 四、设计性实验 方波高度为2，频率为1Hz，分别求其1、3、5次谐波的合成和1、3、5、7、9次谐波的合成，把两个合成结果画在同一个图形窗口的两个子窗口上。 计算这两个合成波的均方误差。 解：其程序如下 t=0:0.01:1; w=2*pi; y=zeros(10,length(t)); for k=1:2:9 y(k,:)=4 *sin(k*w*t)/k/pi; end f135=sum(y(1:5,:)) f19=sum(y(1:9,:)) subplot(2,1,1),plot(t,f135),grid subplot(2,1,2),plot(t,f19),grid half=length(t)/2; e1=f135(1:half)-1; ee1=e1^2; me=sum(ee1)/half 其幅度谱为下图所示： 且这两个合成波的均方误差为0.0679。 grass-roots party and members team construction, anti-corruption construction, thought political work, group mission work, 　　For the development of various undertakings in our team played a political role. A level newly promoted party building,