数学实验选修课第1讲MATLAB入门资料.ppt

* * * 结果:当a=1,b=2,c=3时的数值解是 1. for循环：允许一组命令以固定的和预定的次数重复 for x=array {commands} end 在for和end语句之间的命令串{commands}按数组(array)中的每一列执行一次. 在每一次迭代中，x被指定为数组的下一列，即在第n次循环中，x=array(：，n) 控制流 MATLAB提供三种决策或控制流结构： for循环、while循环、if-else-end结构. 这些结构经常包含大量的MATLAB命令，故经常出现在MATLAB程序中，而不是直接加在MATLAB提示符下. 例 对n=1,2,…,10,求xn= 的值. while (expression) {commands} end 只要在表达式(expression)里的所有元素为真，就执行while和end语句之间的命令串{commands}. 2. While循环 与for循环以固定次数求一组命令相反，while循环以不定的次数求一组语句的值. 例 设银行年利率为11.25%.将10000元钱存入银行，问多长时间会连本带利翻一番？ 3. if-else-end结构 （1）有一个选择的一般形式是： if (expression) {commands} end 如果在表达式(expression)里的所有元素为真，就执行if和end语句之间的命令串{commands}. MATLAB(fun1) MATLAB命令窗口输入fun1(2),fun1(-1)即可. 2) 有三个或更多的选择的一般形式是： if （expression1） {commands1} else if （expression2） {commands2} else if （expression3） {commands3} else if …… ………………………………… else {commands} end end end …… end 先建立M文件fun2.m来定义函数f（x）,再在MATLAB命令窗口输入fun2(2),fun2(0.5), fun2(-1)即可. * * * * * * * * * * * * * * * * * * * 【例8】试用矩阵求逆解法求解矩阵A为系数矩阵的线性代数方程组Ax=b的解。 矩阵求逆及线性代数方程组求解 矩阵求逆及线性代数方程组求解 3. 直接解法 对于线性代数方程组Ax=b，我们可以运用左除运算符“\”像解一元一次方程那样简单地求解： x=A\b 当系数矩阵A为N*N的方阵时，MATLAB会自行用高斯消去法求解线性代数方程组。 矩阵求逆及线性代数方程组求解 解法1：分别解方程组 (1)Ax=b1；(2)Ay=b2 矩阵求逆及线性代数方程组求解 解法2：将两个方程组连在一起求解：Az=b 很明显，这里的解z的两个列向量便是前面分别求得的两组解x和y 矩阵求逆及线性代数方程组求解 习题 1、解方程组Ax＝b，分别用求逆解法与直接解法求其解。 2、编一个m程序，求N阶方阵A的行列式的值。 数据处理方法 1、查取最大值 MAX函数的命令格式有： [Y,I]= max (X)：将max(X)返回矩阵X的各列中的最大元素值及其该元素的位置赋予行向量Y与I；当X为向量时，则Y与I为单变量。 [Y,I]=max(X,[],DIM)：按数组X的第DIM维的方向查取其最大的元素值及其该元素的位置赋予向量Y与I。 【例1】查找下面数列x的最大值。 x=[3 5 9 6 1 8] % 产生数列x x = 3 5 9 6 1 8 y=max(x) % 查出数列x中的最大值赋予y y = 9 [y,I]=max(x) % 查出数列x中的最大值及其该元素的位置赋予y,I y = 9 I = 3 数据处理方法 【例2】分别查找下面3×4的二维数组x中各