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)二重积分symsxyf=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所围成的体积symsxyzRr=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('是闰年')elsey=xforn=1:3x=f*xy=[y,x]endyfigure(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=x1forn=1:3x1=f*x1y1=[y1,x1]endy1figure(4),bar(y1(1,:))figure(5),bar(y1(2,:))figure(6),bar(y1(3,:))杨辉三角functionY=yang(n)ifnargin==0,n=3;endY=eye(n);Y(:,1)=ones(n,1);fork=3:nY(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:9forb=0:9;forc=0:9;ifa^3+b^3+c^3==a*100+b*10+c;sxh=a*100+b*10+cendendendend四、实验结果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)%每月还款金额int