Discussion :

[nox78]

salut,

j'ai le problème suivant :

je fais

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
beta=0:0.1:90;
beta=beta.*pi/180; %conversion en radian
 
ws = 2 * pi * 60;
C = 1/(15*ws);
L = 2.56/ws;
w0 = 1/sqrt(L*C);
nu=wo/ws;
eta = w0 / ws;
S = 1 / ( 1 - 1 / eta^2 );
 
Ceff = 1 ./ (1./C - 4./pi .* ( 1./(2.*C) .* S .* (beta + sin(2.*beta)./2) + ws.^2 .* L .* S.^2 .* (cos(beta)).^2 .* ( tan(beta) - nu.*tan(nu.*beta) ) ) );
plot(beta,Ceff);

Je veux interpoler cette fonction. Donc je vais dans tools, basic fitting, je choisis le degré 9 et je copie le polynome dans une fonction.

Code :

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

1
2
3
function y = poly_Ceff(x)
 
y=-0.0033*x.^9+0.023*x.^8-0.069*x.^7+0.11*x.^6-0.092*x.^5+0.045*x.^4-0.013*x.^3+0.0017*x.^2-(9.8e-5)*x+0.00018;

Et pourtant lorsque je veux vérifier :

Code :

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

1
2
hold on;
plot(beta,poly_Ceff(beta),'r');

je ne trouve pas du tout le même graph :'(

Quelqu'un aurait une explication ?
J'obtiens le meme polynome en utilisant la fonction polyfit...

[salseropom]

Salut, je n'ai pas exécuté ton prgm mais ton résultat ne m'étonne pas du tout. Je pense que tu dois avoir un pb de bord : pour augmenter la précision de l'interpolation, tout le monde pense qu'il faut augmenter l'ordre du polynome (égalité entre f et P, entre f' et P', entre f'' et P'' etc...) mais plus tu montes dans l'ordre et plus les "bosses" de ton polynomes seront grandes.
Une solution consiste donc à ne pas augmenter l'ordre du polynome mais plutôt à choisir judicieusement ses points d'interpolation. Pour cela, utilise les polynomes de Chebychef qui place les points sur les extrémités de ton segment (si tu interpoles sur [-1 1] je crois que ces polynomes ont 90 % des points d'interpolations proches de -1 et 1). Je ne sais pas si j'ai été clair...
Ton pb est plus d'ordre numérique que de programation.

[nox78]

oui je connaissais le pb des effets de bord. J'ai choisi le degré de mon polynôme en comparant les restes. Quand les restes cessent de décroître je considère que j'ai atteint le degré optimal.

C'est une bonne facon de faire?

[nox78]

bon en fait j'ai trouvé le probleme.

Il faut utiliser polyfit et un affichage en format long pour avoir plus de précision.

Normal si je vais à x^9 il faut des coefficients tres précis et l'affichage par defaut ne va pas en dessous de 10e-4.

avec des coefficients en format long ca passe ^^


merci du coup de main en tout cas

[salseropom]

OK, je ne suis pas sûr de t'avoir été d'une grande aide, mais bon...ton pb est résolu.