数值分析部分算法的可视化界面设计.docVIP

数值分析部分算法的可视化界面设计 数值计算实践—课程设计报告 课题名称 数值分析部分算法的可视化界面设计 完成时间 2013年 01 月 17 日 姓名 班级 学号 101611 111584 成绩 一、实验目的，内容 二、相关背景知识介绍 三、程序代码 四、运行结果 五、制作过程中出现的问题，解决方法及体会 实验目的及内容 1.1实验目的：运用MATLAB将数值分析部分算法的程序制作成可视化文件，以之方便一些不熟悉编程或是不熟悉计算所用函数的人员对相关问题进行数值求解。 1.2实验内容： MATLAB 软件可以创建人机互动可视化界面，然后在互动界面上实现处理数据及显示相关图形的功能。 1.2.1设计的思想: （1）创建可视化界面（图1.1） （2）为可视化界面添加控件 （3）设置控件的的相应属性 （4）添加控件的事件代码实现处理数据及显示相关结果的功能 图1.1 二．相关背景知识介绍 2.1插值： 插值基函数为，，(1 i=j; 0 ij)为n个线性无关的n次多项式,可作为的基底，任意，，根据插值条件：，可得：为次插值多项式。 2.2插值： 由于如果要增加一个节点，则插值基函数就都要改变，计算量大大增加。插值具有递推性。假定为在个节点上对函数所做的次插值多项式，，，有，所以为的零点，，为待定系数，根据插值条件及解方程组可得：n次差商型插值公式为，进而可得到节点的插值公式为。 2.3最小2乘多项式拟合： （1） 当拟合函数为多项式时，称为多项式拟合，满足式（1）的称为最小二乘拟合多项式。 特别地，当时，称为线性拟合或直线拟合。显然 为的多元函数，因此上述问题即为求的极值问题。由多元函数求极值的必要条件，得 （2） 即 2.4线性方程组迭代解法： 将系数矩阵A分为D,L,U三部分，其中A=D-L-U;D为A的对角元素部分。L为下三角矩阵，U为上三角矩阵。 迭代格式为：= （3）雅克比迭代法：, （4）高斯-塞德尔迭代法：, 2.5数值积分 2.5.1复合梯形： 将所给区间[a,b]n等分，求步长 由公式 求得点 通过 即将 中求得积分数值 2.5.2复合辛普森： 将所给区间[a,b]n等分，求步长 由公式，，求得点 3）求得积分数值 2.5.3龙贝格算法 引入记号，，，，根据递推公式进行递推，二分后得到梯形值，进行次加速后得到 ，取，求得，按递推公式计算，得到表格，得出结果。 2.6常微分数值解 2.6.1欧拉法： 一阶常微分的方程及，其中， 给出步长h，用下列公式进行计算 进行迭代。 2.6.2改进欧拉法： 给定求解的初值问题及，其中，给出步长h 用欧拉公式求其值进行迭代。 2.6.3四阶龙格库塔： 给定求解的初值问题及，其中，以及步长h， 用 和计算出；用for循环。在for循环中使用公式 三．程序代码 3.1退出键程序(通用): function pushbuttonEXIT_Callback(hObject, eventdata, handles) % hObject handle to pushbuttonEXIT (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) ss=questdlg(你真的要退出,退出信息窗口！,不，我还想再看看！,是的，我要退出！,是的，我要退出！); switch ss case 是的我要退出！ delete(handles.figure1); end 3.2返回键程序（通用）: function pushbutton5_Callback(hObject, eventdata, handles) % hObject handle to pushbutton5 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) fopen(SZFX); delete(handles.figure1); 3.3计算键程序： 3.3.1