Tài liệu hạn chế xem trước, để xem đầy đủ mời bạn chọn Tải xuống
1
/ 24 trang
THÔNG TIN TÀI LIỆU
Thông tin cơ bản
Định dạng
Số trang
24
Dung lượng
143,09 KB
Nội dung
180
Chơng 11 : nộisuyvàxấpxỉhàm
Đ1.Nội suy Lagrange
Trong thực tế nhiều khi phải phục hồi một hàm y = f(x) tại mọi giá trị x trong một
đoạn [ a,b ] nào đó mà chỉ biết một số nhất định các giá trị của hàm tại một số điểm cho
trớc.Các giá trị này đợc cung cấp qua thực nghiệm hay tính toán.Vì vậy nảy sinh vấn đề
toán học là trên đoạn a x b cho một loạt các điểm x
i
( i= 0,1,2 ) và tại các điểm x
i
này
giá trị của hàm là y
i
= f(x
i
) đã biết.Bây giờ ta cần tìm đa thức :
P
n
(x) = a
o
x
n
+ a
1
x
n-1
+ +a
n-1
x
+ a
n
sao cho P
n
(x
i
) = f(x
i
) = y
i
.Đa thức P
n
(x) đợc gọi là đa thức nộisuy của hàm y = f(x).Ta chọn
đa thức để nộisuyhàm y = f(x) vì đa thức là loại hàm đơn giản,luôn có đạo hàmvà nguyên
hàm.Việc tính giá trị của nó theo thuật toán Horner cũng đơn giản.
Bây giờ ta xây dựng đa thức nộisuy kiểu Lagrange.Gọi L
i
là đa thức :
)xx) (xx)(xx) (xx(
)xx) (xx)(xx) (xx(
L
ni1ii1ii0i
n1i1i0
i
=
+
+
Rõ ràng là L
i
(x) là một đa thức bậc n và :
=
=
ij
0
ij1
)x(
L
j
i
Ta gọi đa thức này là đa thức Lagrange cơ bản.
Bây giờ ta xét biểu thức :
=
=
n
0i
ii
n
)x(L)x(f)x(
P
Ta thấy P
n
(x) là một đa thức bậc n vì các L
i
(x) là các đa thức bậc n và thoả mãn điều
kiện P
n
(x
i
) = f(x
i
) = y
i
.Ta gọi nó là đa thức nộisuy Lagrange.
Với n = 1 ta có bảng
x x
0
x
1
y y
0
y
1
Đa thức nộisuy sẽ là :
P
1
(x) = y
o
L
0
(x) + y
1
L
1
(x
1
)
10
1
0
xx
xx
L
=
01
0
1
xx
xx
L
=
nên
01
0
1
10
1
01
xx
xx
y
xx
xx
y)x(P
+
=
Nh vậy P
1
(x) là một đa thức bậc nhất đối với x
Với n = 2 ta có bảng
x x
0
x
1
x
2
y y
0
y
1
y
2
Đa thức nộisuy sẽ là :
P
2
(x) = y
o
L
0
(x) + y
1
L
1
(x
1
) + y
2
L
2
(x
2
)
)xx)(xx(
)xx)(xx(
L
2010
21
0
=
)xx)(xx(
)xx)(xx(
L
2101
20
1
=
181
)xx)(xx(
)xx)(xx(
L
1202
10
2
=
Nh vậy P
1
(x) là một đa thức bậc hai đối với x
Trên cơ sở thuật toán trên ta có chơng trình tìm đa thức nộisuy của một hàm khi
cho trớc các điểm và sau đó tính trị số của nó tại một giá trị nào đó nh sau :
Chơng trình 11-1
#include <conio.h>
#include <stdio.h>
#include <ctype.h>
#define max 21
int maxkq,n;
float x[max],y[max],a[max],xx[max],yy[max];
float x0,p0;
void main()
{
int i,k;
char ok ;
void vaosolieu(void);
float lagrange(int,float [],float [],float);
void inkq(void);
clrscr();
printf("%24cNOI SUY DA THUC LAGRANGE\n",' ');
vaosolieu();
k=0;
ok='c';
while (ok=='c')
{
printf("Tinh gia tri cua y voi x la x0 = ");
scanf("%f",&x0);
p0=lagrange(n,x,y,x0);
printf("Gia tri cua y = %15.5f\n",p0);
printf("\n");
k=k+1;
maxkq=k;
xx[k]=x0;
yy[k]=p0;
flushall();
printf("Tinh tiep khong(c/k)?");
scanf("%c",&ok);
}
inkq();
182
}
void vaosolieu()
{
int i,t;
char ok;
printf("\n");
printf("Ham y = f(x)\n");
printf("So cap (x,y) nhieu nhat la max = 20\n");
printf("So diem da cho truoc n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
printf("\n");
printf(" SO LIEU BAN VUA NHAP\n");
printf(" x y\n");
for (i=1;i<=n;i++)
printf("%8.4f %8.4f\n",x[i],y[i]);
ok=' ';
t=1;
flushall();
while (t)
{
printf("\nCo sua so lieu khong(c/k):?");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
}
float lagrange(int n,float x[max],float y[max],float x0)
{
int i,k;
183
float g0;
p0=0.0;
for (k=1;k<=n;k++)
{
g0=1.0;
for (i=1;i<=n;i++)
if (i!=k)
g0=g0*(x0-x[i])/(x[k]-x[i]);
p0=p0+y[k]*g0;
}
return(p0);
}
void inkq()
{
int i,j,k;
printf("\n");
printf("%24cBANG SO LIEU\n",' ');
printf("%18cx %24cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%20.4f %25.4f\n",x[i],y[i]);
printf("\n");
printf("%24cKET QUA TINH TOAN\n",' ');
printf("%14cx %10cy\n",' ',' ');
for (k=1;k<=maxkq;k++)
printf("%15.5f %15.5f\n",xx[k],yy[k]);
getch();
}
Giả sử ta có bảng các giá trị x,y :
x 0 3 -2 2 4
y 0 -3.75 10 -2 4
vậy theo chơng trình tại x = 2.5 y = -3.3549.
Đ2.Nội suy Newton
Bây giờ ta xét một cách khác để xây dựng đa thức nộisuy gọi là phơng pháp
Newton.Trớc hết ts đa vào một khái niệm mới là tỉ hiệu
Giả sử hàm y = y(x) có giá trị cho trong bảng sau :
x x
0
x
1
x
2
x
n-1
x
n
y y
0
y
1
y
2
y
n-1
y
n
Tỉ hiệu cấp 1 của y tại x
i
,x
j
là :
184
ji
ji
ji
xx
yy
]x,x[y
=
Tỉ hiệu cấp hai của y tại x
i
,x
j
,x
k
là :
ki
kjji
kji
xx
]x,x[y]x,x[y
]x,x,x[y
=
v.v.
Với y(x) = P
n
(x) là một đa thức bậc n thì tỉ hiệu cấp 1 tại x,x
0
:
0
0nn
0n
xx
)x(P)x(P
]x,x[P
=
là một đa thức bậc (n-1).Tỉ hiệu cấp 2 tại x,x
0
,x
1
:
1
10n0n
10n
xx
]x,x[P]x,x[P
]x,x,x[P
=
là một đa thức bậc (n-2) v.v và tới tỉ hiệu cấp (n+1) thì :
P
n
[ x,x
o
, ,x
n
] = 0
Từ các định nghĩa tỉ hiệu ta suy ra :
P
n
(x) = P
n
(x
0
) + ( x- x
0
)P
n
[x,x
o
]
P
n
[x,x
0
] = P
n
[x
0
,x
1
] + ( x- x
1
) P
n
[x,x
o
,x
1
]
P
n
[x,x
o
,x
1
] = P
n
[x
0
,x
1
,x
2
] + ( x- x
2
) P
n
[x,x
o
,x
1
,x
2
]
P
n
[x,x
o
, ,x
n-1
] = P
n
[x
0
,x
1
, ,x
n
] + ( x- x
n
) P
n
[x,x
o
, ,x
n
]
Do P
n
[ x,x
o
, ,x
n
] = 0 nên từ đó ta có :
P
n
(x) = P
n
(x
0
) + (x - x
0
)P
n
[x
o
,x
1
] + (x - x
0
)(x - x
1
)P
n
[x
0
,x
1
,x
2
] +
+(x - x
0
)(x - x
n-1
)P
n
[x
0
,,x
n
]
Nếu P
n
(x) là đa thức nộisuy của hàm y=f(x) thì :
P
n
(x
i
) = f(x
i
) = y
i
với i = 0 ữ n
Do đó các tỉ hiệu từ cấp 1 đến cấp n của P
n
và của y là trùng nhau và nh vậy ta có :
P
n
(x) = y
0
+ (x - x
0
)y[x
0
,x
1
] + (x - x
0
)(x - x
1
)y[x
0
,x
1
,x
2
] + +
(x - x
0
)(x - x
1
) (x - x
n-1
)y[x
0
, ,x
n
]
Đa thức này gọi là đa thức nộisuy Newton tiến xuất phát từ nút x
0
của hàm y =
f(x).Ngoài đa thức tiến còn có đa thức nộisuy Newton lùi xuất phát từ điểm x
n
có dạng nh
sau :
P
n
(x) = y
n
+ (x - x
n
)y[x
n
,x
n-1
] + (x - x
n
)(x - x
n-1
)y[x
n
,x
n-1
,x
n-2
] + +
(x - x
n
)(x - x
n-1
) (x - x
1
)y[x
n
, ,x
0
]
Trờng hợp các nút cách đều thì x
i
= x
0
+ih với i = 0,1, ,n.Ta gọi sai phân tiến cấp 1
tại i là :
y
i
= y
i+1
- y
i
và sai phân tiến cấp hai tại i :
2
y
i
= (y
i
) = y
i+2
- 2y
i+1
+ y
i
và sai phân tiến cấp n là :
n
y
i
= (
n-1
y
i
)
Khi đó ta có :
h
y
]
x,
x[y
0
1
0
=
h2
y
]
x,x,
x[y
2
0
2
21
0
=
185
)h!n(
y
]
x.,.,.x[y
n
0
n
n0
=
Bây giờ đặt x = x
0
+ ht trong đa thức Newton tiến ta đợc :
y
!n
)1nt.(.).1t(t
y
!2
)1t(t
yty
)htx(
P
0
n
0
2
00
0
n
+
++++=+
thì ta nhận đợc đa thức Newton tiến xuất phát từ x
0
trong trờng hợp nút cách đều.Với n =1
ta có :
P
1
(x
0
+ht) = y
0
+ y
0
Với n =2 ta có :
y
2
)1t(t
yty
)htx(
P
0
2
00
0
2
++=+
Một cách tơng tự ta có khái niệm các sai phân lùi tại i :
y
i
= y
i
- y
i-1
2
y
i
= (y
i
) = y
i
- 2y
i-1
+ y
i-2
n
y
i
= (
n-1
y
i
)
và đa thức nộisuy Newton lùi khi các điểm nộisuy cách đều :
y
!n
)1nt.(.).1t(t
y
!2
)1t(t
yty
)htx(
P
n
n
n
2
nn
0
n
+
+
+
++++=+
Ví dụ : Cho hàm nh bảng sau :
x 0.1 0.2 0.3 0.4
y 0.09983 0.19867 0.29552 0.38942
Ta tính giá trị của hàm tại 0.14 bằng đa thức nộisuy Newton vì các mốc cách đều
h = 0.1.Ta có bảng sai phân sau :
i x y
y
2
y
3
y
0 0.1 0.09983
0.09884
1 0.2 0.19867 -
0.00199
0.09685 -0.00096
2 0.3 0.29552 -
0.00295
0.09390
3 0.4 0.38942
Ta dùng công thức Newton tiến với điểm gốc là x
0
= 0.1.h = 0.1.Với x = 0.14 ta có
0.14 = 0.1 + 0.1t nên t = 0.4 và kết quả là :
1395433
6
.000096.0
!3
)2t)(1t(t
00199.0
!2
)1t(t
099884.0.t09983.0)t1.01.0(P =
++=+ Chơng
trình nộisuy Newton nh sau :
Chơng trình 11-2
186
//Noi suy Newton
#include <conio.h>
#include <stdio.h>
#include <ctype.h>
#define max 11
void main()
{
int i,j,k,n,t;
float a[max],b[max],x[max],y[max];
char ok;
float x0,p;
clrscr();
printf("So diem da cho n = ");
scanf("%d",&n);
for (i=1;i<=n;i++)
{
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
}
printf("%10cBANG SO LIEU\n",' ');
printf("%8cx%30cy\n",' ',' ');
for (i=1;i<=n;i++)
printf("%4c%8.4f%23c%8.4f\n",' ',x[i],' ',y[i]);
ok=' ';
t=0;
flushall();
while (t)
{
printf("Co sua so lieu khong(c/k): ");
scanf("%c",&ok);
if (toupper(ok)=='C')
{
printf("Chi so cua phan tu can sua i = ");
scanf("%d",&i);
printf("Gia tri moi : ");
printf("x[%d] = ",i);
scanf("%f",&x[i]);
printf("y[%d] = ",i);
scanf("%f",&y[i]);
flushall();
}
if (toupper(ok)!='C')
t=0;
}
a[1]=y[1];
for (j=1;j<=n-1;j++)
187
{
for (i=1;i<=n-j;i++)
y[i]=(y[i+1]-y[i])/(x[i+j]-x[i]);
a[j+1]=y[1];
}
b[n]=a[n];
for (k=n-1;k>=1;k )
{
for (j=n-1;j>=1;j )
b[j]=a[j] ;
for (i=n-1;i>=k;i )
a[i]=a[i]-b[i+1]*x[k];
}
for (i=n;i>=1;i )
printf("He so bac %d la :%8.4f\n",i-1,a[i]);
printf("\n");
k=0;
ok='c';
flushall();
while (ok=='c')
{
printf("Tinh gia tri cua y tai x = ");
scanf("%f",&x0);
p=0;
for (k=n;k>=1;k )
p=p*x0+a[k];
printf("Tri so noisuy tai x0 = %4.2f la : %10.5f\n",x0,p);
getch();
printf("Ban co muon tinh tiep cac diem khac khong(c/k)");
do
scanf("%c",&ok);
while ((ok!='c')&&(ok!='k'));
}
}
Dùng chơng trình này nộisuy các giá trị cho trong bảng sau
0 0.2 0.4 0.6 0.8 1.0
1 1.2214027
6
1.4918247 1.8221188 2.2255409
3
2.7182818
3
ta có các hệ số của đa thức nộisuy : 0.0139(bậc 5),0.0349(bậc 4),0.1704(bậc3),0.4991(bậc
2),1.0001(bậc 1) và 1.0000(bậc 0).
Đ3.Nội suy Aitken
Một dạng khác của đa thức nộisuyđợc xác định bằng thuật toán Aitken.Giả sử ta có
n điểm đã cho của hàm f(x).Nh vậy qua hai điểm x
0
và x
1
ta có đa thức nộisuy Lagrange
của hàm f(x) đợc viết dới dạng :
188
01
11
00
01
xx
xxy
xxy
)x(P
=
là một đa thức bậc 1 :
01
0
1
10
1
001
xx
xx
y
xx
xx
y)x(P
+
=
.Khi x = x
0
thì :
0
01
011
000
001
y
xx
xxy
xxy
)x(P =
=
Khi x = x
1
thì :
1
01
111
100
101
y
xx
xxy
xxy
)x(P =
=
Đa thức nộisuy Lagrange của f(x) qua 3 điểm x
0
,x
1
,x
2
có dạng :
02
212
001
012
xx
xx)x(P
xx)x(P
)x(P
=
và là một đa thức bậc 2:
)xx)(xx(
)xx)(xx(
y
)xx)(xx(
)xx)(xx(
y
)xx)(xx(
)xx)(xx(
y)x(P
1202
10
2
2101
20
1
2010
21
0012
+
+
=
Khi x = x
0
thì :
0
02
0212
000
0012
y
xx
xx)x(P
xxy
)x(P =
=
Khi x = x
1
thì :
1
02
121
101
1012
y
xx
xxy
xxy
)x(P =
=
Khi x = x
2
thì :
2
02
222
20201
2012
y
xx
xxy
xx)x(P
)x(P =
=
Tổng quát đa thức nộisuy Lagrange qua n điểm là :
02
nn 12
0)1n (01
n 012
xx
xx)x(P
xx)x(P
)x(P
=
Nh vậy ta có thể dùng phép lặp để xác định lần lợt các đa thức Lagrange.Sơ đồ tính
toán nh vậy gọi là sơ đồ Neville-Aitken.
Ví dụ : Cho các cặp điểm (0,0.4),(1.4,1.5),(2.6,1.8),(3.9,2.6),tính y tại x=2
189
97143.1
04.1
6.05.1
24.0
xx
xxy
xxy
)x(P
01
11
00
01
=
−
−
−
=
−
−
−
=
65.1
4.16.2
6.08.1
6.05.1
xx
xxy
xxy
)x(P
12
22
11
12
=
−
−
=
−
−
−
=
7242.1
06.2
6.065.1
297143.1
xx
xx)x(P
xx)x(P
)x(P
02
212
001
012
=
−
−
=
−
−
−
=
4308.1
6.29.3
9.16.2
6.08.1
xx
xxy
xxy
)x(P
23
33
22
23
=
−
=
−
−
−
=
5974.1
4.19.3
9.14308.1
6.065.1
xx
xx)x(P
xx)x(P
)x(P
13
323
112
123
=
−
−
=
−
−
−
=
6592.1
09.3
9.15974.1
27242.1
xx
xx)x(P
xx)x(P
)x(P
03
3123
0012
0123
=
−
−
=
−
−
−
=
Ch−¬ng tr×nh ®−îc viÕt nh− sau
Ch−¬ng tr×nh 11-3
//Noi suy Aitken
#include <conio.h>
#include <stdio.h>
#include <ctype.h>
#define max 11
void main()
{
float x[max],y[max],yd[max];
float x1;
int j,k,n,n1;
clrscr();
printf("Cho so diem da co n = ");
scanf("%d",&n1);
n=n1-1 ;
for (k=0;k<=n;k++)
{
printf("x[%d] = ",k+1);
scanf("%f",&x[k]);
printf("y[%d] = ",k+1);
scanf("%f",&y[k]);
}
printf("Cho diem can tinh gia tri cua ham x1 = ");
[...]... 8,4 9 9,1 10 9,4 11 9 ,5 12 9 ,5 13 9,4 ta có n = 7 và chọn m = 2 và tính đợc theo chơng trình các hệ số : a0 = -0 .1119 05 ; a1 = 2 .54 5238 ; a2 = -4 . 857 143 vàhàmxấpxỉ sẽ là : f(x) = -0 .1119 05 + 2 .54 5238x -4 . 857 143x2 2 .Hàm dạng Aecx : Khi các số liệu thể hiện một sự biến đổi đơn điệu ta dùng hàmxấpxỉ là y = Aecx.Lấy logarit hai vế ta có : lny = lnA + cxlne 193 Theo điều kiện đạo hàm S = 0 ta có hệ phơng...scanf("%f",&x1); for (k=0;k=0;j ) yd[j]=(yd[j]*(x1-x[k+1])-yd[j+1]*(x1-x[j]))/(x[j]-x[k+1]); } printf("Gia tri ham tai x = %6.3f la y = %8.4f\n",x1,yd[0]); getch(); } Dùng chơng trình này để nộisuy các cặp số (1,3),(2 ,5) ,(3,7),(4,9) và (5, 11) tại x = 2 .5 ta có y = 6 Đ4 .Xấp xỉhàm bằng phơng pháp bình phơng... trình các hệ số : A = 7.1641 và q = 1. 953 1 vàhàmxấpxỉ sẽ là : f(x) = 12 85. 44x1. 953 1 4 .Hàm lợng giác : Khi quan hệ y=f(x) có dạng tuần hoàn ta dùng hàm xấp xỉ là tổ hợp tuyến tính của các hàm sin và cosin dạng : n n i =1 i =1 f (x) = a 0 + a i cos(ix) + b i sin(ix) Để đơn giản trớc hết ta xét hàm chỉ có một số hạng sin-cos,nghĩa là : f (x) = a 0 + a 1 cos x + b1 sin x Hàm S sẽ có dạng : n S = [y... 10,chu kì T = 15 ta nhận đợc kết quả tính a0 = 1.7 ; a1 = 0 .5 ; b1 = -0 .8661 và = 4.18879.Nh vậy hàm xấp xỉ có dạng : f(x) = 1.7 + 0.5cos(4.18879x) - 0.8661sin(4.18879x) 5.Hàm hữu tỉ : Khi quan hệ y = f(x) có dạng đờng cong bão hoà hay dạng arctan,tan v.v ta dùng hàm xấp xỉ là hàm hữu tỉ dạng đơn giản : y= ax b+x Lấy nghịch đảo của nó ta có : 1 b1 1 = + y ax a Đặt 1/y = Y,1/x = X,b/a = B và 1/a = A phơng... trớc ta đã nộisuy giá trị của hàm. Bài toán đó là cho một hàm dới dạng bảng số và phải tìm giá trị của hàm tại một giá trị của đối số không nằm trong bảng Trong thực tế,bên cạnh bài toán nộisuy ta còn gặp một dạng bài toán khác.Đó là tìm công thức thực nghiệm của một hàm .Nội dung bài toán là từ một loạt các điểm cho trớc (có thể là các giá trị của một phép đo nào đó) ta phải tìm một hàm xấp xỉ các giá... bảng x y 0 128 0 2 6 35 4 324 6 162 8 76 10 43 12 19 ta có n = 7 và tính đợc theo chơng trình các hệ số : A = 12 85. 44 va c = -0 .3476 vàhàm xấp xỉ sẽ là : f(x) = 12 85. 44 3 .Hàm dạng Axq : Khi các số liệu thể hiện một sự biến đổi đơn điệu ta cũng có thể dùng hàmxấpxỉ là y = Axq.Lấy logarit hai vế ta có : lny = lnA + qlnx Theo điều kiện đạo hàm triệt tiêu ta có hệ phơng trình : n n q ln x i + n ln A... SO GOC OMEGA = %10.5f\n",omg); 200 printf("HE SO HANG\n"); printf("a[0] = %8.4f\n",a[0]); printf("CAC HE SO BAC CAO\n"); printf("%5ccos%25csin\n",' ',' '); for (i=1;i . theo chơng trình các hệ số : a 0 = -0 .1119 05 ; a 1 = 2 .54 5238 ; a 2 = -4 . 857 143 và hàm xấp xỉ sẽ là : f(x) = -0 .1119 05 + 2 .54 5238x -4 . 857 143x 2 2 .Hàm dạng Ae cx : Khi các số liệu thể. printf("% 15. 5f % 15. 5f
",xx[k],yy[k]); getch(); } Giả sử ta có bảng các giá trị x,y : x 0 3 -2 2 4 y 0 -3 . 75 10 -2 4 vậy theo chơng trình tại x = 2 .5 y = -3 . 354 9. Đ2 .Nội suy. này để nội suy các cặp số (1,3),(2 ,5) ,(3,7),(4,9) và (5, 11) tại x = 2 .5 ta có y = 6. Đ4 .Xấp xỉ hàm bằng phơng pháp bình phơng bé nhất Trong các mục trớc ta đã nội suy giá trị của hàm. Bài