《Matlab语言》实验四.docx

《Matlab语言》实验四

第一篇：《Matlab语言》实验四实验四、Matlab的程序设计（一）实验目的1.掌握利用if语句实现选择结构的方法。2.掌握利用switch语句实现多分支选择结构的方法。3.掌握利用for语句实现循环结构的方法。4.掌握利用while语句实现循环结构的方法。5.掌握Matlab函数的编写及调试方法。(二）实验环境1．计算机2．MATLAB7.0集成环境（三）实验内容1、编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。functiony=myfun1(x)sinx,y(x)x,x6,x00x3x3选择一些数据测试你编写的函数。function[y]=myfun1(x)ifxy=sin(x);elseifx>0&x3y=-x+6;End结果:>>x=5;>>[y]=myfun1(x)y=12、编写一个函数求向量x中元素的平均值、最大值、最小值、均方根值。function[m_x,max_x,min_x,rms_x]=myfun2(x)方均根值(RootMeanSquare)的计算公式为：rms1NNi1xi2用下面数据测试你写的函数：function[m_x,max_x,min_x,rms_x]=myfun2(x)%求平均值sum_x=sum(x);%向量元素求和[m,n]=size(x);%最好用n=length(x);m_x=sum_x/n;%求最大值采用逐个比较方式ifx(1)>x(2)max_x=x(1);elsemax_x=x(2);endfork=3:nifmax_xelsemax_x=max_x;%可省略endend%求最小值ifx(1)min_x=x(2);endfork=3:nifmin_x>x(k)min_x=x(k);elsemin_x=min_x;%可省略endend%求均方根值sum_x2=0;fork=1:nsum_x2=sum_x2+x(k).^2;rms_x=sqrt(sum_x2/n);endm_x;max_x;min_x;rms_x;%按照函数值行（1）x=sin(0:0.01:6*pi)运行结果:>>x=sin(0:0.01:6*pi);>>[m_x,max_x,min_x,rms_x]=myfun2(x)m_x=-1.1256e-007max_x=1.0000min_x=-1.0000rms_x=0.7071（2）x=rand(1,200)，得到的x为200个（0，1）之间均匀分布的随机数。>>x=rand(1,200);>>[m_x,max_x,min_x,rms_x]=myfun2(x)m_x=0.5231max_x=0.9994min_x=0.0056rms_x=0.58693、编写求解方程axbxc0的根的函数（这个方程不一定为一元二次方程，因a、b、c的不同取值而定），这里应根据a、b、c的不同取值分别处理，有输入参数提示，当a0,b0,c~0时应提示“为恒不等式!”。并输入几组典型值加以检验。2M文件:disp('关于方程ax^2+bx+c=0')a=input('请输入系数a=');b=input('b=');c=input('c=');ifa==0&&b==0ifc==0disp('解为全体实数');elsedisp('为恒不等式!');endelseP=[a,b,c];x=roots(P)enddisp('关于方程ax^2+bx+c=0')a=input('请输入系数a=');b=input('b=');c=input('c=');ifa==0&&b==0ifc==0disp('解为全体实数');elsedisp('为恒不等式!');endelseP=[a,b,c];x=roots(P)end工作空间测试:关于方程ax^2+bx+c=0请输入系数a=2b=4c=6x=-1.0000+1.4142i-1.0000-1.4142i4、输入一个百分制成绩，要求输出成绩等级A+、A、B、C、D、E。其中100分为A+，90分～99分为A，80分～89分为B，70分～79分为C，60分～69分为D，60分以下为E。要求：（1）用switch语句实现。（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。g=input('请输入成绩');ifg100disp('错误，输入的成绩应为0~100间的数')elseg=fix(g/10);switchgcase{10}disp('A+');case{9}disp('A');case{8}disp('B');case{7}disp('C');case{6}disp('D')otherwisedisp('E');end5、利用f