2023年实验二MATLAB程序设计含实验报告

实验二MATLAB程序设计一、实验目的.掌握运用if语句实现选择结构的方法。.掌握运用switch语句实现多分支选择结构的方法。.掌握运用for语句实现循环结构的方法。.掌握运用while语句实现循环结构的方法。.掌握MATLAB函数的编写及调试方法。二、实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。M文献的编写：启动MATLAB后，点击FileINew|M-File,启动MATLAB的程序编辑及调试器（Editor/Debugger）,编辑以下程序，点击FileISave保存程序，注意文献名最佳用英文字符。点击DebuglRun运营程序，在命令窗口查看运营结果,程序如有错误则改正三、实验内容.编写求解方程o?+〃x+c=0的根的函数（这个方程不一定为一元二次方程，因4、b、C的不同取值而定），这里应根据4、b、C的不同取值分别解决，有输入参数提醒，当。=0力=0,c〜=（）时应提醒“为恒不等式！、并输入儿组典型值加以检查。（提醒:提醒输入使用了〃夕〃才函数）.输入一个百分制成绩，规定输出成绩等级A+、A、B、C、D、E。其中100分为A+,90分〜99分为A,80分〜89分为B,70分〜79分为C,60分〜69分为D,60分以下为E。

b=0c=»2此为恒不等式第二题成绩等级划分一一请输入一个百分制成绩,查询划分等级:111输出犯错成绩等级划分一一请输入一个百分制成绩,查询划分等级：55E第三题%c1assic〃3n+1”problemfromnumbertheory.while1n=inputCEntern,negativequits：z;先输入一个非负整数ifn<=0ifn<=0breakifn<=0break%假如输入的数为负数,则不做操作ifn<=0break%假如输入的数为负数,则不做操作enda=n;whilen>1ifrem(n,2)==0n=n/2;eIsen=3*n+l；enda=[a,n];enda=n;whilen>1ifrem(n,2)==0n=n/2;enda=n;whilen>1ifrem(n,2)==0n=n/2;eIsen=3*n+l；enda=[a,n];%将。的值赋值给a%满足条件n>l时循环下列语句%当。能被2整除则n除以2%假如不能被整除，则n乘以3加1%输出数组[a,n],并对比n是否>1,大于则继续上诉循环endend第四题请输入x的值：-1y=-0.8497请输入x的值：2y二2请输入X的值:9ans=1y二-3第五题请输入x的值:sin(0:0.01:6*pi)平均值为mx=-1.1256e-007最大值为max_x=1.0000最小值为minx=-1.0000均方根值为rms_x=0.7071请输入x的值:rand(l,200)平均值为m_x=0.4969最大值为maxx=0.9943最小值为min_x=0.0099均方根值为rms_x=0.5707第六题(pi八2)/6=(1/1*2)+(1/2-2)+(1/3()+……+(1/n，2),求pi的近似值请输出n值：100approximation=3.1321(Pi-2)/6=(1/厂2)+(1/2八2)+(l/3-2)+……+(l/rf2),求pi的近似值请输出n值：1000approximation:3.1406(p「2)/6=(1/T2)+(1/2人2)+(1/3-2)+……+(l/rT2),求pi的近似值请输出n值：10000approximation=3.1415.实验中出现的问题及解决方法对函数的嵌套关系理解混乱解决：反复尝试，得出对的的嵌套关系。'”没有运营成功。解决：需要加符号不然无法得出预期结果.思考题的回答a=3;b=2;s=1/2；fori=l:l5s=s+a/b；c=a；a=a+b;b=c；end规定：(1)用switch语句实现。(2)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出犯错信息。(提醒:注意单元矩阵的用法)3.数论中一个有趣的题目：任意一个正整数,若为偶数，则用2除之,若为奇数，则与3相乘再加上1。反复此过程，最终得到的结果为1。如：2-1—>10—>5―>16—>8—>4―>2―>16-3—10—5—16.8—4—2—1运营下面的程序，按程序提醒输入n=l,2,3,5,7等数来验证这一结论。请为关键的Matlab语句填写上相关注释,说明其含义或功能。%classic"3n+1"problemfromnumbertheory.while1nsinputfEntern,negativequits:');ifn<=0breakenda=n;whilen>1ifrem(n,2)==0n=n/2;.编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。function[y]=myfun1(x)

