作业帮 > 综合 > 作业

基于加窗FFT的频谱分析 27号之前 急用

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/10/01 19:49:59
基于加窗FFT的频谱分析 27号之前 急用
设信号为x(t)=24sin(wt)+12sin(3wt+40pi/180)+6sin(5wt+60pi/180)+3sin(7wt+80pi/180)+1.5sin(9wt+120pi/180),w=100pi.
1、分别用矩形窗、汉宁窗、哈明窗、布莱克曼窗+FFT对信号进行频谱分析,分析各窗函数对频谱分析的影响;
2、在采样频率一定时,增加截断时间长度,分析截断时间长度对频谱分析的影响;
3、在截断时间长度一定时,修改采样频率,分析采样频率对频谱分析的影响;
4、给出所有MATLAB程序及相关频谱图;
程序:
Fs=1000;
T=1/Fs;
Tp=0.1;
N=Tp*Fs;
w=100*pi;
n=1:N;
Xn=24*sin(w*n*T)+12*sin(3*w*n*T+40*pi/180)+6*sin(5*w*n*T+60*pi/180)+3*sin(7*w*n*T+(80*pi)/180)+1.5*sin(9*w*n*T+120*pi/180);
Xn=Xn/max(abs(Xn));
Xk=fft(Xn,4096);
fk=Fs*(0:4095);
plot(fk,abs(Xk)/max(abs(Xk)));
xlabel('Hz');
ylabel('幅值');
title('矩形窗');
2)
Fs=1000;
T=1/Fs;
Tp=0.1;
N=Tp*Fs;
w=100*pi;
n=1:N;
Xn=24*sin(w*n*T)+12*sin(3*w*n*T+40*pi/180)+6*sin(5*w*n*T+60*pi/180)+3*sin(7*w*n*T+(80*pi)/180)+1.5*sin(9*w*n*T+120*pi/180);
Xn=Xn/max(abs(Xn)),wn=hanning(N);
Xn1=Xn.*wn';
Xk=fft(Xn1,4096);
fk=Fs*(0:4095);
plot(fk,abs(Xk)/max(abs(Xk)));
xlabel('Hz');
ylabel('幅值');
title('汉宁窗');
剩下的根据这个改改就行了