今天做信号与系统实验,有用到这部分的知识,感觉这里面道道很多,因此拿出来与大家分享一下。

首先,我们知道matlab是一个很好的编程环境(或许这种解释性语言算不上编程?)。总之用matlab对系统的仿真尤其是对数字信号的处理简单便捷。

今天做了一下傅里叶变换求系统响应这部分的实验。matlab工具箱提供了傅里叶变换的符号求解函数,而数值方法求解则通过傅里叶变换的基本定义来实现,详细如下:

dt = 0.01;

t = 0:dt:8;

et = exp(-3*t).*uCT(t); % 激励信号,其中uCT(t)是我自己写的单位阶跃函数

N = 2000;k = -N:N;w = pi*k/(N*dt);

fet = dt*et*exp(-1i*t'*w);

下面卖点来了,怎样用频域分析法求系统的零输入响应。

我左思右想怎么都觉得用符号法来实现不怎么靠谱。我们知道频域分析法主要体现了系统函数的功能,即幅度加权和相位修正。用符号法的话,怎么能改变’sym’里的幅度和相位呢?

这里就体现了数值法比符号法的优势了,数字信号处理大都使用数值法也是这样的原因:编程人员自由度高,可以做任意的合法操作。

既然确定用数值法来完成频域分析求零输入响应,那么,我们要实现的还是那两句话:幅度加权,相位修正。

先得到系统函数,matlab提供了函数freqs,详细参数自己查吧,这里不啰嗦了。

假设我们已经求出来了系统函数记为H。幅度加权就很简单了。

设激励信号为e(t),那么加权后即为e(t) * abs(H)。

难点还在相位修正。

我们知道有这样的公式:F (w) = |F(w)| * exp(j * Ψ(w))。于是我们可以先求出Ψ(w),在对相位函数进行相位修正Ψ1(w) = Ψ(w) + angle(H(w)),最后计算exp(j * Ψ1(w))。则经过幅度加权和相位修正后的F1(w) = |F1(w)| * Ψ1(w)。在此基础上进行傅里叶逆变换得到零状态响应r(t)。

fai = log(fet./abs(fet))/1i;

fai = fai + angle(H);

fet = abs(fet).*exp(1i.*fai);

frt = abs(H).*fet;

dw = pi/20;

rt = dw*frt*(exp(-1i*t'*w))'/2/pi;

希望本文对正在研究傅里叶变换的同学有帮助吧。

matlab零状态响应幅度频谱,频域分析法求解零状态响应的matlab过程相关推荐

  1. 简单迭代法、牛顿法、弦割法、布洛依登法求解方程或方程组【Matlab】

    利用迭代法求解定非线性方程及方程组,使得误差不超过10^(-8).同时应用迭代加速技术,提交迭代运算效率. 此题需要用到的MATLAB代码及附录: 附录6 二分法作根的隔离 % 附录6 二分法作根的隔 ...

  2. matlab零状态响应幅度频谱,matlab零状态、零输入响应

    <matlab零状态.零输入响应>由会员分享,可在线阅读,更多相关<matlab零状态.零输入响应(9页珍藏版)>请在人人文库网上搜索. 1.1.已知离散时间系统的差分方程为: ...

  3. 卷积法求解系统的零状态响应_连续LTI系统零状态响应求解方法的分析

    连续 LTI 系统零状态响应求解方法的分析 张淑敏 [摘 要] [摘 要]零状态响应是电子技术相关课程中的一个重要概念,本文将 通过对时域分析法和(复)频域分析法求解连续 LTI 系统零状态响应的分析 ...

  4. Chapter7.1:频域分析法理论基础

    该系列博客主要讲述Matlab软件在自动控制方面的应用,如无自动控制理论基础,请先学习自动控制系列博文,该系列博客不再详细讲解自动控制理论知识. 自动控制理论基础相关链接:https://blog.c ...

  5. 卷积法求解系统的零状态响应_【信号与系统】课件16

    信号与系统课件16 1. 连续时间傅里叶变换的性质回顾 傅里叶变换的公式: F(jω) = int(f(t)·exp(-jωt),t,-inf,inf)=FT[f(t)] 傅里叶反变换的公式: f(t ...

  6. 卷积法求解系统的零状态响应_因果系统的零状态响应的一种简易计算方法

    年 月 宜春师专学报 年第 期 . 因果系统的零状态响应的一种简易计算方法 王 卫林 连续线性系统的研究方法 , 可以按系统的数学模型的求解方式分为时域法和变换域法两大类时域法是直接处理系统的数学模型 ...

  7. P7 频域分析法-《Matlab/Simulink与控制系统仿真》程序指令总结

    上一篇 回到目录 下一篇 <Matlab/Simulink与控制系统仿真>程序指令总结 Matlab_Simulink_BookExample 7. 频域分析法 7.1 Matlab 函数 ...

  8. matlab二阶阻尼震荡衰减,二阶欠阻尼电路的零响应输入-MATLAB课程设计.doc

    二阶欠阻尼电路的零响应输入-MATLAB课程设计.doc 摘要 MATLAB具有强大的运算和作图功能,为电路分析引入计算机数值方法提供了技术保证.MATLAB中的变量与常量都是矩阵,其元素可以使复数和 ...

  9. matlab 复频谱图,基于MATLAB实现连续信号与系统复频域分析

    题目:基于MATLAB实现连续信号与系统复频域分析 摘要 拉普拉斯变换(Laplace Transform),是工程数学中常用的一种积分变换.它是为简化计算而建立的实变量函数和复变量函数间的一种函数变 ...

最新文章

  1. 2 CO配置-企业结构-定义-维护成本控制范围(Controlling Area)
  2. Windows10 环境下Jupyter Notebook的安装与使用
  3. 张志华:机器学习的发展历程及启示
  4. 海龟交易策略要点总结
  5. win11升级不满足最低系统要求怎么办 windows11升级不满足最低系统要求的解决方法
  6. [转载] JAVA的throw用法及扩展
  7. 利用PHP的VLD查询OPCODE
  8. 百度网盘客户端使用的一些小窍门
  9. Windows服务器IE浏览器无法下载文件解决方法
  10. 宋朝人物第一,朱熹都说他是“天地间第一流人物”
  11. php自我介绍50字,简短自我介绍50字
  12. 十大流行开源软件推荐
  13. android 获取IP地址
  14. 【CUDA 基础】3.4 避免分支分化
  15. python可视化之matplotlib绘图--蜘蛛侠钢铁侠(组合分区绘图)
  16. 函数型数据主成分分析(FPCA)
  17. pixfllow光流传感器数据获取及悬停程序编写
  18. php filter 文件包含,文件包含lfi
  19. Linux云主机开启IPv6服务
  20. 三菱PLC与第三方设备TCP通讯_PLC | 8台FX3U三菱PLC的无线网络通讯方案

热门文章

  1. linux开发板添加qt库文件_QT程序怎么移植到开发板并运行
  2. iOS获取网页上数据(图片、文字、视频)
  3. 『HarmonyOS』Ability基础(类比Android中Activity学习)
  4. 使用迅雷下载百度网盘2G以上文件
  5. Java新生管理系统的设计与实现
  6. 数论:求一个数的因子专题(因子数,因子和,质因子)
  7. 2021-07-14莫得感情的流水账日记之Day1
  8. 谈谈开发者眼中的编程语言是怎样的?
  9. Processing - 练习(3)飞碟追踪 - (渐进、延迟效果)
  10. tiktok无法关注和喜欢_微软收购Tiktok的出价是该公司最终让您喜欢它的辉煌计划