kopia lustrzana https://github.com/Dsplib/libdspl-2.0
110 wiersze
40 KiB
HTML
110 wiersze
40 KiB
HTML
<!-- HTML header for doxygen 1.8.13
|
|
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
|
<html xmlns="http://www.w3.org/1999/xhtml">
|
|
<head>
|
|
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
|
|
-->
|
|
<html>
|
|
<head>
|
|
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
|
|
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
|
|
<meta name="generator" content="Doxygen 1.8.15"/>
|
|
<meta name="viewport" content="width=device-width, initial-scale=1"/>
|
|
<title>libdspl-2.0: Исходный файл F:/dsplib.org/libdspl-2.0/dspl/src/filter_iir.c</title>
|
|
<link href="tabs.css" rel="stylesheet" type="text/css"/>
|
|
<script type="text/javascript" src="jquery.js"></script>
|
|
<script type="text/javascript" src="dynsections.js"></script>
|
|
<link rel="stylesheet" href="https://maxcdn.bootstrapcdn.com/font-awesome/4.4.0/css/font-awesome.min.css">
|
|
<meta name="viewport" content="width=device-width, initial-scale=1.0">
|
|
<link href="doxy_stylesheet.css" rel="stylesheet" type="text/css" />
|
|
<link href="http://dsplib.org/dsplib-stylesheet.css".css" rel="stylesheet" type="text/css">
|
|
</head>
|
|
<body>
|
|
<!-- Yandex.Metrika counter -->
|
|
<script type="text/javascript" >
|
|
(function (d, w, c) {
|
|
(w[c] = w[c] || []).push(function() {
|
|
try {
|
|
w.yaCounter32971694 = new Ya.Metrika({
|
|
id:32971694,
|
|
clickmap:true,
|
|
trackLinks:true,
|
|
accurateTrackBounce:true
|
|
});
|
|
} catch(e) { }
|
|
});
|
|
var n = d.getElementsByTagName("script")[0],
|
|
s = d.createElement("script"),
|
|
f = function () { n.parentNode.insertBefore(s, n); };
|
|
s.type = "text/javascript";
|
|
s.async = true;
|
|
s.src = "https://mc.yandex.ru/metrika/watch.js";
|
|
if (w.opera == "[object Opera]") {
|
|
d.addEventListener("DOMContentLoaded", f, false);
|
|
} else { f(); }
|
|
})(document, window, "yandex_metrika_callbacks");
|
|
</script>
|
|
<noscript><div><img src="https://mc.yandex.ru/watch/32971694" style="position:absolute; left:-9999px;" alt="" /></div></noscript>
|
|
<!-- /Yandex.Metrika counter -->
|
|
<div class = "header-bar">
|
|
<div class = "menu-bar">
|
|
<nav>
|
|
<ul>
|
|
<div class = "dsplib-logo"><img src="http://dsplib.org/logo.png" /> </div>
|
|
<li><a href="http://ru.dsplib.org">Содержание</a></li>
|
|
<li><a href="http://ru.dsplib.org/dspl">DSPL</a></li>
|
|
<li><a href="http://ru.dsplib.org/forum">Форум</a></li>
|
|
<li>
|
|
<div class="search-bar">
|
|
<div class="ya-site-form ya-site-form_inited_no" onclick="return {'action':'http://ru.dsplib.org/search_results.html','arrow':false,'bg':'transparent','fontsize':12,'fg':'#000000','language':'ru','logo':'rb','publicname':'ru.dsplib.org поиск','suggest':true,'target':'_self','tld':'ru','type':2,'usebigdictionary':true,'searchid':2332185,'input_fg':'#000000','input_bg':'#ffffff','input_fontStyle':'normal','input_fontWeight':'normal','input_placeholder':'поиск','input_placeholderColor':'#000000','input_borderColor':'#7f9db9'}"><form action="https://yandex.ru/search/site/" method="get" target="_self" accept-charset="utf-8"><input type="hidden" name="searchid" value="2332185"/><input type="hidden" name="l10n" value="ru"/><input type="hidden" name="reqenc" value=""/><input type="search" name="text" value="" style = "height: 24px; font-family: verdana,arial; font-size: 12px"/><input type="submit" value="Найти" style = "display: none;"/></form></div><style type="text/css">.ya-page_js_yes .ya-site-form_inited_no { display: none; }</style><script type="text/javascript">(function(w,d,c){var s=d.createElement('script'),h=d.getElementsByTagName('script')[0],e=d.documentElement;if((' '+e.className+' ').indexOf(' ya-page_js_yes ')===-1){e.className+=' ya-page_js_yes';}s.type='text/javascript';s.async=true;s.charset='utf-8';s.src=(d.location.protocol==='https:'?'https:':'http:')+'//site.yandex.net/v2.0/js/all.js';h.parentNode.insertBefore(s,h);(w[c]||(w[c]=[])).push(function(){Ya.Site.Form.init()})})(window,document,'yandex_site_callbacks');</script>
|
|
</div>
|
|
</li>
|
|
</ul>
|
|
</nav>
|
|
</div>
|
|
</div>
|
|
<div class="main-dsplib">
|
|
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
|
|
<div id="titlearea">
|
|
<table cellspacing="0" cellpadding="0">
|
|
<tbody>
|
|
<tr style="height: 56px;">
|
|
<td id="projectalign" style="padding-left: 0.5em;">
|
|
<div id="projectname">libdspl-2.0
|
|
</div>
|
|
<div id="projectbrief">Библиотека алгоритмов цифровой обработки сигналов</div>
|
|
</td>
|
|
</tr>
|
|
</tbody>
|
|
</table>
|
|
</div>
|
|
<!-- end header part -->
|
|
<!-- Создано системой Doxygen 1.8.15 -->
|
|
<div id="nav-path" class="navpath">
|
|
<ul>
|
|
<li class="navelem"><a class="el" href="dir_12d635cfd40a13a86dc83e7a5cc309f3.html">src</a></li> </ul>
|
|
</div>
|
|
</div><!-- top -->
|
|
<div class="header">
|
|
<div class="headertitle">
|
|
<div class="title">filter_iir.c</div> </div>
|
|
</div><!--header-->
|
|
<div class="contents">
|
|
<div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">/*</span></div><div class="line"><a name="l00002"></a><span class="lineno"> 2</span> <span class="comment">* Copyright (c) 2015-2019 Sergey Bakhurin</span></div><div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">* Digital Signal Processing Library [http://dsplib.org]</span></div><div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment">*</span></div><div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment">* This file is part of libdspl-2.0.</span></div><div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment">*</span></div><div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment">* is free software: you can redistribute it and/or modify</span></div><div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment">* it under the terms of the GNU Lesser General Public License as published by</span></div><div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment">* the Free Software Foundation, either version 3 of the License, or</span></div><div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment">* (at your option) any later version.</span></div><div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment">*</span></div><div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment">* DSPL is distributed in the hope that it will be useful,</span></div><div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment">* but WITHOUT ANY WARRANTY; without even the implied warranty of</span></div><div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment">* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span></div><div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment">* GNU General Public License for more details.</span></div><div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment">*</span></div><div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment">* You should have received a copy of the GNU Lesser General Public License</span></div><div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment">* along with Foobar. If not, see <http://www.gnu.org/licenses/>.</span></div><div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment">*/</span></div><div class="line"><a name="l00020"></a><span class="lineno"> 20</span> </div><div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="preprocessor">#include <stdlib.h></span></div><div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="preprocessor">#include <stdio.h></span></div><div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="preprocessor">#include <string.h></span></div><div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="preprocessor">#include "dspl.h"</span></div><div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="preprocessor">#include "dspl_internal.h"</span></div><div class="line"><a name="l00026"></a><span class="lineno"> 26</span> </div><div class="line"><a name="l00027"></a><span class="lineno"> 27</span> </div><div class="line"><a name="l00028"></a><span class="lineno"> 28</span> </div><div class="line"><a name="l00029"></a><span class="lineno"> 29</span> </div><div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="comment">/*******************************************************************************</span></div><div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="comment">\ingroup IIR_FILTER_DESIGN_GROUP</span></div><div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="comment">\fn int bilinear(double* bs, double* as, int ord, double* bz, double* az)</span></div><div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="comment">\brief Analog filter transfer function H(s) bilinear transform to the</span></div><div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment"> digital filter transfer function H(z).</span></div><div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="comment"> </span></div><div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment">Function calculates digital filter coefficients by rational substitution</span></div><div class="line"><a name="l00037"></a><span class="lineno"> 37</span> <span class="comment"></span></div><div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment">\f[</span></div><div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment">s \leftarrow \frac{1 - z^{-1}}{1 - z^{-1}}.</span></div><div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="comment">\f]</span></div><div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="comment"></span></div><div class="line"><a name="l00042"></a><span class="lineno"> 42</span> <span class="comment">Digital filter order still the same as analog prototype order. </span></div><div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="comment">Analog prototype frequency \f$\Omega\f$ related with the digital filter </span></div><div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="comment">normalized frequency \f$\omega\f$ as:</span></div><div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="comment"></span></div><div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="comment">\f[</span></div><div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="comment">\Omega = \tan(\omega / 2).</span></div><div class="line"><a name="l00048"></a><span class="lineno"> 48</span> <span class="comment">\f]</span></div><div class="line"><a name="l00049"></a><span class="lineno"> 49</span> <span class="comment"></span></div><div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="comment">\param[in] bs Pointer to the numerator coefficients of an </span></div><div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="comment"> analog prototype transfer function \f$H(s)\f$. \n</span></div><div class="line"><a name="l00052"></a><span class="lineno"> 52</span> <span class="comment"> Array size is `[ord+1 x 1]`. \n</span></div><div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="comment"> Memory must be allocated. \n \n</span></div><div class="line"><a name="l00054"></a><span class="lineno"> 54</span> <span class="comment"></span></div><div class="line"><a name="l00055"></a><span class="lineno"> 55</span> <span class="comment">\param[in] as Pointer to the denominator coefficients of an </span></div><div class="line"><a name="l00056"></a><span class="lineno"> 56</span> <span class="comment"> analog prototype transfer function \f$H(s)\f$. \n</span></div><div class="line"><a name="l00057"></a><span class="lineno"> 57</span> <span class="comment"> Array size is `[ord+1 x 1]`. \n</span></div><div class="line"><a name="l00058"></a><span class="lineno"> 58</span> <span class="comment"> Memory must be allocated. \n \n</span></div><div class="line"><a name="l00059"></a><span class="lineno"> 59</span> <span class="comment"> </span></div><div class="line"><a name="l00060"></a><span class="lineno"> 60</span> <span class="comment">\param[in] ord Filter order. \n</span></div><div class="line"><a name="l00061"></a><span class="lineno"> 61</span> <span class="comment"> Number of coefficients \f$H(s)\f$ and \f$H(z)\f$</span></div><div class="line"><a name="l00062"></a><span class="lineno"> 62</span> <span class="comment"> numerator and denominator equals `ord+1`. \n \n</span></div><div class="line"><a name="l00063"></a><span class="lineno"> 63</span> <span class="comment"> </span></div><div class="line"><a name="l00064"></a><span class="lineno"> 64</span> <span class="comment">\param[out] bz Pointer to the numerator coefficients of a</span></div><div class="line"><a name="l00065"></a><span class="lineno"> 65</span> <span class="comment"> digital filter transfer function \f$H(z)\f$. \n</span></div><div class="line"><a name="l00066"></a><span class="lineno"> 66</span> <span class="comment"> Array size is `[ord+1 x 1]`. \n</span></div><div class="line"><a name="l00067"></a><span class="lineno"> 67</span> <span class="comment"> Memory must be allocated. \n \n</span></div><div class="line"><a name="l00068"></a><span class="lineno"> 68</span> <span class="comment"></span></div><div class="line"><a name="l00069"></a><span class="lineno"> 69</span> <span class="comment">\param[out] az Pointer to the numerator coefficients of a</span></div><div class="line"><a name="l00070"></a><span class="lineno"> 70</span> <span class="comment"> digital filter transfer function \f$H(z)\f$. \n</span></div><div class="line"><a name="l00071"></a><span class="lineno"> 71</span> <span class="comment"> Array size is `[ord+1 x 1]`. \n</span></div><div class="line"><a name="l00072"></a><span class="lineno"> 72</span> <span class="comment"> Memory must be allocated. \n \n</span></div><div class="line"><a name="l00073"></a><span class="lineno"> 73</span> <span class="comment"></span></div><div class="line"><a name="l00074"></a><span class="lineno"> 74</span> <span class="comment">\return</span></div><div class="line"><a name="l00075"></a><span class="lineno"> 75</span> <span class="comment"> `RES_OK` if filter is calculated successfully. \n \n</span></div><div class="line"><a name="l00076"></a><span class="lineno"> 76</span> <span class="comment"> Else \ref ERROR_CODE_GROUP "code error". \n</span></div><div class="line"><a name="l00077"></a><span class="lineno"> 77</span> <span class="comment"></span></div><div class="line"><a name="l00078"></a><span class="lineno"> 78</span> <span class="comment">\author Sergey Bakhurin www.dsplib.org </span></div><div class="line"><a name="l00079"></a><span class="lineno"> 79</span> <span class="comment">*******************************************************************************/</span></div><div class="line"><a name="l00080"></a><span class="lineno"><a class="line" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga1730d572761655ef7104b1b2cd1f9c11"> 80</a></span> <span class="keywordtype">int</span> DSPL_API <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga1730d572761655ef7104b1b2cd1f9c11">bilinear</a>(<span class="keywordtype">double</span>* bs, <span class="keywordtype">double</span>* as, <span class="keywordtype">int</span> ord, <span class="keywordtype">double</span>* bz, <span class="keywordtype">double</span>* az)</div><div class="line"><a name="l00081"></a><span class="lineno"> 81</span> {</div><div class="line"><a name="l00082"></a><span class="lineno"> 82</span>  <span class="keywordtype">double</span> c[2] = {1.0, -1.0};</div><div class="line"><a name="l00083"></a><span class="lineno"> 83</span>  <span class="keywordtype">double</span> d[2] = {1.0, 1.0};</div><div class="line"><a name="l00084"></a><span class="lineno"> 84</span>  <span class="keywordflow">return</span> <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#gadab58ef5d80c29882e663ffbe10f9a80">ratcompos</a>(bs, as, ord, c, d, 1, bz, az);</div><div class="line"><a name="l00085"></a><span class="lineno"> 85</span> }</div><div class="line"><a name="l00086"></a><span class="lineno"> 86</span> </div><div class="line"><a name="l00087"></a><span class="lineno"> 87</span> </div><div class="line"><a name="l00088"></a><span class="lineno"> 88</span> </div><div class="line"><a name="l00089"></a><span class="lineno"> 89</span> </div><div class="line"><a name="l00090"></a><span class="lineno"> 90</span> <span class="comment">/******************************************************************************</span></div><div class="line"><a name="l00091"></a><span class="lineno"> 91</span> <span class="comment">Digital IIR filter coefficients calculation</span></div><div class="line"><a name="l00092"></a><span class="lineno"> 92</span> <span class="comment">*******************************************************************************/</span></div><div class="line"><a name="l00093"></a><span class="lineno"> 93</span> <span class="keywordtype">int</span> DSPL_API iir(<span class="keywordtype">double</span> rp, <span class="keywordtype">double</span> rs, <span class="keywordtype">int</span> ord, <span class="keywordtype">double</span> w0, <span class="keywordtype">double</span> w1,</div><div class="line"><a name="l00094"></a><span class="lineno"> 94</span>  <span class="keywordtype">int</span> type, <span class="keywordtype">double</span>* b, <span class="keywordtype">double</span>* a)</div><div class="line"><a name="l00095"></a><span class="lineno"> 95</span> {</div><div class="line"><a name="l00096"></a><span class="lineno"> 96</span>  <span class="keywordtype">double</span> *bs = NULL;</div><div class="line"><a name="l00097"></a><span class="lineno"> 97</span>  <span class="keywordtype">double</span> *as = NULL;</div><div class="line"><a name="l00098"></a><span class="lineno"> 98</span>  <span class="keywordtype">double</span> *bt = NULL;</div><div class="line"><a name="l00099"></a><span class="lineno"> 99</span>  <span class="keywordtype">double</span> *at = NULL;</div><div class="line"><a name="l00100"></a><span class="lineno"> 100</span>  <span class="keywordtype">double</span> wa0, wa1, ws;</div><div class="line"><a name="l00101"></a><span class="lineno"> 101</span>  <span class="keywordtype">int</span> err, ord_ap = ord;</div><div class="line"><a name="l00102"></a><span class="lineno"> 102</span>  </div><div class="line"><a name="l00103"></a><span class="lineno"> 103</span>  <span class="keywordflow">if</span>(((type & DSPL_FILTER_TYPE_MASK) == DSPL_FILTER_LPF) ||</div><div class="line"><a name="l00104"></a><span class="lineno"> 104</span>  ((type & DSPL_FILTER_TYPE_MASK) == DSPL_FILTER_HPF))</div><div class="line"><a name="l00105"></a><span class="lineno"> 105</span>  {</div><div class="line"><a name="l00106"></a><span class="lineno"> 106</span>  bs = (<span class="keywordtype">double</span>*)malloc((ord_ap+1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00107"></a><span class="lineno"> 107</span>  as = (<span class="keywordtype">double</span>*)malloc((ord_ap+1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00108"></a><span class="lineno"> 108</span>  bt = (<span class="keywordtype">double</span>*)malloc((ord_ap+1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00109"></a><span class="lineno"> 109</span>  at = (<span class="keywordtype">double</span>*)malloc((ord_ap+1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00110"></a><span class="lineno"> 110</span>  }</div><div class="line"><a name="l00111"></a><span class="lineno"> 111</span>  </div><div class="line"><a name="l00112"></a><span class="lineno"> 112</span>  </div><div class="line"><a name="l00113"></a><span class="lineno"> 113</span>  <span class="keywordflow">if</span>(((type & DSPL_FILTER_TYPE_MASK) == DSPL_FILTER_BPASS) ||</div><div class="line"><a name="l00114"></a><span class="lineno"> 114</span>  ((type & DSPL_FILTER_TYPE_MASK) == DSPL_FILTER_BSTOP))</div><div class="line"><a name="l00115"></a><span class="lineno"> 115</span>  {</div><div class="line"><a name="l00116"></a><span class="lineno"> 116</span>  <span class="keywordflow">if</span>(ord % 2)</div><div class="line"><a name="l00117"></a><span class="lineno"> 117</span>  <span class="keywordflow">return</span> ERROR_FILTER_ORD_BP;</div><div class="line"><a name="l00118"></a><span class="lineno"> 118</span>  <span class="keywordflow">else</span></div><div class="line"><a name="l00119"></a><span class="lineno"> 119</span>  {</div><div class="line"><a name="l00120"></a><span class="lineno"> 120</span>  ord_ap = ord / 2;</div><div class="line"><a name="l00121"></a><span class="lineno"> 121</span>  bs = (<span class="keywordtype">double</span>*)malloc((ord_ap + 1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00122"></a><span class="lineno"> 122</span>  as = (<span class="keywordtype">double</span>*)malloc((ord_ap + 1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00123"></a><span class="lineno"> 123</span>  bt = (<span class="keywordtype">double</span>*)malloc((ord + 1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00124"></a><span class="lineno"> 124</span>  at = (<span class="keywordtype">double</span>*)malloc((ord + 1)*<span class="keyword">sizeof</span>(double));</div><div class="line"><a name="l00125"></a><span class="lineno"> 125</span>  }</div><div class="line"><a name="l00126"></a><span class="lineno"> 126</span>  }</div><div class="line"><a name="l00127"></a><span class="lineno"> 127</span>  err = iir_ap(rp, rs, ord_ap, type, bs, as);</div><div class="line"><a name="l00128"></a><span class="lineno"> 128</span>  <span class="keywordflow">if</span>(err != <a class="code" href="group___e_r_r_o_r___c_o_d_e___g_r_o_u_p.html#ga312c04ce416e0afdbd653394e36a117d">RES_OK</a>)</div><div class="line"><a name="l00129"></a><span class="lineno"> 129</span>  <span class="keywordflow">goto</span> error_proc;</div><div class="line"><a name="l00130"></a><span class="lineno"> 130</span>  </div><div class="line"><a name="l00131"></a><span class="lineno"> 131</span>  <span class="comment">/* frequency transformation */</span></div><div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  wa0 = tan(w0 * M_PI * 0.5);</div><div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  wa1 = tan(w1 * M_PI * 0.5);</div><div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  </div><div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  <span class="keywordflow">switch</span>(type & DSPL_FILTER_TYPE_MASK)</div><div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  {</div><div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  </div><div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  <span class="keywordflow">case</span> DSPL_FILTER_LPF:</div><div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#gacac2324392316d6e27c05bc0f839e81c">low2low</a>(bs, as, ord_ap, 1.0, wa0, bt, at);</div><div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  </div><div class="line"><a name="l00142"></a><span class="lineno"> 142</span>  <span class="keywordflow">case</span> DSPL_FILTER_HPF:</div><div class="line"><a name="l00143"></a><span class="lineno"> 143</span>  ws = filter_ws1(ord_ap, rp, rs, type);</div><div class="line"><a name="l00144"></a><span class="lineno"> 144</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#gacac2324392316d6e27c05bc0f839e81c">low2low</a>( bs, as, ord_ap, 1.0, 1.0 / ws, bs, as);</div><div class="line"><a name="l00145"></a><span class="lineno"> 145</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga33bbbb274dc08cda15273963ae919187">low2high</a>(bs, as, ord_ap, 1.0, wa0, bt, at);</div><div class="line"><a name="l00146"></a><span class="lineno"> 146</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00147"></a><span class="lineno"> 147</span>  </div><div class="line"><a name="l00148"></a><span class="lineno"> 148</span>  <span class="keywordflow">case</span> DSPL_FILTER_BPASS:</div><div class="line"><a name="l00149"></a><span class="lineno"> 149</span>  err = low2bp(bs, as, ord_ap, 1.0, wa0, wa1, bt, at);</div><div class="line"><a name="l00150"></a><span class="lineno"> 150</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00151"></a><span class="lineno"> 151</span>  </div><div class="line"><a name="l00152"></a><span class="lineno"> 152</span>  <span class="keywordflow">case</span> DSPL_FILTER_BSTOP:</div><div class="line"><a name="l00153"></a><span class="lineno"> 153</span>  <span class="comment">/* need frequency transform ws -> 1 rad/s */</span></div><div class="line"><a name="l00154"></a><span class="lineno"> 154</span>  </div><div class="line"><a name="l00155"></a><span class="lineno"> 155</span>  ws = filter_ws1(ord_ap, rp, rs, type);</div><div class="line"><a name="l00156"></a><span class="lineno"> 156</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#gacac2324392316d6e27c05bc0f839e81c">low2low</a>( bs, as, ord_ap, 1.0, 1.0 / ws, bs, as);</div><div class="line"><a name="l00157"></a><span class="lineno"> 157</span>  err = low2bs(bs, as, ord_ap, 1.0, wa0, wa1, bt, at);</div><div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  </div><div class="line"><a name="l00160"></a><span class="lineno"> 160</span>  <span class="keywordflow">default</span>:</div><div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  err = ERROR_FILTER_TYPE;</div><div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  }</div><div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  <span class="keywordflow">if</span>(err != <a class="code" href="group___e_r_r_o_r___c_o_d_e___g_r_o_u_p.html#ga312c04ce416e0afdbd653394e36a117d">RES_OK</a>)</div><div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="keywordflow">goto</span> error_proc;</div><div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  </div><div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  </div><div class="line"><a name="l00168"></a><span class="lineno"> 168</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga1730d572761655ef7104b1b2cd1f9c11">bilinear</a>(bt, at, ord, b, a);</div><div class="line"><a name="l00169"></a><span class="lineno"> 169</span> </div><div class="line"><a name="l00170"></a><span class="lineno"> 170</span> </div><div class="line"><a name="l00171"></a><span class="lineno"> 171</span> error_proc:</div><div class="line"><a name="l00172"></a><span class="lineno"> 172</span> </div><div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  <span class="keywordflow">if</span>(bs)</div><div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  free(bs);</div><div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  <span class="keywordflow">if</span>(as)</div><div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  free(as);</div><div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  <span class="keywordflow">if</span>(bt)</div><div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  free(bt);</div><div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  <span class="keywordflow">if</span>(at)</div><div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  free(at);</div><div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  </div><div class="line"><a name="l00182"></a><span class="lineno"> 182</span>  <span class="keywordflow">return</span> err;</div><div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  </div><div class="line"><a name="l00184"></a><span class="lineno"> 184</span> }</div><div class="line"><a name="l00185"></a><span class="lineno"> 185</span> </div><div class="line"><a name="l00186"></a><span class="lineno"> 186</span> </div><div class="line"><a name="l00187"></a><span class="lineno"> 187</span> </div><div class="line"><a name="l00188"></a><span class="lineno"> 188</span> </div><div class="line"><a name="l00189"></a><span class="lineno"> 189</span> </div><div class="line"><a name="l00190"></a><span class="lineno"> 190</span> </div><div class="line"><a name="l00191"></a><span class="lineno"> 191</span> <span class="comment">/******************************************************************************</span></div><div class="line"><a name="l00192"></a><span class="lineno"> 192</span> <span class="comment">Analog prototype for IIR </span></div><div class="line"><a name="l00193"></a><span class="lineno"> 193</span> <span class="comment">*******************************************************************************/</span></div><div class="line"><a name="l00194"></a><span class="lineno"> 194</span> <span class="keywordtype">int</span> iir_ap(<span class="keywordtype">double</span> rp, <span class="keywordtype">double</span> rs, <span class="keywordtype">int</span> ord, <span class="keywordtype">int</span> type, <span class="keywordtype">double</span>* b, <span class="keywordtype">double</span>* a)</div><div class="line"><a name="l00195"></a><span class="lineno"> 195</span> {</div><div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  <span class="keywordtype">int</span> err;</div><div class="line"><a name="l00197"></a><span class="lineno"> 197</span>  <span class="keywordflow">switch</span>(type & DSPL_FILTER_APPROX_MASK)</div><div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  {</div><div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  <span class="keywordflow">case</span> DSPL_FILTER_BUTTER:</div><div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga65df99359f211937e4dc1141f57ef571">butter_ap</a>(rp, ord, b, a);</div><div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  <span class="keywordflow">case</span> DSPL_FILTER_CHEBY1:</div><div class="line"><a name="l00203"></a><span class="lineno"> 203</span>  err = <a class="code" href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga6215e6c0a94c617370c94855cc31c2fe">cheby1_ap</a>(rp, ord, b, a);</div><div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  <span class="keywordflow">case</span> DSPL_FILTER_CHEBY2:</div><div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  err = cheby2_ap_wp1(rp, rs, ord, b, a);</div><div class="line"><a name="l00207"></a><span class="lineno"> 207</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00208"></a><span class="lineno"> 208</span>  <span class="keywordflow">case</span> DSPL_FILTER_ELLIP:</div><div class="line"><a name="l00209"></a><span class="lineno"> 209</span>  err = ellip_ap(rp, rs, ord, b, a);</div><div class="line"><a name="l00210"></a><span class="lineno"> 210</span>  <span class="keywordflow">break</span>;</div><div class="line"><a name="l00211"></a><span class="lineno"> 211</span>  <span class="keywordflow">default</span>:</div><div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  err = ERROR_FILTER_APPROX; </div><div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  }</div><div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  </div><div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  </div><div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  <span class="keywordflow">return</span> err;</div><div class="line"><a name="l00217"></a><span class="lineno"> 217</span> }</div><div class="line"><a name="l00218"></a><span class="lineno"> 218</span> </div><div class="line"><a name="l00219"></a><span class="lineno"> 219</span> </div><div class="line"><a name="l00220"></a><span class="lineno"> 220</span> </div><div class="line"><a name="l00221"></a><span class="lineno"> 221</span> </div><div class="ttc" id="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p_html_gacac2324392316d6e27c05bc0f839e81c"><div class="ttname"><a href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#gacac2324392316d6e27c05bc0f839e81c">low2low</a></div><div class="ttdeci">int low2low(double *b, double *a, int ord, double w0, double w1, double *beta, double *alpha)</div><div class="ttdoc">Частотное преобразование ФНЧ-ФНЧ</div><div class="ttdef"><b>Definition:</b> <a href="filter__ft_8c_source.html#l00179">filter_ft.c:179</a></div></div>
|
|
<div class="ttc" id="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p_html_ga65df99359f211937e4dc1141f57ef571"><div class="ttname"><a href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga65df99359f211937e4dc1141f57ef571">butter_ap</a></div><div class="ttdeci">int butter_ap(double rp, int ord, double *b, double *a)</div><div class="ttdoc">Расчет передаточной характеристики аналогового нормированного ФНЧ Баттерворта.</div><div class="ttdef"><b>Definition:</b> <a href="filter__ap_8c_source.html#l00033">filter_ap.c:33</a></div></div>
|
|
<div class="ttc" id="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p_html_gadab58ef5d80c29882e663ffbe10f9a80"><div class="ttname"><a href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#gadab58ef5d80c29882e663ffbe10f9a80">ratcompos</a></div><div class="ttdeci">int ratcompos(double *b, double *a, int n, double *c, double *d, int p, double *beta, double *alpha)</div><div class="ttdoc">Рациональная композиця</div><div class="ttdef"><b>Definition:</b> <a href="filter__ft_8c_source.html#l00204">filter_ft.c:204</a></div></div>
|
|
<div class="ttc" id="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p_html_ga33bbbb274dc08cda15273963ae919187"><div class="ttname"><a href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga33bbbb274dc08cda15273963ae919187">low2high</a></div><div class="ttdeci">int low2high(double *b, double *a, int ord, double w0, double w1, double *beta, double *alpha)</div><div class="ttdoc">Частотное преобразование ФНЧ-ФВЧ</div><div class="ttdef"><b>Definition:</b> <a href="filter__ft_8c_source.html#l00151">filter_ft.c:151</a></div></div>
|
|
<div class="ttc" id="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p_html_ga1730d572761655ef7104b1b2cd1f9c11"><div class="ttname"><a href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga1730d572761655ef7104b1b2cd1f9c11">bilinear</a></div><div class="ttdeci">int bilinear(double *bs, double *as, int ord, double *bz, double *az)</div><div class="ttdoc">Билинейное преобразование передаточной характеристики аналогового фильтра , в передаточную характерис...</div><div class="ttdef"><b>Definition:</b> <a href="filter__iir_8c_source.html#l00080">filter_iir.c:80</a></div></div>
|
|
<div class="ttc" id="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p_html_ga6215e6c0a94c617370c94855cc31c2fe"><div class="ttname"><a href="group___i_i_r___f_i_l_t_e_r___d_e_s_i_g_n___g_r_o_u_p.html#ga6215e6c0a94c617370c94855cc31c2fe">cheby1_ap</a></div><div class="ttdeci">int cheby1_ap(double rp, int ord, double *b, double *a)</div><div class="ttdoc">Расчет передаточной характеристики аналогового нормированного ФНЧ Чебышёва первого рода.</div><div class="ttdef"><b>Definition:</b> <a href="filter__ap_8c_source.html#l00129">filter_ap.c:129</a></div></div>
|
|
<div class="ttc" id="group___e_r_r_o_r___c_o_d_e___g_r_o_u_p_html_ga312c04ce416e0afdbd653394e36a117d"><div class="ttname"><a href="group___e_r_r_o_r___c_o_d_e___g_r_o_u_p.html#ga312c04ce416e0afdbd653394e36a117d">RES_OK</a></div><div class="ttdeci">#define RES_OK</div><div class="ttdoc">Функция завершилась корректно. Ошибки отсутствуют.</div><div class="ttdef"><b>Definition:</b> <a href="dspl_8h_source.html#l00094">dspl.h:94</a></div></div>
|
|
</div><!-- fragment --></div><!-- contents -->
|
|
<!-- HTML footer for doxygen 1.8.13-->
|
|
<!-- start footer part -->
|
|
<hr class="footer"/><address class="footer"><small>
|
|
Документация по libdspl-2.0. Последние изменения: Пн 4 Ноя 2019 16:56:51. Создано системой  <a href="http://www.doxygen.org/index.html">
|
|
<img class="footer" src="doxygen.png" alt="doxygen"/>
|
|
</a> 1.8.15
|
|
</small></address>
|
|
</div>
|
|
</body>
|
|
</html> |