MATLAB简单程序大全.docx
- 文档编号:865105
- 上传时间:2022-10-13
- 格式:DOCX
- 页数:7
- 大小:17.31KB
MATLAB简单程序大全.docx
《MATLAB简单程序大全.docx》由会员分享,可在线阅读,更多相关《MATLAB简单程序大全.docx(7页珍藏版)》请在冰豆网上搜索。
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=(aA2-cA2)A(l/2);
e=c/a;
f=sqrt(1-exp
(2).*cos(t)A2);
l=int(f£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)
二重积分
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)
dispC是闰年J
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;
End
数学实验
symsx;x=l:
100;f(x)=exp(・0・2*x)・*sin(0・5*x);symsx;x=l:
100;
f(x)=exp(-0.2*x).*sin(0.5*x);plot(x/(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・l:
2*pi;
f(x)=exp(-0.2*x).*sin(0.5*x);plot(x/(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)
fl=sym(,(l+a/x)Ax,)
limit(fl)
fl=sym(,(l+a/x)Ax,),x=mf;
limit(fl)
clc
fl=sym(,(l+l/x)Ax,);
fl=sym(t(l+l/x)Ax,);
clc
fl=sym(,(l+l/x)Ax,);
limit(fl/x\inf)
Clc
相遇问丿
a=o;b=ioo;
va=10;vb=8;vc=60;
f=l;k=O;
plot(A,OJn/,B,OJgoJ,holdon
while(B-A)>0.2
iff==l
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・'),paiise(l)f=-f;k=k+l;
end
数学实验
f=[0913.5;0.100;00.20]x=[100;100;100]
y=x
forn=l:
3
x=f*x
y=[y,x]
end
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
forn=l:
3
xl=f*xl
yl=[yl,xl]
end
yi
figure(4),bar(yl(l,:
))
figure⑸,bar(yl(2,:
))
figure⑹,bar(y1(3,:
))
杨辉三角
functionY=yang(n)
ifnargin==0,n=3;end
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);
End
实例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:
9;
ifaA3+bA3+cA3==a*100+b*10+c;
sxh=a*100+b*10+c
end
end
end
end
四、实验结果
sxh=153,370,371,407。
实例4住房贷款的等额本息还款法计算
1.实验内容
当前在银行放贷项目中,住房按揭贷款占相当大的比重。
这一贷款促进了房地产事业的迅速发展,它不但改善了居民的居住条件,也促进了建筑行业、钢铁行业和其他相关行业的发展,与此同时,也创造了大量的就业机会。
在这中间,银行利率起了巨大的杠杆作用,过高的利率会使按揭购房减少,房价下跌;过低的利率会使按揭贷款购房者数量增加,当房源供不应求时,受供求关系的影响,房价将上涨。
假设某购房者向银行贷款的金额为po,银行的月利率为a,贷款期限为n个月,求每月的还款金额为
2.实验分析
我们用等比级数求解:
第一月末贷款的本息之和为
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元。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 简单 程序 大全
![提示](https://static.bdocx.com/images/bang_tan.gif)