Добавил:
korayakov
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:Мат Моделирование 2002 (ЛисовЕЦ) / лабы / задания / lab_05 / kw_0804
.txt Протокол Курсовой работы
v0=1000;g=9.8;a=pi/4;T1=v0*sin(a)/g;T2=2*T1;L=v0^2/g;v0x=v0*cos(a);v0y=v0*sin(a);u0=2000;
» x0=.6*L;x1=.8*L;t=0:1:T2;
» x=v0x*t;y=v0y*t-g*t.^2/2;
» whos
Name Size Elements Bytes Density Complex
L 1 by 1 1 8 Full No
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
u0 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
x 1 by 145 145 1160 Full No
x0 1 by 1 1 8 Full No
x1 1 by 1 1 8 Full No
y 1 by 145 145 1160 Full No
Grand total is 446 elements using 3568 bytes
x=x';y=y';
» plot(x,y)
alfa=.005;xi=x(15:40);for i=1:26,yi(i)=y(i+14)*(1+alfa*randn);end,yi=yi';
size(xi)=26 1
hold
» plot(xi,yi,'r+')
» hold
Current plot released
» plot(x(15:40),y(15:40),xi,yi,'r+')
» AI=[xi.^2 xi ones(26,1)];
» AAI=AI'*AI;
» cond(AAI)
ans =
3.2945e+019
» A=[x.^2 x ones(length(x),1)];
» AA=A'*A;
cond(AA)
ans =
1.9021e+020
» koefi=inv(AAI)*(AI'*yi)
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 3.035223e-020
koefi =
0.0000
1.0016
-3.5444
» koefi(1)
ans =
-9.8771e-006
» koefi(2)
ans =
1.0016
» koefi(3)
ans =
-3.5444
koef=inv(AA)*(A'*y)
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 5.257195e-021
koef =
0.0000
1.0000
0.0000
» koef(1)
ans =
-9.8000e-006
» koef(2)
ans =
1.0000
» koef(3)
ans =
-1.7010e-011
» max(abs(koef-koefi))
ans =
3.5444
» max(abs(A*koef-A*koefi))
ans =
636.7089
»
» plot(x,A*koef,x,A*koefi,'r+')
»
» cond(AA)
ans =
1.9021e+020
size(t)
ans =
1 145
» t0=.6*145
t0 =
87
» t0=.8*145
t0 =
116
» L
L =
1.0204e+005
»
» LL=.6*L;
» LL-x(87)
ans =
413.3066
» .6*L-x(87);
» .6*L-x(87)
ans =
413.3066
» x(145)
ans =
1.0182e+005
» y(145)
ans =
216.9765
» x0
x0 =
6.1224e+004
» x1
x1 =
8.1633e+004
for i=1:26,xar(i)=x(87)+u0x*(t(87+i)-t(87));yar(i)=u0y*(t(i+87)-t(87))-g*(t(i+87)-t(87)).^2/2;end
plot(x,y,xar,yar,'*')
»
v0=1000;g=9.8;a=pi/4;T1=v0*sin(a)/g;T2=2*T1;L=v0^2/g;v0x=v0*cos(a);v0y=v0*sin(a);u0=2000;
» x0=.6*L;x1=.8*L;t=0:1:T2;
» x=v0x*t;y=v0y*t-g*t.^2/2;
» whos
Name Size Elements Bytes Density Complex
L 1 by 1 1 8 Full No
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
u0 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
x 1 by 145 145 1160 Full No
x0 1 by 1 1 8 Full No
x1 1 by 1 1 8 Full No
y 1 by 145 145 1160 Full No
Grand total is 446 elements using 3568 bytes
x=x';y=y';
» plot(x,y)
alfa=.005;xi=x(15:40);for i=1:26,yi(i)=y(i+14)*(1+alfa*randn);end,yi=yi';
size(xi)=26 1
hold
» plot(xi,yi,'r+')
» hold
Current plot released
» plot(x(15:40),y(15:40),xi,yi,'r+')
» AI=[xi.^2 xi ones(26,1)];
» AAI=AI'*AI;
» cond(AAI)
ans =
3.2945e+019
» A=[x.^2 x ones(length(x),1)];
» AA=A'*A;
cond(AA)
ans =
1.9021e+020
» koefi=inv(AAI)*(AI'*yi)
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 3.035223e-020
koefi =
0.0000
1.0016
-3.5444
» koefi(1)
ans =
-9.8771e-006
» koefi(2)
ans =
1.0016
» koefi(3)
ans =
-3.5444
koef=inv(AA)*(A'*y)
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 5.257195e-021
koef =
0.0000
1.0000
0.0000
» koef(1)
ans =
-9.8000e-006
» koef(2)
ans =
1.0000
» koef(3)
ans =
-1.7010e-011
» max(abs(koef-koefi))
ans =
3.5444
» max(abs(A*koef-A*koefi))
ans =
636.7089
»
» plot(x,A*koef,x,A*koefi,'r+')
»
» cond(AA)
ans =
1.9021e+020
size(t)
ans =
1 145
» t0=.6*145
t0 =
87
» t0=.8*145
t0 =
116
» L
L =
1.0204e+005
»
» LL=.6*L;
» LL-x(87)
ans =
413.3066
» .6*L-x(87);
» .6*L-x(87)
ans =
413.3066
» x(145)
ans =
1.0182e+005
» y(145)
ans =
216.9765
» x0
x0 =
6.1224e+004
» x1
x1 =
8.1633e+004
for i=1:26,xar(i)=x(87)+u0x*(t(87+i)-t(87));yar(i)=u0y*(t(i+87)-t(87))-g*(t(i+87)-t(87)).^2/2;end
plot(x,y,xar,yar,'*')
»
Соседние файлы в папке lab_05