chouyang.m

离散傅里叶变换的程序

clear all;                   %清空matlab工作空间所有变量
close all;                   %关闭所有绘图对话框 

fs = 4;                      %采样频率fs = 4Hz    
dt = 0.001;                  %定义增量（步长）为0.001，
t  = 0:dt:(50-1)/fs;         %计算出需要绘制的模拟信号上点的X轴坐标值 t值为0,0.001,0.002 .....12.248，12.249，12.25     
xt = 1.5*sin(2.5*pi*t);      %原模拟信号为1.5sin(2.5*pi*t)，计算出需要绘制的模拟信号上点的Y轴坐标值

Ts = 1/fs;                   %采样时间间隔
n = 0:50-1;                  %采样50个点  
nTs = n*Ts;                 
st = 1.5*sin(2.5*pi*nTs);    %将原模拟信号1.5sin(2.5*pi*t)中的t用nT取代完成抽样

figure(1);                   %创建一个新的图形对话框 
subplot(3,1,1);              %将已创建的图形对话框的分割成3行1列的子区域，把从上到下第一个子窗口作为当前活动窗口
plot(t,xt);                  %画原模拟信号的波形图
title('原始信号X(t)');       %把当前绘图坐标的上方的标题设置为“ 原始信号X(t)” 

subplot(3,1,2);              %将已创建的图形对话框的分割成3行1列的子区域，把从上到下第二个子窗口作为当前活动窗口 
stem(nTs,st,'.');            %画原模拟信号经过采样过的波形图
title('信号X(nTs)');         

subplot(3,1,3);              %将已创建的图形对话框的分割成3行1列的子区域，把从上到下第三个子窗口作为当前活动窗口
stem(0:length(nTs)-1,st,'.');%画原模拟信号经过采样过的数字信号X(n)的离散波形图，length(nTs)函数得到nTs的长度
title('信号X(n)');