經(jīng)模03：程序設(shè)計

1、主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件MatlabMatlab程序設(shè)計程序設(shè)計經(jīng)濟模型經(jīng)濟模型與與Matlab應(yīng)用應(yīng)用主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件1、基礎(chǔ)、基礎(chǔ)一、一、Matlab程序程序vMatlab m-文件v命令文件運行命令窗口：輸入M文件的文件名M文件窗口：F5v函數(shù)文件function out1, out2, . = myfun(in1, in2, .)函數(shù)表達式使用：命令窗口：文件名（）主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件簡單條件語句多條件條件語句 if (條件式) 語句組 end if (條件式1) 語句組1 elseif (

2、條件式2) 語句組2 elseif (條件式3) 語句組3 else . . end多選擇條件語句 if (條件式) 語句組 1 else 語句組 2 endvif -else-end分支結(jié)構(gòu)分支結(jié)構(gòu)2、條件語句、條件語句主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件v例1函數(shù)100010 xyxxif (x0) y=-1elseif (x=0) y=0else y=1endv編程1l01.mv編程2函數(shù)？?另：x？輸入一個x的值,輸出符號函數(shù)y的值主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件varray 是一個矩陣：循環(huán)列一般：數(shù)組一般：循環(huán)變量=初值:步長:終值默認(rèn)：步長=1

3、第一類循環(huán)語句第二類循環(huán)語句結(jié)構(gòu)for 循環(huán)變量= array 循環(huán)體語句組endwhile (條件式) 循環(huán)體語句組 end3、循環(huán)語句、循環(huán)語句主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件 生成一個6階矩陣，使其主對角線上元素皆為1，與主對角線相鄰元素皆為2，其余皆為0。120000212000021200002120000212000021av例3for 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結(jié)果：v編程l04.m主講人：孫云龍主講人：孫云龍數(shù)

4、學(xué)建模課件數(shù)學(xué)建模課件v例4 求自然數(shù)的前n項和v編程結(jié)果：n=100 Sum=5050n=input(n=)sum=0; k=1;while k=n sum=sum+k; k=k+1;endsuml05.m好不好？好不好？怎么改？怎么改？1234主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件v指令4、程序流控制指令、程序流控制指令 break 跳出循環(huán) 一般和 if結(jié)構(gòu) 結(jié)合使用 continue 結(jié)束本次循環(huán) 繼續(xù)進行下次循環(huán) return 終止當(dāng)前指令 返回上一級或等待鍵盤輸入命令 pause 暫停 按任意鍵繼續(xù) input 輸入 提示鍵盤輸入主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課

5、件數(shù)學(xué)建模課件v編程v例6：連續(xù)奇數(shù)求和clearsum1=0;for i=1:100 n=2*i-1; if sum11000 sum1=sum1+n; else break endendsum1,nl07.mv程序?qū)Σ粚Γ繌?開始一直到和達到1000為止問：加了多少項主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件v工資、薪金所得適用級年終獎稅率1(0,18000)32(18000,54000)103(54000,108000)204(108000,420000)255(420000,660000)306(660000,960000)357(960000,)45二、初等模型二、初等模型

6、1、個人所得稅問題、個人所得稅問題主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件問題問題v1、納稅額計算函數(shù)條件判斷程序程序t1.m t2.mfunction 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級數(shù)月應(yīng)納稅所得額（基數(shù)3500元）稅率(%)10-1500321500-45001034500-90002049000-3500025535000-5500030655000-8000035780000-4

7、5主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件問題問題v2、分配月收入、年獎金合理避稅為什么如何分程序程序v搜索：定步長for-endwhile -endv紀(jì)錄最優(yōu)點if-endt3.m算法算法function min,k1 =tax3(x)k=0;min=t1(x);k1=0;while kd1 min=d1; k1=k; end k=k+1;end另：程序有問題？mint,k1=t3(5000)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件v小李夫婦買房需向銀行貸款6萬元，按月分期等額償還房屋抵押貸款，月利率是0.01，貸款期為25年。小李夫婦每月能有900元的結(jié)余。v（1）

8、小李夫婦是否有無能力買房？月供多少？v（2）有一則廣告：本公司能幫您提前三年還清貸款，只要每半月還錢一次，但由于文書工作多了，要求您先付三個月的錢。是否劃算？ v（3）小李夫婦若將結(jié)余全部用來還貸，多長時間還清房貸？ 2、現(xiàn)金流計算、現(xiàn)金流計算主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件v規(guī)則現(xiàn)金流v基本公式 1(1)(1)TtTtPFVPVraterate1(1)(1)(1)TTPFVrateraterate現(xiàn)值現(xiàn)金流終值利率主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件v月供？ 631.9345 v提前三年，半月還錢一次，先付三個月，是否劃算？現(xiàn)值60550v貸款時間2問題*1(1)(1)(1)TT