《Matlab语言》实验四.doc

实验四、 Matlab的程序设计 （一）实验目的 1.掌握利用if语句实现选择结构的方法。 2.掌握利用switch语句实现多分支选择结构的方法。 3.掌握利用for语句实现循环结构的方法。 4.掌握利用while语句实现循环结构的方法。 5. 掌握Matlab函数的编写及调试方法。 (二实验环境 1．2． 1、编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。 function y=myfun1(x) 选择一些数据测试你编写的函数。 function [y]=myfun1(x) if x=0 y=sin(x); elseif x0x=3 y=x; elseif x3 y=-x+6; End 结果: x=5; [y]=myfun1(x) y = 1 2、编写一个函数求向量x中元素的平均值、最大值、最小值、均方根值。 function [m_x,max_x,min_x,rms_x]=myfun2(x) 方均根值(Root Mean Square)的计算公式为： 用下面数据测试你写的函数： 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; %求最大值 采用逐个比较方式 if x(1)x(2) max_x=x(1); else max_x=x(2); end for k=3:n if max_xx(k) max_x=x(k); else max_x=max_x; %可省略 end end %求最小值 if x(1)x(2) min_x=x(1); else min_x=x(2); end for k=3:n if min_xx(k) min_x=x(k); else min_x=min_x; %可省略 end end %求均方根值 sum_x2=0; for k=1:n sum_x2=sum_x2+x(k).^2; rms_x=sqrt(sum_x2/n); end m_x; max_x; min_x; rms_x; %按照函数值行 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-007 max_x = 1.0000 min_x = -1.0000 rms_x = 0.7071 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.5231 max_x = 0.9994 min_x = 0.0056 rms_x = 0.5869 3、编写求解方程的根的函数（这个方程不一定为一元二次方程，因的不同取值而定），这里应根据的不同取值分别处理，有输入参数提示，当时应提示“为恒不等式!”。并输入几组典型值加以检验。 M文件: disp(关于方程ax^2+bx+c=0) a=input(请输入系数a=); b=input( b=); c=input( c=); if a==0b==0 if c==0 disp(解为全体实数); else disp(为恒不等式!); end else P=[a,b,c]; x=roots(P) end disp(关于方程ax^2+bx+c=0) a=input(请输入系数a=); b=input( b=); c=input( c=); if a==0b==0 if c==0 disp(解为全体实数); else disp(为恒不等式!); end else P=[a,b,c]; x=roots(P) end 工作空间测试: 关于方程ax^2+bx+c=0 请输入系数a=2 b=4 c=6 x = -1.0000 + 1.4142i -1.0000 - 1.4142i 4、输入一个百分制成绩，要求输出成绩等级A+、A、