没有合适的资源?快使用搜索试试~ 我知道了~
华工信号与系统实验3 实验报告
资源推荐
资源详情
资源评论
实 验 三 基本信号的产生和实现
一、实验目的
应 用 离 散 傅 里 叶 变 换 (DFT), 分 析 模 拟 信 号
x
(t)的 频 谱 。 深 刻 理 解 利 用 DFT
分 析 模 拟 信 号频 谱 的 原 理 , 分 析 过 程 中 出 现 的 现 象 及 解 决 方 法 。
二、实验原理
连 续 周 期 信 号 相 对 于 离 散 周 期 信 号 ,连 续 非 周 期 信 号 相 对 于 离 散 非 周 期 信 号 ,
都 可 以 通 过 时 域 抽 样 定 理 建 立 相 互 关 系 。因 此 ,在 离 散 信 号 的 DFT 分析 方 法 基 础
上 , 增 加 时 域 抽 样 的 步 骤 , 就 可 以 实 现 连 续 信 号 的 DFT 分 析 。
连 续 周 期 信 号 相 对 于 离 散 周 期 信 号 ,连 续 非 周 期 信 号 相 对 于 离 散 非 周 期 信 号 ,都
可 以 通 过 时 域 抽 样 定 理 建 立 相 互 关 系 。因 此 ,在 离 散 信 号 的 DFT 分 析 方 法 基 础 上 ,
增 加 时 域 抽 样 的 步 骤 , 就 可 以 实 现 连 续 信 号 的 DFT 分 析 。
1. 利 用 DFT 分 析 连 续 周 期 信 号 的 频 谱
分 析 步 骤 为 :
(1) 确 定 周 期 信 号 的 基 本 周 期 ;
(2) 计 算 一 个 周 期 内 的 抽 样 点 数 N。 若 周 期 信 号 的 最 高 次 谐 频 为 p 次 谐 波 ,
则 频 谱 中 有 2p+1 根 谱 线;若 周 期 信 号 的 频 谱 无 限 宽 ,则 认 为 集 中 信 号 90%以 上(或
根 据 工 程 允 许 而 定)能 量 的 前(p+1)次 谐 波 为 近 似 的 频 谱 范 围 ,其 余 谐 波 忽 略 不 计 。
取 N>=2p+1;
(3) 对 连 续 周 期 信 号 以 抽 样 间 隔 T= T0 /N 进 行 抽 样 , 得 到 x[k] ;
(4) 利 用 FFT 函 数 对 x[k]作 N 点 FFT 运 算 , 得 到 X[m];
(5) 最 后 求 得 连 续 周 期 信 号 的 频 谱 为 X(nw0)=X[m]/N。
2. 利 用 DFT 计 算 连 续 非 周 期 信 号 x(t) 的 频 谱
分 析 步 骤 为 :
(1)根 据 时 域 抽 样 定 理 , 确 定 时 域 抽 样 间 隔 T, 得 到 离 散 序 列 x[k];
(2) 确 定 信 号 截 短 的 长 度 M 及 窗 函 数 的 类 型 , 得 到 有 限 长 M 点
离 散 序 列 xM[k]=x[k]w[k];
(3) 确 定 频 域 抽 样 点 数 N, 要 求 N>=M;
(4) 利 用 FFT 函 数 进 行 N 点 FFT 计 算 得 到 N 点 的 X[m];
(5) 由 X[m]可 得 连 续 信 号 频 谱 X(jw)样 点 的 近 似 值 。
三、实验内容(题目+代码+结果截图+实验过程中碰到的问题和如何解决的)
1. 利 用 FFT 分 析 信 号 的 频 谱 。
(1) 确 定 DFT 计 算 的 各 参 数 ( 抽 样 间 隔 , 截 短 长 度 , 频 谱 分 辨 率 等 ) ;
(2) 比 较 理 论 值 与 计 算 值 , 分 析 误 差 原 因 , 提 出 改 善 误 差 的 措 施 。
2. 分 析 周 期 信 号 的 频 谱 时 ,如 果 分 析 长 度 不 为 整
周 期 ,利 用 fft 函 数 计 算 并 绘 出 其 频 谱 ,总 结 对 周 期 信 号 进 行 频 谱 分 析 时 ,如
何 选 取 信 号 的 分 析 长 度 。
3. 假 设 一 实 际 测 得 的 一 段 信 号 的 长 度 为 0.4s, 其 表 达 式 为
x(t)=cos(2pf1t)+0.75 cos(2pf2t)
其 中 f1=100Hz, f2=110Hz。 当 利 用 FFT 近 似 分 析 该 信 号 的 频 谱 时 , 需 要 对 信 号
进 行 时 域 抽 样 。 试 确 定 一 合 适 抽 样 频 率 , 利 用 DFT 分 析 信 号 x(t)的 频 谱 。
若 在 信 号 截 短 时 使 用 Hamming 窗 , 由 实 验 确 定 能 够 分 辨 最 小 谱 峰 间 隔
和 信 号 长 度 的 关 系 。 若 采 用 不 同 参 数 的 Kaiser 窗 , 重 新 确 定 能 够 分 辨 最
小 谱 峰 间 隔 和 信 号 长 度 的 关 系 。
4. 产 生 一 个 淹 没 在 噪 声 中 的 信 号 x(t),例如 由 50Hz 和 120Hz 的 正 弦 信 号 以 及
一 个 零 均 值 的 随 机 噪 声 叠 加 而 成 。确 定 抽 样 间 隔 和 信 号 截 短 长 度 ,分 析 信 号 的 频
谱 , 指 出 50Hz 和 120Hz 的 正 弦 成 分 对 应 的 谱 峰 位 置 , 详 细 写 出 检 测 信 号 的 步 骤
和 原 理 。
1、代码:
fsam=50;
Tp=6;
N=512;
T=1/fsam;
t=0:T:Tp;
x=exp(-2*t);
X=T*fft(x,N); %频域幅度 去归一化,乘以采样间隔 T
subplot(3,1,1);
plot(t,x);
xlabel('t');
ylabel('|X(jw)|')
title('时域波形');
w=(-N/2:N/2-1)*(2*pi/N)*fsam; %频域频率 去归一化,除以 T
y=1./(j*w+2);
subplot(3,1,2);
plot(w,abs(fftshift(X)),w,abs(y),'r-.');
title('幅度谱');
xlabel('w');
ylabel('y=|1/(j*w+1)|')
legend('理论值','计算值');
axis([-10,10,0,1.4])
subplot(3,1,3);
error = abs(abs(fftshift(X))-abs(y));
plot(w,error)
title('理论与计算的误差');
剩余19页未读,继续阅读
资源评论
夜雨聆风400
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功