第2章MATLAB程序设计1.pptVIP

例2-2 利用M文件建立MYMAT矩阵。 （9）eye 函数和 size 函数 作用：size函数用来返回指定矩阵的行数和列数。 eye 函数产生指定行数和列数的单位矩阵。 两个函数共同使用可产生与指定矩阵同阶的单位矩阵。 格式：Y＝eye（size（D）） 三、矩阵的拆分 1．矩阵元素通过下标引用矩阵的元素，例如A(3,2)=200采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在MATLAB中，矩阵元素按列存储，先第一列，再第二列，依次类推。例如A=[1,2,3;4,5,6];A(3)ans = 2显然，序号(Index)与下标(Subscript )是一一对应的，以m×n矩阵A为例，矩阵元素A(i,j)的序号为(j-1)*m+i。其相互转换关系也可利用sub2ind和ind2sub函数求得。 三、矩阵的拆分 1．矩阵元素 （a）通过下标引用矩阵的元素，例如A(3,2)=200 例：A=[1,2,3;4,5,6]; sub2ind(size(A),1,2) [i,j]=ind2sub(size(A),3) 三、矩阵的拆分 2．矩阵拆分 (1) 利用冒号表达式获得子矩阵 三、矩阵的拆分 三、矩阵的拆分 例求一元二次方程ax2+bx+c=0的根 程序: a=input(a=); b=input(b=); c=input(c=); d=b^2-4*a*c; x1=(-b+sqrt(d))/(2*a); x2 =(-b-sqrt(d))/(2*a); disp([x1,x2]); 运行结果 a=2 b=7 c=3 -0.5000 -3.0000 pause 程序运行时，为了查看程序的中间结果，有时需要暂停程序的执行，可使用pause函数 格式：pause(暂停秒数) 如果省略时间，则知道按任意键后程序继续运行． 强行终止程序运行可使用Ctrl+C 2.3.3选择结构 单分支if语句 格式:if 条件 语句组 end 条件成立时执行if与end之间的语句,执行完后继续执行后边的语句. 条件不成立时直接执行后边的语句. 单分支 x=input(‘x=‘) 　if x=3 　 y=sqrt(3-x); 　end 　if x3 　 y=sqrt(x-3); 　end 双分支if结构 双分支if语句 格式:if 条件 语句组1 else 语句组２ end 例 双分支写： 　x=input(‘x=‘) 　if x=3 　 y=sqrt(3-x); 　else 　 y=sqrt(x-3); 　end 多分支if语句 if 条件１ 　　语句组１ 　elseif 条件２ 　　　　语句组２ 　．．． 　else 条件n 　　　语句组n end x=input(‘x=‘) if x=2 y=x+x-2; elseif x2 x=0 y=x-(x-2); else y=-x-(x-2); end switch语句 根据表达式的取值不同分别执行不同的语句 格式：switch 表达式 　　　　　　case 表达式1 语句组１ 　　　　　　case 表达式２ 语句组２ 　　　　　　．．． 　　　　　　otherwise 　　　　　　　语句组n end try语句 这是试探性结构语句 格式：try 　　　　　语句组１ 　　　　catch 语句组２ 　　　　end 先试探语句组１，如果出现错误，则将错误信息符给系统变量lasterr并转去执行语句组２． 例 程序：a=[1,2,3;4,5,6]; b=[6,5,4;3,2,1]; try c=a*b; catch c=a.*b; end disp(lasterr); disp(c); 运行结果 Error using == * Inner matrix dimension