MATLAB简单程序大全.docx

MATLAB简单程序大全.docx

《MATLAB简单程序大全.docx》由会员分享，可在线阅读，更多相关《MATLAB简单程序大全.docx（10页珍藏版）》请在冰豆网上搜索。

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'）

求x^2+y^2=1和x^2+z^2=1所围成的体积

symsxyzR

r=1;

Z=sqrt（1-x^2）;

y0=Z;

V=8*int（int（Z,y,0,y0）,x,0,1）

求导数及图像

f='1/（5+4*cos（x））';

subplot（1,2,1）;ezplot（f）

f1=diff（f）

subplot（1,2,2）;ezplot（f1）

绕x轴旋转

t=（0:

20）*pi/10;

r=exp（-.2*t）.*sin（.5*t）;

theta=t;

x=t'*ones（size（t））;

y=r'*cos（theta）;

z=r'*sin（theta）;

mesh（x,y,z）

colormap（[000]）

某年是否闰年

year=input（'inputyear:

='）;

n1=year/4;

n2=year/100;

n3=year/400;

ifn1==fix（n1）&n2~=fix（n2）

disp（'是闰年'）

elseifn1==fix（n1）&n3==fix（n3）

disp（'是闰年'）

else

disp（'不是闰年'）

End

玫瑰线的绘制

theta=0:

0.001:

2*pi;

rho=2*cos（3*theta）;

figure

（1）;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;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元。