sin羽y(x)=\x,—x+6,x<Q0<x<3x>3x<Q0<x<3x>3.编写一个函数求向量x中元素的平均值、最大值、最小值、均方根值。tunction[m_x,max_x,min_x,rms_x]=myfun2(x)方均根值(RootMeanSquare)的计算公式为：用下面数据测试你写的函数：x=sin(0:0.01:6*pi)x=rand(1,200),得到的x为200个(0,1)之间均匀分布的随机数。216.根据3=1+3+1+......+!，求乃的近似值。当n分别取100、1000、61~2-3-n~10000时，结果是多少？思考题：有一分数序列：j_358B212，2，3?5?TTT…编写一段程序，求前16项的和。四、实验报告规定(包含预习报告规定和最终报告规定).实验名称[..实验目的葫习报告.实验设备及条件口J最终报告.实验内容及规定.实验程序设计指程序代码。.实验结果及结果分析实验结果规定必须客观,有数据的可以记录数据，没有数据的简朴描述实验现象。结果分析是对实验结果的理论评判。.实验中出现的问题及解决方法.思考题的回答四、实验报告的提交方式Word文档，命名方式:实验号—你的学号—姓名!!!例如本次实验:实验—一张三.doc（信息101提交报告邮箱）:E_mail：（网络工程101提交作业邮箱）:E_mail：（注意网络班的M是大写的）下一次课前提交，过期不收！五、参考文献参考教材和Matlab帮助文献。.实验名称MATLAB程序设计.实验目的.掌握运用if语句实现选择结构的方法。.掌握运用switch语句实现多分支选择结构的方法。.掌握运用f。r语句实现循环结构的方法。.掌握运用while语句实现循环结构的方法。.掌握MATLAB函数的编写及调试方法。.实验设备及条件计算机一台(带有MATLAB7.0以上的软件环境)。.实验内容及规定把实验内容的应用题，用MATLAB的语法编写出来,并运营成功，注意题与题之间用相应注释分割。.实验程序设计%1dispd一元二次方程计算器，请输入下列数值进行计算：’)；a=input('a=，);b=input('b=')；c=input('c=')；ifa=0&b==0&c〜=0disp('此为恒不等式‘)；elsed=b*b-4*a*c；x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];disp(「x1=',num2str(x(1)),z,x2=，,num2str(x⑵)]);end%2SCOres=input('成绩等级划分请输入一个百分制成绩，查询划分等级：')；switchscorescase{100}rate=，A+’;casenum2cell(90:99)rate='A';casenurn2cel1(80:89)rate='B'；casenum2ce11(70:79)rate=，C';casenum2cell(60:69)rate=，D';casenum2cel1(0:59)rate='E';otherwisedisp(z输出犯错');enddisp(rate)%3%classic〃3n+l"problemfromnumbertheory.while1n=input('Entern,negativequits:f);%输入一个非负整数ifn<=0%假如输入的数为负数，则不做操作breakenda二n;%将口的值赋值给awhi1en>l%满足条件n>1时循环下列语句ifrem(n,2)==0%当n能被2整除则n除以2n=n/2;e1se与假如不能被整除，则n乘以3加1n=3*n+l;enda=[a,n];蟒俞出数组[a,n],并对比n是否＞1,大于则继续上诉循环enda%输出最后的结果aend%4myfun1.m%c1cformatcompactformatlonggx=input('请输入x的值:‘);y=tranl(x);disp('y=')disp(y)%4tran.m%functiony=myfun1(x)%当x<=0y=sin(x);%当3>x>0y=x;%当x>3y=-x+6;ifx<=0y=sin(x);elseif3>x>0y=x;elsex>3y=-x+6;end%5tran.m%function[m_x,max_x,minx,rms_x]=myfun2(x)%分别是平均值，最大值，最小值，均方根值inx=mean(mean(x))max_x=max(max(x))min_x=min(min(x))a=sqrt(mean(x.八2));rms_x=sqrt(mean(a.2))%5myfun2.mx二input。平均值、最大值、最小值、均方根值的计算，请输入x的值，;[m_x,max_x,min_x,rms_x]=tran（x）;dispC平均值为'）m_xdisp（'最大值为'）maxxdisp（'最小值为'）min_xdisp（'均方根值为‘）rms_x%6dispC(pi^2)/6=(l/r2)+(l/2八2)+(1/3-2)+……+(l/n*2),求pi的近似值’)appr