模型03：Matlab程序.pptxVIP

数 学 建 模Matlab 程序设计一、 Matlab程序1、基础命令行 变量m文件编辑[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;endmovie(ft)网格点矩阵 绘图取坐标绘图去网格设坐标获取影片的帧画面 播放l01.m2、条件语句if -else-end分支结构简单条件语句多条件条件语句 if (条件式) 语句组 end if (条件式1) 语句组1 elseif (条件式2) 语句组2 elseif (条件式3) 语句组3 else ... ... end多选择条件语句 if (条件式) 语句组 1 else 语句组 2 end例1函数输入一个x的值,输出符号函数y的值编程1if (x0) y=-1elseif (x=0) y=0else y=1end?另：x？编程2函数？l02.mswitch-case结构例2month=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 endswith 开关表达式 case表达式1 语句段1 case表达式2 语句段2 …… otherwise 语句段nendl03.m3、循环语句第一类循环语句第二类循环语句结构for 循环变量= array 循环体语句组endwhile (条件式) 循环体语句组 endarray 是一个矩阵：循环列一般：数组一般：循环变量=初值:步长:终值默认：步长=1例3 生成一个6阶矩阵，使其主对角线上元素皆为1，与主对角线相邻元素皆为2，其余皆为0。编程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 endend a结果：l04.m例4求自然数的前n项和编程n=input(‘n=‘)sum=0; k=1;while k=n sum=sum+k; k=k+1;endsum结果：n=100 Sum=5050好不好？怎么改？1234注：赋指语句——自相关应用：求和l05.m例5sum=zeros(6,1);for n=eye(6,6) sum=sum+n;endsum结果是多少？l06.m4、程序流控制指令 break 跳出循环 一般和 if结构 结合使用 continue 结束本次循环 继续进行下次循环 return 终止当前指令 返回上一级或等待键盘输入命令 pause 暂停 按任意键继续 input 输入 提示键盘输入指令编程例6clearsum1=0;for i=1:100 n=2*i-1; if sum11000 sum1=sum1+n; else break endendsum1,n连续奇数求和从1开始一直到和达到1000为止问：加了多少项程序对不对？l07.m二、程序设计实例Matlab1、 Fibonacci数列问题：求前n项？Fibonacci数列元素为：1，1，2，3，5，8，…，满足关系F(1)=1F(2)=1F(n)=F(n-1)+F(n-2)程序pf.mpf1.mfunction f=pf(n)n=?if n==1 f=1;elseif n==2 f=[1 2];else f=[1 2]; for i=3:n f(i)=f(i-1)+f(i-2); endend黄金分割2、个人所得税问题工资、薪金所得适用级数全月应纳税所得额（基数3500元）年终奖税率(%)10-1500(0,18000)321500-4500(18000,54000)1034500-9000(54000,108000)2049000-35000(108000,42000055000(420000,66000080000(660000,960000(960000,?)45年奖金有什么问题？程序问题t1.m t2.m1、纳税额计算函数条件判断function f=t1(x)t=x-3500;if t=0 f=0;elseif t=1500 f=t*0.03;elseif t=4500 f=45+(t-1500)*0.1;……else f=29625+(t-80000)*0.45;end级数月应纳税所得额（基数3500元）税率(%)10-15