Prev: inverse z transform
Next: GUI: create a pushbutton to restart from the beginning
From: Bruno Luong on 12 Apr 2010 15:58
Here is the "Conv" way:

a = 0.9;

P=1;
for k=1:10
P = conv(P,[1 -a^k]);
end
disp(P)

% Bruno
From: Walter Roberson on 12 Apr 2010 16:12
Julian Geiger wrote:
> Great, thx. Is there a way to order the result according to coefficients
> of x^k?
> (example: x*(-a^3 - a^5) + x^2*(a^5 + a^6) ...)

In Maple, the call would be named collect(); the MuPad name is probably fairly
similar.

Note: in Maple, sometimes the coefficients corresponding to x^0 occur at the
beginning and sometimes they occur at the end.

If you want to find the coefficients of a particular power of x, see coeff()
From: Bruno Luong on 12 Apr 2010 17:08
I have simplified the Vandermond method to a single line code:

a=0.9;
k = 10;
u=a.^(1:k);

% Engine
P=[1; vander(u)\(u(:).^k)]

% Bruno
First  |  Prev  | 
Pages: 1 2
Prev: inverse z transform
Next: GUI: create a pushbutton to restart from the beginning