Functions |
1 | figure () |
| plot (x0(1), x0(2),'rx') |
| plot (x1(1), x1(2),'rx') plot(x2(1) |
| x2 (2) |
rx | plot (x3(1), x3(2),'rx') plot(x4(1) |
rx | x4 (2) |
rx rx | axis ([0 11.5 0 7]) xlabel('x') = A0(1)*t.^3 + B0(1)*t.^2 + C0(1)*t + D0(1) |
| ylabel ('y') |
| title ('Kubicky spline z Matlabu') |
| plot (X2(i), Y2(i)) |
| plot (X3(i), Y3(i)) |
| plot (X4(i), Y4(i)) |
stejne jako v prikladu | konst () |
axes, cmd, par | coordtst () |
vykresleni splinu z naseho
interpolatoru Coordmv | figure (2) plot(x0(1) |
vykresleni splinu z naseho
interpolatoru Coordmv | x0 (2) |
| title ('Kubicky spline z Coordmv') |
| plot (pos(:, 1)/konst, pos(:, 2)/konst) |
Variables |
clear | all |
| clc |
zadani uzlovych bodu | x0 = [0,0]' |
| x1 = [3,4]' |
| x2 = [4,2]' |
| x3 = [9,5]' |
| x4 = [10.6,6]' |
vypocet tecnych vektoru v
krajnich uzlovych bodech x0 a
x4 pomoci hyperboly | v0 = 1/2*( - x2 + 4*x1 - 3*x0) |
| v4 = -1/2*( - x2 + 4*x3 - 3*x4) |
vypocet tecnych vektoru ve
vnitrnich uzlovych bodech x2 a
x3 | M = [ 4 1 0 |
| Q = [ [-3*x0 + 3*x2 - v0]' |
| v |
| v2 = v(2,:)' |
| v3 = v(3,:)' |
hold | on |
| t = 0:0.005:1 |
polynom | X1 = A0(1)*t.^3 + B0(1)*t.^2 + C0(1)*t + D0(1) |
| Y1 = A0(2)*t.^3 + B0(2)*t.^2 + C0(2)*t + D0(2) |
polynom | X2 = A1(1)*t.^3 + B1(1)*t.^2 + C1(1)*t + D1(1) |
| Y2 = A1(2)*t.^3 + B1(2)*t.^2 + C1(2)*t + D1(2) |
polynom | X3 = A2(1)*t.^3 + B2(1)*t.^2 + C2(1)*t + D2(1) |
| Y3 = A2(2)*t.^3 + B2(2)*t.^2 + C2(2)*t + D2(2) |
polynom | X4 = A3(1)*t.^3 + B3(1)*t.^2 + C3(1)*t + D3(1) |
| Y4 = A3(2)*t.^3 + B3(2)*t.^2 + C3(2)*t + D3(2) |
for | i |
vykresleni splinu z naseho
interpolatoru Coordmv | rx |