用matlab编写拉格朗日插值算法的程序10[ 标签：matlab,插值算法,程序 ]用matlab编写拉格朗日插值算法的程序，并以下面给出的函数表为数据基础，在整个插值区间上采用拉格朗日插值法计算f（0.6），写出程序源代码，输出计算结果x -2.15 -1.00 0.01 1.02 2.03 3.25y 17.03 7.24 1.05 2.03 17.06 23.05匿名回答:1人气:6解决时间:2011-05-2419:58满意
答案
八年级地理上册填图题岩土工程勘察试题省略号的作用及举例应急救援安全知识车间5s试题及答案
好评率：83%做了一个测试，希望有所帮助。代码：%用matlab编写拉格朗日插值算法的程序，并以下面给出的函数表为数据基础，%在整个插值区间上采用拉格朗日插值法计算f（0.6），写出程序源代码，输出计算结果%x -2.15 -1.00 0.01 1.02 2.03 3.25%y 17.03 7.24 1.05 2.03 17.06 23.05functionmain()clc;x=[-2.15 -1.00 0.01 1.02 2.03 3.25];y=[17.03 7.24 1.05 2.03 17.06 23.05];x0=0.6;f=Language(x,y,x0)functionf=Language(x,y,x0)%求已知数据点的拉格朗日插值多项式%已知数据点的x坐标向量:x%已知数据点的y坐标向量:y%插值点的x坐标:x0%求得的拉格朗日插值多项式或在x0处的插值:fsymstl;if(length(x)==length(y))n=length(x);elsedisp('x和y的维数不相等！');return; %检错endh=sym(0);for(i=1:n)l=sym(y(i));for(j=1:i-1)l=l*(t-x(j))/(x(i)-x(j));end;for(j=i1:n)l=l*(t-x(j))/(x(i)-x(j));end;h=hl;endsimplify(h);if(nargin==3)f=subs(h,'t',x0);%计算插值点的函数值elsef=collect(h);f=vpa(f,6);%将插值多项式的系数化成6位精度的小数end结果：f=0.0201>>如何用MATLAB编写的拉格朗日插值算法的程序、二阶龙格-库塔方法的程序和SOR迭代法的程序，要能运行的∮初夏戀雨￠回答:2人气:29解决时间:2009-12-0819:04满意答案好评率：100%拉格朗日functiony=lagrange(x0,y0,x)n=length(x0);m=length(x);fori=1:mz=x(i);s=0.0;fork=1:np=1.0;forj=1:nifj~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)s;endy(i)=s;endSOR迭代法的Matlab程序function[x]=SOR_iterative(A,b)%用SOR迭代求解线性方程组,矩阵A是方阵x0=zeros(1,length(b));%赋初值tol=10^(-2);%给定误差界N=1000;%给定最大迭代次数[n,n]=size(A);%确定矩阵A的阶w=1;%给定松弛因子k=1;%迭代过程whilek<=Nx(1)=(b(1)-A(1,2:n)*x0(2:n)')/A(1,1);fori=2:nx(i)=(1-w)*x0(i)w*(b(i)-A(i,1:i-1)*x(1:i-1)'-A(i,i1:n)*x0(i1:n)')/A(i,i);endifmax(abs(x-x0))<=tolfid=fopen('SOR_iter_result.txt','wt');fprintf(fid,'\n********用SOR迭代求解线性方程组的输出结果********\n\n');fprintf(fid,'迭代次数:%d次\n\n',k);fprintf(fid,'x的值\n\n');fprintf(fid,'.8f\n',x);break;endk=k1;x0=x;endifk==N1fid=fopen('SOR_iter_result.txt','wt');fprintf(fid,'\n********用SOR迭代求解线性方程组的输出结果********\n\n');fprintf(fid,'迭代次数:%d次\n\n',k);fprintf(fid,'超过最大迭代次数，求解失败！');fclose(fid);endMatlab中龙格-库塔(Runge-Kutta)方法原理及实现龙格-库塔(Runge-Kutta)方法是一种在工程上应用广泛的高精度单步算法。由于此算法精度高，采取措施对误差进行抑制，所以其实现原理也较复杂。该算法是构建在数学支持的基础之上的。龙格库塔方法的理论基础来源于泰勒MATCH_
word
word文档格式规范word作业纸小票打印word模板word简历模板免费word简历
_1714189395175_0和使用斜率近似表达微分，它在积分区间多预计算出几个点的斜率，然后进行加权平均，用做下一点的依据，从而构造出了精度更高的数值积分计算方法。如果预先求两个点的斜率就是二阶龙格库塔法，如果预先取四个点就是四阶龙格库塔法。一阶常微分方程可以写作：y'=f(x,y),使用差分概念。(Yn1-Yn)/h=f(Xn,Yn)推出（近似等于，极限为Yn'）Yn1=Ynh*f(Xn,Yn)另外根据微分中值定理，存在0
>y1='3*a(2)2*a(3)-(2*a(1)^21)*exp(2*a(1))'; 解释：a(1)为x，a(2)、a(3)……为y1，y2……>>y2='4*a(2)a(3)(a(1)^22*a(1)-4)*exp(2*a(1))';>>y3='2*a(2)-a(3)-a(1)*exp(3*a(1))';>>y4='a(2)a(1)^2*exp(a(1))';>>y5='a(3)-exp(2*a(1))';>>A=[{y1},{y2},{y3},{y4},{y5}]A=[1x38char] [1x41char] [1x28char] [1x21char] 'a(3)-exp(2*a(1))'>>rk(A,0:0.1:1,0.1,ones(1,5))y=1.0000 1.0000 1.0000 1.0000 1.00001.4692 1.1649 1.1312 1.1226 0.99632.1246 1.5112 1.3309 1.3031 0.99303.0680 2.1507 1.6105 1.5679 1.00784.4629 3.2638 1.9873 1.9573 1.07186.5725 5.1403 2.4884 2.5335 1.23769.8237 8.2476 3.1587 3.3940 1.593014.9117 13.3402 4.0745 4.6925 2.283722.9721 21.6384 5.3681 6.6746 3.549735.8640 35.1212 7.2699 9.7346 5.784356.6365 57.0045 10.1818 14.5095 9.6315matlab解微分方程组dx/dt=xydy/dt=x-y2010-10-615:56提问者：刘の鱼|浏览次数：409次推荐答案2010-10-620:50不知道解得对不对程序：dsolve('Dx=xy','Dy=x-y','t')解得：x=C1*exp(2^(1/2)*t)C2*exp(-2^(1/2)*t)y=C1*2^(1/2)*exp(2^(1/2)*t)-C2*2^(1/2)*exp(-2^(1/2)*t)-C1*exp(2^(1/2)*t)-C2*exp(-2^(1/2)*t)>>dsolve('Dx=xy','Dy=x-y')ans=y:[1x1sym]x:[1x1sym]>>disp(ans.x)C1/exp(2^(1/2)*t)C2*exp(2^(1/2)*t)-(2^(1/2)*C1)/exp(2^(1/2)*t)2^(1/2)*C2*exp(2^(1/2)*t)>>disp(ans.y)C1/exp(2^(1/2)*t)C2*exp(2^(1/2)*t)参考资料：
浙公网安备 33021202002483