1
.TH mrfit 1 "April 1993" "Scilab Group" "Scilab Function"
4
mrfit - frequency response fit
8
[num,den]=mrfit(w,mag,order)
9
sys=mrfit(w,mag,order,weight)
10
[num,den]=mrfit(w,mag,order,weight)
15
: positive real vector of frequencies (Hz)
18
: real vector of frequency responses magnitude (same size as \fVw\fR)
21
: integer (required order, degree of \fVden\fR)
24
: positive real vector (default value \fVones(w)\fR).
29
\fVsys=mrfit(w,mag,order,weight)\fR returns a bi-stable transfer function
30
\fVG(s)=sys=num/den\fR, of of given \fVorder\fR such that
31
its frequency response magnitude \fVabs(G(w(i)))\fR
32
matches \fVmag(i)\fR i.e. \fVabs(freq(num,den,%i*w))\fR should be
34
\fVweight(i)\fR is the weigth given to \fVw(i)\fR.
37
w=0.01:0.01:2;s=poly(0,'s');
38
G=syslin('c',2*(s^2+0.1*s+2),(s^2+s+1)*(s^2+0.3*s+1)); // syslin('c',Num,Den);
42
frespfit=repfreq(Gid,w);
43
plot2d([w',w'],[mag(:),abs(frespfit(:))])
46
cepstrum, frfit, freq, calfrq