MATLAB_02_程序设计.pptVIP

例题：求排列的逆序数 求排列32514的逆序数 用M函数实现： 输入：一个排列，如 k = [3 2 1 4 5] 输出：排列的逆序数 t = 5 斐波纳契数列 斐波纳契数列还暗含着许多有趣的数字规律，如从第3个数开始每隔两个必是2的倍数，从第4个数开始每隔3个必是3的倍数，从第5个数开始每隔4个必是5的倍数……另外，这个数列最具有和谐之美的地方是，越往后，相邻两项的比值会无限趋向于黄金比1.61803…… 斐波纳契数列在现代物理、准晶体结构、化学等领域都有直接的应用，为此，美国数学会从1960年代起出版了《斐波纳契数列》季刊，专门刊载这方面的研究成果。 例题：用for循环指令来寻求Fibonacc数组中第一个大于10000的元素。 n=100; a=ones(1,n); for i=3:n a(i)=a(i-1)+a(i-2); if a(i)=10000 a(i) break; end; end i 例题：水仙花数 　　1. 找出所有的“水仙花数”及水仙花数总的个数。所谓“水仙花数”是一个3位数，其各位数字的立方和等于该数本身。例如： 。 　　2. 找出所有这样的4位正整数，其各位数字的4次方的和等于该数本身，例如： 并请给出满足这样条件的4位正整数一个有多少个。 用蒙特卡罗 ( Monte Carlo ) 投点法计算 ? 的值 n=100000; a=2; m=0; for i=1:n x=rand(1)*a/2;y=rand(1)*a/2; if ( x^2 + y^2 = (a/2)^2 ) m=m+1; end end fprintf(计算出来的 pi 为：%f\n,4*m/n); 例题：蒙特卡罗投点法 循环结构：while while expression (commands) end expression称为条件表达式，(commands)为要执行的循环代码，一般情况下，expression的计算结果为一个标量，但也可以是一个数组表达式。当expression的结果为一个标量，且该标量为true时，(commands)就会被一直执行下去；当expression的结果为一个数组时，只有当数组中的所有元素均为True时，(commands)才会被一直执行下去。 例题：用while循环求Fibonacci数组中第一个大于10000的元素。 a(1)=1; a(2)=1; i=2; while a(i)=10000 a(i+1)=a(i-1)+a(i); i=i+1; end; i,a(i) 例题：用while循环计算特殊变量eps n＝0； EPS=1； while(1+EPS) 1 EPS＝EPS／2； n＝n＋l； end n, EPS＝EPS * 2 逻辑数组 A=zeros(2,5); A(:)=-4:5 L=abs(A)3 X=A(L) 找出数组 中所有绝对值大于3的元素 L与A数组同样大小，L的元素是0或1，它是“逻辑数组”。 通过逻辑数组L中“逻辑值1”所在的位置，找出A中满足条件的元素。 逻辑数组常用来代替if和loop结构 Find命令 find是matlab特有的命令，也常用来代替if和loop结构 find(A)：找出矩阵非零元素所在的下标 例： A=[0,4,0;-1,0,0]; [I,J]=find(A) find(条件)：找出符合条件的元素所在的位置 a=[4,5,78,121,3,65,24,2]; b=find(a10) A=[0,4,0,-1,0,0]; b=find(A) 控制程序流的其他常用指令 keyboard 将“控制权”交给键盘，用户可以输入各种matlab指令，仅当用户输入return指令后，“控制权”才交还给程序 break 终止while、for循环，或在if-end，switch-case，try-catch中导致中断 continue 跳过位于其后的循环中的其他指令，执行循环的下一次迭代 pause pause(n) 使程序暂停，等待用户按任意键继续 使程序暂停n秒后，在继续执行 return 结束return指令所在函数的执行 error(‘message’) 显示出错信息message，终止程序 lasterr 显示最新出错原因，并终止程序 wlastwarn 显示matlab自动给出的最新警告，程序继续运行 warning(‘message’) 显示警告信息message，程序继续运行 问题：下面的加法等式中：两个真分数之和等于第三个真分数；其中每个方框表示一个非零数字，且任意两个数字互不相同，比如：