Добавил:
korayakov
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Мат Моделирование 2002 (ЛисовЕЦ) / лабы / labs2004 / Задания 2004 / LAB_05 / KW_03
.TXT » 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.
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.
Соседние файлы в папке LAB_05