[信息与通信]第3章 MATLAB的程序设计基础.ppt

例3.1 电阻电路的计算 例3.2 戴维南定理 例3.2 戴维南定理（续） 写成 例3.3 含受控源的电阻电路 例3.3的电路图 续例3.3 含受控源的电阻电路 3.5.4 程序举例 1．编制M文件实现分段函数： （1）用函数文件来实现 1．编制M文件实现分段函数： （2）用程序文件来实现 2. 编制程序，判断输入数据的奇偶性。 （1）只考虑输入数字的情况 2. 编制程序，判断输入数据的奇偶性。 （2）考虑输入是空格或直接按回车键的情况 3. 编制程序，寻找输入数组中的最大数。 4. 求任意度数的三角函数值，并将正弦、余弦、正切及余切同时按表格形式显示出来。 函数的递归调用 调用函数的过程中出现直接或间接地调用该函数本身的现象称为函数的递归调用。在递归调用的函数中一般要有条件递归调用的语句，否则函数会一直循环下去。 例：求一个正整数的阶乘 function y=fun(n) %compute n! if n0 error(n should be positive integral or );%display the error information return; end if n==0|n==1 y=1; else y=n*fun(n-1); end 把方程组写成矩阵形式为 用基尔霍夫定理列方程组 (1)如us＝10v，求i3，u4，u7 (2)如u4＝6v，求us，i3，u7 0 ) ( 0 ) ( ) ( c 7 6 5 b 5 c 5 b 5 4 3 a 3 s b 3 a 3 2 1 = + + + - = - + + + - = - + + i R R R i R i R i R R R i R u i R i R R R clear,format compact R1=2;R2=4;R3=12;R4=4;R5=12;R6=4;R7=2;%为给定元件赋值 % 解问题（1） display(解问题（1）) a11=R1+R2+R3;a12=-R3;a13=0; %将系数矩阵各元素赋值 a21=-R3;a22=R3+R4+R5; a23=-R5; a31=0;a32=-R5;a33=R5+R6+R7; b1=1;b2=0;b3=0; us=input(us=); % 输入解(1)的已知条件A=[a11,a12,a13;a21,a22,a23;a31,a32,a33] % 列出系数矩阵A B=[b1;0;0]; I=A\B*us; % I=[ia;ib;ic] ia=I(1);ib=I(2);ic=I(3); i3=ia-ib,u4=R4*ib,u7=R7*ic% 解出所需变量 % 利用电路的线性性质及问题（1）的解 display(解问题（2）) u42=input(给定u42= ); k1=i3/us;k2=u4/us;k3=u7/us; % 由问题（1）得出待求量与us的比例系数 us2=u42/k2,i32=k1/k2*u42,u72=k3/k2*u42 % 按比例方法求出所需变量 全局变量是由MATLAB提供的global命令来设置，一般在M函数的开头定义。命令形式为： global a b c “global”命令应当在工作空间和M函数中都出现，如果只在一方出现，则不被承认为全局变量。另外，在 MATLAB 中对变量名是区分大小写的，因此，在程序中为了不与普通变量相混淆，习惯上我们将全局变量用大写字母表示。 3.5 程序流程控制 计算机程序通常都是从前到后逐条执行的，但有时也会根据实际情况，中途改变执行次序，称为流程控制。 3.5.1 循环结构 在循环结构中，一组被重复执行的语句称为循环体，每循环一次，都必须作出是否继续重复的决定，这个决定所依据的条件称为循环的终止条件。MATLAB提供了两种循环结构：for - end 循环和 while - end 循环。 1．for - end 循环结构 for - end循环语句从for开始，用end结束，其基本格式为： for 循环变量 = 表达式 循环体 end for和end是必需的，不可省略，且必须配对使用。 表达式是一个矩阵，用来表示循环的次数。表达式通常的形式为“m:s:n”，m是循环初值，n是循环终值，s为步长，s可以缺省，缺省值为 1。 for 循环变量 = 表达式 循环体 end 循环体被循环执行，执行的次数由表达式控制。循环变量依次取表达式矩阵的各列，每取一次，循环体执行一次。 循环不会因为在循环体内对循环变量重新设置值而中断。 for - end 循环结构的执行过程是：从表达式矩阵的第一列开始，依