Matlab观察高斯脉冲频谱代码

最新推荐文章于 2023-01-12 14:46:14 发布

JMAVI

最新推荐文章于 2023-01-12 14:46:14 发布

阅读量4.8k

点赞数 1

分类专栏：计算电磁学文章标签： matlab plot fft

本文链接：https://blog.csdn.net/wangjiannuaa/article/details/6016631

版权

计算电磁学专栏收录该内容

41 篇文章 11 订阅

订阅专栏

TsampleNums = 1000;
FsampleNums = 10 ;

cc=2.99792458e8;
mur=1.0;
epsr=1.0;

delta=0.001 ;
dt=delta*sqrt(epsr*mur)/(2.0*cc);

ndelay = 30;
t0 = ndelay*dt;
tau = 30*dt;

Taxis = zeros(TsampleNums);
Faxis = zeros(FsampleNums);
J = zeros(TsampleNums);

n=[1:1:TsampleNums];

J = exp(-4*pi*(n*dt-t0).^2 /tau^2);

NFFT=2^nextpow2(TsampleNums);

subplot(2,1,1);

plot((1:1000),J,'-b','LineWidth',2,...
                'MarkerEdgeColor','k',...
                'MarkerFaceColor','g',...
                'MarkerSize',10);
        xlabel('Time   s','fontsize',12);
        ylabel('Amplitude of E/H V/m,A/m','fontsize',12);
        title('TimeDomain','fontsize',12);
grid;

Y = fft(J,NFFT);
fs = 1/dt;
f=fs*(1:400)/(NFFT);
subplot(2,1,2);
plot(f,abs(Y(1:400)),'-b','LineWidth',2,...
                'MarkerEdgeColor','k',...
                'MarkerFaceColor','g',...
                'MarkerSize',10);
        xlabel('Frequency   Hz','fontsize',12);
        ylabel('Amplitude','fontsize',12);
        title('FrequencyDomain','fontsize',12);
grid;