Discussion :

[diana82]

voila,
mais j'ai probleme,j'essai de programmer un filtre recursif.
la fonction du genre:

y(k)=a*y(k-1)+b*y(k-1)+c*x(k)+d*x(k-1)

les x(k):sont les vecteurs d'entrees
j'ai un probleme a utiliser les indices.

[tug83]

C'est quoi excatement ton problème, ton erreur, ton code?
La section 9 de cette note qui utilise la fonction FILTER pourrait peut-être t'aider pour en tous cas vectoriser ton code et donc gagner en terme de performance

[diana82]

voila le code que j'ai fait:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
I=imread('test.bmp')
m=size(I,1);
n=size(I,2);
n00=1.68;
n11=0.6318;
n22=3.735;
n33=1.783;
d11=1.997;
d22=-0.6803;
d33=0.2598;
for i=1:m
        x=I(i,:);
        if i==1
       y(i)=n00*x(i)
        if i==2
            y(i)=-d11*y(i-1)+n00*x(i)+n11*x(i-1)
            if i==3
            y(i)=-d11*y(i-1)-d22*y(i-2)+n00*x(i)+n11*x(i-1)+n22*x(i-2)
               if i==4
            y(i)=-d11*y(i-1)-d22*y(i-2)-d33*y(i-3)+n00*x(i)+n11*x(i-1)+n22*x(i-2)+n33*x(i-3)
             if i>4
            y(k)=-d11*y(i-1)-d22*y(i-2)-d33*y(i-3)+n00*x(i)+n11*x(i-1)+n22*x(i-2)+n33*x(i-3)
               end
            end
        end
        end
        end
end

je suis sense avoir pour chaque vecteur d'entree,un vecteur de sortie.
mais ce n'est pas le cas.(je travail sous matlab7.0)