MATLAB简单程序大全备课讲稿.doc

MATLAB简单程序大全

精品文档
MATLAB简单程序大全

求特征值特征向量
A=[2 3 4;1 5 9;8 5 2]
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)

二重积分
syms x y
f=x^2*sin(y);
int(int(f,x,0,1),y,0,pi)
ezmesh(f,[0,1,0,pi])

函数画图
syms x;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所围成的体积
syms x y z R
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([0 0 0])

某年是否闰年
year=input('input year:=');
n1=year/4;
n2=year/100;
n3=year/400;
if n1==fix(n1)&n2~=fix(n2)
disp('是闰年')
elseif n1==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)

相遇问题
function k=moto(A,B)
if nargin==0,A=0;B=100;end
va=10;vb=8;vc=60;
f=1;k=0;
while(B-A)>0.2
if f==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

数学实验
syms x; x=1:100;
f(x)=exp(-0.2*x).*sin(0.5*x);
syms x; x=1:100;
f(x)=exp(-0.2*x).*sin(0.5*x);
plot(x,f(x))
syms x; x=0:0.1:2*pi;
f(x)=exp(-0.2*x)*sin(0.5*x);
plot(x,f(x))
syms x; x=0:0.1:2*pi;
f(x)=exp(-0.2*x).*sin(0.5*x);
plot(x,f(x))
syms x; 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'),hold on
while(B-A)>0.2
if f==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=[0 9 13.5;0.1 0 0;0 0.2 0]
x=[100;100;100]
y=x
for n=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
for n=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,:))

杨辉三角
function Y=yang(n)
if nargin==0,n=3;end
Y=eye(n);Y(:,1)=ones(n,1);
for k=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。
三、实验程序（寻求水仙花数经典的程序）
for a=1:9
for b=0:9;
for c=0:9;
if a^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元。

收集于网络，如有侵权请联系管理员删除