Код программы в среде MATLAB. Результат программы

⇐ ПредыдущаяСтр 2 из 2

Код программы в среде MATLAB

Файл Vandermond. m

function z=Vandermond(x)

M=length(x);

z=ones(M, M);

for i=1: M

for j=1: M

z(i, j)=x(i). ^(j-1);

end;

Файл pol. m

function z=pol(a, x1)

M=length(a);

s=0;

for i=1: M

s=s+a(i)*x1. ^(i-1);

end;

z=s;

> > x=[0. 43; 0. 48; 0. 55; 0. 62; 0. 70; 0. 75]

x =

0. 4300

0. 4800

0. 5500

0. 6200

0. 7000

0. 7500

> > y=[1. 6359; 1. 7323; 1. 8768; 2. 0304; 2. 2284; 2. 3597]

y =

1. 6359

1. 7323

1. 8768

2. 0304

2. 2284

2. 3597

> > M=Vandermond(x)

M =

1. 0000 0. 4300 0. 1849 0. 0795 0. 0342 0. 0147

1. 0000 0. 4800 0. 2304 0. 1106 0. 0531 0. 0255

1. 0000 0. 5500 0. 3025 0. 1664 0. 0915 0. 0503

1. 0000 0. 6200 0. 3844 0. 2383 0. 1478 0. 0916

1. 0000 0. 7000 0. 4900 0. 3430 0. 2401 0. 1681

1. 0000 0. 7500 0. 5625 0. 4219 0. 3164 0. 2373

> > a=M^-1*y

a =

10. 0997

-80. 6863

291. 7494

-511. 6367

444. 9904

-152. 9063

> > x=0. 53;

> > y1=pol(a, x)

y1 =

1. 8349

> > x1=0. 43: 0. 01: 0. 75;

> > y1=pol(a, x1);

> > plot(x1, y1)

Результат программы

⇐ Предыдущая 12