如何将时域转模型化为频域模型

lx2008

大家好！ 有几个问题请教高手
在使用函数oe做参数辨识的时候，返回了多项式后，想通过这个多项式推出传递函数，但是根据帮助转化的公式，
与直接使用pem得到的传函作bode图比较，总是不一样！
用帮助中的例子也是这个问题。 请问是不是单位有问题？请高手指点。
还有个细节问题不太明白，按oe的帮助中，转化过程是这样的：
y(t)=B(q)/F(q)*u(t-nk)+e(t)
B(q)=b1 + b2*q^(-1) + ... + bnb*q^(-nb+1)
F(q)=1 + f1 q^(-1) + ... +fnf*q^(-nf)

               B(s)                bnb*s^(nb-1)+bnb-1*s^(nb-1)+...+b1
G(s)=------------- = ---------------------------------------------------                                          
               F(s)                  s^nf + fn* s^(nf-1) + .... + f1

在帮助中这个例子中（下面黑体），nb=2, nf=2,
按上述公式，应该推出：B(q)最高阶为-1，F(q)最高阶为-2， 而用matlab运行后的结果却都是最高-2，这是第二个问题。
因为有第二个问题，在代入G(s)中时就更迷惑了。

                                          0.4974*s^2+0.9982*s                                
不知应该是：G(s)= -------------------------------------------，                     
                                           s^2 + 0.6998* s - 1.499

                                         0.4974*s+0.9982
还应该是：G(s)= ------------------------------------------                                 
                                        s^2 + 0.6998* s - 1.499

但不管是哪个，好像都与下面pem得到的传函不同！请问可以从oe计算结果推出传递函数么？应该怎么做？谢谢！

====================================================
下面是帮助中例子的一部分，对同一个频率特性数据用oe和用pem的结果：
m1 = oe(gfr,[2 2 1]) % Discrete-time Output error (transfer function) model

Discrete-time IDPOLY model: y(t) = [B(q)/F(q)]u(t) + e(t)
B(q) = 0.9982 q^-1 + 0.4974 q^-2                                                                                 
F(q) = 1 - 1.499 q^-1 + 0.6998 q^-2

lx2008

自己顶一下

ljelly

从算法上来讲，oe和pem和armax用的是同样的算法
具体的区别应该在它们各自实现的各项参数指标上，侧重点不同，结果有些微的差别
要不然也就不用做这么多函数，只做一个不就完了

其次，从传递函数的角度讲，Z变换的形式有几种方法，以及系统是否是最小相位系统来说
bode图是有可能不一样的

具体的应该仔细研究一下，这几个函数的使用，找出他们的侧重点，针对什么样的问题，是否能
解决同样的问题