MATLAB简单程序大全Word格式文档下载.doc
《MATLAB简单程序大全Word格式文档下载.doc》由会员分享,可在线阅读,更多相关《MATLAB简单程序大全Word格式文档下载.doc(7页珍藏版)》请在冰豆网上搜索。
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)
else
disp('
不是闰年'
End
玫瑰线的绘制
0.001:
rho=2*cos(3*theta);
figure
(1);
polar(theta,rho)
x=rho.*cos(theta);
y=rho.*sin(theta);
figure
(2);
相遇问题
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;
数学实验
x=1:
100;
f(x)=exp(-0.2*x).*sin(0.5*x);
plot(x,f(x))
x=0:
0.1:
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)
),x=inf;
(1+1/x)^x'
limit(f1,'
inf'
x'
inf)
Clc
A=0;
plot(A,0,'
ro'
B,0,'
go'
),holdon
plot(A,0,'
R.'
g.'
),pause
(1)
k=k+1;
k,tk,A,B
f=[0913.5;
0.100;
00.20]
x=[100;
100]
y=x
forn=1:
3
x=f*x
y=[y,x]
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
x1=f*x1
y1=[y1,x1]
y1
figure(4),bar(y1(1,:
figure(5),bar(y1(2,:
figure(6),bar(y1(3,:
杨辉三角
functionY=yang(n)
ifnargin==0,n=3;
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);
实例2寻找水仙花数
一、实验内容
数153有一个特殊的性质,即153=1^3+5^3+3^3,那么有没有其他这样的数呢?
我们怎样利用MATLAB程序设计快速寻求这些水仙花数呢?
二、实验分析
把任一数各位数字的立方和求出,再求所得的各位数字的立方和,依次进行下去,总可以得到一个水仙花数.如把76按这样的运算进行下去:
结果得到一个水仙花数370。
三、实验程序(寻求水仙花数经典的程序)
fora=1:
9
forb=0:
9;
forc=0:
ifa^3+b^3+c^3==a*100+b*10+c;
sxh=a*100+b*10+c
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元。