【信号与系统】【北京航空航天大学】实验四、幅频、相频响应和傅里叶变换

最新推荐文章于 2024-07-24 14:54:00 发布
最新推荐文章于 2024-07-24 14:54:00 发布
阅读量678 收藏 8
点赞数 7
分类专栏： MATLAB 信号与系统 文章标签： 算法 matlab 开发语言
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43031313/article/details/135617859
版权

一、实验目的

1、 掌握利用MATLAB计算系统幅频、相频响应的方法；
2、 掌握使用MATLAB进行傅里叶变换的方法；
3、 掌握使用MATLAB验证傅里叶变换的性质的方法。

二、实验内容

1、
在这里插入图片描述
MATLAB代码：

>> clear all;
>> a = [1 3 2];
>> b = [1];
>> w = -2 * pi : 0.01 : 2 * pi;
>> [H, w] = freqs(b, a, w);
>> mag = abs(H);
>> pha = angle(H);
>> rea = real(H);
>> ima = imag(H);
>> figure
>> subplot(2, 1, 1);
>> plot(w, mag);
>> grid on;
>> xlabel('\omega/\omega_0');
>> ylabel('|H(jw)|');
>> title('幅频响应曲线');
>> axis([-2 * pi, 2 * pi, 0, 0.6]);
>> subplot(2, 1, 2);
>> plot(w, pha);
>> grid on;
>> xlabel('\omega/\omega_0');
>> ylabel('∠H(jw) (rad)');
>> title('相频响应曲线');
>> axis([-2 * pi, 2 * pi, 0, 0.6]);
>> axis([-2 * pi, 2 * pi, -3, 3]);
>> figure
>> subplot(2, 1, 1);
>> plot(w, rea);
>> axis([-2 * pi 2 * pi -0.1 0.6]);
>> grid on;
>> xlabel('\omega/\omega_0');
>> ylabel('Re{H(jw)}');
>> title('频率响应的实部');
>> subplot(2, 1, 2);
>> plot(w, ima);
>> xlabel('\omega/\omega_0');
>> ylabel('Im{H(jw)}');
>> axis([-2 * pi 2 * pi -0.4 0.4]);
>> grid on;
>> title('频率响应的虚部');

幅频响应、相频响应曲线：

在这里插入图片描述

频率响应的实部和虚部：

在这里插入图片描述

2、

在这里插入图片描述

MATLAB代码：

（1）X1(t) 信号

