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