信息与计算张椿木.doc

文档引见:

传达与计算张椿木.doc数值辨析实验公报

全体教职员:算学与计算理科

专业:传达与计算理科

一级:学分一级

学号: 2012750217

姓名: 张椿木

训诫者:曹雪叶

2014年6月14日

湘潭综合性大学算学与计算理科全体教职员

先生实验公报

专业

算学与计算理科

很好的

1

姓名

张椿木

学号

2012750217

实验日期

2014 年 3 月 20 日至 2014 年 6 月 20 日

实验项主语

齐式插值说得中肯振荡

实验愿意的及请求允许

实验愿意的:

思索交替[-1],等容的1分],春分是

x(i)=-1+2i/n,i=0,1,2…,n

zeragrangian插值齐式i

L(x)=∑L(i)(x)/(1 25x(j)^2 ) i=0,1,…n

带着l(i)(x), i=0,1,…n,n是n次添写的基功能

实验请求允许:

(1)选择递加的点数n,3…,画f(x)与插值齐式功能

L(x)在[-1,1]上的图象,对实验结实举行了比拟和辨析

(2)选择等等功用,比如,在区间[-5中精确地解释,功用5]

h(x)=x/(1+x^4) , g(x)=arctanx

反复在上文中实验,看一眼结实方法

(3)交替[a],b]上切比雪夫点精确地解释为

xk=(b+a)/2+((b-a)/2)cos((2k-1)π/(2(n+1))),k=1,2,^,n+1,

带x1,x 2^x(n 1)解释上述的添写齐式

式,比拟结实。

注:实验结实、附参考文献及相符合次

指导教授的提议

成就: 署名:

年月日

凡例

注:实验结实、附参考文献及相符合次

实验课题与设计

插补

实验快速地流动:

次:

齐式插值说得中肯振荡(实验)

for m=1:6

子一批(2,3,m) 将窗口陷于2*3胶料的窗口

largrang(6*m) 运转largrang功能

if m==1

title(”longn=6”)

elseif m==2

title(”longn=12”)

elseif m==3

title(”longn=18”)

elseif m==4

title(”longn=24”)

elseif m==5

title(”longn=30”)

elseif m==6

title(”longn=36”)

end 每个窗口的插值点数

end

另存为:chazhi.m.

function largrang(longn)

mm=input(”please input mm(要运转的功能数是输出mm的音量):mm='

if mm==1 d表现域的局限于值

d=1;

elseif mm==2||mm==3

d=5;

end

x0=linspace(-d,d,朗恩) x的包装

if mm==1

y0=1./(1.+25.*x0.^2);

elseif mm==2

y0=x0./(1.+x0.^4);

elseif mm==3

y0=atan(x0);

end

x=sym(”x”);n=length(x0); s=0.0;

for k=1:n

p=1.0;

for j=1:n

if j~=k

p=p*(x-x0(j))/(x0(k)-x0(j));

end

end

s=p*y0(k)+s;

end

y=s;

if mm==1

ezplot(”1/(1+25*x^2)”)

elseif mm==2

ezplot(”x/(1+x^4)”)

elseif mm==3

Ezplot('atan(x)')

end

hold on

ezplot(y,[-d,d])

hold off

使守恒为:

数值实验结实与辨析

为第一点钟功能f(x)=1/(1 25×2

为另外的个功能h(x)=x/(1 x4

为第三个功能g(x)=arctan(x

议论:

由三个功能获得的largrang插值齐式及其运算,获取功能图像,它暗示一点钟功能的中心点指责尽量多,但功能两端的中心点越多,拉格朗格插值齐式不近似额使移近功能,相反,它离功用更远,功能两端的悄悄溜走更骗子,argrang插值齐式不收敛于功能

实验总结

功能添写齐式成绩的实验看重,尽管不愿意结实是笔误的,即使补充中心点的音量并做很多实验,功能的拉格朗插值齐式的普通规则是,中心点增大时,拉格朗格插值齐式不收敛于功能的两端,不再近似,它离得更远。,剧跳性更强.因而为功能的largrang插值齐式成绩可以用MATLAB来举行成绩的辨析,获得了正确的实验结实

湘潭综合性大学算学与计算理科全体教职员

先生实验公报

专业

算学与计算理科

很好的

1

姓名

张椿木

学号

2012750217

实验日期

2014 年 3 月 20 日至 2014 年 6 月 20 日

实验项主语

线性方程组的导演求解过程

实验愿意的及请求允许

实验:思索线性方程

编辑可以自动行为选择的主部件,选择prin解线性方程组的高斯被取消的事物法
实验请求允许:
(1)矩阵选择,和方程有一点钟解。取n=10计算必要的数,结实方法?
(2)现选择次中手工操作拔取主元的功用.每步被取消的事物快速地流动总拔取按模最小或按模尽量小的元素作认为优先元,值夜并记载计算结实。以防最大的元素,结实方法?辨析实验结实
(3)取矩阵阶数n=20或更大,反复上述的实验快速地流动,值夜、记载、辨析确切的成绩及计算的不符合,阐明在开除PR中选择次要元素的功能
(4)选择您感兴趣的等等成绩或生产随机矩阵,计算必要的的音量。反复实验,对实验结实举行了值夜、记载和辨析
注:实验结实、附参考文献及相符合次
指导教授的提议
成就: 署名:
年月日
凡例
注:实验结实、附参考文献及相符合次
另外的实验课题与设计
线性方程组的导演求解过程
次:
M-记录的成立
function x=gauss(n,r)
n=input('请输出矩阵a的次:n='
A=diag(6*ones(1,n)) 判断(1个,n-1页),1) 判断(8*1,n-1页),-1)
b=A*ones(n,1)
P=INPUT('与准则对应的必要的号为P准则:P='
pp=cond(A,p)
pause
[m,N]=胶料(A)
nb=n+1;Ab=[A b]
R=INPUT('请输出如果手工操作,手工操作输出1,自动行为
愿意的转载自桃豆。请转位本源

发表评论

电子邮件地址不会被公开。 必填项已用*标注