Мат Моделирование 2002 (ЛисовЕЦ) / лабы / labs2004 / Задания 2004 / LAB_05 / KW_0804

Протокол Курсовой работы

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,'*')
»