clc,clear;
fs=100e9;
T=2e-6;
t=0:1/fs:T;
tao=20e-12;xita=3.18e-9;
beita=4;
T1=5e-9;
fai=pi/4;
h=beita*(exp(-t/tao)/tao-exp(-t/xita)/xita);
x=ones(1,length(t));
x=conv(h,(cos(x+fai)).^2);
for i=2:length(t)
end
% for i=1:1000
% tmp(i+1)=t(i+1)-T1;
% x(i)=conv(beita*(exp(-tmp(i+1)/tao)/tao-exp(-tmp(i+1)/xita)/xita),(cos(x(i)+fai)).^2);
% x(i+1)=conv(beita*(exp(-t(i+1)/tao)/tao-exp(-t(i+1)/xita)/xita),(cos(x(i)+fai)).^2);
% y(i+1)=x(i+1);
% % c=xm(i);
% % x(i+1)=beita*(exp(-t(i+1)/tao)/tao-exp(-t(i+1)/xita)/xita).*(cos(c+fai)).^2;
% % tm(i+1)=t(i+2)-T1;
% % xm(i+1)=beita*(exp(-tm(i+1)/tao)/tao-exp(-tm(i+1)/xita)/xita).*(cos(c+fai)).^2;
% end
% plot(x);
% title('输入信号时域');
% ylabel('强度'); xlabel('时间/ms');
% X=fft(x);
% X=abs(X);%求振幅大小
% Y=10*log(fftshift(X));
% figure(1);
% subplot(1,2,1);
% subplot(1,2,2);
% k=fs*linspace(0,1,length(t))-fs/2;
% plot(k,Y);
% title('输入信号频域');
% ylabel('强度'); xlabel('频率/Hz');