MATLAB简单程序大全.docx
《MATLAB简单程序大全.docx》由会员分享，可在线阅读，更多相关《MATLAB简单程序大全.docx（10页珍藏版）》请在冰点文库上搜索。
MATLAB简单程序大全
MATLAB简单程序大全
MATLAB简单程序大全
求特征值特征向量
A=[234;159;852]
det（A）
A'
rank（A）
inv（A）
rref（A）
eig（A）%求特征值和特征向量
卫星运行问题
h=200,H=51000,R=6378;
a=（h+H+2*R）/2;
c=（H-h）/2;
b=（a^2-c^2）^（1/2）;
e=c/a;
f=sqrt（1-exp
（2）.*cos（t）^2）;
l=int（f,t,0,pi/2）
L=4*a.*l
动态玫瑰线
n=3;N=10000;
theta=2*pi*（0:
N）/N;
r=cos（n*theta）;
x=r.*cos（theta）;
y=r.*sin（theta）;
comet（x,y）
二重积分
symsxy
f=x^2*sin（y）;
int（int（f,x,0,1）,y,0,pi）
ezmesh（f,[0,1,0,pi]）
函数画图
symsx;f=exp（-0.2*x）*sin（0.5*x）;
ezplot（f,[0,8*pi]）
玫瑰线
theta=0:
0.01:
2*pi;
r=cos（3*theta）;
polar（theta,r,'r'）
相遇问题
functionk=moto（A,B）
ifnargin==0,A=0;B=100;end
va=10;vb=8;vc=60;
f=1;k=0;
while（B-A）>0.2
iff==1
tk=（B-A）/（vb+vc）;
else
tk=（B-A）/（vc+va）;
end
A=A+va*tk;
B=B-vb*tk;
f=-f;
k=k+1;
End
数学实验
symsx;x=1:
100;
f（x）=exp（-0.2*x）.*sin（0.5*x）;
symsx;x=1:
100;
f（x）=exp（-0.2*x）.*sin（0.5*x）;
plot（x,f（x））
symsx;x=0:
0.1:
2*pi;
f（x）=exp（-0.2*x）*sin（0.5*x）;
plot（x,f（x））
symsx;x=0:
0.1:
2*pi;
f（x）=exp（-0.2*x）.*sin（0.5*x）;
plot（x,f（x））
symsx;x=0:
0.1:
2*pi;
f（x）=exp（-0.2*x）.*sin（0.5*x）;
zplot（x,f（x））
clc
f=sym（'sin（x）/x'）;
limit（f）
f1=sym（'（1+a/x）^x'）
limit（f1）
f1=sym（'（1+a/x）^x'）,x=inf;
limit（f1）
clc
f1=sym（'（1+1/x）^x'）;
limit（f1,'inf'）
f1=sym（'（1+1/x）^x'）;
limit（f1,'inf'）
clc
f1=sym（'（1+1/x）^x'）;
limit（f1,'x',inf）
Clc
相遇问题
A=0;B=100;
va=10;vb=8;vc=60;
f=1;k=0;
plot（A,0,'ro',B,0,'go'）,holdon
while（B-A）>0.2
iff==1
tk=（B-A）/（vb+vc）;
else
tk=（B-A）/（vc+va）;
end
A=A+va*tk;
B=B-vb*tk;
plot（A,0,'R.',B,0,'g.'）,pause
（1）
f=-f;k=k+1;
end
k,tk,A,B
数学实验
f=[0913.5;0.100;00.20]
x=[100;100;100]
y=x
forn=1:
3
x=f*x
y=[y,x]
end
y
figure
（1）,bar（y（1,:
））
figure
（2）,bar（y（2,:
））
figure（3）,bar（y（3,:
））
[p,d]=eig（f）
a=p（:
1）
b=a/sum（a）
x1=b*300
y1=x1
forn=1:
3
x1=f*x1
y1=[y1,x1]
end
y1
figure（4）,bar（y1（1,:
））
figure（5）,bar（y1（2,:
））
figure（6）,bar（y1（3,:
））
杨辉三角
functionY=yang（n）
ifnargin==0,n=3;end
Y=eye（n）;Y（:
1）=ones（n,1）;
fork=3:
n
Y（k,2:
k-1）=Y（k-1,1:
k-2）+Y（k-1,2:
k-1）;
End
实例2寻找水仙花数
一、实验内容
数153有一个特殊的性质,即153=1^3+5^3+3^3,那么有没有其他这样的数呢？
我们怎样利用MATLAB程序设计快速寻求这些水仙花数呢？
二、实验分析
把任一数各位数字的立方和求出,再求所得的各位数字的立方和,依次进行下去,总可以得到一个水仙花数.如把76按这样的运算进行下去:
结果得到一个水仙花数370。
三、实验程序（寻求水仙花数经典的程序）
fora=1:
9
forb=0:
9;
forc=0:
9;
ifa^3+b^3+c^3==a*100+b*10+c;
sxh=a*100+b*10+c
end
end
end
end
四、实验结果
sxh=153，370,371,407。
实例4住房贷款的等额本息还款法计算
一、实验内容
当前在银行放贷项目中，住房按揭贷款占相当大的比重。
这一贷款促进了房地产事业的迅速发展，它不但改善了居民的居住条件，也促进了建筑行业、钢铁行业和其他相关行业的发展，与此同时，也创造了大量的就业机会。
在这中间，银行利率起了巨大的杠杆作用，过高的利率会使按揭购房减少，房价下跌；过低的利率会使按揭贷款购房者数量增加，当房源供不应求时，受供求关系的影响，房价将上涨。
假设某购房者向银行贷款的金额为p0，银行的月利率为a，贷款期限为n个月，求每月的还款金额为r.
二、实验分析
我们用等比级数求解：
第一月末贷款的本息之和为
，
第二月末贷款的本息之和为：
，
第三月末贷款的本息之和为
，
第n月末贷款的本息之和为
，
考虑第n个月还清贷款，则p（n）=0,求得
.
假设银行贷款20万，月利率0.465%，借款期限为10年，即为120个月，问每月应还金额和10年内共计支付多少利息。
三、实验程序
n=120;
p0=2e+6;
a=0.00465;
r=p0*a*（1+a）^n/（（1+a）^n-1）%每月还款金额
interest=r*n-p0%10年内应支付的利息
四、实验结果
r=2.1785e+004,即每个月还款2178.5元。
interest=6.1415e+005，即10应支付的利息为61415元。
升级会员 