% dm10406
% 傅里葉變換的頻域卷積定理
dt = 0.005;                                     %設(shè)置時間間隔
t = -0.5:dt:0.5;
f = Heaviside(t+0.5)-Heaviside(t-0.5);          %脈寬為1的門信號
f1 = dt*conv(f,f);                              %通過兩個脈寬相等的門信號卷積產(chǎn)生三角波信號
n = -1:dt:1;                                    %三角波信號的取值范圍為門信號取值范圍的兩倍(-1~+1)
subplot(231);
plot(n,f1);                                     %繪制三角波信號
title('三角波信號');
xlabel('t');
f2 = cos(10*pi*n);                               %定義余弦信號
subplot(232);
plot(n,f2);
title('余弦信號');
xlabel('t');
f = f1.*f2;                                     %計算余弦調(diào)幅信號
subplot(233);
plot(n,f);
title('余弦調(diào)幅信號');
xlabel('t');
W1 = 10*pi*2;                                   %設(shè)定采樣角頻率
N = 200;                                        %設(shè)置頻率采樣點數(shù)
k = -N:N;
W = k*W1/N;
F1 = f1*exp(-j*n'*W)*dt;                        %求f1(t)的傅里葉變換F1(ω)
F2 = f2*exp(-j*n'*W)*dt;                        %求f2(t)的傅里葉變換F2(ω)
F = conv(F1,F2)/(2*pi);                         %求F1(ω)和F2(ω)的卷積
F1 = abs(F1);                                   %求f1(t)幅度頻譜
subplot(234);
plot(W/pi,F1);
axis([-20 20 0 1.1]);
title('三角波信號幅度頻譜');
xlabel('\omega    單位pi');
F2 = abs(F2);      %求f2(t)幅度頻譜
subplot(235);
plot(W/pi,F2);
axis([-20 20 0 1.1]);
title('余弦信號幅度頻譜');
xlabel('\omega   單位:pi');
F = abs(F);                 %求f(t)幅度頻譜
subplot(236);
k1 = -2*N:2*N;
W11 = 2*pi*k1/(N*0.1);
plot(W11/pi,F);
axis([-20 20 0 2]);
title('余弦調(diào)幅信號的幅度頻譜');
xlabel('\omega    單位pi');
%End