2) Результат работы программы

0

0

-0,30945

0,073447

-0,00543

0,309451095037108

0,0734474907336602

0,00543335786926846

0,2

0,201336002541094

0,037356

0,190515

0,206291

0,163980249619322

0,0108208153110148

0,00495535440957337

0,4

0,410752325802816

0,384163

0,358636

0,416482

0,0265897249221636

0,0521162973068814

0,00572965833526906

0,6

0,636653582148241

0,730969

0,57781

0,638285

0,0943158666912908

0,0588435676170165

0,00163174964523061

0,8

0,888105982187623

1,077776

0,848037

0,884848

0,189670314610789

0,0400688368516718

0,00325777417030526

1

1,17520119364380

1,424583

1,169317

1,169317

0,249381951113491

0,00588377273368801

0,00588377273369023

1,2

1,50946135541217

1,77139

1,541651

1,50484

0,261928637304000

0,0321894858415388

0,00462157683983144

1,4

1,90430150145153

2,118197

1,965037

1,904562

0,213895339223519

0,0607359049152112

0,000260587652961775

1,6

2,37556795320023

2,465004

2,439477

2,381631

0,0894357354337023

0,0639091630489852

0,00606320740683453

1,8

2,94217428809568

2,811811

2,96497

2,949194

0,130363751502868

0,0227956828054405

0,00701951308485560

2

3,62686040784702

3,158617

3,541516

3,620397

0,468243023295327

0,0853444375245576

0,00646358892162180

= 0.589020998976464

= 0.0299431387393548

=0.00028404128889588



3) Текст программы

clc

clear all

hold on

format long g;

a=0;

b=3;

n=10;

x=linspace(a,b,n+1);

y=exp(2*x)-cos(2*x);

x=x';

y=y';

A=[x,y];

plot(x,y,'o');

for s=1:3

for p=0:s

for q=0:s

C(p+1)=0;

B(p+1,q+1)=0;

end;

end;

for p=0:s

for i=0:10

C(p+1)=C(p+1)+(((x(i+1))^(p))*y(i+1));

end;

end;

for p=0:s

for q=0:s

for i=0:10

B(p+1,q+1)=B(p+1,q+1)+((x(i+1))^(p+q));

end;

end;

end;

D=B;

for i=1:s+1

D(i,s+2)=C(i);

end;

N=size(D,1);

for i=1:N

amax=D(i,i);

imax=i;

for j=(i+1):N

if abs(D(j,i))>abs(amax)

amax=D(j,i);

imax=j;

end;

end;

for j=i:(N+1)

h=D(i,j);

D(i,j)=D(imax,j);

D(imax,j)=h;

end;

h=D(i,i);

for j=i:(N+1)

D(i,j)=D(i,j)/h;

end;

for k=i+1:N

h=D(k,i);

for j=i:N+1

D(k,j)=D(k,j)-D(i,j)*h;

end;

end;

end;

v(N)=D(N,N+1);

for i=1:(N-1)

i1=N-i;

sum=0;

for j=(i1+1):N

sum=sum+v(j)*D(i1,j);

end;

v(i1)=D(i1,N+1)-sum;

end;

k='g-.';

if s<3

v(4)=0;

k='m-';

end;

if s<2

v(3)=0;

k='r:';

end;

f=v(1)+v(2)*x+v(3)*x.^2+v(4)*x.^3;

pr1=f'

for i=1:11

eps(i)=abs(y(i)-f(i));

end

pr2=eps'

sumeps=0;

for i=1:11

sumeps=sumeps+(eps(i)).^2;

end;

sumeps

plot(x,f,k)

B=0;

C=0;

D=0;

end;

legend ('Рассчитанные данные', 'linear', 'parabolic', 'cubic')




Случайные файлы

Файл
77550-1.rtf
13109.doc
131923.rtf
94661.rtf
11.doc




Чтобы не видеть здесь видео-рекламу достаточно стать зарегистрированным пользователем.
Чтобы не видеть никакую рекламу на сайте, нужно стать VIP-пользователем.
Это можно сделать совершенно бесплатно. Читайте подробности тут.