clear all;t=0:0.1:50;yp=sin(2*pi*2*t);yq=sin(2*pi*2.0125*t);y=yp+yq; Y1=fft(y,512);Y2=fft(y,4096);Y1=fftshift(Y1);Y2=fftshift(Y2);c1=[0:255]/51.2;c2=[0:2047]/409;subplot(2,1,1);plot(c1,abs(Y1(257:512)));axis([1.2,2.2,-30,300]);title('512点FFT的频
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/17 19:29:11
clear all;t=0:0.1:50;yp=sin(2*pi*2*t);yq=sin(2*pi*2.0125*t);y=yp+yq; Y1=fft(y,512);Y2=fft(y,4096);Y1=fftshift(Y1);Y2=fftshift(Y2);c1=[0:255]/51.2;c2=[0:2047]/409;subplot(2,1,1);plot(c1,abs(Y1(257:512)));axis([1.2,2.2,-30,300]);title('512点FFT的频
clear all;
t=0:0.1:50;
yp=sin(2*pi*2*t);
yq=sin(2*pi*2.0125*t);
y=yp+yq;
Y1=fft(y,512);
Y2=fft(y,4096);
Y1=fftshift(Y1);
Y2=fftshift(Y2);
c1=[0:255]/51.2;
c2=[0:2047]/409;
subplot(2,1,1);plot(c1,abs(Y1(257:512)));
axis([1.2,2.2,-30,300]);
title('512点FFT的频谱');
grid on;
subplot(2,1,2);plot(c2,abs(Y2(2049:4096)));
axis([1.2,2.2,-30,300]);
title('4096点FFT的频谱');
grid on;
clear all;t=0:0.1:50;yp=sin(2*pi*2*t);yq=sin(2*pi*2.0125*t);y=yp+yq; Y1=fft(y,512);Y2=fft(y,4096);Y1=fftshift(Y1);Y2=fftshift(Y2);c1=[0:255]/51.2;c2=[0:2047]/409;subplot(2,1,1);plot(c1,abs(Y1(257:512)));axis([1.2,2.2,-30,300]);title('512点FFT的频
只是个MATLAB的程序吧.
t=0:0.1:50;t取值 0-50 间隔是0.1
yp=sin(2*pi*2*t); 三角函数 SIN(2π2t)
yq=sin(2*pi*2.0125*t);同上
y=yp+yq; 相加 也就是在yp的图形上加上yq的图形合成
Y1=fft(y,512);用FFT函数进行频谱分析
Y2=fft(y,4096);同上
Y1=fftshift(Y1);让正半轴部分和负半轴部分的图像分别关于各自的中心对称,将fft计算输出的零频移到输出的中心位置
Y2=fftshift(Y2);同上
c1=[0:255]/51.2;C1分别取值除51.2
c2=[0:2047]/409;
subplot(2,1,1);plot(c1,abs(Y1(257:512)));画图函数subplot(2,1,1)表示画2个图 此为第一个
axis([1.2,2.2,-30,300]);横纵坐标起始和结束位置(此为纵坐标1.2-2.2,横坐标-30到300)
title('512点FFT的频谱');给上图命名
grid on;
subplot(2,1,2);plot(c2,abs(Y2(2049:4096)));同上
axis([1.2,2.2,-30,300]);
title('4096点FFT的频谱');
grid on;