Polinoame' PDF
Polinoame' PDF
Polinoame' PDF
L8/1. Să se calculeze rădăcinile polinomului h(x) = 62x9 – 43x8 + 45x6 – 34x3 + 61x – 96 şi să se determine eroarea de
calcul pentru o soluţie reală. Câte soluţii reale are ecuaţia h(x)=0.
L8/2. Sǎ se afle rădăcinile polinomului t(x) = x11 – 10x4 + 10x3 – 5x2 - 5x + 150
L8/3. Care sunt rădăcinile reale are polinomului f(x) = 345x8 – 4x6 + 45x5 – 2x2 + 15x – 56 ?
L8/4. Sǎ se determine rezultatul înmulţirii polinoamelor g(x) = 22x4 + 15x3 – 56 şi f(x).
L8/5. Sǎ se împartǎ h(x) / g(x) şi să se găsească restul şi câtul împărţirii.
L8/6. Sǎ se evalueze polinomul h(x) pentru x1 = 10, x2 = 25, x3 = -54.6 şi polinomul g(x) pentru x = -12.
L8/7. Să se aproximeze în sensul celor mai mici pătrate setul de date de mai jos şi să se traseze graficele polinomului de grad
3, grad 5 şi liniei care aproximează datele. Să se determine (prin calcul) care este cea mai bună aproximare.
X={0.1; 0.3; 0.4; 0.6; 0.7; 0.8; 0.9; 1; 1.2; 2.22; 5.6; 7.8}
Y={0.7; 0.43; 0.2; 0.02; -0.04; -0.09; -0.07; 0; 0.5; 0.65; 0.77; 0.86}
L8/8. Să se reprezinte grafic linia poligonală şi polinomul de interpolare determinate de punctele din planul Oxy având
coordonatele date în tabelul de mai jos:
X 0,25 0,5 0,75 1 1,25 1,5 1,75 2
Y 3 2,5 0,5 1,5 0,5 -1 3 0
Să se determine valoarea polinomului în x={1,35;1,99}
L8/9. Se dă tabelul următor ce reprezintă înregistrările efectuate în cuptorul de tratament termic pentru un set de roţi dinţate
supuse revenirii:
Timp [min] 1 5 12 25 30 34 60 120
Temperatură 630 575 506 420 300 287 115 64
[°C]
Să se interpoleze liniar, spline şi cubic pentru a se afla valorile temperaturii piesei pentru t1 = 20 minute, t2 = 50 minute şi t3 =
110 minute. Să se reprezinte grafic.
L8/10. Se dă tabelul următor care reprezintă valorile cotei unei suprafeţe aerodinamice(extradosul unei aripi) pentru
x,y cele două coordonate în plan. Să se interpoleze liniar, spline şi cubic pentru a determina cota corespunzătoare punctelor
(12, 25), (55, 32) şi (30, 45).
x [mm] 5 15 25 35 45 55 65
y [mm]
10 625 637 699 705 752 684 623
20 685 705 721 750 763 675 595
30 705 764 795 821 652 632 587
40 745 780 845 865 601 598 555
50 790 825 888 842 746 647 520
L8/11. Să se determine polinomul de interpolare Lagrange de gradul 2 pentru vectorul y={2,1,3,7} pe diviziunea
x=0/0,25/0,5/1 şi să se determine valoarea acelui polinom în x= 0,75.
L8/12. Experimentările efectuate pe un aparat ELITRON17, utilizând epruvete din 205Cr115 cu electrozi din WCo8
în formă de bare cu diametrul de 2,5 mm, la o tensiune de regim de 50 V şi tensiune de 10V au condus la obţinerea
următorelor rezultate:
Timp, Microduritatea Microduritatea Creşterea Cantitatea de
min înainte de după de microdurităţ material
durificare, durificare, ii, depusă pe
daN/mm2 daN/mm2 daN/mm2 piesă, g
1 854 1160 3.5
2 946 1375 3.7
3 946 1389 6.7
4 868 1390 8.8
5 829 1366 7.6
6 843 1640 9.9
7 837 1616 8.6
8 854 1583 13.0
9 844 1554 11.2
10 946 1669 10.1
Să se găsescă relaţiile matematice (modelele matematice) care reprezintă creşterea microdurităţii funcţie de timp şi respectiv
cantitatea de material depusă în funcţie de timp.
L8/13. Să se determine polinomul de interpolare pentru funcţia f:[-0,5;0,5], f(x)=sin(πx) pe o diviziune uniformă a
intervalului [-0,5;0,5]cu pasul k=0,05. Să se compare apoi valoarea acestui polinom cu valoarea funcţiei date în punctul 0,25.
Calcul Numeric Asistat de Calculator
S8/1.
h=[62 -43 0 45 0 0 -34 0 61 -96]
sol=roots(h)
polyval(h,sol(5))
Ecuaţia are doar o singură soluţie reală.
S8/2.
t=[1 0 0 0 0 0 0 -10 10 -5 -5 150]
roots(t)
ans = -1.46
S8/5.
» [C R]=deconv(h,g)
Câtul împărţirii polinoamelor este:
C =
Columns 1 through 4
2.81818181818182 -3.87603305785124 2.64274981217130 0.24357967351957
Columns 5 through 6
7.00747666888129 -14.64409096694949
Restul împărţirii polinoamelor este:
R =
1.0e+002 *
Columns 1 through 4
0 0.00000000000000 -0.00000000000000 0.00000000000000
Columns 5 through 8
0 0 3.33655353985835 0.13640461717096
Columns 9 through 10
4.53418693457352 -9.16069094149171
S8/6.
» polyval(h,10)
ans = 5.774496651400000e+010 5
» polyval(h,25)
ans = 2.299609369701790e+014
4
» polyval(h,-54.6)
ans = -2.707730144931591e+017
» polyval(g,-12) 3
ans = 430216
2
S8/7.
x=[0.1 0.3 0.4 0.6 0.7 0.8 0.9 1 1.2 2.22 1
5.6 7.8]
y=[0.7 0.43 0.2 0.02 -0.04 -0.09 -0.07 0 0
0.5 0.65 0.77 0.86]
Determinarea celei mai bune aproximări:
c1=polyfit(x,y,1) -1
c3=polyfit(x,y,3)
c5=polyfit(x,y,5) -2
y1=polyval(c1,x) 0 1 2 3 4 5 6 7 8
y3=polyval(c3,x)
y5=polyval(c5,x)
R1=sum((y-y1).^2)
R3=sum((y-y3).^2)
2
CNAC - LABORATOR 2
R5=sum((y-y5).^2)
R1 = 0.8568
R3 = 0.7509
R5 = 0.1269
Cea mai mică sumă este cea a polinomului de gradul 5.
S8/8.
x=[0.25 0.5 0.75 1 1.25 1.5 1.75 2]
y=[3 2.5 .5 1.5 0.5 -1 3 0]
xi=0.25:0.01:2
yy1=interp1(x,y,xi, 'linear')
yyc=interp1(x,y,xi,'cubic')
yysp=interp1(x,y,xi,'spline')
plot(x,y, 'ro',xi,yy1,'b',xi,yyc,'k',xi,yysp, 'g');grid on
xcalc=[1.35 1.99]
yy1=interp1(x,y,xcalc,'linear')
yyc=interp1(x,y,xcalc,'cubic') 4
yysp=interp1(x,y,xcalc,'spline')
3
Valorile obţinute sunt:
xcalc = 1.3500 1.9900
yy1 = -0.1000 0.1200 2
yyc = -0.3280 0.2544
yysp = -0.5380 0.4753
1
S8/9.
t=[1 5 12 25 30 34 60 120] 0
T=[630 575 506 420 300 287 115 64]
interp1(t,T,[20 50 110])
T20L=interp1(t,T,20,’linear’) -1
T20SP=interp1(t,T,20,’spline’)
T20C=interp1(t,T,20,’cubic’)
-2
0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
S8/10.
z=[ 625 637 699 705 752 684 623; 685
705 721 750 763 675 595;705 764 795 821 652 632 587;745 780 845 865 601 598 555;790 825 888 842
746 647 520]
[x y]=meshgrid(5:10:65, 10:10:50)
» interp2(x,y,z,12,25)
ans = 722.6500
» interp2(x,y,z,55,32)
ans = 625.2000
» interp2(x,y,z,30,45)
ans = 860