经模03：程序设计.pptxVIP

Matlab程序设计 经济模型与Matlab应用 1、基础 一、Matlab程序 Matlab m-文件 命令文件 运行 命令窗口：输入M文件的文件名 M文件窗口：F5 函数文件 function [out1, out2, ...] = myfun(in1, in2, ...) 函数表达式 使用：命令窗口：文件名（……） 简单条件语句 多条件条件语句 if (条件式) 语句组 end if (条件式1) 语句组1 elseif (条件式2) 语句组2 elseif (条件式3) 语句组3 else ... ... end 多选择条件语句 if (条件式) 语句组 1 else 语句组 2 end if -else-end分支结构 2、条件语句 例1 函数 if (x0) y=-1 elseif (x=0) y=0 else y=1 end 编程1 l01.m 编程2函数？ ? 另：x？ 输入一个x的值,输出符号函数y的值 array 是一个矩阵：循环列 一般：数组 一般：循环变量=初值:步长:终值 默认：步长=1 第一类循环语句 第二类循环语句结构 for 循环变量= array 循环体语句组 end while (条件式) 循环体语句组 end 3、循环语句 生成一个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 例4 求自然数的前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、程序流控制指令 break 跳出循环 一般和 if结构 结合使用 continue 结束本次循环 继续进行下次循环 return 终止当前指令 返回上一级或等待键盘输入命令 pause 暂停 按任意键继续 input 输入 提示键盘输入 编程 例6：连续奇数求和 clear sum1=0; for i=1:100 n=2*i-1; if sum11000 sum1=sum1+n; else break end end sum1,n l07.m 程序对不对？ 从1开始 一直到和达到1000为止 问：加了多少项 工资、薪金所得适用 级 年终奖 税率 1 (0,18000) 3 2 (18000,54000) 10 3 (54000,108000) 20 4 (108000,420000) 25 5 (420000,660000) 30 6 (660000,960000) 35 7 (960000,?) 45 二、初等模型 1、个人所得税问题 问题 1、纳税额计算函数 条件判断 程序 t1.m t2.m 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元） 税率(%) 1 0-1500 3 2 1500-4500 10 3 4500-9000 20 4 9000-35000 25 5 35000-55000 30 6 55000-80000 35 7 80000- 45 问题 2、分配月收入、年奖金 合理避税 为什么 如何分 程序 搜索：定步长 for-end while -end 纪