matlab数学实验报告.pdf

该【matlab数学实验报告 】是由【青山代下】上传分享，文档一共【29】页，该文档可以免费在线阅读，需要了解更多关于【matlab数学实验报告 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。:..MATLAB数学实验报告指导老师:班级:小组成员:时间:_/_/_:..Matlab第二次实验报告小组成员:题目:实验四;MATLAB选择结构与应用实验目的:掌握if选择结构与程序流程控制;重点掌握break;return;pause语句的应用..问题:问题:验证“哥德巴赫猜想”;即:任何一个正偶数n>=均可表示为两个质数的和..要求编制一个函数程序;输入一个正偶数;返回两个质数的和..问题分析:由用户输入一个大于的偶数;由input语句实现..由if判断语句判断是否输入的数据符合条件..再引用质数判断函数来找出两个质数;再向屏幕输出两个质数即可..编程:functionz;z=geden;n=input'pleaseinputn'ifn<:..disp'dataerror';returnendifmodn;==fori=:n/k=;forj=:sqrtiifmodi;j==k=k+;endendforj=:sqrtn-iifmodn-i;j==k=k+;endendifk==fprintf'twonumbersare'fprintf'%.f;%.f';i;n-ibreakendend:..end结果分析:如上图;用户输入了大于的偶数返回两个质数和;通过不断试验;即可验证哥德巴赫猜想..纪录:if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰;更快的解决问题..题目:实验四;MATLAB选择结构与应用实验目的:用matlab联系生活实际;解决一些生活中常见的实际问题..:..问题:问题四:在一边长为的四个顶点上各站有一个人;他们同时开始以等速顺时针沿跑道追逐下一人;在追击过程中;每个人时刻对准目标;试模拟追击路线;并讨论..()四个人能否追到一起()若能追到一起;每个人跑过多少路程()追到一起所需要的时间设速率为问题分析:由正方形的几何对称性和四个人运动的对称性可知;只需研究个人的运动即可解决此问题..编程:holdonaxis;a=;;b=;;k=;dt=.;v=;whilek<d=norma-b;k=k+;:..plota;a;'r.';'markersize';;plotb;b;'b.';'markersize';;fprintf'k=%.fb%.f;%.fa%.f;%.fd=%.f';k;b;b;a;a;da=a+b-a/d*dt;b-a/d*dt;b=b+b-a/d*dt;-b-a/d*dt;ifd<=.breakendendfprintf'每个人所走的路程为:%.f';k*v*dtfprintf'追到一起所需要的时间为%.f';k*dt结果分析::..上图为人的模拟运动路线;有对称性可解决所提问题..-上图为运算过程和运算结果..四个人可以追到一起;走过的路程为.;时间也为..纪录:此题利用正方形和运动的对称性可以简便运算..题目:实验八;河流流量估计与数据插值目的:由一些测量数据经过计算处理;解决一些生活实际问题..:..问题:实验八上机练****题第三题:瑞士地图如图所示;为了算出他的国土面积;做以下测量;由西向东为x轴;由南向北为y轴;从西边界点到东边界点划分为若干区域;测出每个分点的南北边界点y和y;得到以下数据mm..已知比例尺:;计算瑞士国土面积;精确值为平方公里..测量数据如下:x=............;y=;y=;问题分析:先由题目给定的数据作出瑞士地图的草图;再根据梯形法;使用trapz语句;来估算瑞士国土的面积..编程:x=............;:..y=;y=;plotx;y;'r.';'markersize';;plotx;y;'r.';'markersize';;axisgrid;holdont=:;u=splinex;y;t;u=splinex;y;t;plott;uplott;us=trapzt;u;s=trapzt;u;s=s-s**/;fprintf'S=%.f';s结果分析::....上图为运算结果;计算出瑞士的国土面积为平方公里;与准确值较为接近..:..;trapz语句可以方便计算不规则图形面积;但存在一定误差..题目:实验七:圆周率的计算与数值积分目的:将数值积分最基本的原理应用于matlab之中;解决一些与积分有关的问题..问题:实验七上机练****题第一题:;为了计算最大排洪量;需要计算其断面积;试根据所给数据m用梯形法计算其断面积..问题分析:河床断面可近似分割成若干曲边梯形;近似处理把它们当做梯形来计算面积可使问题得到简化..编程:clc;clear;x=;y=;:..plotx;y;'k.';'markersize';;axis;grid;holdont=:;u=splinex;y;t;plott;u;s=*-trapzt;u;fprintf's=%.f';s果分析:上图为河床的断面图..:...平方米..纪录:使用梯形法计算不规则图形面积十分简便易行..题目:实验七:圆周率的计算与数值积分目的:使用matlab计算解决一些有关积分的问题..问题:实验七上机练****题第三题:从地面发射一枚火箭;在最初秒内记录其加速度如下;试求火箭在秒时的速度..Ts=;:.....;加速度为速度的微分;已知微分求积分;类似于面积问题;可使用梯形法来计算..编程:clc;clear;x=;y=...........;plotx;y;'k.';'markersize';;axis;grid;holdons=::;z=splinex;y;s;plots;y;v=trapzx;y;fprintf'v=%.f';v结果分析::....上图为计算结果;..纪录:梯形法可以推广解决许多已知微分求积分的其他问题..:..题目:实验七:圆周率的计算与数值积分目的:计算曲线弧长闭曲线周长可使用微元法;ds=sqrtdx^+dy^;在转化微积分问题;累加即可得到结果..问题:实验七上机练****题第三题:计算椭圆想x^/+y^=的周长;使结果具有五位有效数字..问题分析:编程:s=;dx=.;forx=:.:.dy=.-x+..^/-.-x.^/;ds=sqrtdx.^+dy.^;s=s+ds;ends=*s;fprintf'thelengthis'fprintf'%.f';s结果分析::..上图为计算结果;:计算不规则曲线弧长;可使用微元法;划分为若干小的看做直角三角形;利用勾股定理解决..题目:实验九人口预测与数据拟合目的:掌握一些曲线拟合的方法;了解曲线拟合常用函数..问题:用电压U=v的电池给电容器充电;t时刻的电压Vt=U-U-Vexp-t/τ;其中V是电容器的初始电压;τ是充电常数;由所给数据确定V和τ..t=.;V=........;:..问题分析:题中已给出函数关系式;为指数函数曲线拟合;将所给函数式整理可得标准的exp形函数曲线;从而便于解决..编程:t=.;V=........;plott;V;'k.';'markersize';;axis;grid;=;a=sumt:n;b=sumt:n.*t:n;c=sumlogV:n;d=sumt:n.*logV:n;A=na;ab;B=c;d;p=invA*Bx=:;y=expp+p*x;plotx;y;'r-';'linewidth';:..结果分析:上图为电压与时间关系图..上图为计算结果;即U-V=.;所以V=.;-/τ=-.;所以τ=.纪录:曲线拟合的一个重难点是选择合适的曲线函数;才能提高拟合度..:..题目:实验七圆周率的计算与数值积分目的:拓展圆周率的各种计算方法;掌握其他数值的近似计算方法..问题:实验七练****计算ln的近似值精确到的-次方()利用级数展开的方法来计算()利用梯形法计算()利用抛物线法问题分析:级数展开;梯形法;抛物线法是常见的近似运算方法..编程:级数展开的方法clc;clear;n=;r=;p=;k=-;whiler>=.e-n=n+;:..k=k*-;p=p+k/n;r=absp-p;fprintf'n=%.f;p=%.f';n;p;p=p;end梯形法clc;clear;f=inline'./x';x=:.:;y=fx;p=trapzx;y;fprintf'p=%.f';p抛物线法clc;clear;f=inline'./x';a=;b=;n=;z=quadf;a;b;fprintf'z=%.f';z:..结果分析:()级数展开的方法()梯形法:..抛物线法纪录:级数展开法;梯形法;抛物线法;计算近似值时应合理利用..梯形法和抛物线法不易提高精确度;级数展开法可以提高精确度..题目:实验八河流流量估计与数据插值目的:掌握求插值多项式的方法;并利用此计算近似值..问题:已知y=fx的函数表如下x=......;y=......;求四次拉格朗日插值多项式;并由此求f.:..问题分析:利用所给函数表可计算拉格朗日插值多项式..编程:functionp=lagrangex;yL=lengthx;a=onesL;forj=:La:;j=a:;j-.*x';endx=inva*y';fori=:Lpi=xL-i+;endx=......;y=......;plotx;y;'k.';'markersize';axisgrid;holdon;p=lagrangex;y;t=:.:.;u=polyvalp;t;plott;u;'r-':..a=polyvalp;.结果分析:上图为所求结果;估算值和插值多项式..纪录:插值多项式是一项十分实用的方法..:..求正整数n的阶乘:p=***…*n=n;并求出n=时的结果目的:练****使用循环变量解决数学问题问题:对程序:Clear;clc;n=;p=;fori=:np=p*i;fprintfi=%.f;p=%.f’;i;pend进行修改使它:利用input命令对n惊醒赋值问题分析:题中给出程序中“n=”修改;使用input命令;讲题中的输出命令放出循环之外..编程:clear;clc;n=input'n=';p=;fori=:np=p*i;end:..n=i=;p=>>结果分析:使用input命令可以实现人机对话;使用户自由赋值;输出语句在程序中的位置对输出的结果有很大的影响;在循环内部可以在计算过城中不断输出结果;在循环之外则可以控制只输出最后结果..题目:对于数列{√};n=;;…;求当其前n项和不超过时的n的值及合的大小..目的:运用条件循环解决文帝个项数的循环程序求解;问题:对程序:clear;clc;n=;s=;whiles<=n=n+;s=s+sqrtn;fprintfn=%.f;s=%.f’;n;s:..ends;导致s超过上限后仍有一次的循环;若把循环变量改为这次的s;则可以避免这种情况的发生..编程:clear;clc;n=;s=;whiles+sqrtn<=n=n+;s=s+sqrtnfprintf'n=%.f;s=%.f';n;send结果:……s=.n=;s=.s=:...n=;s=.s=.n=;s=.>>从结果中可以看出;最后一步为我们需要的答案;从这道题我们可以得出循环变量对一道编程的重要性..