kopia lustrzana https://github.com/Dsplib/libdspl-2.0
modified: dspl/dox/ru/fft.dox
modified: dspl/dox/ru/groups_define.dox modified: dspl/dox/ru/randgen.dox modified: dspl/dox/ru/win.dox modified: include/dspl.hpull/6/merge
Dsplib
rodzic
5e49724aee
commit
2acda0aeea
Plik diff jest za duży
Load Diff
|
|
@ -0,0 +1 @@
|
|||
\relax
|
||||
Plik binarny nie jest wyświetlany.
|
|
@ -0,0 +1,126 @@
|
|||
This is pdfTeX, Version 3.14159265-2.6-1.40.20 (TeX Live 2019/W32TeX) (preloaded format=latex 2020.2.28) 25 APR 2020 16:58
|
||||
entering extended mode
|
||||
restricted \write18 enabled.
|
||||
%&-line parsing enabled.
|
||||
**_formulas.tex
|
||||
(./_formulas.tex
|
||||
LaTeX2e <2020-02-02> patch level 4
|
||||
L3 programming layer <2020-02-25>
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/base/article.cls
|
||||
Document Class: article 2019/12/20 v1.4l Standard LaTeX document class
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/base/size10.clo
|
||||
File: size10.clo 2019/12/20 v1.4l Standard LaTeX file (size option)
|
||||
)
|
||||
\c@part=\count167
|
||||
\c@section=\count168
|
||||
\c@subsection=\count169
|
||||
\c@subsubsection=\count170
|
||||
\c@paragraph=\count171
|
||||
\c@subparagraph=\count172
|
||||
\c@figure=\count173
|
||||
\c@table=\count174
|
||||
\abovecaptionskip=\skip47
|
||||
\belowcaptionskip=\skip48
|
||||
\bibindent=\dimen134
|
||||
)
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/base/ifthen.sty
|
||||
Package: ifthen 2014/09/29 v1.1c Standard LaTeX ifthen package (DPC)
|
||||
)
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics/epsfig.sty
|
||||
Package: epsfig 2017/06/25 v1.7b (e)psfig emulation (SPQR)
|
||||
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics/graphicx.sty
|
||||
Package: graphicx 2019/11/30 v1.2a Enhanced LaTeX Graphics (DPC,SPQR)
|
||||
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics/keyval.sty
|
||||
Package: keyval 2014/10/28 v1.15 key=value parser (DPC)
|
||||
\KV@toks@=\toks14
|
||||
)
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics/graphics.sty
|
||||
Package: graphics 2019/11/30 v1.4a Standard LaTeX Graphics (DPC,SPQR)
|
||||
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics/trig.sty
|
||||
Package: trig 2016/01/03 v1.10 sin cos tan (DPC)
|
||||
)
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics-cfg/graphics.cfg
|
||||
File: graphics.cfg 2016/06/04 v1.11 sample graphics configuration
|
||||
)
|
||||
Package graphics Info: Driver file: dvips.def on input line 105.
|
||||
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/graphics-def/dvips.def
|
||||
File: dvips.def 2017/06/20 v3.1d Graphics/color driver for dvips
|
||||
))
|
||||
\Gin@req@height=\dimen135
|
||||
\Gin@req@width=\dimen136
|
||||
)
|
||||
\epsfxsize=\dimen137
|
||||
\epsfysize=\dimen138
|
||||
)
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/base/inputenc.sty
|
||||
Package: inputenc 2018/08/11 v1.3c Input encoding file
|
||||
\inpenc@prehook=\toks15
|
||||
\inpenc@posthook=\toks16
|
||||
)
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/newunicodechar/newunicodechar.sty
|
||||
Package: newunicodechar 2018/04/08 v1.2 Defining Unicode characters
|
||||
)
|
||||
|
||||
Package newunicodechar Warning: Redefining Unicode character on input line 7.
|
||||
|
||||
|
||||
Package newunicodechar Warning: Redefining Unicode character on input line 8.
|
||||
|
||||
(c:/texlive/2019/texmf-dist/tex/latex/l3backend/l3backend-dvips.def
|
||||
File: l3backend-dvips.def 2020-02-23 L3 backend support: dvips
|
||||
\l__pdf_internal_box=\box45
|
||||
\g__pdf_backend_object_int=\count175
|
||||
\l__pdf_backend_content_box=\box46
|
||||
\l__pdf_backend_model_box=\box47
|
||||
\g__pdf_backend_annotation_int=\count176
|
||||
\g__pdf_backend_link_int=\count177
|
||||
\g__pdf_backend_link_sf_int=\count178
|
||||
)
|
||||
(./_formulas.aux)
|
||||
\openout1 = `_formulas.aux'.
|
||||
|
||||
LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: Checking defaults for TS1/cmr/m/n on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 11.
|
||||
LaTeX Font Info: ... okay on input line 11.
|
||||
LaTeX Font Info: External font `cmex10' loaded for size
|
||||
(Font) <7> on input line 12.
|
||||
LaTeX Font Info: External font `cmex10' loaded for size
|
||||
(Font) <5> on input line 12.
|
||||
[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] [100] [101] [102]
|
||||
[103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114]
|
||||
[115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126]
|
||||
[127] [128] [129] [130] [131] [132] [133] [134] [135] [136] [137]
|
||||
(./_formulas.aux) )
|
||||
Here is how much of TeX's memory you used:
|
||||
847 strings out of 480998
|
||||
12712 string characters out of 5916110
|
||||
246038 words of memory out of 5000000
|
||||
16407 multiletter control sequences out of 15000+600000
|
||||
532338 words of font info for 24 fonts, out of 8000000 for 9000
|
||||
1141 hyphenation exceptions out of 8191
|
||||
38i,6n,29p,293b,90s stack positions out of 5000i,500n,10000p,200000b,80000s
|
||||
|
||||
Output written on _formulas.dvi (137 pages, 19024 bytes).
|
||||
|
|
@ -0,0 +1,423 @@
|
|||
\documentclass{article}
|
||||
\usepackage{ifthen}
|
||||
\usepackage{epsfig}
|
||||
\usepackage[utf8]{inputenc}
|
||||
\usepackage{newunicodechar}
|
||||
\newunicodechar{⁻}{${}^{-}$}% Superscript minus
|
||||
\newunicodechar{²}{${}^{2}$}% Superscript two
|
||||
\newunicodechar{³}{${}^{3}$}% Superscript three
|
||||
|
||||
\pagestyle{empty}
|
||||
\begin{document}
|
||||
\[ e = \max \left( \frac{|x(k) - y(k)| }{ |x(k)|} \right), \quad \quad |x(k)| > 0, \]
|
||||
\pagebreak
|
||||
|
||||
\[ e = \max(|x(k) - y(k)| ), ~\qquad \quad~|x(k)| = 0, \]
|
||||
\pagebreak
|
||||
|
||||
$ e$
|
||||
\pagebreak
|
||||
|
||||
\[ e = \max \left( \frac{|x(k) - y(k)|}{|x(k)|} \right), \quad \quad |x(k)| > 0, \]
|
||||
\pagebreak
|
||||
|
||||
$ C_ord(x)$
|
||||
\pagebreak
|
||||
|
||||
\[ C_ord(x) = 2 x C_{ord-1}(x) - C_{ord-2}(x), \]
|
||||
\pagebreak
|
||||
|
||||
$ C_0(x) = 1 $
|
||||
\pagebreak
|
||||
|
||||
$ C_1(x) = x$
|
||||
\pagebreak
|
||||
|
||||
$ U_{ord}(x)$
|
||||
\pagebreak
|
||||
|
||||
\[ U_{ord}(x) = 2 x U_{ord-1}(x) - U_{ord-2}(x), \]
|
||||
\pagebreak
|
||||
|
||||
$ U_0(x) = 1 $
|
||||
\pagebreak
|
||||
|
||||
$ U_1(x) = 2x$
|
||||
\pagebreak
|
||||
|
||||
$ x = a + j b $
|
||||
\pagebreak
|
||||
|
||||
\[ |x|^2 = x x^* = a^2 + b^2. \]
|
||||
\pagebreak
|
||||
|
||||
\[ \textrm{Arccos}(x) = \frac{\pi}{2} - \textrm{Arcsin}(x) = \frac{\pi}{2} -j \textrm{Ln}\left( j x + \sqrt{1 - x^2} \right) \]
|
||||
\pagebreak
|
||||
|
||||
\[ \textrm{Arcsin}(x) = j \textrm{Ln}\left( j x + \sqrt{1 - x^2} \right) \]
|
||||
\pagebreak
|
||||
|
||||
\[ \textrm{cos}(x) = \frac{\exp(jx) + \exp(-jx)}{2} \]
|
||||
\pagebreak
|
||||
|
||||
\[ \textrm{Ln}(x) = j \varphi + \ln(|x|), \]
|
||||
\pagebreak
|
||||
|
||||
$\varphi$
|
||||
\pagebreak
|
||||
|
||||
\[ \textrm{sin}(x) = \frac{\exp(jx) - \exp(-jx)}{2j} \]
|
||||
\pagebreak
|
||||
|
||||
\[ y(k) = \sqrt{x(k)}, \qquad k = 0 \ldots n-1. \]
|
||||
\pagebreak
|
||||
|
||||
$ c = a * b$
|
||||
\pagebreak
|
||||
|
||||
$a$
|
||||
\pagebreak
|
||||
|
||||
$b$
|
||||
\pagebreak
|
||||
|
||||
$n = n_0 \times n_1 \times n_2 \times \ldots \times n_p \times m$
|
||||
\pagebreak
|
||||
|
||||
$n_i = 2,3,5,7$
|
||||
\pagebreak
|
||||
|
||||
$m $
|
||||
\pagebreak
|
||||
|
||||
\[ H(z) = \frac{\sum_{n = 0}^{N} b_n z^{-n}} {1+{\frac{1}{a_0}}\sum_{m = 1}^{M} a_m z^{-m}}, \]
|
||||
\pagebreak
|
||||
|
||||
$a_0$
|
||||
\pagebreak
|
||||
|
||||
$N=M=$
|
||||
\pagebreak
|
||||
|
||||
$H(z)$
|
||||
\pagebreak
|
||||
|
||||
$s(t) = \sin(2\pi \cdot 0.05 t) + n(t)$
|
||||
\pagebreak
|
||||
|
||||
$n(t)$
|
||||
\pagebreak
|
||||
|
||||
$ n $
|
||||
\pagebreak
|
||||
|
||||
$ x(m) $
|
||||
\pagebreak
|
||||
|
||||
$ m = 0 \ldots n-1 $
|
||||
\pagebreak
|
||||
|
||||
\[ Y(k) = \sum_{m = 0}^{n-1} x(m) \exp \left( -j \frac{2\pi}{n} m k \right), \]
|
||||
\pagebreak
|
||||
|
||||
$ k = 0 \ldots n-1 $
|
||||
\pagebreak
|
||||
|
||||
$x(m)$
|
||||
\pagebreak
|
||||
|
||||
$n$
|
||||
\pagebreak
|
||||
|
||||
$Y(k)$
|
||||
\pagebreak
|
||||
|
||||
$ n^2 $
|
||||
\pagebreak
|
||||
|
||||
\[ y(k) = \sum_{m = 0}^{n-1} x(m) \exp \left( j \frac{2\pi}{n} m k \right), \]
|
||||
\pagebreak
|
||||
|
||||
$y(k)$
|
||||
\pagebreak
|
||||
|
||||
$ u = \textrm{cd}^{-1}(w, k)$
|
||||
\pagebreak
|
||||
|
||||
$ w $
|
||||
\pagebreak
|
||||
|
||||
$ k $
|
||||
\pagebreak
|
||||
|
||||
$ u = \textrm{sn}^{-1}(w, k)$
|
||||
\pagebreak
|
||||
|
||||
$ y = \textrm{cd}(u K(k), k)$
|
||||
\pagebreak
|
||||
|
||||
$ u $
|
||||
\pagebreak
|
||||
|
||||
$ k_i $
|
||||
\pagebreak
|
||||
|
||||
$ K(k) $
|
||||
\pagebreak
|
||||
|
||||
\[ K(k) = \frac{\pi}{2} \prod_{i = 1}^{\infty}(1+k_i), \]
|
||||
\pagebreak
|
||||
|
||||
$ k_0 = k$
|
||||
\pagebreak
|
||||
|
||||
\[ k_i = \left( \frac{k_{i-1}} { 1+\sqrt{1-k_{i-1}^2} } \right)^2 \]
|
||||
\pagebreak
|
||||
|
||||
$ k<1 $
|
||||
\pagebreak
|
||||
|
||||
$ y = \textrm{sn}(u K(k), k)$
|
||||
\pagebreak
|
||||
|
||||
\[ Y(k) = \frac{1}{N} \sum_{m = 0}^{n-1} x(m) \exp \left( j \frac{2\pi}{n} m k \right), \]
|
||||
\pagebreak
|
||||
|
||||
$n = n_0 \times n_1 \times n_2 \times n_3 \times \ldots \times n_p \times m$
|
||||
\pagebreak
|
||||
|
||||
$n = n_0 \times n_1 \times n_2 \ldots \times n_p \times m$
|
||||
\pagebreak
|
||||
|
||||
$ n = 725760 $
|
||||
\pagebreak
|
||||
|
||||
$725760 = 2 \cdot 3 \cdot 4 \cdot 5 \cdot 6 \cdot 7 \cdot 9 \cdot 16 $
|
||||
\pagebreak
|
||||
|
||||
$ n = 172804 = 43201 \cdot 4 $
|
||||
\pagebreak
|
||||
|
||||
$ n = 13 \cdot 17 \cdot 23 \cdot 13 = 66079 $
|
||||
\pagebreak
|
||||
|
||||
$\sqrt{2^{31}} = 46340.95$
|
||||
\pagebreak
|
||||
|
||||
$x_0$
|
||||
\pagebreak
|
||||
|
||||
$x_1$
|
||||
\pagebreak
|
||||
|
||||
$x(k) = x_0 + k \cdot dx$
|
||||
\pagebreak
|
||||
|
||||
$dx = \frac{x_1 - x_0}{n-1}$
|
||||
\pagebreak
|
||||
|
||||
$k = 0 \ldots n-1.$
|
||||
\pagebreak
|
||||
|
||||
$dx = \frac{x_1 - x_0}{n}$
|
||||
\pagebreak
|
||||
|
||||
$10^{x_0}$
|
||||
\pagebreak
|
||||
|
||||
$10^{x_1}$
|
||||
\pagebreak
|
||||
|
||||
$x(k) = 10^{x_0} \cdot dx^k$
|
||||
\pagebreak
|
||||
|
||||
$dx = \sqrt[n-1]{10^{x_1 - x_0}}$
|
||||
\pagebreak
|
||||
|
||||
$dx = \sqrt[n]{10^{x_1 - x_0}}$
|
||||
\pagebreak
|
||||
|
||||
$ H(j \omega) $
|
||||
\pagebreak
|
||||
|
||||
$ H(j \omega)$
|
||||
\pagebreak
|
||||
|
||||
$ H(s) $
|
||||
\pagebreak
|
||||
|
||||
\[ H(s) = \frac {\sum_{k = 0}^{N} b_k s^k} {\sum_{m = 0}^{N} a_m s^m}, \]
|
||||
\pagebreak
|
||||
|
||||
$ N $
|
||||
\pagebreak
|
||||
|
||||
$ s = j \omega $
|
||||
\pagebreak
|
||||
|
||||
$ \omega $
|
||||
\pagebreak
|
||||
|
||||
$H(s)$
|
||||
\pagebreak
|
||||
|
||||
$ H \left(\mathrm{e}^{j\omega} \right) $
|
||||
\pagebreak
|
||||
|
||||
$ 2\pi $
|
||||
\pagebreak
|
||||
|
||||
$ \pi $
|
||||
\pagebreak
|
||||
|
||||
$ -\pi $
|
||||
\pagebreak
|
||||
|
||||
$ H \left(e^{j \omega} \right)$
|
||||
\pagebreak
|
||||
|
||||
\[ H(z) = \frac {\sum_{k = 0}^{N} b_k z^{-k}} {\sum_{m = 0}^{N} a_m z^{-m}}, \]
|
||||
\pagebreak
|
||||
|
||||
$N$
|
||||
\pagebreak
|
||||
|
||||
$z = e^{j \omega} $
|
||||
\pagebreak
|
||||
|
||||
$\omega$
|
||||
\pagebreak
|
||||
|
||||
$ 2 \pi-$
|
||||
\pagebreak
|
||||
|
||||
$2 \pi$
|
||||
\pagebreak
|
||||
|
||||
$-\pi$
|
||||
\pagebreak
|
||||
|
||||
$ \pi$
|
||||
\pagebreak
|
||||
|
||||
$ H \left(e^{j \omega} \right) = H^* \left(e^{-j \omega} \right)$
|
||||
\pagebreak
|
||||
|
||||
$\pi$
|
||||
\pagebreak
|
||||
|
||||
$ -R_p $
|
||||
\pagebreak
|
||||
|
||||
$ H(s)$
|
||||
\pagebreak
|
||||
|
||||
$-R_p$
|
||||
\pagebreak
|
||||
|
||||
$ R_p $
|
||||
\pagebreak
|
||||
|
||||
$-R_s$
|
||||
\pagebreak
|
||||
|
||||
$H(j\cdot 1) = -R_s$
|
||||
\pagebreak
|
||||
|
||||
\[ H(s) = \frac{\sum_{n = 0}^{N_z} b_n \cdot s^n}{\sum_{m = 0}^{N_p} a_m \cdot s^m} = \frac{\prod_{n = 0}^{N_z}(s-z_n)}{\prod_{m = 0}^{N_p} (s-p_m)} \]
|
||||
\pagebreak
|
||||
|
||||
\[ H(z) = \sum_{n = 0}^{ord} h_n z^{-n} \]
|
||||
\pagebreak
|
||||
|
||||
$ F(s) $
|
||||
\pagebreak
|
||||
|
||||
$F(s)$
|
||||
\pagebreak
|
||||
|
||||
$Y(s) = (H \circ F)(s) = H(F(s))$
|
||||
\pagebreak
|
||||
|
||||
\[ H(s) = \frac{\sum\limits_{m = 0}^{n} b_m s^m} {\sum\limits_{k = 0}^{n} a_k s^k}, \quad F(s) = \frac{\sum\limits_{m = 0}^{p} d_m s^m} {\sum\limits_{k = 0}^{p} c_k s^k}, \quad Y(s) = \frac{\sum\limits_{m = 0}^{n p} \beta_m s^m} {\sum\limits_{k = 0}^{n p} \alpha_k s^k} \]
|
||||
\pagebreak
|
||||
|
||||
$Y(s) = (H \circ F)(s)$
|
||||
\pagebreak
|
||||
|
||||
\[ s \leftarrow \frac{1 - z^{-1}}{1 - z^{-1}}. \]
|
||||
\pagebreak
|
||||
|
||||
$\Omega$
|
||||
\pagebreak
|
||||
|
||||
\[ \Omega = \tan(\omega / 2). \]
|
||||
\pagebreak
|
||||
|
||||
\[ s(t) = \sum\limits_{n = 0}^{n_{\omega}-1} S(\omega_n) \exp(j\omega_n t) \]
|
||||
\pagebreak
|
||||
|
||||
$\omega_n$
|
||||
\pagebreak
|
||||
|
||||
$S(\omega_n)$
|
||||
\pagebreak
|
||||
|
||||
$ I_0(x)$
|
||||
\pagebreak
|
||||
|
||||
$ x $
|
||||
\pagebreak
|
||||
|
||||
$[0 \ 3]$
|
||||
\pagebreak
|
||||
|
||||
$ \textrm{sinc}(x,a) = \frac{\sin(ax)}{ax}$
|
||||
\pagebreak
|
||||
|
||||
\[ \textrm{Si}(x) = \int_{0}^{x} \frac{\sin(x)}{x} \, dx\]
|
||||
\pagebreak
|
||||
|
||||
$[-6\pi \ 6\pi]$
|
||||
\pagebreak
|
||||
|
||||
$P_N(x)$
|
||||
\pagebreak
|
||||
|
||||
$N-$
|
||||
\pagebreak
|
||||
|
||||
\[ P_N(x) = a_0 + a_1 \cdot x + a_2 \cdot x^2 + a_3 \cdot x^3 + ... a_N \cdot x^N. \]
|
||||
\pagebreak
|
||||
|
||||
\[ P_N(x) = a_0 + x \cdot (a_1 + x \cdot (a_2 + \cdot ( \ldots x \cdot (a_{N-1} + x\cdot a_N) \ldots ))) \]
|
||||
\pagebreak
|
||||
|
||||
$10^{56}$
|
||||
\pagebreak
|
||||
|
||||
$x(i)$
|
||||
\pagebreak
|
||||
|
||||
$i = 0,1,\ldots n$
|
||||
\pagebreak
|
||||
|
||||
$y(i)$
|
||||
\pagebreak
|
||||
|
||||
\[ y(i) = k_x x(i) + d_x, \qquad k_x = \frac{h}{x_{\textrm{max}} - x_{\textrm{min}}}. \]
|
||||
\pagebreak
|
||||
|
||||
$x_{\textrm{min}}$
|
||||
\pagebreak
|
||||
|
||||
$x_{\textrm{max}}$
|
||||
\pagebreak
|
||||
|
||||
$d_x$
|
||||
\pagebreak
|
||||
|
||||
$h + d_x$
|
||||
\pagebreak
|
||||
|
||||
\end{document}
|
||||
|
|
@ -35,7 +35,6 @@
|
|||
\param t1
|
||||
Указатель на вектор промежуточных вычислений алгоритма БПФ. \n
|
||||
Размер вектора `[n x 1]`. \n
|
||||
|
||||
Память должна быть выделена функцией \ref fft_create. \n \n
|
||||
Структура заполняется функцией \ref fft_create один раз
|
||||
до использования алгоритма БПФ. \n
|
||||
|
|
@ -44,7 +43,6 @@
|
|||
Перед выходом из программы выделенную память под поворотные
|
||||
коэффициенты и массивы промежуточных данных
|
||||
необходимо очистить функцией \ref fft_free. Например:
|
||||
|
||||
\code
|
||||
fft_t pfft; /* объявляем объект fft_t */
|
||||
int n = 64; /* Размер БПФ */
|
||||
|
|
@ -70,14 +68,14 @@ err = fft_create(&pfft, n);
|
|||
fft_free(&pfft);
|
||||
\endcode
|
||||
|
||||
\note
|
||||
Важно отметить, что если объект `fft_t` был создан для размера БПФ равного `n`,
|
||||
то он может быть использован только для БПФ размера `n`. \n
|
||||
|
||||
то он может быть использован только для БПФ размера `n`. \n\n
|
||||
Также необходимо заметить, что функции БПФ самостоятельно контролируют размер,
|
||||
и самостоятельно выделяют память объекта БПФ при необходимости.
|
||||
Так если вызвать любую функцию использующую структуру `fft_t` с заполненными
|
||||
данными для длины БПФ `k` для расчета БПФ длины `n`,
|
||||
то массивы стуктуры будут автоматически пересозданы для длины `n`.
|
||||
то массивы структуры будут автоматически пересозданы для длины `n`.
|
||||
|
||||
\author
|
||||
Бахурин Сергей.
|
||||
|
|
|
|||
|
|
@ -1,4 +1,4 @@
|
|||
/*!
|
||||
/*! ****************************************************************************
|
||||
|
||||
\defgroup TYPES_GROUP Типы данных DSPL
|
||||
В данной группе описаны типы данных библиотеки DSPL и методы их преобразования.
|
||||
|
|
@ -37,7 +37,11 @@
|
|||
\defgroup SPEC_MATH_ELLIP_GROUP Эллиптические функции Якоби
|
||||
Эллиптические функции Якоби вещественного и комплексного аргумента
|
||||
|
||||
\defgroup SPEC_MATH_RAND_GEN_GROUP Генерация случайных чисел
|
||||
|
||||
\defgroup SPEC_MATH_RAND_GEN_GROUP Генерация псевдослучайных чисел
|
||||
В данном разделе описаны функции генерации псевдослучайных чисел с
|
||||
использованием различных датчиков.
|
||||
|
||||
|
||||
\defgroup SPEC_MATH_STAT_GROUP Функции математической статистики
|
||||
|
||||
|
|
@ -56,4 +60,4 @@
|
|||
В данной группе приведены возможные коды ошибок,
|
||||
возвращаемые функциями библиотеки в процессе работы.
|
||||
|
||||
*/
|
||||
**************************************************************************** */
|
||||
|
|
|
|||
|
|
@ -1,3 +1,105 @@
|
|||
/*! ****************************************************************************
|
||||
\ingroup SPEC_MATH_RAND_GEN_GROUP
|
||||
\struct random_t
|
||||
\brief Структура параметров датчиков псевдослучайных чисел.
|
||||
|
||||
Структура хранит инициализацию и текущие регистры различных датчиков
|
||||
псевдослучайных чисел. В библиотеке используются следующие датчики:
|
||||
\li MRG32K3A -- 32 битный датчик разработан Пьером Лекуэром [1].
|
||||
\li MT19937-64 -- 64-битный датчик
|
||||
<a href = "https://en.wikipedia.org/wiki/Mersenne_Twister">
|
||||
Вихрь Мерсенна
|
||||
</a> [2, 3].
|
||||
|
||||
\note
|
||||
[1] Pierre L'Ecuyer, (1999) Good Parameters and Implementations for Combined
|
||||
Multiple Recursive Random Number Generators. Operations Research
|
||||
47(1):159-164. https://doi.org/10.1287/opre.47.1.159 \n\n
|
||||
[2] T. Nishimura, ``Tables of 64-bit Mersenne Twisters // ACM Transactions
|
||||
on Modeling and Computer Simulation 10. (2000) 348--357. \n\n
|
||||
[3] M. Matsumoto and T. Nishimura Mersenne Twister: a 623-dimensionally
|
||||
equidistributed uniform pseudorandom number generator // ACM Transactions
|
||||
on Modeling and Computer Simulation 8. (Jan. 1998) 3--30. \n\n
|
||||
|
||||
\param mrg32k3a_seed
|
||||
Начальная инициализация датчика MRG32K3A. \n \n
|
||||
|
||||
|
||||
\param mrg32k3a_x
|
||||
Первый вектор состояния рекурсивного датчика MRG32K3A. \n \n
|
||||
|
||||
\param mrg32k3a_y
|
||||
Второй вектор состояния рекурсивного датчика MRG32K3A. \n \n
|
||||
|
||||
\param mt19937_mt
|
||||
Первый вектор состояния рекурсивного датчика MT19937-64. \n \n
|
||||
|
||||
\param mt19937_mti
|
||||
Текущий индекс в векторе состояния датчика MT19937-64. \n \n
|
||||
|
||||
Параметры данной структуры заполняются автоматически функцией `random_init`
|
||||
и используются функциями генерации псевдослучайных векторов.
|
||||
|
||||
\author
|
||||
Бахурин Сергей.
|
||||
www.dsplib.org
|
||||
***************************************************************************** */
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/*! ****************************************************************************
|
||||
\ingroup SPEC_MATH_RAND_GEN_GROUP
|
||||
\fn int random_init(random_t* prnd, int type, void* seed)
|
||||
\brief
|
||||
Инициализация датчиков псевдослучайных чисел.
|
||||
|
||||
\param[in,out] prnd
|
||||
Указатель на структуру параметров и векторов состояния
|
||||
датчиков псевдослучайных чисел, которая будет инициализирована. \n\n
|
||||
|
||||
\param[in] type
|
||||
Тип датчика псевдослучайных чисел:
|
||||
\verbatim
|
||||
RAND_TYPE_MRG32K3A - 32-битный датчик MRG32K3A
|
||||
RAND_TYPE_MT19937 - 64-битный датчик MT19937-64
|
||||
\endverbatim
|
||||
|
||||
\param[in] seed
|
||||
Указатель на начальную инициализацию датчика. \n
|
||||
Данный указатель имеет тип `void*`, поскольку параметр инициализации
|
||||
зависит от типа датчика. Например если инициализируем датчик MRG32K3A,
|
||||
т.е. параметр `type` задан как `RAND_TYPE_MRG32K3A`, то данный указатель
|
||||
приводится к типу `double`:
|
||||
\code
|
||||
random_t rnd = {0};
|
||||
double seed = 1234.0;
|
||||
random_init(&rnd, RAND_TYPE_MRG32K3A, (void*)&seed);
|
||||
\endcode
|
||||
Если же используется 64-битный датчик Вихрь Мерсенна
|
||||
(`type` задан как `RAND_TYPE_MT19937`), то `seed` приводится к типу
|
||||
`unsigned long long`:
|
||||
\code
|
||||
random_t rnd = {0};
|
||||
unsigned long long = 1234353456;
|
||||
random_init(&rnd, RAND_TYPE_MT19937, (void*)&seed);
|
||||
\endcode
|
||||
При фиксированном начальном значении датчика, псевдослучайные числа будут
|
||||
повторяться при каждом запуске программы. \n
|
||||
Указатель `seed` может быть `NULL`. В этом случае начальная инициализация
|
||||
датчиков будет задаваться случайными значениями и генерируемые псевдослучайные
|
||||
числа будут различными при каждом запуске программы.
|
||||
|
||||
\author
|
||||
Бахурин Сергей.
|
||||
www.dsplib.org
|
||||
**************************************************************************** */
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
/*! ****************************************************************************
|
||||
|
|
@ -6,7 +108,7 @@
|
|||
|
||||
\brief
|
||||
Генерация вектора нормально распределенных псевдослучайных чисел с
|
||||
заданным математическим ожиданием и среднеквадратическим отклонением.
|
||||
заданным математическим ожиданием и среднеквадратичным отклонением.
|
||||
|
||||
Генерация производится при помощи преобразования Бокса — Мюллера равномерно-распределенной
|
||||
случайной величины в нормально распределенную. \n
|
||||
|
|
@ -28,11 +130,11 @@
|
|||
\n
|
||||
|
||||
\param[in] sigma
|
||||
Среднеквадратическое отклонение (СКО). \n
|
||||
Среднеквадратичное отклонение (СКО). \n
|
||||
\n
|
||||
|
||||
\return
|
||||
`RES_OK` --- вектор случайных чисел сгенерирован успешно. \n
|
||||
`RES_OK` --- вектор случайных чисел рассчитан успешно. \n
|
||||
В противном случае \ref ERROR_CODE_GROUP "код ошибки".
|
||||
|
||||
\author
|
||||
|
|
@ -46,24 +148,32 @@ www.dsplib.org
|
|||
|
||||
/*! ****************************************************************************
|
||||
\ingroup SPEC_MATH_RAND_GEN_GROUP
|
||||
\fn int randu(double* x, int n);
|
||||
\fn int randu(double* x, int n, random_t* prnd)
|
||||
\brief
|
||||
Генерация вектора равномерно-распределенных в интервале
|
||||
от 0 до 1 псевдослучайных чисел.
|
||||
|
||||
Генерация производится при помощи рекурсивного алгоритма L'Ecluyer.
|
||||
Период датчика порядка \f$10^{56}\f$. \n
|
||||
|
||||
|
||||
\param[in,out] x Указатель на вектор случайных чисел. \n
|
||||
\param[in,out] x
|
||||
Указатель на вектор случайных чисел. \n
|
||||
Размер вектора `[n x 1]`. \n
|
||||
Память должна быть выделена. \n \n
|
||||
Память должна быть выделена. \n\n
|
||||
|
||||
\param[in] n
|
||||
Размер вектора случайных чисел. \n \n
|
||||
Размер вектора случайных чисел. \n\n
|
||||
|
||||
\param[in] prnd
|
||||
Указатель на структуру `random_t` параметров датчиков
|
||||
псевдослучайных чисел. \n
|
||||
Структура должна быть предварительно заполнена функцией \ref random_init. \n
|
||||
Данный указатель может быть `NULL`, тогда будет использоваться
|
||||
встроенный датчик, определенный стандартом языка Си. Однако для серьезных нужд,
|
||||
например в криптографии, данный режим использовать не рекомендуется.
|
||||
Нет гарантии в качестве произведенной случайной последовательности если
|
||||
параметр `prnd` задан как `NULL`. \n\n
|
||||
|
||||
\return
|
||||
`RES_OK` --- вектор случайных чисел сгенерирован успешно. \n
|
||||
`RES_OK` --- вектор равномерно-распределенных
|
||||
псевдослучайных чисел рассчитан успешно. \n
|
||||
В противном случае \ref ERROR_CODE_GROUP "код ошибки".
|
||||
|
||||
\author
|
||||
|
|
|
|||
|
|
@ -68,8 +68,7 @@
|
|||
DSPL_WIN_RECT | Непараметрическое прямоугольное окно
|
||||
-------------------------------------------------------------------------
|
||||
\endverbatim
|
||||
\n
|
||||
|
||||
\n
|
||||
Маска `DSPL_WIN_SYM_MASK` задает симметричное
|
||||
или периодическое окно: \n
|
||||
\verbatim
|
||||
|
|
|
|||
|
|
@ -54,8 +54,8 @@ typedef struct
|
|||
|
||||
#define RAND_TYPE_MRG32K3A 0x00000001
|
||||
#define RAND_TYPE_MT19937 0x00000002
|
||||
|
||||
#define RAND_MT19937_NN 312
|
||||
|
||||
typedef struct
|
||||
{
|
||||
|
||||
|
|
|
|||
Ładowanie…
Reference in New Issue