信号的频谱分析.docVIP

南 京 工 程 学 院 通信工程学院 实 验 报 告 题 目： 信号频谱分析 课 程 名 称 信号与系统实验 专 业 电子信息工程（传感网） 班 级 传感121 学 生 姓 名 史少毅 学 号 208120723 设 计 地 点 通信原理实验室 指 导 教 师 徐伟业 设计起止时间： 年 月 日至 年 月 日 一、实验目的： 1、掌握傅立叶级数（FT），学会分析连续时间周期信号的频谱及MATLAB实现； 2、掌握傅立叶变换F（jw），了解傅立叶变换的性质以及MATLAB实现。 3、掌握信号抽样与恢复的原理，能够用MATLAB实现一般信号的采样与恢复。 二、课前验证实验内容： 1、练习并验证实验指导书上实验十----十四的内容。 2、验证下列实验， 用matlab演示周期矩形脉冲信号的周期、脉宽与频谱宽度之间的变化关系。 以课本P133周期矩形脉冲信号为例，画出该信号的频谱图。改变周期及其脉宽的值，察看频谱的变化，与书中的内容进行比较，并由不同谱图分析出周期，脉宽、谱密度、谱幅度以及带宽之间的关系。 function x=time_fun_e %定义的函数 T=5; t=-2*T:0.01:2*T; tao=T/5; x=rectpuls(t,1); subplot(211) plot(t,x) hold on x=rectpuls(t-5,1); plot(t,x) hold on x=rectpuls(t+5,1); plot(t,x) hold on axis([-10,10,0,1.2]); function y=time_fun_s %定义的函数 syms a a1 T=input(please input 周期T=); M=input(please input 周期与脉宽之比M=); A=1; tao=T/M; a=tao/2; y1=sym(Heaviside(t+a1))*A; y=y1-sym(Heaviside(t-a1))*A; y=subs(y,a1,a); y=simple(y); % ezplot(y) %Nn 输出数据的准确位数 % A_sym cos展开系数An % B_sym sin展开系数Bn % T 周期 % Nf 谐波阶数 % m 周期与脉宽比 % tao 脉宽 syms t n y Nf=60; %T=input(please input 周期T=); Nn=32;%输出数据的准确解取32位数 y=time_fun_s %调用符号变量写出的周期矩形脉冲 A0=2*int(y,t,0,T)/T; %计算直流系数A0 f=1/T; As=int(2*y*cos(2*pi*f*n*t)/T,t,0,T); %cos展开系数An Bs=int(2*y*sin(2*pi*f*n*t)/T,t,0,T); %sin展开系数Bn A_sym(1)=double(vpa(A0,Nn)); %获取串数组A0所对应的ASC2码数值数组 for k=1:Nf A_sym(k+1)=double(vpa(subs(As,n,k),Nn));%获取串数组A所对应的ASC2码数值数组 B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn));%获取串数组B所对应的ASC2码数值数组 end s1=fliplr(A_sym);%对A_sym阵左右对称交换 s1(1,k+1)=A_sym(1);%对A_sym的1*K阵扩展为1*（K+1）阵 s2=fliplr(1/2*s1);%对扩展后的S1阵左右对称交换回原位置 s3(1,k+1)=0; %对B_sym阵左右对称交换 s3=fliplr(1/2*B_sym);%对B_sym的1*K阵扩展为1*（K+1）阵 s4=fliplr(s3); %对扩展后的S3阵左右对称交换回原位置 s5=s2-s4*i; %用三角函数展开系数A、B合成傅立叶指数系数 s6=fliplr(s5);%对傅立叶指数系数阵左右对称交换 N=Nf*2*pi*f; %指出谐波点数 k2=-N:2*pi*f:N; s7=[s6,s5(2:end)];%形成-N：N的傅立叶复指数系数 x=timefun_e; %调用时间函