MATLAB简单程序大全

1、百度文库-让每个人平等地提升自我MATLAB简单程序大全求特征值特征向量A=2 3 4;1 5 9;8 5 2det(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(1/2);e=c/a;f=sqrt(1-exp(2).*cos(t)A2);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

2、.*sin(theta);comet(x,y)二重积分syms x yf=xA2*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')求xA2+yA2=1和乂八2+2八2=1所围成的体积syms x y z Rr=1;Z=sqrt(1-xA2);y0=Z;V=8*int(int(Z,y,0,y

3、0),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

4、;n3=year/400;if n1=fix(n1)&n2-=fix(n2)disp('是闰年')elseif n1=fix(n1)&n3=fix(n3) disp('是闰年') else百度文库-让每个人平等地提升自我disp'不是I闰年,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

5、,B) if nargin=0,A=0;B=100;end va=10;vb=8;vc=60;f=1;k=0;while(B-A)>0.2if f=1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);endA=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

6、(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)clcf=sym('sin(x)/x');limit(f)f1=sym('(1+a/x)Ax')limit(fl)f1=sym('(1+a/x)Ax'),x=inf;limit(fl)clcf1=sym('(1+1/x)Ax&

7、#39;);limit(f1,'inf)f1=sym('(1+1/x)Ax');limit(f1,'inf)clcf1=sym('(1+1/x)Ax');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 onwhile(B-A)>0.2if f=1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);endA=A+va*tk;B=B-vb*tk;p

8、lot(A,0,'R.',B,0,'g.'),pause(1) f=-f;k=k+1;endk,tk,A,B数学实验f=0 9 13.5;0.1 0 0;0 0.2 0x=100;100;100百度文库-让每个人平等地提升自我y=xfor n=1:3x=f*xy=y,xendyfigure(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*300y1=x1for n=1:3x1=f*x1y1=y1,x1endy1figure(4),bar(y

9、1(1,:)figure(5),bar(y1(2,:)figure(6),bar(y1(3,:)杨辉三角function Y=yang(n)if nargin=0,n=3;endY=eye(n);Y(:,1)=ones(n,1);for k=3:nY(k2k-1)=Y(k-1,1:k-2)+Y(k-1,2:k-1);End实例2寻找水仙花数一、实验内容数153有一个特殊的性质,即153=1A3+5A3+3A3,那么有没有其他这样的数呢？我们怎样利用MATLAB程序设计快速寻求这些水仙花数呢？二、实验分析把任一数各位数字的立方和求出,再求所得的各位数字的立方和,依次进行下去,总可以得到一个水仙花

10、数.如把76按这样的运算进行下去:百度文库-让每个人平等地提升自我3333333333336559,5 59979,9 791801,1 801514,333 3_3_3 3_3_3 _ _ _33_3_ _5o14 190,190 730,730 370,370370结果得到一个水仙花数370.三、实验程序寻求水仙花数经典的程序for a=1:9for b=0:9;for c=0:9;if aA3+bA3+cA3=a*100+b*10+c;sxh=a*100+b*10+cendendendend四、实验结果sxh =153, 370, 371,407.实例4住房贷款的等额本息还款法计算一、实

11、验内容当前在银行放贷工程中,住房按揭贷款占相当大的比重.这一贷款促进了房 地产事业的迅速开展,它不但改善了居民的居住条件, 也促进了建筑行业、钢铁 行业和其他相关行业的开展,与此同时,也创造了大量的就业时机.在这中间, 银行利率起了巨大的杠杆作用,过高的利率会使按揭购房减少,房价下跌；过低 的利率会使按揭贷款购房者数量增加,当房源供不应求时,受供求关系的影响, 房价将上涨.假设某购房者向银彳T贷款的金额为p0,银行的月利率为a,贷款期限为n个月,求每月的还款金额为r.二、实验分析我们用等比级数求解：第一月末贷款的本息之和为p(1) p0(1 a) r ,第二月末贷款的本息之和为：_2p(2)p(1)(1 a) r p0(1 a) r(1 a) r,第三月末贷款的本息之和为_32p(3) p(2)(1 a)r p0(1 a) r (1a) r (1a) r ,第n月末贷款的本息之和为p(n) p(n 1)(1 a)r p0(1 a)nr(1a)n1 .(1 a)1,考虑第n个月还清贷款,那么p(n)=0,求得百度文库-让每个人平等地提升自我r p0*a(1 a)n/(1 a)n 1.假设银彳T贷款20万,月利率0.465%,借款期限为10年,即为120个月,问每