[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: trivial digital filter
From: |
A S Hodel |
Subject: |
Re: trivial digital filter |
Date: |
Fri, 5 Mar 2004 14:34:41 -0600 |
Sorry for the delay in response; I've been updating my octave
installation on both of my machines.
Try either
sys = tf([1 1][1,0],1e-3);
or, equivalently,
sys = fir2sys([1, 1],1e-3);
The octave control systems toolbox requires "proper" transfer
functions, i.e., at least as many poles
as zeros. A causal FIR system with N zeros has at least N poles, all
of them at z = 0.
On Wednesday, March 3, 2004, at 09:21 AM, Hugo Coolens wrote:
I'm trying to simulate the following digital filter in octave:
sys=tf([1 1],[1],1e-3)
However this gives an error message:
# of poles (0) < # of zeros (1)
Which is of course clear enough but I wonder whether there is a way to
use
a "similar method" as the one above in octave for this kind of filter?
regards,
Hugo
To demonstrate the filter I already simulated its behaviour in the time
domain as follows:
a=[1];
b=[1 1];
fs=1000;
time=0:1/fs:40e-3;
x=sin(2*pi*200*time);
y=filter(b,a,x);
stem(time,y)
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------
A. S. Hodel Dept. ECE, 200 Broun Hall, Auburn University AL 36849-5201
(334) 844-1854/fax(334) 844-1809,
http://www.eng.auburn.edu/users/hodelas , address@hidden
-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.
Octave's home on the web: http://www.octave.org
How to fund new projects: http://www.octave.org/funding.html
Subscription information: http://www.octave.org/archive.html
-------------------------------------------------------------