ch5555作业.docx
- 文档编号:7590064
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:15
- 大小:198.30KB
ch5555作业.docx
《ch5555作业.docx》由会员分享,可在线阅读,更多相关《ch5555作业.docx(15页珍藏版)》请在冰豆网上搜索。
ch5555作业
成绩:
《数字信号处理》
作业与上机实验
(第三、四章)
班级:
N电信12-1F
学号:
24122201801
姓名:
聂诚
任课老师:
李宏民
完成时间:
2014.10.27
信息与通信工程学院
2013—2014学年第2学期
第3-4章离散傅里叶变换及快速算法与应用
1、对信号在x(n)=0.01n,n=0,1,2....N-1,应用FFT程序求N=128点的信号的FFT。
画出其幅度谱与相位谱。
n=0:
128;N=128;
x=0.01*n;
stem(n,x);
figure
(2)
x=0.01*n;N=128;
Y=fft(x,N);
n=0:
128;stem(abs(Y))
figure(3)
x=0.01*n;N=128;
Y=fft(x,N);
n=0:
128;stem(angle(Y))
原函数为图1:
图1
幅度谱为图1.1:
图1.1
相位谱为图1.2:
图1.2
2、一个由40Hz和100Hz正弦信号构成的信号,受零均值随机噪声的干扰(噪声服从标准正态分布,由randn函数产生,n(t)=randn(m,n)返回一个m*n的随机矩阵),即
。
数据采样率为500Hz,试用FFT函数来分析其信号频率成分。
①求其幅度频谱,从频谱图(横坐标以HZ为单位,用plot函数画图)中能否观察出信号的2个频率分量?
n=0:
500;
t=1/500;
s=sin(2*pi*100*t)+sin(2*pi*40*t*n)+randn(1,501);
figure
(1)
plot(n,s);
xlabel('HZ');
Y=fft(s,501)
A=abs(Y);
figure
(2);
plot(n/(t*501),A);
频谱为图2:
图2
幅度谱为图2.1:
图2.1
②提高采样点数,再求该信号的幅度频谱图,此时幅度频谱发生了什么变化?
信号的2个模拟频率和数字频率各为多少?
FFT频谱分析结果与理论上是否一致?
n=0:
1000;
t=1/500;
s=sin(2*pi*100*t)+sin(2*pi*40*t*n)+randn(1,1001);
figure
(1)
plot(n,s);
xlabel('HZ');
Y=fft(s,1001)
A=abs(Y);
figure
(2);
plot(n/(t*1001),A);
频谱为图2.2:
图2.2
幅度谱为图2.3:
图2.3
3、研究高密度频谱与高分辨率频谱。
频率分辨率是指所用的算法能将信号中两个靠得很近的谱峰分开的能力。
信号末尾补零由于没有对原信号增加任何新的信息,因此不能提高频率分辨率,但可以减小栅栏效应,所得到的频谱称为高密度频谱。
在维持采样频率不变的情况下,为提高分辨率只能增加采样点数N,此时所得到的频谱称为高分辨率频谱。
设有连续信号
以采样频率
对信号
采样,分析下列几种情况的幅频特性。
1采集数据长度N=16点,做N=16的FFT;采集数据长度N=16点,补零到256点,做256点的FFT;
(1)fs=32000;%采样频率
N=16;%采集16点
n=0:
N-1;%做16点fft
xa=cos(2*pi*6.5*10^3*n/fs)+cos(2*pi*7*10^3*n/fs)+cos(2*pi*9*10^3*n/fs);
Y=fft(xa);
n=[0:
(length(xa)-1)];
w=(2*pi/length(xa))*n;
figure
(1)
Am=abs(xa);
plot(w,Am,'b');
title('Magnitudepart');
xlabel('frequencyinradians');
ylabel('|X(exp(jw))|');
运行结果:
N=16点DFT幅频特性:
图3
(2)x=[xa(1:
16),zeros(1,240)];%补零到256
Y=fft(x)%做256点DFT
n=[0:
(length(x)-1)];
w=(2*pi/length(x))*n;
figure
(2)
Am=abs(x);
plot(w,Am,'b');grid;
title('Magnitudepart');
xlabel('frequencyinradians');
ylabel('|X(exp(jw))|');
2采集数据长度N=64点,做N=64的FFT;采集数据长度N=64点,补零到256点,做256点的FFT;
3采集数据长度N=256点,做N=16点的FFT。
N=input('N=');
N1=input('N1=');
n=0:
N-1;
n1=0:
N1-1;
t=1/32000;
xa=cos(2*pi*6.5*10^3*n*t)+cos(2*pi*7*10^3*n*t)+cos(2*pi*9*10^3*n*t);
a=fft(xa);b=abs(a);
w=2*pi*n/16;
subplot(211);plot(w,b);
x=[xa(1:
16),zeros(1,240)];
a1=fft(x);b1=abs(a1);
w=2*pi*n1/16;
subplot(212);plot(w,b1);
4观察以上几幅不同的幅频特性曲线,分析和比较它们的特点,并说明形成的原因。
采样数据为16的时候图像由于点数太少,很难反映出频谱的细节特征,只能分辨出两个频率分量,随着采样点的不段增加频率间隔缩小,连续曲线越来越平滑,并且由16点经过补零所得到的,频率分辨率并未提高,仍然只能看出两个频率分量。
5注:
在MATLAB中,可用zeros函数实现填零运算,例如x1=[x(1:
1:
16),zeros(1,240)]。
所有频谱图的横坐标以HZ为单位,用plot函数画图。
6一个连续信号含三个频率分量,经采样得以下序列
。
已知N=16,
分别为1/16,1/64,观察其频谱;当N=64,128,
不变,其结果有何不同,为什么?
N=16
分别为1/16时的源程序:
N=16;
n=0:
N-1;
df=1/16;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
plot(n,x);
title('F=16的时域')
dd=fft(x);
dd=abs(dd);
figure
(2)
stem(n,x)
title('F=16的频域')
N=16
分别为1/64时的源程序:
N=16;
n=0:
N-1;
df=1/64;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
plot(n,x);
title('F=64的时域')
dd=fft(x);
dd=abs(dd);
figure
(2)
stem(n,x)
title('F=64的频域')
N=64
分别为1/64时的源程序:
N=64;
n=0:
N-1;
df=1/64;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
plot(n,x);
title('F=64的时域')
dd=fft(x);
dd=abs(dd);
figure
(2)
stem(n,x)
title('F=64的频域')
N=128
分别为1/64时的源程序:
N=128;
n=0:
N-1;
df=1/64;
x=sin(2*pi*0.15*n)+cos(2*pi*(0.15+df)*n)+cos(2*pi*(0.15+2*df)*n);
plot(n,x);
title('F=64的时域')
dd=fft(x);
dd=abs(dd);
figure
(2)
stem(n,x)
title('F=64的频域')
4、语音信号频谱分析
首先应用录音软件(见下图),录音“我是XXX(学生姓名)”,得到语音数据“录音.WAV”并将“录音.WAV”文件复制到MATLABWORK目录下。
其次在MATLAB环境下,X(n)=wavread('录音.WAV',2048),得到2048点的语音数据X(n),对X(n)进行FFT频谱分析,画出频谱图(横坐标以HZ为单位,用plot函数画图)。
源程序为:
fs=25600;
[x,fs,bits]=wavread('gg.wav');
sound(x,fs,bits);
y1=fft(x,4096);
f=fs*(0:
2047)/4096;
figure
(1)
plot(f,abs(y1(1:
2048)));
figure
(1)
plot(xn);x=fft(xn,2048);a=abs(x);
xlabel('X(n)');
figure
(2)
plot(w,a);
xlabel('Hz');
频谱图为:
图4
图4
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ch5555 作业