← Back to branch summary

~ubuntu-branches/ubuntu/trusty/gnuradio/trusty-updates

« back to all changes in this revision

Viewing changes to gnuradio-core/src/utils/plot_cic_decimator_response.m

  • Committer: Package Import Robot
  • Author(s): A. Maitland Bottoms
  • Date: 2012-02-26 21:26:16 UTC
  • mfrom: (1.1.4)
  • Revision ID: package-import@ubuntu.com-20120226212616-vsfkbi1158xshdql
Tags: 3.5.1-1
http://bugs.debian.org/642716
http://bugs.debian.org/645332
http://bugs.debian.org/394849
http://bugs.debian.org/616832
http://bugs.debian.org/590048
http://bugs.debian.org/642580
* new upstream version, re-packaged from scratch with modern tools
    closes: #642716, #645332, #394849, #616832, #590048, #642580,
    #647018, #557050, #559640, #631863
* CMake build

Show diffs side-by-side

added added

removed removed

Lines of Context:
gnuradio-core/src/utils/plot_cic_decimator_response.m
 
1
#
 
2
# Copyright 2004 Free Software Foundation, Inc.
 
3
# 
 
4
# This file is part of GNU Radio
 
5
# 
 
6
# GNU Radio is free software; you can redistribute it and/or modify
 
7
# it under the terms of the GNU General Public License as published by
 
8
# the Free Software Foundation; either version 3, or (at your option)
 
9
# any later version.
 
10
# 
 
11
# GNU Radio is distributed in the hope that it will be useful,
 
12
# but WITHOUT ANY WARRANTY; without even the implied warranty of
 
13
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 
14
# GNU General Public License for more details.
 
15
# 
 
16
# You should have received a copy of the GNU General Public License
 
17
# along with GNU Radio; see the file COPYING.  If not, write to
 
18
# the Free Software Foundation, Inc., 51 Franklin Street,
 
19
# Boston, MA 02110-1301, USA.
 
20
# 
 
21
 
 
22
function plot_cic_decimator_response (R, N, M)
 
23
  ## R = decimation rate
 
24
  ## N = number of stages (4)
 
25
  ## M = 1
 
26
  gain = (R*M)^N
 
27
  npoints = 1024;
 
28
  w = 0:1/npoints:1-1/npoints;
 
29
  w = w * 1 * pi;
 
30
  ## w = w * R;
 
31
  length(w);
 
32
  num = sin (w*R*M/2);
 
33
  length (num);
 
34
  ## H = sin (w*R*M/2) ./ sin (w/2)
 
35
  denom = sin(w/2);
 
36
  length (denom);
 
37
  H = (num ./ denom) .^ N;
 
38
  H(1) = gain;
 
39
  H = H ./ gain;
 
40
  plot (R*w/(2*pi), 10 * log10 (H));
 
41
  ## plot (R*w/(2*pi), H);
 
42
endfunction
 
43
 
 
44
 
 
45
  
 
 
b'\\ No newline at end of file'