MATLAB实验报告3

《MATLAB实验报告3》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、年级、专业09级数学与应用数学姓名卢海宽学号150142009040名单序号117实验时间2011年11月8日使用设备、软件PC,MATLAB注：实验报告的最后一部分是实验小结与收获实验（3）：差分方程（组）（DifferenceEquation(s)）1.Yourparentsareconsideringa30-year$100,000mortgagethatcharges1%interesteachmonth.Ifthemortgage(loan)istobepaidoffafter360pay

2、ments，whatisthenmonthlypaymentsp?Hint:Ifanrepresentstheamountowedafternmonths,whatarea0anda360?Experimentbybuildingnumericalsolutionstofindavalueofpthatworks.设a(0)代表贷款总额，a(k)为第k个月的欠款数，p为月还款数，r为月利率。于是有a(k)=(1+r)*a(k-1)-p=(1+r)^2*a(k-2)-(1+r)p-p=…=(1+r)^

3、k*a(0)-p*[(1+r)^k-1]/r所以有p=[（1+r）^360*r*a(0)]/[(1+r)^360-1]此处k=0,1,2…360a(0)=100000r=0.01p=((1.01)^360*0.01*100000)/((1.01)^360-1)p=1.0286e+003求第k个月欠款数（图）k=0:1:360;y=(1.01).^k*100000-1028.6*((1.01).^k-1)/0.01;plot(k,y,'ro');gridon;2.Readthecoursewareon

4、thebattleofTrafalgarin1805.a)Formulateasystemofdifferenceequationstomodelthenumberofshipspossessedbyeachforce.AssumetheFrench-Spanishstartwith33shipsandtheBritishstartwith272011秋数学实验基础实验报告（3）：差分方程（组）Page7of7年级、专业09级数学与应用数学姓名卢海宽学号150142009040名单序号117实验时间

5、2011年11月8日使用设备、软件PC,MATLAB注：实验报告的最后一部分是实验小结与收获实验（3）：差分方程（组）（DifferenceEquation(s)）shipsandeachsidelost10%ofthenumberofshipsoftheopposingforce.设f(i)表示第i阶段西法联军的船只数b(i)表示第i阶段英军的船只数f(1)=33;b(1)=27;i=1;whilef(i)>=1&&b(i)>=1i=i+1;f(i)=f(i-1)-0.1*b(i-1);b(i)=

6、b(i-1)-0.1*f(i-1);endplot(f,'ro-');holdonplot(b,'b*-');holdoff;gridon;legend('联军船只数','英军船只数')a)Analternativestrategyforovercomingasuperiorforceistoincreasethetechnologyemployedbytheinferiorforce.SupposethattheBritishshipswereequippedwithsuperiorweaponr

7、y,andthattheFrench-Spanishsufferedlossesequalto15%ofthenumberofshipsoftheopposingforce,whereastheBritishsufferedcasualtiesequalto5%oftheopposingforce.Buildanumericalsolutiontodeterminewhowinsundernewassumptioninahead-onengagement.f(1)=33;b(1)=27;i=1;wh

8、ilef(i)>=1&&b(i)>=1i=i+1;2011秋数学实验基础实验报告（3）：差分方程（组）Page7of7年级、专业09级数学与应用数学姓名卢海宽学号150142009040名单序号117实验时间2011年11月8日使用设备、软件PC,MATLAB注：实验报告的最后一部分是实验小结与收获实验（3）：差分方程（组）（DifferenceEquation(s)）f(i)=f(i-1)-0.15*b(i-1);b(i)=b(i-1)-0.05*f(i-1);en