1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
|
<html><head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<title>invfreqs</title>
<style type="text/css" media="all">
@import url("scilab_code.css");
@import url("xml_code.css");
@import url("c_code.css");
@import url("style.css");
</style>
</head>
<body>
<div class="manualnavbar">
<table width="100%"><tr>
<td width="30%">
<span class="previous"><a href="invfreq.html"><< invfreq</a></span>
</td>
<td width="40%" class="center">
<span class="top"><a href="section_e54aa8aac34aa55341e8b4b782fe1a74.html">FOSSEE Signal Processing Toolbox</a></span>
</td>
<td width="30%" class="next">
<span class="next"><a href="invfreqz.html">invfreqz >></a></span>
</td>
</tr></table>
<hr />
</div>
<span class="path"><a href="index.html">FOSSEE Signal Processing Toolbox</a> >> <a href="section_e54aa8aac34aa55341e8b4b782fe1a74.html">FOSSEE Signal Processing Toolbox</a> > invfreqs</span>
<br /><br />
<div class="refnamediv"><h1 class="refname">invfreqs</h1>
<p class="refpurpose">Fit filter B(s)/A(s)to the complex frequency response H at frequency points F. A and B are real polynomial coefficients of order nA and nB.</p></div>
<div class="refsynopsisdiv"><h3 class="title">Calling Sequence</h3>
<div class="synopsis"><pre><span class="default">[</span><span class="default">B</span><span class="default">,</span><span class="default">A</span><span class="default">,</span><span class="default">C</span><span class="default">] = </span><span class="functionid">invfreqs</span><span class="default">(</span><span class="default">H</span><span class="default">,</span><span class="default">F</span><span class="default">,</span><span class="default">nB</span><span class="default">,</span><span class="default">nA</span><span class="default">,</span><span class="default">W</span><span class="default">,</span><span class="default">iter</span><span class="default">,</span><span class="default">tol</span><span class="default">,</span><span class="default">trace</span><span class="default">)</span>
<span class="default">[</span><span class="default">B</span><span class="default">,</span><span class="default">A</span><span class="default">,</span><span class="default">C</span><span class="default">] = </span><span class="functionid">invfreqs</span><span class="default">(</span><span class="default">H</span><span class="default">,</span><span class="default">F</span><span class="default">,</span><span class="default">nB</span><span class="default">,</span><span class="default">nA</span><span class="default">,</span><span class="default">W</span><span class="default">)</span>
<span class="default">[</span><span class="default">B</span><span class="default">,</span><span class="default">A</span><span class="default">,</span><span class="default">C</span><span class="default">] = </span><span class="functionid">invfreqs</span><span class="default">(</span><span class="default">H</span><span class="default">,</span><span class="default">F</span><span class="default">,</span><span class="default">nB</span><span class="default">,</span><span class="default">nA</span><span class="default">)</span></pre></div></div>
<div class="refsection"><h3 class="title">Parameters</h3>
<dl><dt><span class="term">H:</span>
<dd><p class="para">desired complex frequency response.</p></dd></dt>
<dt><span class="term">F:</span>
<dd><p class="para">frequency (must be same length as H).</p></dd></dt>
<dt><span class="term">nB:</span>
<dd><p class="para">order of the numerator polynomial B.</p></dd></dt>
<dt><span class="term">nA:</span>
<dd><p class="para">order of the denominator polynomial A.</p></dd></dt>
<dt><span class="term">W:</span>
<dd><p class="para">vector of weights (must be same length as F).</p></dd></dt></dl></div>
<div class="refsection"><h3 class="title">Description</h3>
<p class="para">This is an Octave function.
Fit filter B(s)/A(s)to the complex frequency response H at frequency points F. A and B are real polynomial coefficients of order nA and nB.
Optionally, the fit-errors can be weighted vs frequency according to the weights W.
Note: all the guts are in invfreq.m</p></div>
<div class="refsection"><h3 class="title">Examples</h3>
<div class="programlisting"><table border="0" width="100%"><tr><td width="98%"><pre class="scilabcode"><span class="scilabid">B</span> <span class="scilaboperator">=</span> <span class="scilabopenclose">[</span><span class="scilabnumber">1</span><span class="scilaboperator">/</span><span class="scilabnumber">2</span> <span class="scilabnumber">1</span><span class="scilabopenclose">]</span><span class="scilabdefault">;</span>
<span class="scilabid">A</span> <span class="scilaboperator">=</span> <span class="scilabopenclose">[</span><span class="scilabnumber">1</span> <span class="scilabnumber">1</span><span class="scilabopenclose">]</span><span class="scilabdefault">;</span>
<span class="scilabid">w</span> <span class="scilaboperator">=</span> <a class="scilabmacro" href="scilab://linspace">linspace</a><span class="scilabopenclose">(</span><span class="scilabnumber">0</span><span class="scilabdefault">,</span><span class="scilabnumber">4</span><span class="scilabdefault">,</span><span class="scilabnumber">128</span><span class="scilabopenclose">)</span><span class="scilabdefault">;</span>
<span class="scilabid">H</span> <span class="scilaboperator">=</span> <span class="scilabid">freqs</span><span class="scilabopenclose">(</span><span class="scilabid">B</span><span class="scilabdefault">,</span><span class="scilabid">A</span><span class="scilabdefault">,</span><span class="scilabid">w</span><span class="scilabopenclose">)</span><span class="scilabdefault">;</span>
<span class="scilabopenclose">[</span><span class="scilabid">Bh</span><span class="scilabdefault">,</span><span class="scilabid">Ah</span><span class="scilabdefault">,</span> <span class="scilabid">C</span><span class="scilabopenclose">]</span> <span class="scilaboperator">=</span> <span class="scilabid">invfreqs</span><span class="scilabopenclose">(</span><span class="scilabid">H</span><span class="scilabdefault">,</span><span class="scilabid">w</span><span class="scilabdefault">,</span><span class="scilabnumber">1</span><span class="scilabdefault">,</span><span class="scilabnumber">1</span><span class="scilabopenclose">)</span><span class="scilabdefault">;</span>
<span class="scilabid">Bh</span> <span class="scilaboperator">=</span>
<span class="scilabnumber">0.50000</span> <span class="scilabnumber">1.00000</span>
<span class="scilabid">Ah</span> <span class="scilaboperator">=</span>
<span class="scilabnumber">1.0000</span> <span class="scilabnumber">1.0000</span>
<span class="scilabid">C</span> <span class="scilaboperator">=</span> <span class="scilaboperator">-</span><span class="scilabnumber">3.0964e-15</span></pre></td><td valign="top"><a href="scilab://scilab.execexample/"><img src="ScilabExecute.png" border="0"/></a></td><td valign="top"><a href="scilab://scilab.editexample/"><img src="ScilabEdit.png" border="0"/></a></td><td></td></tr></table></div></div>
<br />
<div class="manualnavbar">
<table width="100%">
<tr><td colspan="3" class="next"><a href="http://bugzilla.scilab.org/enter_bug.cgi?product=Scilab%20software&component=Documentation%20pages" class="ulink">Report an issue</a></td></tr>
<tr>
<td width="30%">
<span class="previous"><a href="invfreq.html"><< invfreq</a></span>
</td>
<td width="40%" class="center">
<span class="top"><a href="section_e54aa8aac34aa55341e8b4b782fe1a74.html">FOSSEE Signal Processing Toolbox</a></span>
</td>
<td width="30%" class="next">
<span class="next"><a href="invfreqz.html">invfreqz >></a></span>
</td>
</tr></table>
<hr />
</div>
</body>
</html>
|
