MATLAB简单程序大全Word文档格式.docx
《MATLAB简单程序大全Word文档格式.docx》由会员分享，可在线阅读，更多相关《MATLAB简单程序大全Word文档格式.docx（7页珍藏版）》请在冰豆网上搜索。
theta=2*pi*（0:
N）/N;
r=cos（n*theta）;
x=r.*cos（theta）;
y=r.*sin（theta）;
comet（x^y）
二重积分
symsxyf=xA2*sin（y）;
int（int（f,x90,l）9y,09pi）ezmesh（fJ05lAp>
l）
函数画图
symsx;
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所围成的体积
symsxyzR
r=l;
Z=sqrt（l-xA2）;
yo=z；
V=8*int（intZy，0，y0），x，0，l）
求导数及图像
f=T/（5+4*cos（x））‘；
subplot（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（[000]）
某年是否闰年
year=input（inputyear:
=f）;
nl=year/4;
n2=year/100;
n3=year/400;
ifnl==fix（nl）&
n2~=fix（n2）
dispC是闰年J
elseifnl==fix（nl）&
n3==ftx（n3）
else
disp（'
不是闰年J
End
玫瑰线的绘制
theta=0:
0.001:
2^pi;
rho=2*cos（3*theta）;
figure（l）;
polar（theta,rho）x=rho.*cos（theta）;
y=rho.*sin（theta）;
figure
（2）;
comet（x,y）
相遇问丿
functionk=moto（A,B）
ifnargin==0,A=0;
B=100;
endva=10;
vb=8;
vc=60;
f=l;
k=O;
while（B-A）>
0.2iff==l
tk=（B-A）/（vb+vc）;
else
tk=（B-A）/（vc+va）;
end
A=A+va*tk;
B=B-vb*tk;
f=-f;
k=k+l;
数学实验
x=l:
100;
f（x）=exp（・0・2*x）・*sin（0・5*x）;
symsx;
f（x）=exp（-0.2*x）.*sin（0.5*x）;
plot（x/（x））
x=0:
0.1:
f（x）=exp（-0・2*x）*sin（0・5*x）;
plot（x,f（x））
0・l:
zplot（x,f（x））
clc
f=sym（,sin（x）/x,）;
limit（f）
fl=sym（,（l+a/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/x\inf）
Clc
a=o；
b=ioo；
va=10;
plot（A，OJn/，B，OJgoJ，holdon
while（B-A）>
0.2
iff==l
end
plot（A，0，'
R・‘，B，0，'
g・'
），paiise（l）f=-f;
k=k+l;
f=[0913.5;
0.100;
00.20]x=[100;
100]
y=x
forn=l:
3
x=f*x
y=[y，x]
y
figure⑴,bar（y（l,:
））
figure⑵,bar（y（2,:
figure⑶,bar（y（3,:
[p,d]=eig（f）
a=p（：
l）
b=a/sum（a）
xl=b*300
yl=xl
xl=f*xl
yl=[yl,xl]
yi
figure（4）,bar（yl（l,:
figure⑸,bar（yl（2,:
figure⑹,bar（y1（3,:
杨辉三角
functionY=yang（n）
ifnargin==0,n=3;
Y=eye（n）;
Y（:
l）=ones（n,l）;
fork=3:
n
Y（k,2:
k-l）=Y（k-14：
k-2）+Y（k-l,2:
k-l）;
实例2寻找水仙花数
一、实验内容
数153有一个特殊的性质，即153二"
3+5人3+3八3,那么有没有其他这样的数呢？
我们怎样利用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。
三、实验程序（寻求水仙花数经典的程序）
fora=l:
9
forb=0:
9;
forc=0:
ifaA3+bA3+cA3==a*100+b*10+c;
sxh=a*100+b*10+c
四、实验结果
sxh=153,370,371,407。
实例4住房贷款的等额本息还款法计算
1.实验内容
当前在银行放贷项目中，住房按揭贷款占相当大的比重。
这一贷款促进了房地产事业的迅速发展，它不但改善了居民的居住条件，也促进了建筑行业、钢铁行业和其他相关行业的发展，与此同时，也创造了大量的就业机会。
在这中间，银行利率起了巨大的杠杆作用，过高的利率会使按揭购房减少，房价下跌；
过低的利率会使按揭贷款购房者数量增加，当房源供不应求时，受供求关系的影响，房价将上涨。
假设某购房者向银行贷款的金额为po,银行的月利率为a,贷款期限为n个月，求每月的还款金额为
我们用等比级数求解：
第一月末贷款的本息之和为
p⑴=〃0（l+a）-厂，
第二月末贷款的本息之和为：
p
（2）=p（l）（l+a）—r=/?
0（1+a）2一厂（1+ci）-rt
第三月末贷款的本息之和为
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+a}1-4（1+a）n~}+...+（1+a）+1],考虑第n个月还清贷款，则p（n）=0,求得
r=p0*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元。
升级会员 