龙格—库塔法分析lorenz方程课程设计(编辑修改稿)

龙格—库塔法分析lorenz方程课程设计(编辑修改稿)内容摘要：

(1,i)。 end for r=2:n L(r,1)=A(r,1)/U(1,1)。 end for i=2:n for j=i:n z=0。 for r=1:i1 z=z+L(i,r)*U(r,j)。 end U(i,j)=A(i,j)z。 end if abs(U(i,i))eps flag=39。 failure39。 return。 end for k=i+1:n m=0。 for q=1:i1 m=m+L(k,q)*U(q,i)。 end L(k,i)=(A(k,i)m)/U(i,i)。 end end L U end 4. 结果 LU_de(A) m = 4 n = 4 L = 1 0 0 0 2 1 0 0 1 2 1 0 3 3 2 1 U = 2 4 2 6 0 1 2 3 0 0 3 6 0 0 0 1 在编写 程序过程中由于角标较多因此在运行过程中出现了不少角标不对的错误 题目；给出函数 f(x)=1/(1+25x^2)，求 f（ x 在 [1,1]上取 5 个和 9个等距节点，做最小二乘拟合，得出均方误差 方误差。 五个节点时， matlab编码为： 首先建立 M文件，并保存 function y=f(x) y=1/(1+25*x^2)。 end x=[1 0 1]。 for i=1:5 y(i)=f(x(i))。 end a=polyfit(x,y,3) syms x f1=a(1)*x^3+a(2)*x^2+a(3)*x+a(4) x=[1 0 1]。 for i=1:5 E(i)=(f(x(i))(a(1)*x(i)^3+a(2)*x(i)^2+a(3)*x(i)+a(4)))^2。 end sum(E) 输出结果为 a = f1 = 4878849915647781/1298074214633706907132624082305024*x^31600/2639*x^25348847520430703/649037107316853453566312041152512*x+1514/2639 （拟合的多项式） ans = （均方误差） 九个点的时候， matlab 编码为： x=[1 0 1]。 for i=1:9 y(i)=f(x(i))。 end。