![]()
|
|||||||
Код программы в среде 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; 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) Результат программы
|
|||||||
|