ptrans.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 26 行

%MATLAB PROGRAM ****
%CSD*********
%Create input signal
clf
N=1024;
Nfft=256;
window=hanning(256);
noverlap=128;
dflag='none';
Fs=1000;
n=0:N-1;
t=n/Fs;
randn('state',0)
xn=sin(2*pi*100*t);
yn=2*sin(2*pi*50*t)+0.5*rand(1,N);

[pxy,f]=psd(yn,Nfft,Fs,window,noverlap,dflag);

subplot(211)
plot(f,10*log10(abs(pxy)));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Frequency Response')
grid