授之以渔： 卡尔曼滤波器 ...大泻蜜 ...

aaaalook

最小二乘法获取系统传递函数：

1. function [numd, dend] = LeastSquare(x, y, N)
   
2. count = length(y);
   
3. M = count - 1;
   
4. ai = zeros(N*2, M);
   
5. for i=1:N
   
6.     ai(i, i:M) = x(1:(count-i));
   
7.     ai(i+N, i:M) = -y(1:(count-i));
   
8. end
   
9. bi = y(2:(count));
   
10. xd = (inv(ai*ai')*ai*bi)';
    
11. numd = [0 xd(1:N)];
    
12. dend = [1 xd(N+1: N+N)];

复制代码

输入参数：
x  ---  系统输入阵列
y  ---  系统输出阵列
N  ---  系统传递函数的阶数

函数输出：
系统 Z 传递函数分子与分母的系数，即： h(z) = numd(z) / dend(z)

这个函数实在是居家旅行 ..., 啊，是建模仿真必备良器。

篪骋

无与伦比119

以后肯定有用的，mark了

虫总

无图无真相

wenhaoyuan

强贴留名

邛于

学习中。。。。。。

断翅at雄鹰

:victory::Q:Q

津生有你

fanyu19900706

楼主，大神啊

我の小车

顶下，谢了啊