Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Polinoame' PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 3

LABORATOR 8

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

Indicaţii şi soluţii pentru Laboratorul 8:

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/3. f=[345 0 -4 45 0 0 -2 15 -56]


S8/4.
» conv(g,f)
ans =
Columns 1 through 6
7590 5175 -88 930 -18645 0
Columns 7 through 12
180 -2220 -1007 -840 112 -840
Column 13
3136
Rezultatul este polinomul următor:
7590·x12 +5175·x11– 88·x10 + 930·x9 –18645·x8 +180·x6 –2220·x5 –1007·x4 –840·x3 +112·x2 –840·x +3136

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.

Pentru trasarea graficului:


xi=0.1:0.1:7.8
y1i=polyval(c1,xi)
y3i=polyval(c3,xi)
y5i=polyval(c5,xi)
plot(x,y,'ro',xi,y1i,xi,y3i,xi,y5i)

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

You might also like