Advertisement
Guest User

Untitled

a guest
Feb 22nd, 2019
54
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.38 KB | None | 0 0
  1. function[vektor]=proba4(p)
  2.  
  3. elso_ertek=0;
  4. masodik_ertek=0;
  5.  
  6. gyokei=roots(p);
  7.  
  8. ertek=roots(polyder(p));
  9.  
  10. derivalt=polyder(p);
  11. masodik_derivalt=polyder(derivalt);
  12.  
  13. x=ertek-1; %teszt értékek
  14. y=ertek+1;
  15.  
  16. if polyval(derivalt,x)<0
  17. elso_ertek=-1;
  18. else
  19. elso_ertek=1;
  20. end
  21.  
  22. if polyval(derivalt,y)<0
  23. masodik_ertek=-1;
  24. else
  25. masodik_ertek=1;
  26. end
  27.  
  28. if elso_ertek==-1 && masodik_ertek==1
  29. loc_min=polyval(p,ertek);
  30. vektor(1)=loc_min;
  31. vektor(2)=0;
  32. elseif elso_ertek==1 && masodik_ertek==-1
  33. loc_max=polyval(p,ertek);
  34. vektor(1)=0;
  35. vektor(2)=loc_max;
  36. else
  37. vektor(1)=0;
  38. vektor(2)=0;
  39. end
  40.  
  41.  
  42. %itt az inflexiós pontot szeretném meghatározni.
  43.  
  44. inf_pont_ertek=roots(polyder(derivalt));
  45. %vizsgálom a nála kisebb és nagyobb elemeket, hogy ismerjem a konvexitást,
  46. %mert ezekből határozható meg az inflexiós pont
  47. if polyval(masodik_derivalt,inf_pont_ertek-1)<0
  48. masodik_elso=-1;
  49. else
  50. masodik_elso=1;
  51. end
  52. if polyval(masodik_derivalt,inf_pont_ertek+1)<0
  53. masodik_masodik=-1;
  54. else
  55. masodik_masodik=1;
  56. end
  57. %ugye ott van inflexiós pont ahol - ból + ba vált vagy fordítva, ha nem
  58. %vált nincs inflexiós pont
  59. if masodik_elso==-1 && masodik_masodik==1 || masodik_elso==1 && masodik_masodik==-1
  60. vektor(3)=polyval(p,inf_pont_ertek);
  61. else
  62. vektor(3)=0;
  63. end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement