数值分析的MATLAB程序设计.doc

数值分析作业 题目 利用MATLAB编写数值算法程序 班级 / 学号20100414028 学 生 姓 名 张 阳 指 导 教 师 刘勇进 编写lagrange插值法 计算拉氏基函数，再求和。 图 1 拉氏基函数求解流程图 Lagrange_base_function.m文件 %Copyright(c) 2012 ShenYang Areospace University %All rights reserved % %file name :Lagrange_base_function.m %summary : % %current edition : %programmer :Yung Zhung %accomplish time :2012-04-15 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %function name :Lagrange_base_function.m %function describe : %function in detailed : %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function Y = Lagrange_base_function(x, k) n = length(x); y = ones(n, 2); Y = 1; for ii = 1 : n if ii == k continue; end y(ii, 2) = -x(ii); Y = conv(Y, y(ii, :)); end end 图 2 计算lagrange插值函数及差值值 Lagrange.m文件 %Copyright(c) 2012 ShenYang Areospace University %All rights reserved % %file name :Lagrange.m %summary : % %current edition : %programmer :Yung Zhung %accomplish time :2012-04-15 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %function name :Lagrange.m %function describe : %function in detailed : %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function Y = Lagrange(x, y, point) if min(size(x)) 1 || min(size(y)) 1 error(the muxt be two vectors!!!); end n = length(x); if n ~= length(y) error(the dimention of the two vectors must aggree!!!); end Y = 0; for ii = 1 : n y_1 = Lagrange_base_function(x, ii); Y = polyadd(Y, y(ii) * y_1 / polyval(y_1, x(ii))); end if nargin == 2 disp(Y); else Y = polyval(Y, point); end 编写复合二重simpson求积公式 图 3 二重复合simpson公式流程图 转化为一重simpson求积公式。 %Copyright(c) 2012 ShenYang Areospace University %All rights reserved % %file name :simpson_con.m %summary : % %current edition : %programmer :Yung Zhung %accomplish tim