信息与计算张椿木.doc

文档绍介:

通讯与计算张椿木.doc数值辨析试验传闻

算学与电脑学科学院

专业:通讯与计算学科

出色的:出色的函件

学号: 2012750217

姓名: 张椿木

指导教授:曹雪念

2014年6月14日

湘潭大学人员算学与电脑学科学院

先生试验传闻

专业

算学与计算学科

出色的

1

姓名

张椿木

学号

2012750217

试验日期

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

试验以协议约束

Polynomial interpolation of the shock phenomenon

试验物质及销路

试验物质:

思索区间[ 1,1 ]等容重新分配,分点为:

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

Ze Lagrange插值同次多项式

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

里面,L(我)(X), i=0,1,…n,n是n次插补功用。

试验销路:

对n = 2点越来越多的选择,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,X2 ^ x(n 1)插补混合物建筑物的福

式,喻为的终于。

注:试验终于、商议和相符合的顺序衔接

教员的暗示

成果: 署名:

年月日

凡例

注:试验终于、商议和相符合的顺序衔接

试验课题设计

插补

试验换异:

顺序:

Polynomial interpolation of the shock phenomenon(试验)

for m=1:6

subplot(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(运转输出几个的毫米水银柱高的几种功用):毫米水银柱高=)

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功用插值同次多项式收敛。

试验总结

插补同次多项式功用的MATLAB试验成绩,仍然他们的终于是认不出的。,但加法运算的总共供养多种试验,插补同次多项式功用的普通圣职授任可以找到,当使绕枢轴旋转加法运算,largrang插值同次多项式收敛功用的两端,归咎于更将近,但此外远离,紧张不安性更强.因而关闭功用的largrang插值同次多项式成绩可以依靠MATLAB来举行成绩的辨析,让法度更精确。

湘潭大学人员算学与电脑学科学院

先生试验传闻

专业

算学与计算学科

出色的

1

姓名

张椿木

学号

2012750217

试验日期

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

试验以协议约束

求求解线性方程组的直线方式

试验物质及销路

试验物质:思索线性方程组

对首要元件的自动行为选择制剂,可以人工操作选择求解线性方程的高斯删去换异
试验销路:
(1)矩阵,则方程有解。对n = 10计算的条件数,终于方式?
(2)现选择顺序中人工操作拔取主元的功用.每步删去换异总拔取按模最小或按模尽量小的元素作认为优先元,守候并记载每一步的终于。假使全体数量距离换异,试验终于与辨析。终于呢?
(3)矩阵阶数为n=20或更多。,反复试验的换异,守候、记载并辨析差别计算和距离,解说在距离PR首要元素的选择说话中肯功能
(4)选择你感兴趣的停止成绩或随机制造的矩阵,计算它的条件数。反复上述的试验,对试验终于举行记载和辨析。
注:试验终于、商议和相符合的顺序衔接
教员的暗示
成果: 署名:
年月日
凡例
注:试验终于、商议和相符合的顺序衔接
两试验课题设计
求求解线性方程组的直线方式
顺序:
M寄给报社的使被安排好
function x=gauss(n,r)
N =输出(请输出矩阵的阶数:n = ')
A=diag(6*ones(1,n)) (的诊断法(1,N-1),1) (8×的诊断法(1,N-1),-1)
b=A*ones(n,1)
P =输出(用于条件数旗是p准则:P =)
pp=cond(A,p)
pause
[m,] =重大(一)
nb=n+1;Ab=[A b]
R =输出(请输出手册,人工操作输出1,自动行为
从厦门活栓网转载请划出出处。

This entry was posted in dafabet手机版. Bookmark the <a href="https://www.ebuym.com/dafabetsjb/75.html" title="Permalink to 信息与计算张椿木.doc" rel="bookmark">permalink</a>.

发表评论

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