专题软件开发实践报告_第1页
专题软件开发实践报告_第2页
专题软件开发实践报告_第3页
专题软件开发实践报告_第4页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

教学实践报告“反问题”的数值求解学生姓名专业信息与计算科学班级1123202指导教师2014年9月5日第一周实践报告任务安排时间计划本次实践共三周,又第二周上课时间只有两天,所以总共开发时间12天。经讨论小组将第一周开发最初的界面,第二周整合各组员的界面,第三周修改代码、程序调试、优化界面。任务分工计划小组成员,界面开发程序Filter_fn、Gen_data等六个,每人负责两个程序的界面开发,实现组员之间分工清晰、相互合作。熟悉“反问题”的数值解法研究2.1文献资料的翻译通过指导老师及谷歌翻译得到资料的大致翻译,对研究问题有一个比较清晰的理解。2.2“反问题”数值解法第一周所做:3.1反问题的理解及翻译3.2设置整体的界面及控件callback3.3把图形画到指定的区域3.4设置变量调用函数程序理解:functionvarargout=d(varargin)%---Executesonbuttonpressinpushbutton1.functionpushbutton1_Callback(hObject,eventdata,handles)%设置的控件反应%hObjecthandletopushbutton1(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)fg;%调用函数%---Executesonbuttonpressinpushbutton2.functionpushbutton2_Callback(hObject,eventdata,handles)%hObjecthandletopushbutton2(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)jiemian;%---Executesonbuttonpressinpushbutton3.functionpushbutton3_Callback(hObject,eventdata,handles)%hObjecthandletopushbutton3(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)s1=str2double(get(handles.edit2,'String'))%从文本框中获得数据赋值给s1s2=str2double(get(handles.edit3,'String'))%从文本框中获得数据赋值给s2s3=str2double(get(handles.edit4,'String'))%从文本框中获得数据赋值给s3h=1/s1;x=[h/2:h:1-h/2]';%ComputematrixKcorrespondingtoconvolutionwithGaussiankernel.kernel=(1/sqrt(pi)/s2)*exp(-(x-h/2).^2/s2^2);K=toeplitz(kernel)*h;%Setuptruesolutionf_trueanddatad=K*f_true+error.f_true=.75*(.1<x&x<.25)+.25*(.3<x&x<.32)+(.5<x&x<1).*sin(2*pi*x).^4;Kf=K*f_true;eta=s3/100*norm(Kf)*randn(s1,1)/sqrt(s1);d=Kf+eta;%Displaythedata.%figure(2)plot(x,f_true,'-',x,d,'o',x,Kf,'--')xlabel('xaxis')axes(handles.axes1);%把图形画到指定的位置%---Executesonbuttonpressinpushbutton4.functionpushbutton4_Callback(hObject,eventdata,handles)%hObjecthandletopushbutton4(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)s1=str2double(get(handles.edit2,'String'))s2=str2double(get(handles.edit3,'String'))s3=str2double(get(handles.edit4,'String'))%Setupgrid.h=1/s1;x=[h/2:h:1-h/2]';%ComputematrixKcorrespondingtoconvolutionwithGaussiankernel.kernel=(1/sqrt(pi)/s2)*exp(-(x-h/2).^2/s2^2);K=toeplitz(kernel)*h;%figure(1)mesh(K)title('MeshPlotRepresentationofMatrixK')axes(handles.axes3);%---Executesonbuttonpressinpushbutton5.functionpushbutton5_Callback(hObject,eventdata,handles)%hObjecthandletopushbutton5(seeGCBO)%eventdatareserved-tobedefinedinafutureversionofMATLAB%handlesstructurewithhandlesanduserdata(seeGUIDATA)s1=str2double(get(handles.edit2,'String'))s2=str2double(get(handles.edit3,'String'))s3=str2double(get(handles.edit4,'String'))%Setupgrid.h=1/s1;x=[h/2:h:1-h/2]';%ComputematrixKcorrespondingtoconvolutionwithGaussiankernel.kernel=(1/sqrt(pi)/s2)*exp(-(x-h/2).^2/s2^2);K=toeplitz(kernel)*h;%Setuptruesolutionf_trueanddatad=K*f_true+error.f_true=.75*(.1<x&x<.25)+.25*(.3<x&x<.32)+(.5<x&x<1).*sin(2*pi*x).^4;Kf=K*f_true;eta=s3/100*norm(Kf)*randn(s1,1)/sqrt(s1);d=Kf+eta;%Displaythedata.%ComputeaneigendecompositionofK.Kissymmetric,sothisis%equivalenttoanSVD.[V,svals]=eig(K);[svals,indx]=sort(-diag(svals));%Sort-eigenvaluesindecreasingorder.svals=-svals;%+eigsareinincreasingorde

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论