典型序列频谱分析.docx
- 文档编号:24510005
- 上传时间:2023-05-28
- 格式:DOCX
- 页数:17
- 大小:291.04KB
典型序列频谱分析.docx
《典型序列频谱分析.docx》由会员分享,可在线阅读,更多相关《典型序列频谱分析.docx(17页珍藏版)》请在冰豆网上搜索。
典型序列频谱分析
第1章设计任务与要求1
1.1对于三种典型序列——单位采样序列、实指数序列、矩形序列1
1.2自行设计一个周期序列1
第2章原理及过程2
1设计原理2
第3章设计内容4
1.1单位采样序列4
1.1.1时域波形4
1.2傅里叶变换4
1.3幅度谱及相位谱5
1.4频移6
1.5时移7
2.1时域图形7
2.2傅里叶变换8
2.3幅度谱与相位谱9
2.4频移10
2.5时移10
3.1时域图形错误!
未定义书签。
3.2傅里叶变换错误!
未定义书签。
3.3幅度谱与相位谱错误!
未定义书签。
3.4时移错误!
未定义书签。
3.5频移错误!
未定义书签。
4.1幅度特性曲线11
4.4周期序列的DFS12
4.5傅里叶变换12
第4章心得与体会13
参考文献14
第1章设计任务与要求
1.1对于三种典型序列——单位采样序列、实指数
序列、矩形序列
要求:
(1)画出以上序列的时域波形图;
(2)求出以上序列的傅里叶变换;(3)
画出以上序列的幅度谱及相位谱,并对相关结果予以理论分析;(4)对以上序列分别进行
时移,画出时移后序列的频谱图,验证傅里叶变换的时移性质;(5)对以上序列的频谱分
别进行频移,求出频移后频谱所对应的序列,并画出序列的时域波形图,验证傅里叶变换的
频移性质。
1.2自行设计一个周期序列
要求:
(1)画出周期序列的时域波形图;
(2)求周期序列的DFS并画出幅度特性
曲线;(3)求周期序列的FT,并画出幅频特性曲线;(4)比较DFS和FT的结果,从中
可以得出什么结论。
第2章原理及过程
1设计原理
运用Matlab及相关软件,对三种典型序列进行频谱分析,得到并验证相应的结果。
MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程
序设计的高科技计算环境。
它将数值分析、矩阵计算、科学数据可视化以及非线性
动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学
研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决
方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑
模式,代表了当今国际科学计算软件的先进水平。
[1]
MATLAB和Mathematica、Maple并称为三大数学软件。
它在数学类科技应用软件
中在数值计算方面首屈一指。
MATLAB可以进行矩阵运算、绘制函数和数据、实现算
法、创建用户界面、连
接其他编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图
像处理、信号检测、金融建模设计与分析等领域。
MATLAB的基本数据单位是矩阵,它的指令表达式与数学、工程中常用的形式十
分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得
多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学
软件。
在新的版本中也加入了对C,FORTRANC++,JAVA的支持。
可以直接调用,用
户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许
多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用。
利用傅里叶变换的方法对振动的信号进行分解,并按频率顺序展开,使其成为
频率的函数,进而在频率域中对信号进行研究和处理的一种过程,称为频谱分析。
将信号在时间域中的波形转变为频率域的频谱,进而可以对信号的信息作定量
解释。
对信号进行频谱分析,是对其进行傅里叶变换,得到其振幅谱与相位谱。
分析软件
主要为Matlab。
对于信号来说,分模拟信号与数字信号。
进行频谱分析时,对于模拟信号
来说,首先对其进行抽样,使其离散化,然后利用离散傅里叶变换(DFT)或者快速
傅里叶变换(FFT),然后对其幅度(ABS)和相位(ANGLE的图像进行分析,而对
于数字信号来说,则可直接进行离散傅里叶变换或快速傅里叶变换
第3章设计内容
1.1单位采样序列
1.1.1时域波形
程序:
n=-5:
10;
y=[zeros(1,5),1,zeros(1,10)];
stem(n,y)
axis([-5,10,0,2]);
title('单位取样序列')
1.2傅里叶变换
程序:
clear
clc
%%俞入信号
n=-5:
10;%时域信号的时间范围
x=[zeros(1,5),1,zeros(1,10)];%时域信号x
w=[0:
1e-2:
2];%想要观察的频率范围
%%预定义y=w;
a=w;
j=sqrt(-1);%先定义变量维度,提高运算速度
%%计算频点
fori=1:
length(w)
f=trapz(n,x.*exp(-j*w(i)*n));y(i)=abs(f);
a(i)=angle(f);
end
%%输出
subplot(3,1,1),plot(n,x)subplot(3,1,2),plot(w,y)subplot(3,1,3),plot(w,a)
图2
1.3幅度谱及相位谱
n=1:
50;%定义序列的长度是50
x=zeros(1,50);%注意:
MATLAB中数组下标从1开始x
(1)=1;
closeall;
subplot(3,1,1);stem(x);title('单位冲击信号序列');
k=-25:
25;
X=x*(exp(-j*pi/12.5)).A(n'*k);
magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('单位冲击信号的幅度谱’);
angX=angle(X);%绘制x(n)的相位谱
subplot(3,1,3);stem(angX);title('单位冲击信号的相位谱');
1.4频移
程序:
fs=256;N=256;%采样频率和数据点数
n=0:
N-1;t=n/fs;%时间序列f=20;%定义信号中心频率(hz)x=[zeros(1,5),1,zeros(1,10)];%信号
淋序区
y1=fft(x,N);%对信号进行快速Fourier变换
y2=fftshift(y1);
mag2=abs(y2);
f仁n*fs/N;%频率序列
f2=n*fs/N-fs/2;%shift后频率序列
-Asn--iiDfia并詩non1ah
图4
1.5时移
程序N=50;
n=0:
N-1;M=20;
xn=[1zeros(1,N)];
nm=mod((n-M),N);
xm=xn(nm+1);
subplot(3,1,1),stem(xm);
title('|X(k)|');
k=-25:
25;
X=x*(exp(-j*pi/25))4((nm)'*k);magX=abs(X);%绘制x(m)的幅度谱subplot(3,1,2);stem(magX);title('矩形序列的幅度谱');
angX=angle(X);%绘制x(m)的相位谱
subplot(3,1,3);stem(angX);title('矩形序列相位谱');
|X(k)|
2.1时域图形
程序:
n=-5:
10;
y=[zeros(1,5),ones(1,5),zeros(1,6)];
矩彫信号序列
图6
2.2傅里叶变换
程序:
clear
clc
%%俞入信号
n=-5:
10;%时域信号的时间范围
x=[zeros(1,5),1,zeros(1,10)];%时域信号x
w=[0:
1e-2:
2];%想要观察的频率范围
%%预定义
y=w;
a=w;
j=sqrt(-1);%先定义变量维度,提高运算速度
%%计算频点
fori=1:
length(w)
f=trapz(n,x.*exp(-j*w(i)*n));
y(i)=abs(f);
a(i)=angle(f);
end
%%输出
subplot(3,1,1),
plot(n,x)
subplot(3,1,2),
plot(w,y)
subplot(3,1,3),
plot(w,a)
图7
2.3幅度谱与相位谱
程序:
n=1:
50;
x=sign(sign(10-n)+1);
closeall;
subplot(3,1,1);stem(x);title('矩形信号序列’);
k=-25:
25;
X=x*(exp(-j*pi/25))A(n'*k);
magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('矩形序列的幅度谱');
angX=angle(X);%绘制x(n)的相位谱
subplot(3,1,3);stem(angX);title('矩形序列相位谱');
矩形序列
0.5--
0$0510152025303540
矩形序列的幅度谱
矩形序列的相位谱
图8
2.4频移
程序:
fs=256;N=256;%采样频率和数据点数
n=0:
N-1;t=n/fs;%时间序列
f=20;%定义信号中心频率(hz)
x=[zeros(1,5),ones(1,5),zeros(1,6)];%信号
淋序区
y1=fft(x,N);%对信号进行快速Fourier变换
y2=fftshift(y1);
mag2=abs(y2);
f仁n*fs/N;%频率序列
f2=n*fs/N-fs/2;%shift后频率序列
subplot(4,1,1),plot(f2,mag2,丫');%绘出随频率变化的振幅
-150-100-50050100150
图9
2.5时移
程序:
N=50;
n=0:
N-1;M=20;
xn=sign(sign(10-n)+1);
nm=mod((n-M),N);
xm=xn(nm+1);
subplot(3,1,1),stem(xm);
title('|X(k)|');
k=-25:
25;
X=x*(exp(-j*pi/25))4((nm)'*k);magX=abs(X);%绘制x(m)的幅度谱subplot(3,1,2);stem(magX);title('矩形序列的幅度谱');
angX=angle(X);%绘制x(m)的相位谱
subplot(3,1,3);stem(angX);title('矩形序列相位谱');
|X(k)l
图10
4.1幅度特性曲线
程序:
n=1:
50;%定义序列的长度是50
x=sin(pi*n);%注意:
MATLAB中数组下标从1开始
closeall;
subplot(3,1,1);stem(x);title('单位冲击信号序列');
k=-25:
25;
X=x*(exp(-j*pi/12.5)).A(n'*k);
magX=abs(X);%绘制x(n)的幅度谱subplot(3,1,2);stem(magX);title('单位冲击信号的幅度谱’);
angX=angle(X);%绘制x(n)的相位谱
subplot(3,1,3);stem(angX);title('单位冲击信号的相位谱');
单位冲击信号序列
1■1■1■1P
一""0耳一J¥卅de?
朋
P'o
nTat\
〜込宀才goCM2『也j占
<2AA
jiiiiiii/
6
\1
■2IIIIII||吉|
'0;」1U2030他4U屈SO
单位冲击倍号的相位谱
Q?
弃
pMT
曾
11
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 典型 序列 频谱 分析