Facebook
From Chocolate Eider, 5 Years ago, written in Plain Text.
Embed
Download Paste or View Raw
Hits: 232
  1. function aprox
  2. x=[0:5];
  3. y=[1 2 3 5 9 10];
  4. p1=0;
  5. p2=0;
  6. f=0;
  7. for i=1:length(x)
  8.     p1(i)=exp(x(i));
  9.     p2(i)=1;
  10. end
  11. A=[p1*p1' p1*p2'; p2*p1' p2*p2'];
  12. v=[p1*y'; p2*y'];
  13. Q=inv(A)*v;
  14. a=Q(1)
  15. b=Q(2)
  16. X=[0:0.01:5];
  17. for i=1:length(X)
  18.     f(i)=a*exp(X(i))+b;
  19. end
  20. plot(X,f,'-'),xlabel('x'),ylabel('y');
  21. for i=1:length(x)
  22.     p1(i)=(x(i))^2;
  23.     p2(i)=x(i);
  24.     p3(i)=1;
  25. end
  26. A=[p1*p1' p1*p2' p1*p3'; p2*p1' p2*p2' p2*p3'; p3*p1' p3*p2' p3*p3'];
  27. v=[p1*y'; p2*y'; p3*y'];
  28. Q=inv(A)*v;
  29. a=Q(1)
  30. b=Q(2)
  31. c=Q(3)
  32. X=[0:0.01:5];
  33. for i=1:length(X)
  34.     f(i)=a*(X(i))^2+b*X(i)+c;
  35. end
  36. hold on
  37. plot(X,f,'-r'),xlabel('x'),ylabel('y');
  38. plot(x,y,'*')
  39. hold off