MATLAB简单程序大全Word文档格式.docx

1、 theta=2*pi*（0:N）/N; r=cos（n*theta）; x=r.*cos（theta）; y=r.*sin（theta）;comet（xy）二重积分syms x y f=xA2*sin（y）;int（int（f,x90,l）9y,09pi） ezmesh（fJ05lApl）函数画图syms x;f=exp（-0.2*x）*sin（0.5*x）; ezplot（f,0,8*pi）玫瑰线theta=O:O.Ol :2*pi; r=cos（3*theta）; polar（thet%r/r）求xA2+yA2=l和xA2+zA2=l所围成的体积syms x y z Rr=l;Z=sqr

2、t（l-xA2）;yo=z；V=8*int（intZy，0，y0），x，0，l）求导数及图像f=T/（5+4*cos（x） ； subplo t（ 11）; ezplot （f） fl=diff（f）subplot（ 1，2,2）; ezplot（fl）绕X轴旋转t=（0:20）*pi/10; r=exp（2*t）.*sin（5*t）; theta=t;x=V *ones（size（t）; y=rf*cos（theta）; z=rf*sin（theta）;mesh（x，y，z） colormap（0 0 0）某年是否闰年year=input（input year:=f）; nl=year/4;

3、n2=year/100;n3=year/400;if nl=fix（nl）&n 2=fix（n2）dispC是闰年Jelseif nl=fix（nl）&n3=ftx（n3）elsedisp（不是闰年JEnd玫瑰线的绘制theta=0:0.001:2pi; rho=2*cos（3*theta）; figure（l） ;polar（theta,rho） x=rho.*cos（theta）; y=rho.*sin（theta）;figure（2） ;comet（x,y）相遇问丿function k=moto（A,B）if nargin=0, A=0; B=100; end va=10;vb=8;vc

4、=60;f=l;k=O; while（B-A）0.2 if f=ltk=（B-A）/（vb+vc）; elsetk=（B-A）/（vc+va）; endA=A+va*tk;B=B-vb*tk;f=-f; k=k+l;数学实验 x=l:100; f（x）=exp（02*x）*sin（05*x）; syms x;f（x）=exp（-0.2*x）.*sin（0.5*x）; plot（x/（x） x=0:0.1: f（x）=exp（-02*x）*sin（05*x）; plot（x,f（x）0l: zplot（x,f（x）clcf=sym（,sin（x）/x,）;limit（f）fl=sym（,（l+a

5、/x）Ax,）limit （fl）fl=sym（,（l+a/x）Ax,）,x=mf;fl=sym（,（l+l/x）Ax,）;fl=sym（t（l+l/x）Ax,）;limit（fl/xinf）Clca=o；b=ioo；va=10;plot（A，OJn/，B，OJgoJ，hold onwhile（B-A）0.2iff=lendplot（A，0，R，B，0，g），paiise（l） f=-f;k=k+l;f=0 9 13.5;0.1 0 0;0 0.2 0 x=100;100y=xfor n=l:3x=f*xy=y，xyfigure ,bar（y（l,:）figure ,bar（y（2,:figu

6、re ,bar（y（3,:p,d=eig（f）a=p（：,l）b=a/sum（a）xl=b*300yl=xlxl=f*xlyl=yl,xlyifigure（4）,bar（yl（l,:figure ,bar（yl（2,:figure ,bar （y 1（3,:杨辉三角function Y=yang（n）if nargin=0,n=3;Y=eye（n）;Y （: ,l）=ones（n,l）;for k=3:nY（k,2:k-l）=Y（k-14：k-2）+Y（k-l,2:k-l）;实例2寻找水仙花数一、 实验内容数153有一个特殊的性质，即153二3+5人3+3八3,那么有没有其他这样的数 呢？我们

7、怎样利用MATLAB程序设计快速寻求这些水仙花数呢？2.实验分析把任一数各位数字的立方和求出，再求所得的各位数字的立方和，依次进 行下去，总可以得到一个水仙花数.如把76按这样的运算进行下去: 73 + 63 = 559,53 + 53 + 9? = 979,93 + 7彳 + 9? = 1801, F + 83 + 0彳 + F = 514, 53 +13+43 =190, l3 + 93 + 03 = 730,7? + 3? + O3 = 370,3s + 73 + O3 =370 结果得到一个水仙花数370。三、实验程序（寻求水仙花数经典的程序）for a=l:9for b=0:9;fo

8、r c=0:if aA3+bA3+cA3=a*100+b*10+c;sxh=a*100+b*10+c四、实验结果sxh =153, 370,371,407。实例4住房贷款的等额本息还款法计算1.实验内容当前在银行放贷项目中，住房按揭贷款占相当大的比重。这一贷款促进了 房地产事业的迅速发展，它不但改善了居民的居住条件，也促进了建筑行业、 钢铁行业和其他相关行业的发展，与此同时，也创造了大量的就业机会。在这 中间，银行利率起了巨大的杠杆作用，过高的利率会使按揭购房减少，房价下 跌；过低的利率会使按揭贷款购房者数量增加，当房源供不应求时，受供求关 系的影响，房价将上涨。假设某购房者向银行贷款的金额为

9、po,银行的月利率为a,贷款期限为n 个月，求每月的还款金额为我们用等比级数求解：第一月末贷款的本息之和为p = 0（l + a）-厂，第二月末贷款的本息之和为：p（2） = p（l）（l + a） r= /?0（1 + a）2 一厂（1 + ci）-r t第三月末贷款的本息之和为p（3） = p（2）（l + a） r = p0（l + 一 厂（1 + a）2 r（l + d） 厂，第n月末贷款的本息之和为p（n） = p（n -1）（1 + a） _ 厂=p0（l + a1 - 4（1 + a）n +. + （1 + a） +1, 考虑第n个月还清贷款，则p（n）=0,求得r = p 0* a（l + a）”/（l+a）l假设银行贷款20万，月利率0.465%,借款期限为10年，即为120个月，问每 月应还金额和10年内共计支付多少利息。三、 实验程序n=120;p0=2e+6;a=0.00465;r=pO*a*（l+a）Aii/（l+a）An-l） % 每月还款金额interest=r*n-pO % 10年内应支付的利息四、 实验结果r =2.1785e+004,即每个月还款2178.5元。interest =6.1415e+005,即 10 应支付的利息为 61415 元。