第一次作业Matlab编程实现DFT实验.doc
- 文档编号:2504403
- 上传时间:2023-05-03
- 格式:DOC
- 页数:3
- 大小:59.50KB
第一次作业Matlab编程实现DFT实验.doc
《第一次作业Matlab编程实现DFT实验.doc》由会员分享,可在线阅读,更多相关《第一次作业Matlab编程实现DFT实验.doc(3页珍藏版)》请在冰点文库上搜索。
题目:
用Matlab实现以下3个信号的DFT分析,具体要求:
画出对应信号的时域波形及其振幅、相位谱。
(画出物理坐标)
解答:
代码如下:
(1)N=64;
%生成delata函数,并对其做64点的采样
y1=zeros(1,N);
y1
(1)=1;
i=0:
N-1;
subplot(3,3,1);
stem(i,y1);
xlabel('n');
title('delta函数做64点采样');
Y1=fft(y1,N);%delta函数做64点采样图
mag=abs(Y1);
subplot(3,3,2);
stem(i,mag);%delta函数做64点FFT
xlabel('K');
title('delta函数的64点DFT');
ang=angle(Y1)*180/pi;%求DFT的相位,并把弧度转化为角度
subplot(3,3,3);
plot(i,ang);
xlabel('K');
ylabel('角度');
title('delta函数的64点DFT相位图');
(2)
%令A1=1;
A1=1;
t=0:
1/N:
1-1/N;
f0=50;
y2=A1*sin(2*pi*f0*t);%做64点采样
i=0:
N-1;
subplot(3,3,4);
stem(i,y2);
xlabel('n');
title('A1*sin(2*pi*f0*t)函数做64点采样f0=50');
Y2=fft(y2,N);
mag=abs(Y2);
subplot(3,3,5);
stem(i,mag);
xlabel('K');
title('A1*sin(2*pi*f0*t)函数的64点DFTf0=50');
ang=angle(Y2)*180/pi;
subplot(3,3,6);
plot(i,ang);
xlabel('K');
ylabel('角度');
title('A1*sin(2*pi*f0*t)函数的64点DFT相位图');
(3)
%令A2=1;
A2=1;
t=0:
1/N:
1-1/N;
f0=50;
y3=A2*sin(2*pi*f0*t.*t);%做64点采样
i=0:
N-1;
subplot(3,3,7);
stem(i,y3);
xlabel('n');
title('A2*sin(2*pi*f0*t*t)函数做64点采样f0=50');
Y3=fft(y3,N);
mag=abs(Y3);
subplot(3,3,8);
stem(i,mag);
xlabel('K');
title('A2*sin(2*pi*f0*t*t)函数的64点DFTf0=50');
ang=angle(Y3)*180/pi;
subplot(3,3,9);
plot(i,ang);
xlabel('K');
ylabel('角度');
title('A2*sin(2*pi*f0*t*t)函数的64点DFT相位图');
实验结果
截图:
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第一次 作业 Matlab 编程 实现 DFT 实验