ultrawb.m

用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);
W=fft(wt);
W=10^(-4.13)*abs(W)./max(abs(W)); 
W=log10(W)*10;               
f=[-N/2:1:N/2-1]/T/N;
figure(2);
plot(f/1.0e9,fftshift(W));  
grid on;
axis([0 12 -80 -40]);