Add an image for the remezb

[scilab.git] / scilab / modules / signal_processing / help / en_US / filters / remezb.xml

diff --git a/scilab/modules/signal_processing/help/en_US/filters/remezb.xml b/scilab/modules/signal_processing/help/en_US/filters/remezb.xml
index 8f68f93..38d76fc 100644 (file)
--- a/scilab/modules/signal_processing/help/en_US/filters/remezb.xml
+++ b/scilab/modules/signal_processing/help/en_US/filters/remezb.xml
@@ -57,7 +57,7 @@
             can be obtained from the output of the function
             by using the following commands
         </para>
-        <programlisting role=""><![CDATA[ 
+        <programlisting role="no-scilab-exec"><![CDATA[ 
 hn(1:nc-1)=an(nc:-1:2)/2;
 hn(nc)=an(1);
 hn(nc+1:2*nc-1)=an(2:nc)/2;
@@ -112,6 +112,41 @@ hn(nc+1:2*nc-1)=an(2:nc)/2;
 // Plot the filter's magnitude response
 plot(.5*(0:255)/256,frmag(hn,256));
  ]]></programlisting>
+
+<scilab:image>
+nc=21;ngrid=nc*16;
+fg=.24*(0:ngrid/2-1)/(ngrid/2-1);
+fg(ngrid/2+1:ngrid)=fg(1:ngrid/2)+.26*ones(1:ngrid/2);
+
+ds(1:ngrid/2)=ones(1:ngrid/2);
+ds(ngrid/2+1:ngrid)=zeros(1:ngrid/2);
+
+wt=ones(fg);
+
+an=remezb(nc,fg,ds,wt)
+
+hn(1:nc-1)=an(nc:-1:2)/2;
+hn(nc)=an(1);
+hn(nc+1:2*nc-1)=an(2:nc)/2;
+
+plot(.5*(0:255)/256,frmag(hn,256));
+
+nc=21; ngrid=nc*16;
+fg=.5*(0:(ngrid-1))/ngrid;
+
+ds(1:ngrid/2)=(0:ngrid/2-1)/(ngrid/2-1);
+ds(ngrid/2+1:ngrid)=ds(ngrid/2:-1:1);
+
+wt=ones(fg);
+
+an=remezb(nc,fg,ds,wt)
+
+hn(1:nc-1)=an(nc:-1:2)/2;
+hn(nc)=an(1);
+hn(nc+1:2*nc-1)=an(2:nc)/2;
+
+plot(.5*(0:255)/256,frmag(hn,256));
+</scilab:image>
     </refsection>
     <refsection role="see also">
         <title>See Also</title>