matlab方波频谱分析课程设计报告.docVIP

通信系统建模与仿真课程设计 2010 级 通信工程 专业 1013072 班级 题 目 基于Matlab/Simulink的信号频谱的估计 姓 名 学号 指导教师 胡娟，王丹，王娜，闫利超 2013年6月14日 任务书 用Matlab编程方式产生一个100Hz的方波，画出其波形。并用fft指令计算其频谱，做出幅度谱和相位谱，与理论结果进行对比。 用Simulink方式重做上题，并通过统计模块在时域和频域同时计算信号的功率，看两者计算结果是否一致，验证帕萨瓦尔定理。 理论分析 方波的一个周期可用 依据周期信号傅里叶级数系数的定义，有 因此，方波信号的的傅里叶级数展开式为 根据周期信号傅里叶级数同傅里叶变换之间的关系： 可知，方波信号的傅里叶变换是 显然，当n为偶数时，，因此方波信号中只存在奇次谐波，其功率谱为 化为以频率为自变量表示的功率密度谱，得到 可见，方波在几次谐波处存在冲激谱线，其功率谱谱线冲激强度为数列，n取奇数，C为常数。 离散时间信号的帕斯瓦尔定理：对于N点的离散序列及其离散傅里叶变换，其时域能量等于频域能量，即 时域和频域的平均功率关系为 其中，T为采样时间间隔；N为离散时间序列的点数；为离散时间序列的时间长度。 Matlab代码详述 clear; clc; fs=1e6; t1=0:1/fs:0.1; %计算时间范围 ft=square(2*pi*100*t1,50); subplot(3,1,1);plot(t1,ft); %时域波形 axis([0 0.1 -1.2 1.2]); T1=0.01; %信号周期 w1=2*pi/T1; %信号角频率 n=-59:2:59; %奇次谐波数 W=w1.*n; %数字角频率 F_w=-4*j./n; %频谱理论结果 subplot(3,1,2); stem(W,abs(F_w)); %频域幅度谱 hold on; w_m=3e4; %截断频率 T=pi/w_m; %采样间隔 L=5.9; t=0:T:L; %时域截断 x_t=square(2*pi*100*t,50); %信号序列 N=length(x_t); %序列长度（点数） X_k=fft(x_t); %FFT计算 w0=2*pi/(N*T); %离散频率间隔 kw=2*pi/(N*T).*[0:N-1]; %离散频率样点 X_kw=T.*X_k; %乘以T得到连续傅里叶变换频谱的样值 plot(kw-w_m,abs(fftshift(X_kw)),.,MarkerSize,10); %做出数值计算的幅度谱点 subplot(3,1,3); stem(W,angle(F_w)); %频域相位谱 hold on; plot(kw-w_m,angle(fftshift(X_kw))); %做出数值计算的相位谱点 4 SIMULINK各模块说明 由于Simulink中FFT模块只接受2的整数幂次点数数据，故设计变换数据采样率为2048样值/秒。FFT变换数据长度到2048，对应时间长度为1秒。因此，频率分辨率为1Hz。由式 确定功率谱估计值。其中秒，N=2048。依据教材式（3.34）确定频域、时域