使用matlab，对以上程序中y函数进行加窗校正，数字滤波，时频分析
fs=180;N=1024;
n=0:N-1;
t=n/fs;
y=0.68*sin(2*pi*10*t+pi/2)+0.45*cos(2*pi*20*t+1/4*pi)+0.28*randn(size(t));
figure(1)
subplot(3,1,1)
plot(t,y)
xlabel('time/s')
ylabel('amplitude/1')
axis()
title('0.68*sin(2*pi*10*t+pi/2)+0.45*cos(2*pi*20*t+1/4*pi)+randn')
a=fft(y,N);
b=abs(a/N);
b1=b(1:N/2+1);
b1(2:end-1)=2*b1(2:end-1);
f=(0:(N/2))*fs/N;
subplot(3,1,2)
plot(f,b1)
xlabel('frequency/Hz')
ylabel('amplitude/1')
title('A0.68 0.45 fs10 20')
axis()
xlabel('frequency/Hz')
ylabel('amplitude/1')
title('corrected')
axis()
使用matlab，对以上程序中y函数进行加窗校正，数字滤波，时频分析
页:
[1]