Мат Моделирование 2002 (ЛисовЕЦ) / лабы / задания / lab_05 / kw_03

» whos
Name Size Elements Bytes Density Complex

T1 1 by 1 1 8 Full No
T2 1 by 1 1 8 Full No
a 1 by 1 1 8 Full No
g 1 by 1 1 8 Full No
v0 1 by 1 1 8 Full No
v0x 1 by 1 1 8 Full No
v0y 1 by 1 1 8 Full No

T2 =

144.3075

» t=0:1:T2;
» whos
Name Size Elements Bytes Density Complex

T1 1 by 1 1 8 Full No
T2 1 by 1 1 8 Full No
a 1 by 1 1 8 Full No
g 1 by 1 1 8 Full No
t 1 by 145 145 1160 Full No
v0 1 by 1 1 8 Full No
v0x 1 by 1 1 8 Full No
v0y 1 by 1 1 8 Full No
x 1 by 145 145 1160 Full No
y 1 by 145 145 1160 Full No

» plot(x,y,'.')
» for i=1:25, yy(i)=y(i+25);xx(i)=x(i+25);end
» whos
Name Size Elements Bytes Density Complex

T1 1 by 1 1 8 Full No
T2 1 by 1 1 8 Full No
a 1 by 1 1 8 Full No
g 1 by 1 1 8 Full No
i 1 by 1 1 8 Full No
t 1 by 145 145 1160 Full No
v0 1 by 1 1 8 Full No
v0x 1 by 1 1 8 Full No
v0y 1 by 1 1 8 Full No
x 1 by 145 145 1160 Full No
xx 1 by 25 25 200 Full No
y 1 by 145 145 1160 Full No
yy 1 by 25 25 200 Full No

» xx=xx';yy=yy';

» A=[xx.^2 xx ones(25,1)];AA=A'*A;cond(AA)

ans =

4.8026e+020

» koef=inv(AA)*A'*yy

koef =

0.0000
1.0000
0.0000

» yyy=[(x.^2 )' x' ones(length(x),1)]*koef;
» plot(x,y,'.',x,yyy,'r.')

» yye=yy.*(1+.1*randn(length(xx),1));
» koefe=inv(AA)*A'*yye

koefe =

1.0e+003 *

0.0000
0.0012
-2.0914

» max(abs(koefe-koef))

ans =

2.0914e+003

» yyye=[(x.^2 )' x' ones(length(x),1)]*koef;
» plot(x,y,'.',x,yyye,'r.')
» yyye=[(x.^2 )' x' ones(length(x),1)]*koefe;
» plot(x,y,'.',x,yyye,'r.')
» yye=yy.*(1+.01*randn(length(xx),1));
» koefe=inv(AA)*A'*yye

Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 2.082114e-021


koefe =

1.0e+003 *

0.0000
0.0009
1.2141

» max(abs(koefe-koef))

ans =

1.2141e+003

» yyye=[(x.^2 )' x' ones(length(x),1)]*koefe;
» plot(x,y,'.',x,yyye,'r.')
» yye=yy.*(1+.001*randn(length(xx),1));
» koefe=inv(AA)*A'*yye

Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 2.082114e-021


koefe =

0.0000
0.9985
18.1926

» yyye=[(x.^2 )' x' ones(length(x),1)]*koefe;
» plot(x,y,'.',x,yyye,'r.')

» max(abs(y'-yyye))

ans =

199.1117Ї kw
Ї t=0:1:T2;
Ї dop

ans =

-39.1572
-0.8428

Ї quit

1336 flops.