filtic函数 matlab_matlab filtic 函数应用 filter 解差分方程 dft 函数
matlab filtic 函数应用 filter 解差分方程 dft 函数
一、 解差分方程
说明都在代码注释里面了
%这里要利用filtic函数 为滤波器的直接II型实现选择初始条件
%求解查分方程 y(n) - 0.4y(n-1) - 0.45y(n-2) = 0.45x(n) +0.4x(n-1) - x(n-2)
%y(-1) = 0 y(-2) = 1 x(-1) = 1 x(-2) = 2
%x(n) = 0.8 ^u(n) 状态方程H(z) = (0.45 + 0.4z^-1 - z^-2)/(1- 0.4z^-1-0.45z^-2)
% x 为滤波前序列 y 为输出 序列
% 总结一下 首先把num 和den 写出来
num = [ 0.45 0.4 -1];
den = [1 -0.4 -0.45]
% 其次把初始条件写来
x0 = [1 2] ;
y0= [0 1];
N = 50 ;
n = [1 :N-1]';
x = 0.8.^n;
%生成初始条件
Zi = filtic(num, den , y0 , x0);
[y , Zf] = filter(num , den ,x, Zi);
plot(n , x ,'R-', n, y, 'b--');
xlabel('n'); ylabel('(n)--y(n)');
legend('Input x' , 'Output y', 1);
grid;运行结果:
二、和上面类似简单的解差分方程
%这里要借助filter 和filtic 解差分方程
%假设输入信号为 x = 0.8.^n;
%求输出序列
%差分方程y(n) = x(n) - 5x(n-1)+ 8x(n-3)
num = [1 -5 0 8];
N = 30 ;
n =[0: N-1] ;
x = 0.8.^n;
y = filter(num , 1, x);
stem(n, y) ;
title('Response');
grid
运行结果:
三、dft 函数
%这里写出一个dft函数
% Xk = DFT[xn] = 公式
% 总结一下dft 写法首先确定xn N WN = exp(-j*2*pi/N)
%其次 确定nk = n'* k 再确定WNnk Xk = xn*WNnk
function [Xk] = dft(xn, N)
n = [0 :1 : N-1];
k = n ;
WN = exp(-j*2*pi/N);
nk = n'*k;
WNnk = WN.^nk;
Xk = xn*WNnk;
end
四、dft 应用
% x(n) = cos(0.48*pi*n) + cos(0.52*pi*n)
%绘制x(n) Xk(k) dft 变换
clf;
N = 100;
n = 0: N-1;
xn = cos(0.48*pi*n) + cos(0.52*pi*n);
Xk = dft(xn,N);
%幅度值 相位值
magXk = abs(Xk);
phaXk = angle(Xk);
% 原函数图像
subplot(221);
plot(n,xn);
xlabel('n'); ylabel('x(n)');
title('xn');
%幅度频率特性
subplot(222);
k = 0:length(magXk) - 1;
plot(k, magXk);
xlabel('k');
ylabel('abs(Xk)');
title('Xk');
运行结果:
filtic函数 matlab_matlab filtic 函数应用 filter 解差分方程 dft 函数相关推荐
- matlab 解函数方程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
MATLAB程序设计教程(7)--MATLAB解方程与函数极值 第7章MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7. ...
- matlab求函数极值教程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
MATLAB程序设计教程(7)--MATLAB解方程与函数极值 第7章MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7. ...
- Python中lambda详解(包括内置函数map、reduce、filter、sorted、max)
文章目录 一.lambda是什么? 1.lambda语法 2.语法详解 二.lambda的使用 1.定义 2.调用 3.替换 4.作返回值 三.lambda作参数 1.map函数 2.reduce函数 ...
- 【JavaScript 笔记】— 函数高级(变量作用域、解构赋值、方法、高阶函数、闭包、箭头函数、generator)
JavaScript个人笔记 变量作用域 变量提升 全局对象 windows 命名空间 局部作用域 常量 解构赋值 使用场景 方法 apply 装饰器 高阶函数(Array) map reduce m ...
- map函数的用法python,详解Python map函数及Python map()函数的用法
python map函数 map()函数 map()是 Python 内置的高阶函数,它接收一个函数 f 和一个 list,并通过把函数 f 依次作用在 list 的每个元素上,得到一个新的 list ...
- Python 常用内置函数map、zip、filter、reduce、enumerate
Python 中有许多非常实用的内置函数,通过这些函数我们可以方便的实现某些功能,下面就列举一些常用的内置函数. 1. map() 函数 map() 可以根据提供的函数对指定序列做映射,它接受一个函数 ...
- linux动态库注册函数,Linux动态库函数的详解
linux动态库函数的详解 加载动态库 void *dlopen(const char *filename, int flag); flag的可能值: rtld_lazy rtld_now rtld_ ...
- python while函数_详解python while 函数及while和for的区别
1.while循环(只有在条件表达式成立的时候才会进入while循环) while 条件表达式: pass while 条件表达式: pass else: pass 不知道循环次数,但确定循环条件的时 ...
- php调用C代码的方法详解和zend_parse_parameters函数详解
来源:http://my.oschina.net/Customs/blog/490873 http://blog.csdn.net/super_ufo/article/details/3863731 ...
最新文章
- Silverlight/Windows8/WPF/WP7/HTML5周学习导读(6月25日-7月1日)
- JAVA修饰符类型(public,protected,private,friendly)
- python web为什么不火-Python这么火,为何有人说Python不好找工作?
- 模糊数学笔记大全(模糊集、截集、模糊矩阵、模糊关系、模糊聚类、模糊模型识别、模糊综合决策、模糊数)
- 深度学习核心技术精讲100篇(十九)--GBDT(梯度提升树) 和 Resnet (残差网络)的原理
- 200723学习日报
- scrollview下拉刷新_SwiftUI之View Tree 实战3(下拉刷新)
- 一滴水从高处落下来,会不会砸死人?
- 腐蚀单机怎么进_暖气片堵塞是什么原因?要怎么解决呢?
- linux 编码转换-转
- mysql数据类型范围导致失败
- 信息安全期末复习整理
- 基础物理-机械振动与波
- Excel2013每次打开都弹出配置进度窗口的分步解决办法
- java运维工程师简历模板_系统运维工程师个人个人简历模板.doc
- dataframe排序 pd.rank()
- 服装管理系统c语言(大一)
- android 蓝牙串口自动断开,uni-app 蓝牙连接被自动断开问题
- Liunx磁盘管理——LVM
- Linux命令之chown
热门文章
- 有问有答 | 容器精华问答,如何玩转容器服务?
- 如果有这样一台服务器……
- h5 history
- python二叉树去重_python二叉树学习分享
- mysql 闪回_MySQL数据误删除的快速解决方法(MySQL闪回工具)
- Oracle 表空间常用sql
- JAVA实现AES 解密报错Input length must be multiple of 16 when decrypting with padded cipher
- 第七篇:Spring Boot整合Thymeleaf_入门试炼03
- 计算机专业课程群建设,计算机科学与技术专业课程群建设的研究与实践
- mysql pt_MySQL慢查询之pt-query-digest分析慢查询日志