使用matlab，对以上程序中y函数进行加窗校正，数字滤波，时频分析

Richard0824

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([0,1,-1.5,2])
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('A  0.68   0.45    fs  10   20')
axis([0,50,0,1])
xlabel('frequency/Hz')
ylabel('amplitude/1')
title('corrected')
axis([0,50,0,1])

使用matlab，对以上程序中y函数进行加窗校正，数字滤波，时频分析