模型03：Matlab程序

1、Matlab编程、数学建模、x，y =消息网格(-1:0.1:1 ) z=x.* y； surf(z) ta=axis； fori=1:20 surf (sin (0.5* I.* z ) ) shadinginteraxis (ta ) ft (:i)=getframe； 取得end movie(ft )、网格点矩阵绘制坐标绘制，取得网格设定坐标，取得电影的帧画面播放，l01.m，1，基础，一，Matlab程序，命令行变量m文件编辑，2，条件语句，if -else-end分支什么，还有，x？ 输入x的值，输出符号函数y的值，输出switch-case结构，swith开关式case式1状态段1

2、 case式2状态段2 otherwise状态段n end，month=I switchmonthcase 3，4，5 season=spring case 6，7，8 season=summer case 9，10，11 season=autumnotherwiseason=winter end 示例2 array是矩阵：一般循环列：一般循环变量=初始值：步骤：的最终值默认值：步骤=1，3，循环语句产生6次矩阵，使得主对角线上的所有元素都为1，主对角线上的所有元素都为2，其馀的全部为0。 例3，fori=1:6 forj=1:6 ifi=jaa (I，j)=1； elseif abs(i-j

3、)=1 a(i，j)=2； else a(i，j)=0； 最终结束a，结果：编程，l04.m，自然数的上位n项和，求出编程，结果： n=100 Sum=5050，n=input (“) sum=0； k=1； while k=n和=和k； k=k 1； end sum，l05.m，可以吗？ 12三四，注：指定语句自相关应用：合计，示例4，示例5，l06.m，sum=zeros (6，1，1 ) forn=eye (6，6 ) sum=sumn； 最终和的结果是什么？ 指令，4，程序流控制指令，编程，连续奇数加，clear sum1=0； for i=1:100 n=2*i-1； if和1100

4、0和1=和1 n； elsebreakendsum1，n，l07.m，程序是否正确？ 问题从1到和达到1000 :加上了什么，例子6，1，Fibonacci数列，问题：求前n项？ Fibonacci数列的元素是1，1，2，3，5，8，关系F(1)=1 F(2)=1 F(n)=F(n-1) F(n-2 )，2，编程示例，Matlab，n=？ 程序、pf.m pf1.m、if n=1 f=1； elseif n=2 f=1 2； else f=1 2； for I=3: nf (I )=f (I-1 ) f (I-2 )结束，function f=pf(n )，黄金分割，工资，工资收入的适用，2，

5、个人所得税问题，年度奖金问题是什么，问题，1，纳税额计算函数条件判断，步骤，t1.m 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、问题、2、分配月薪、年度奖金合理逃税的理由、程序、检索：固定步骤for-end while -end记录的好处if-end、t3.m、算法、function min，k1 =tax3(最小=t1 (x ) k1=0； while kd1 min=d1； k1=k； 结束k=k 1； end，还有，程序有问题吗

6、？ mint，k1=t3(5000 )，持续：月薪，年度奖金非专业报告，程序，t4.m，算法，循环计算各种收入分配方式分类，for i=1:n A(1，i)=100*i； A(2，I )，A(3，i)=t5(A(1，I )； endb=a (:1 ) k=1； for j=2:n if A(3，j)=B(3，k) k=k 1； B(:k)=A(:j )最终B(:k1 )=a 判断素数不能被2整除的数吗？ 是否验证循环判定？ 素数不能被2整除的整数，i:整数，for m=2: fix(sqrt(i) if mod(i，m)=0 k=1； break； else k=0； 结束、编程、算法、记录像素数、主程序：合计，n=input(n=)； for i=3:n (判定像素数) (记录像素数)结束前和(前)、if k=0； prime=prime i； 结束，gd.m，prime=2；n:终点，确定素数，记录素数的程序，或： primes isprime，门巴赫预期，算法，循环：素数，验证：满足，p=primes(i )； forj=1:长度=j :长度结束，if p(j) p(k)=i m(2:3，