%Samp3_11
N=100; dt=0.01;   %數據點數和采樣間隔
n=0:N-1; t=n*dt;   %時間序列
f=n/(N*dt);     %頻率序列
h=[1 zeros(1,N-1)];  %脈沖信號第一個值為1，其余為零
fo=2;f1=10;          % 頻率漸增函數鳥鳴信號從2Hz增加到10Hz
x=chirp(t,fo,1,f1);    %在時間序列t上產生頻率漸增信號—鳥鳴信號
xh=conv(x,h);    %將鳥鳴信號和脈沖信號進行卷積，參看上面的用法
XH=dfs(xh(1:N),N);  %由于卷積后數據變長，這里只選取與原數據長度相等的數據個數
subplot(2,2,1),plot(f,real(XH)*2/N); %繪出卷積后Fourier變換的實部
xlabel('頻率/Hz');title('信號卷積后Fourier變換XH的實部');
subplot(2,2,2),plot(f,imag(XH)*2/N); %繪出卷積后Fourier變換的虛部
xlabel('頻率/Hz');title('信號卷積后Fourier變換XH的虛部');
X=dfs(x,N);  %對chirp信號進行Fourier變換
H=dfs(h,N);   %對脈沖信號進行Fourier變換
XH1=X.*H;   %對鳥鳴信號和脈沖信號的Fourier變換相乘
subplot(2,2,3),plot(f,real(XH1)*2/N); %繪出Fourier變換后乘積的實部
xlabel('頻率/Hz');title(' 頻率域乘積XH1的實部')
subplot(2,2,4),plot(f,imag(XH1)*2/N); %繪出Fourier變換后乘積的虛部
xlabel('頻率/Hz');title(' 頻率域乘積XH1的虛部');