Interpolation d'une fonction
salut,
j'ai le problème suivant :
je fais
Code:
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:
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:
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...
