.m文件

function yy=lagrange(x1,y1,xx)

%本程序為Lagrange1插值,其中x1,y1

%為插值節(jié)點和節(jié)點上的函數(shù)值,輸出為插值點xx的函數(shù)值,

%xx可以是向量.

syms x

n=length(x1);

for i=1:n

t=x1;t(i)=[];L(i)=prod((x-t)./(x1(i)-t));% L向量用來存放插值基函數(shù)

end

u=sum(L.*y1);

p=simplify(u) % p是簡化后的Lagrange插值函數(shù)(字符串)

yy=subs(p,x,xx);

clf

plot(x1,y1,'ro',xx,yy,'*')

x=[ 0.4 0.5 0.6 0.7 0.8];

y=[-0.916291;-0.693147;-0.510826;-0.356675;-0.223144]';

yy=lagrange(x,y,0.54)

p =

- (14363668061545223*x^4)/6755399441055744 + (229230406283396627*x^3)/33776997205278720 - (6086876668119665137*x^2)/675539944105574400 + (23595121244981107513*x)/3377699720527872000 - 186390055565518223/70368744177664000

yy =

-0.6161

再問： 我們老師給出了部分程序代碼，要把這個補全，能否把這個補全呢？ 不要你的那個啊 怎么那么復(fù)雜？ x0=[o.4;0.5;0.6;0.7;0.8]; y0=[-0.916291;-0.693147;-0.510826;-0.356675;-0.223144]; a=0; i=1; j=1; for i=1:n p=1; n=5; for j=1:n if j~=i; q=(x-x0(i))/(x0(j)-xo(i)); end p=p*q; end a=a+p; end

再答： x0=[0.4;0.5;0.6;0.7;0.8]; y0=[-0.916291;-0.693147;-0.510826;-0.356675;-0.223144]; a=0; x=0.54 i=1; j=1; n=5; for i=1:n p=1; q=1; for j=1:n if j~=i; q=(x-x0(i))/(x0(j)-x0(i)); end p=p*q; end a=a+p*y0(i); end a x = 0.5400 >> a a = -0.6391

再問： 運行了老說 ？？？undefined fuction or variable 'untitled1' 運行不了。。。

總結(jié)

以上是生活随笔為你收集整理的matlab插值法计算根号,怎么用matlab利用拉格朗日插值计算法的原理编写并计算函数所在节点的近似值....的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。