計(jì)算方法實(shí)驗(yàn)報(bào)告三

2023-03-27 20:51 作者:啊啊啊不不不b 0人讀過(guò) | 我要投稿

一、實(shí)驗(yàn)名稱(chēng)：用Chebyshev點(diǎn)插值

?

二、實(shí)驗(yàn)?zāi)康模?/span>

學(xué)會(huì)運(yùn)用chebyshev點(diǎn)進(jìn)行插值

三、實(shí)驗(yàn)內(nèi)容及要求：

$%E5%AF%B9%E5%87%BD%E6%95%B0f(x)%3D1%2F(1%2B25x%5E2%20)%2Cx%E2%88%88%5B-1%2C1%5D%0A%E5%88%A9%E7%94%A8chebyshev%E7%82%B9%E8%BF%9B%E8%A1%8C%E6%8F%92%E5%80%BC%EF%BC%8C%E6%9E%84%E9%80%A0%E6%8F%92%E5%80%BC%E5%A4%9A%E9%A1%B9%E5%BC%8F%0A%E5%8F%96n%3D3%2C8%2C16%2C20%0A$

四、實(shí)驗(yàn)步驟（或記錄）

clear;
? n=[4,8,16,20];
? for i=1:4
? ??? if i==1
? ??? for j=1:n(i)
? ??????? x1(j)=cos((2*j-1)*pi/(2*n(i)));
? ??? end
? ??? y1=1./(1+25.*x1.^2);
? ??? elseif i==2
? ??????? for j=1:n(i)
? ??????? x2(j)=cos((2*j-1)*pi/(2*n(i)));
? ??? end
? ??????? y2=1./(1+25.*x2.^2);
? ??? elseif i==3
? ??????? for j=1:n(i)
? ??????? x3(j)=cos((2*j-1)*pi/(2*n(i)));
? ??? end
? ??????? y3=1./(1+25.*x3.^2);
? ??? elseif i==4
? ??????? for j=1:n(i)
? ??????? x4(j)=cos((2*j-1)*pi/(2*n(i)));
? ??? end
? ??????? y4=1./(1+25.*x4.^2);
? ??? end
? end
? x5=-1:0.005:1;
? y5=1./(1+25.*x5.^2);
? plot(x1,y1,'g',x2,y2,'b',x3,y3,'r',x4,y4,'b',x5,y5,'y')