麻烦大家谁帮把这个程序注释一下吧,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))
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/18 21:55:32
麻烦大家谁帮把这个程序注释一下吧,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))
麻烦大家谁帮把这个程序注释一下吧,
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))
只是个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;