跪求高手帮忙 。频谱问题。
小弟新手，最近要做一个关于FFT求频谱相位图的课业。给出的正弦波形式为y=10+5*sin(2*pi*50*t+(30/180*pi))+5*sin(2*pi*150*t+(45/180*pi))+5*sin(2*pi*300*t+(60/180*pi)); 求频谱图和相位图。 其中 t=0:1e-5:0.1 然后我根据取样间隔时间，确定FS为100000，因为正弦波形地周期为2000我选择2倍的整周期2096个点。 问题是为了避免频谱泄露，我想加窗函数，应该怎么加，在程序里应该怎么样表现。我自己写了程序，可是频谱图貌似是对的， 可是对应的相位图明显差太多了，请高手指教，下面给出程序。
clear;
fs = 100000;
N = 2048;
t=(0:1e-5:0.1)';
f0=50;
f1=150;
f2=300;
y=10+5*sin(2*pi*f0*t+(30/180*pi))+5*sin(2*pi*f1*t+(45/180*pi))+5*sin(2*pi*f2*t+(60/180*pi));
figure(1);
plot(y)
title('Dat')
ylabel('Amp')
xlabel('N ')
grid
H=fft(y,N);
mag=abs(H)/N;
f=(0:length(H)-1)'*fs/length(H);
figure(2);
plot(f(1:N/2),mag(1:N/2));
xlabel('frequenz(Hz)');
ylabel('Ampli');
grid;
phase = angle(H);
figure(3);
plot(f(1:N/2),phase(1:N/2)*180/pi)
xlabel('Frequency (Hz)')
ylabel('Phase (Degrees)')
grid on
figure(4); clf;
subplot(211), loglog(f(1:N/2),mag(1:N/2)); hold on;
xlabel('frequenz(Hz)');
ylabel('Amplitudgang');
grid;
subplot(212), semilogx(f(1:N/2),phase(1:N/2)*180/pi); hold on;
xlabel('Frequency (Hz)')
ylabel('Phase (Degrees)')
grid on
页:
[1]