genpsd.m

本人编写的一些uwb仿真程序,希望对大家有点用处

% function [PSD,df]=genpsd(x,fc)
%  计算信号的功率谱密度并输出图形
%  输入参数
%  x       输入信号
%  fc      输入信号的采样频率
%  
%  返回参数
%  
%  PSD     功率谱密度
%  df      频率分辨率

function [PSD,df]=genpsd(x,fc)
%step1 evaluation of the psd

dt=1/fc;
n=length(x);
t=n*dt;
df=1/t;
x=fft(x);
x=x/n;
mpsd=abs(x).^2/(df^.2);
PSD=fftshift(mpsd);
PSD=(1/t).*PSD;

%step2 graphical representation
 frequency=linspace(-fc/2,fc/2,length(PSD));
 pf=plot(frequency,PSD);
 set(pf,'linewidth',[2]);
 ax=gca;
 set(ax,'fontsize',12);
 x=xlabel('Frequency [Hz]');
 set(x,'fontsize',14);
 y=ylabel('Power Spectral Density [V^2/Hz]');
 set(y,'fontsize',14);