ultrawb.asv

用matlab编写的数字信号处理程序

clc,clear all,close all;

%%%(1111)画出信号的时域波形
T=0.002*10^(-9);                      %抽样时间间隔为0.001ns
t=0:T:10^(-9);                      %t从0到1ns
E=1;                                  %定义个参数 
tao=0.11*10^(-9);
tc=0.5*10^(-9);                       
wt=sqrt(E)*sqrt(8/3/tao)*(1-4*pi*((t-tc)/tao).^2).*exp(-2*pi*((t-tc)/tao).^2);     %构造wt，注意这里t是矢量，要用有关矢量的运算符号
%figure(1);
%plot(t,wt);
%xlabel('Time (ns)');
%ylabel('Amplitude');
%title('the time domain undee');

%%%(22222)画出信号频谱%%%%%%%%%%%%%%%%%%%%%%%
N=length(wt);
WK=fft(wt);
WK1=10^(-4.13)*abs(WK)/max(abs(WK));
WK2=10*log10(WK1);
f=[-N/2:1:N/2]/T/N;
figure(2);
plot(f/e9,fftshift(wt));
grid;
xlabel('Frequency (GHz)');
ylabel('FIRP Emission Level (dbm)');
title('the frequency domain undee');
axis([0 12*10^11]);