首页/文章/ 详情

Matlab计算频谱、功率谱和功率谱密度

1年前浏览2007

本实例以高斯信号为例,计算幅度谱、相位谱、双边功率谱、双边功率谱密度、单边功率谱、单边功率谱密度。高斯信号的半波全宽FWHM=50ps,中心点位于2.5ns处。


MATLAB程序代码:

%================

%作者:yoyoba

%Email:stuyou@126.com

%================

clc;

clear;

FWHM=50e-12;

%高斯信号FWHM宽度,为50ps

time_window=100*FWHM;

%高斯信号的采样窗口宽度,该值决定了傅里叶变换后的频率分辨率


Ns=2048;

%采样点

dt=time_window/(Ns-1);

%采样时间间隔

t=0:dt:time_window;

%采样时间

gauss_time=exp(-0.5*(2*sqrt(2*log(2))*(t-2.5e-9)/FWHM).^2);

%高斯脉冲,中心位于2.5ns处。


plot(t*1e 9,gauss_time,'linewidth',2.5);

xlabel('Time/ns');

ylabel('Amplitude/V');

title('Gauss pulse');


%================

%以下计算双边谱、双边功率谱、双边功率谱密度

%================

gauss_spec=fftshift(fft(ifftshift(gauss_time))); 

%傅里叶变换,并且进行fftshift移位操作。

gauss_spec=gauss_spec/Ns; 

%求实际的幅度值;

df=1/time_window;

%频率分辨率

k=floor(-(Ns-1)/2:(Ns-1)/2),

double_f=k*df;

%双边频谱对应的频点


figure;%幅度谱

plot(double_f*1e-9,abs(gauss_spec),'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Amplitude/V');

title('double Amplitude spectrum');


figure;%相位谱

plot(double_f*1e-9,angle(gauss_spec),'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Phase/rad');

title('double Phase spectrum');


figure;%功率谱

double_power_spec_W=abs(gauss_spec).^2;

%双边功率谱,单位W;

double_power_spec_mW=double_power_spec_W*1e 3;

%双边功率谱,单位mW;

double_power_spec_dBm=10*log10(double_power_spec_mW);

%双边功率谱,单位dBm;

plot(double_f*1e-9,double_power_spec_dBm,'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Power/dBm');

title('double Power spectrum');


figure;%功率谱密度

double_power_specD_W=abs(gauss_spec).^2/(df);

%双边功率谱密度,单位W/Hz

double_power_specD_mW=double_power_specD_W*1e 3;

%双边功率谱密度,单位mW/Hz

double_power_specD_dBm=10*log10(double_power_specD_mW);

%双边功率谱密度,单位dBm/Hz

plot(double_f*1e-9,double_power_specD_dBm,'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Power/(dBm/Hz)');

title('double power spectrum Density');


%================

%以下计算单边谱、单边功率谱及单边功率谱密度

%================

gauss_spec=fft(ifftshift(gauss_time));

%计算单边谱无需fftshift

gauss_spec=gauss_spec/Ns;

%计算真实的幅度值

single_gauss_spec=gauss_spec(1:floor(Ns/2));

single_f=(0:floor(Ns/2)-1)*df;


figure;%幅度谱

plot(single_f*1e-9,abs(single_gauss_spec),'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Amplitude/V');

title('single Amplitude spectrum');


figure;%相位谱

plot(single_f*1e-9,angle(single_gauss_spec),'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Phase/rad');

title('single Phase spectrum');


figure;%功率谱

double_power_spec_W=abs(gauss_spec).^2;

single_power_spec_W=2*double_power_spec_W(1:floor(Ns/2));

%单边功率谱,单位W

single_power_spec_mW=single_power_spec_W*1e 3;

%单边功率谱,单位mW;

single_power_spec_dBm=10*log10(single_power_spec_mW);

%双边功率谱,单位dBm;

plot(single_f*1e-9,single_power_spec_dBm,'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Power/dBm');

title('single Power spectrum');


figure;%功率谱密度

double_power_specD_W=abs(gauss_spec).^2/(df);

single_power_specD_W=2*double_power_specD_W(1:floor(Ns/2));

%单边功率谱密度,单位W/Hz

single_power_specD_mW=single_power_specD_W*1e 3;

%单边功率谱密度,单位mW/Hz

single_power_specD_dBm=10*log10(single_power_specD_mW);

%单边功率谱密度,单位dBm/Hz

plot(single_f*1e-9,single_power_specD_mW,'linewidth',2.5);

xlabel('Frequency/GHz');

ylabel('Power/(dBm/Hz)');

title('single power spectrum density');



运行结果:

【免责声明】本文来自chainunix网stuyou博客,版权归原作者所有,仅用于学习等,对文中观点判断均保持中立,若您认为文中来源标注与事实不符,若有涉及版权等请告知,将及时修订删除,谢谢大家的关注!  
来源:CAE之家
科普MATLAB
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2022-08-08
最近编辑:1年前
CAE之家
硕士 | CAE仿真负责人 个人著作《汽车NVH一本通》
获赞 1059粉丝 5167文章 773课程 18
点赞
收藏
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