>> clear all;
>> syms x;
>> f = exp(-x) .* heaviside(x);
>> F = fourier(f);
>> Fs = abs(F);
>> Fe = angle(F);
>> figure
>> subplot(2, 1, 1);
>> ezplot(Fs, [-2 * pi, 2 * pi]);
>> title('符号计算X1(jw)的幅度谱', 'FontSize', 14);
>> grid on;
>> xlabel('\omega/\omega_0');
>> ylabel('|X1(jw)|');
>> subplot(2, 1, 2);
>> ezplot(Fe, [-2 * pi, 2 * pi]);
>> title('符号计算X1(jw)的相位谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('∠X1(jw) (rad)');
>> grid on;
>> clear all;
>> T = 0.01;
>> dw = 0.01;
>> t = -5 : T : 5;
>> w = -2 * pi : dw : 2 * pi;
>> f = exp(-t) .* (t >= 0);
>> F = f * exp(-j * t' * w) * T;
>> Fe = angle(F);
>> Fs = abs(F);
>> figure
>> subplot(2, 1, 1);
>> plot(w, Fs);
>> axis([-2 * pi 2 * pi 0.1 1.1]);
>> grid on;
>> title('数值计算X1(jw)的幅度谱', 'FontSize', 14);
>> ylabel('|X1(jw)|');
>> xlabel('\omega/\omega_0');
>> subplot(2, 1, 2);
>> plot(w, Fe);
>> axis([-2 * pi 2 * pi -1.5 1.5]);
>> title('数值计算X1(jw)的相位谱', 'FontSize', 14);
>> grid on;
>> ylabel('∠X1(jw) (rad)');
>> xlabel('\omega/\omega_0');

（2）X2(t) 信号

>> clear all;
>> syms x;
>> f = exp(-3 * abs(x)) .* sin(2 * x);
>> F = fourier(f);
>> figure
>> Fs = abs(F);
>> Fe = angle(F);
>> subplot(2, 1, 1);
>> ezplot(Fs);
>> grid on;
>> title('符号计算X2(jw)的幅度谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('|X2(jw)|');
>> subplot(2, 1, 2);
>> ezplot(Fe);
>> fplot(Fe);
>> axis([-2 * pi 2 * pi -1.6 1.6]);
>> axis([-2 * pi 2 * pi -2 2]);
>> grid on;
>> title('符号计算X2(jw)的相位谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('∠X2(jw) (rad)');
>> clear all;
>> T = 0.01;
>> dw = 0.01;
>> t = -10 : T : 10;
>> w = -2 * pi : dw : 2 * pi;
>> f = exp(-3 * abs(t)) .* sin(2 * t);
>> F = f * exp(-j * t' * w) * T;
>> figure
>> Fs = abs(F);
>> Fe = angle(F);
>> subplot(2, 1, 1);
>> plot(w, Fs);
>> grid on;
>> axis([-2 * pi 2 * pi -0.01 0.23]);
>> title('数值计算X2(jw)的幅度谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('|X2(jw)|');
>> subplot(2, 1, 2);
>> plot(w, Fe);
>> grid on;
>> axis([-2 * pi 2 * pi -2 2]);
>> title('数值计算X2(jw)的相位谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('∠X2(jw) (rad)');

（3）X3(t) 信号

>> clear all;
>> syms x;
>> f = heaviside(x + 0.5) - heaviside(x - 0.5);
>> F = fourier(f);
>> Fs = abs(F);
>> Fe = angle(F);
>> figure
>> subplot(2, 1, 1);
>> ezplot(Fs, [-5 * pi, 5 * pi]);
>> grid on;
>> title('符号计算X3(jw)的幅度谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('|X3(jw)|');
>> subplot(2, 1, 2);
>> ezplot(Fe, [-5 * pi, 5 * pi]);
>> grid on;
>> title('符号计算X3(jw)的相位谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('∠X3(jw) (rad)');
>> clear all;
>> T = 0.01;
>> dw = 0.01;
>> t = -5 : T : 5;
>> w = -5 * pi : dw : 5 * pi;
>> f = Heaviside(t + 0.5) - Heaviside(t - 0.5);
>> F = f * exp(-j * t' * w) * T;
>> Fs = abs(F);
>> Fe = angle(F);
>> figure;
>> subplot(2, 1, 1);
>> plot(w, Fs);
>> axis([-5 * pi 5 * pi -0.1 1.1]);
>> grid on;
>> title('数值计算X3(jw)幅度谱', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('|X3(jw)|');
>> subplot(2, 1, 2);
>> plot(w, Fe);
>> plot(w, abs(Fe));
>> axis([-5 * pi 5 * pi -0.2 3.3]);
>> grid on;
>> xlabel('\omega/omega_0');
>> xlabel('\omega/\omega_0');
>> ylabel('∠X3(jw) (rad)');
>> title('数值计算X3(jw)相位谱', 'FontSize', 14);

生成的谱图：

（1）X1(t) 信号

符号计算方法：

在这里插入图片描述

数值计算方法：

在这里插入图片描述
（2）X2(t) 信号

符号计算方法：

在这里插入图片描述

数值计算方法：

在这里插入图片描述
（3）X3(t) 信号

符号计算方法：

在这里插入图片描述
数值计算方法：

在这里插入图片描述

3、

在这里插入图片描述

（1）、验证线性性质

MATLAB代码：

>> clear all;
>> syms t;
>> x1 = exp(-t) .* heaviside(t);
>> x2 = exp(-3 * abs(t)) .* sin(2 * t);
>> x = 2 * x1 + 5 * x2;
>> X = fourier(x);
>> figure
>> subplot(2, 1, 1);
>> aX = abs(X);
>> pX = angle(X);
>> ezplot(aX);
>> grid on;
>> title('直接计算傅里叶变换的实部', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('|X(jw)|');
>> subplot(2, 1, 2);
>> ezplot(pX);
>> grid on;
>> title('直接计算傅里叶变换的虚部', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('∠X(jw) (rad)');
>> X1 = fourier(x1);
>> X2 = fourier(x2);
>> Y = 2 * X1 + 5 * X2;
>> aY = abs(Y);
>> pY = angle(Y);
>> figure
>> subplot(2, 1, 1);
>> ezplot(aY);
>> grid on;
>> title('使用性质计算傅里叶变换的实部', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('|X(jw)|');
>> subplot(2, 1, 2);
>> ezplot(pY);
>> grid on;
>> title('使用性质计算傅里叶变换的虚部', 'FontSize', 14);
>> xlabel('\omega/\omega_0');
>> ylabel('∠X(jw) (rad)');

直接计算：

在这里插入图片描述

使用性质计算：

在这里插入图片描述

（2）验证卷积性质：

// 待更新。

4、

在这里插入图片描述

MATLAB代码：

signal_DTFT函数：

function X = signal_DTFT(xn, n, w)
X = xn * (exp(-1i) .^ (n' * w));

其余代码：

>> clear all;
>> w = -4 * pi : 0.01 : 4 * pi;
>> n = -30 : 30;
>> xn = 0.5 .^ (n - 1) .* (n >= 1);
>> X = signal_DTFT(xn, n, w);
>> magX = abs(X);
>> figure
>> subplot(2, 1, 1);
>> plot(w, magX);
>> grid on;
>> xlabel('\omega/\omega_0');
>> ylabel('|X(e^{jw})|');
>> title('傅里叶变换的模值', 'FontSize', 14);
>> axis([-4 * pi, 4 * pi, 0.5, 2]);
>> subplot(2, 1, 2);
>> plot(w, phaX);
>> axis([-4 * pi, 4 * pi, -4, 4]);
>> grid on;
>> xlabel('\omega/\omega_0');
>> ylabel('∠X(e^{jw}) (rad)');
>> title('傅里叶变换的相位', 'FontSize', 14);

傅里叶变换结果：

在这里插入图片描述
5、思考题

在这里插入图片描述

答： 数值计算和符号计算的精度不同。符号计算的精度取决于MATLAB中 ezplot 函数的精度，而数值计算的精度取决于自己定义的 w 变量的精度。

三、实验收获与感想：

1、 单位阶跃函数 u(t) 可以用 heaviside(t) 函数近似表示。在 t = 0 时刻的函数值为 0.5. 举例：

>> clear all;
>> syms x;
>> u = heaviside(x);
>> ezplot(u, [-5, 5]);
>> xlabel('t');
>> ylabel('u(t)');
>> title('单位阶跃函数', 'FontSize', 14);
>> heaviside(0)

ans =

0.5000

函数图像：

在这里插入图片描述

不是AI
关注
  • 7
    点赞
  • 踩
  • 8
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
matlab幅频图
05-01
这是关于matlab幅频率图的代码，希望能给同学们一下心得体会
数字 一阶低通滤波器 matlab 幅频响应代码
12-15
数字 一阶低通滤波器 matlab 幅频响应代码，详情见博客：数字 一阶低通滤波器 详细分析 冰三点水
信号与系统matlab学习（四）
m0_58970439的博客
03-13 607
线性时不变系统的频域分析法是一种变换域分析法，它把时域中求解响应的问题通过傅立叶变换转换成频域中的问题。整个过程在频域内进行，因此它主要是研究信号频谱通过系统后产生的变化，利用频域分析法可分析系统的频率响应、波形失真、物理可实现等实际问题。在时域分析中我们知道，系统的零状态响应等于输入信号与系统单位冲激响应的卷积积分，即若e(t)、h(t)的傅立叶变换均存在，由时域卷积定理可知称为系统函数，由于它是频率的函数，故又称为频率响应函数。一般为频率的复函数，又可写为幅值与相位的形式。
MATLAB 离散系统变换域分析
vizio.blog.csdn.net
11-12 1044
MATLAB 离散系统变换域分析（华南理工大学DSP实验）
【数字信号处理】序列傅里叶变换 ( 傅里叶变换实例 | 矩形窗函数 | 傅里叶变换 | 傅里叶变换幅频特性 | 傅里叶变换相频特性 )
让 学习 成为一种 习惯 ( 韩曙亮 の 技术博客 )
03-07 1万+
一、序列傅里叶变换实例、 1、傅里叶变换、 2、傅里叶变换幅频特性、 3、傅里叶变换相频特性、
matlab计算数字滤波器的幅频响应
坚持
03-11 1万+
1 freqz uses an FFT-basedalgorithm to calculate the Z-transform frequency response of a digitalfilter. Specifically, the statement [h,w] = freqz(b,a,p) returns the p-point complex fr
MATLAB画频率响应曲线（幅频特性和相频特性）并将横坐标转换为赫兹hz单位
热门推荐
Mr.horse的博客
11-26 10万+
matlab画频率响应曲线的函数为： [h,w] = freqz(b,a,n) b,a：传递函数系数 h：频率响应 w：角频率，0~π 更多参数解释参考官方链接https://ww2.mathworks.cn/help/signal/ref/freqz.html?requesteddomain=true 用freqz画频率响应曲线的一个例子：绘制如下系统的频响曲线： H(z)=(1-0.5z^-1) B=[1 -0.5]; A =[1]; [H,w]=freqz(B,A); Hf=a
【信号与系统实验】实验四 傅里叶变换、系统的频域分析
06-14
实验四的主题聚焦于傅里叶变换以及系统的频域分析，这是信号处理和控制系统理论中的核心概念。傅里叶变换是一种将时域信号转化为频域表示的方法，它在通信工程、图像处理、音频处理等多个领域有着广泛应用。理解...
信号与系统实验傅里叶变换.docx
03-22
信号与系统实验傅里叶变换 在本实验中，我们将学习如何使用 MATLAB 实现傅里叶变换，并绘制信号的幅度频谱和相位频谱。 知识点1：傅里叶变换 傅里叶变换是一种将信号从时域转换到频域的技术。它可以将信号表示...
北京航空航天大学《复变函数与积分变换》期末复习资料.pdf
08-17
最后，试卷的三、四、五题是具体的计算题，包括闭曲线上的积分、幂级数展开以及傅立叶变换和傅立叶积分的计算，这些都需要应用到复变函数的基本定理和方法。 总的来说，《复变函数与积分变换》这门课程涵盖了许多...
信号幅频相频特性的画法(频率响应法)
07-18
信号幅频相频特性的画法(频率响应法)。频率响应法和信号幅频相频特性的画法。
傅里叶变换求解幅频相频
10-12
总的来说，这个项目涉及到使用MATLAB进行傅里叶变换，提取信号的幅频和相频特性，可能还包括系统辨识、噪声注入和频带分析等复杂操作。通过这些步骤，我们可以深入理解信号的频域特性，并对系统的行为进行建模和分析...
北京航空航天大学复变函数与积分变换的一套习题.7z
04-09
北京航空航天大学作为中国顶尖的工科院校，其在复变函数与积分变换的教学上具有较高的水平，这套习题集就是对这一课程知识的深入探讨和实践训练。 复变函数是研究复数作为自变量的函数，它涉及到复数的解析性、奇异...
Matlab 求 LTI 系统的系统函数的幅频特性和相频特性
洪源兄
12-01 7229
Matlab 求 LTI 系统的系统函数的幅频特性和相频特性，信号与系统的频谱分析就是将信号与系统的时域表征经过傅里叶变换转换到频域表征，通过幅度频谱和相位频谱，可以获得信号与系统在频域的分布特性。
MATLAB控制系统仿真8.4：控制系统校正的频率响应法
FUXI_Willard
03-27 1254
该系列博客主要讲述Matlab软件在自动控制方面的应用，如无自动控制理论基础，请先学习自动控制系列博文，该系列博客不再详细讲解自动控制理论知识。
【MATLAB】FFT算法计算幅频谱和相位谱
2201_75598815的博客
04-05 3318
以有限长序列xn＝［1，0.5，0，0.5，1，1，0］为例计算DFT和IDFT、幅度谱和相位谱。
幅频响应 matlab画法,matlab绘制数字滤波器的幅频响应曲线
weixin_39841709的博客
03-16 5914
FREQZ 是计算数字滤波器的频率响应的函数 [H,W] = FREQZ(B,A,N)可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，即可绘制出系统在或 范围内的频响曲线。In its simplest form,freqzaccepts the filtercoefficient vectorsbanda,and an i...
【STM32F407的DSP教程】第28章 FFT和IFFT的Matlab实现（幅频响应和相频响应）
Simon223的博客
04-28 8966
完整版教程下载地址：http://www.armbbs.cn/forum.php?mod=viewthread&tid=94547 第28章 FFT和IFFT的Matlab实现（幅频响应和相频响应） 本章主要讲解fft，ifft和fftshift在matlab上的实现。 目录 第28章 FFT和IFFT的Matlab实现（幅频响应和相频响应） 28.1 初学者重要提示 28.2 Matlab的FFT函数 28.2.1 函数语法 28.2.2...
代码随想录算法训练营第22天-leetcode-二叉树08：669. 修剪二叉搜索树；108.将有序数组转换为二叉搜索树；538.把二叉搜索树转换为累加树
最新发布
weixin_65180740的博客
07-24 936
你可能需要改变树的根节点，所以结果应当返回修剪好的二叉搜索树的新的根节点。给出二叉 搜索 树的根节点，该树的节点值各不相同，请你将其转换为累加树（Greater Sum Tree），使每个节点 node 的新值等于原树中大于或等于 node.val 的值之和。题目要求，实际上是每个节点加上中序排列后面的那个节点，故而可以反向中序排序，即RNL，用prior记录前一个节点的值，加上前一个节点的值即可。左右子树也必须是二叉搜索树。后序遍历每个节点，对于在范围外的节点，考虑保留其左子树/右子树（只可能保留一个）
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

不是AI

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值