随机信号的随机处理，以及波形

clf; Fs=2000; N=1024;Nfft=256;n=0:N-1;t=n/Fs; randn('state',0); xn=cos(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N); [Pxx1,f]=pmtm(xn,4,Nfft,Fs);%Thomson多维度法功率谱密度估计 %此处有问题 subplot(121),plot(f,10*log10(Pxx1)); xlabel('频率/HZ');ylabel('功率谱/dB'); title('多窗口法(MTM)NW=4'); grid on; [Pxx,f]=pmtm(xn,2,Nfft, Fs); subplot(122),plot(f,10*log10(Pxx)); xlabel('频率/HZ');ylabel('功率谱/dB'); title('多窗口法(MTM)NW=2'); grid on; Fs); subplot(122),plot(f,10*log10(Pxx)); xlabel('频率/HZ');ylabel('功率谱/dB'); title('多窗口法(MTM)NW=2'); grid on;