模型09：与Matlab程序 .pptVIP

主讲人：孙云龙 数学建模课件 主讲人：孙云龙 数学建模课件 Email:sunyl@swufe.edu.cn 数学建模 第九讲 程序设计 数学实验 与 [x,y]=meshgrid(-1:0.1:1); z=x.*y; surf(z); ta=axis; for i=1:20 surf(sin(0.5*i.*z)) shading interp axis(ta) ft(:,i)=getframe; end movie(ft) 程序设计 网格点矩阵 绘图 取坐标 绘图 去网格 设坐标 获取影片的帧画面 播放 l01.m 引例 一、 Matlab程序 简单条件语句 多条件条件语句 if (条件式) 语句组 end if (条件式1) 语句组1 elseif (条件式2) 语句组2 elseif (条件式3) 语句组3 else ... ... end 多选择条件语句 if (条件式) 语句组 1 else 语句组 2 end 1、条件语句 if -else-end分支结构 例1 函数 if (x0) y=-1 elseif (x=0) y=0 else y=1 end 编程1 l02.m 编程2函数？ ? 另：x？ 输入一个x的值,输出符号函数y的值 switch-case结构 swith 开关表达式 case 表达式1 语句段1 case 表达式2 语句段2 …… otherwise 语句段n end month=input(‘month=’); switch month case{3,4,5} season=spring case{6,7,8} season=summer case{9,10,11} season=autumn otherwise season=winter end l03.m 例2 array 是一个矩阵：循环列 一般：数组 一般：循环变量=初值:步长:终值 默认：步长=1 第一类循环语句 第二类循环语句结构 for 循环变量= array 循环体语句组 end while (条件式) 循环体语句组 end 2、循环语句 生成一个6阶矩阵，使其主对角线上元素皆为1，与主对角线相邻元素皆为2，其余皆为0。 例3 for i=1:6 for j=1:6 if i==j a(i,j)=1; elseif abs(i-j)==1 a(i,j)=2; else a(i,j)=0; end end end a 结果： 编程 l04.m 求自然数的前n项和 编程 结果： n=100 Sum=5050 n=input(‘n=‘) sum=0; k=1; while k=n sum=sum+k; k=k+1; end sum l05.m 好不好？ 怎么改？ 1 2 3 4 注：赋指语句——自相关 应用：求和 例4 例5 l06.m sum=zeros(6,1); for n=eye(6,6) sum=sum+n; end sum 结果是多少？ 指令 3、程序流控制指令 break 跳出循环 一般和 if结构 结合使用 continue 结束本次循环 继续进行下次循环 return 终止当前指令 返回上一级或等待键盘输入命令 pause 暂停 按任意键继续 input 输入 提示键盘输入 编程 连续奇数求和 clear sum1=0; for i=1:100 n=2*i-1; i