clear ,close all;
%矩形函數(shù)
L=65;%矩形的長度
W=300;%矩形序列的長度
figure(1);%圖表1
n=0:W-1;%分隔數(shù)
x=[(n-L)<0];%n小于L時(shí)x為1
subplot(3,1,1);%3行1列第一個(gè)圖
plot(x);%繪制矩形
grid on;%打網(wǎng)格
title('矩形圖');%標(biāo)題
xlabel('n');%x橫軸
ylabel('x');%縱軸
hold on;%保持

RTG=fft(x,W);%RTG=rectangle
RTGAM=abs(RTG);%取幅值
RTGANG=angle(RTG);%取角度
subplot(3,1,2);%3行1列第2個(gè)圖
plot(n,RTGAM);%繪矩形幅頻特性圖
title('矩形框的幅度');%標(biāo)題
xlabel('頻率');%x橫軸
ylabel('幅度');%縱軸
grid on;%打網(wǎng)格
subplot(3,1,3);%3行1列第3個(gè)圖
plot(RTGANG);%繪矩形相頻特性圖
grid on;%打網(wǎng)格
title('矩形框的相位');%標(biāo)題
xlabel('相位');%x橫軸
ylabel('幅度');%縱軸

%第一個(gè)正弦函數(shù)
f0=20;   %正弦頻率1#
%-------------------------
%前后兩次采樣頻率不一樣就可以分開了
fs=300;   %采樣頻率
%-------------------------
figure(2);%圖表2
x1=[0:L-1];%矩形窗長度
y1=4*sin(2*pi*f0.*x1/fs);%第一個(gè)正弦函數(shù)的表達(dá)式
subplot(2,1,1);%繪圖表2中2行1列第1個(gè)圖
plot(y1);%時(shí)域正弦波形
title('時(shí)域正弦波形');%標(biāo)題
xlabel('n');%x橫軸
ylabel('幅度');%縱軸
grid on;%打網(wǎng)格
yfft1=fft(y1,W);%對(duì)Y1作FFT變換
subplot(2,1,2);%繪圖表2中2行1列第2個(gè)圖
plot(n,abs(yfft1));%繪正弦1的幅頻特性圖
title('FFT變換后的正弦');%標(biāo)題
xlabel('頻率');%x橫軸
ylabel('幅度');%縱軸
grid on;%打網(wǎng)格

%第二個(gè)正弦函數(shù)
f2=21;%正弦頻率2#
figure(3);%圖表3
%-------------------------
%前后兩次采樣頻率不一樣就可以分開了
fs=300; %采樣頻率
%-------------------------
x2=0:L-1;%矩形窗長度
y2=4*sin(2*pi*f2*x2/fs);%第二個(gè)正弦函數(shù)的表達(dá)式
yfft2=fft(y2,W);%對(duì)Y2作FFT變換
yfft2am=abs(yfft2);%取幅值
yfft2angle=angle(yfft2);%取相位
subplot(3,1,1);%繪圖表3中3行1列第1個(gè)圖
plot(n,abs(yfft1),'-r');%繪正弦1的幅頻特性圖
hold on;%保持
plot(n,yfft2am,'-g');%繪正弦2的幅頻特性圖
title('the AM of sin1 and sin2');
xlabel('頻率');
ylabel('幅度');

subplot(3,1,2);%相位信息
plot(angle(yfft1),'-g');%繪正弦1的相頻特性圖
hold on;%保持
plot(yfft2angle,'--r');%繪正弦2的相頻特性圖

%正弦1和正弦2的幅頻特性疊加
y3=abs(yfft1)+abs(yfft2);
subplot(3,1,3);
plot(n,y3);
title('幅頻特性相加');
xlabel('頻率');
ylabel('幅度');
grid on;%打網(wǎng)格