资源描述：
《MATLAB简单程序大全》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、.word可编辑.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动态玫瑰线.专业.专注..word可编辑.n=3;N=10000;the
2、ta=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').专业.
3、专注..word可编辑.求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))
4、;y=r'*cos(theta);z=r'*sin(theta);mesh(x,y,z)colormap([000]).专业.专注..word可编辑.某年是否闰年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('是闰年')elsedisp('不是闰年')End玫瑰线的绘制theta=0:0
5、.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).专业.专注..word可编辑.相遇问题functionk=moto(A,B)ifnargin==0,A=0;B=100;endva=10;vb=8;vc=60;f=1;k=0;while(B-A)>0.2iff==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(
6、vc+va);endA=A+va*tk;B=B-vb*tk;f=-f;k=k+1;End数学实验symsx;x=1:100;f(x)=exp(-0.2*x).*sin(0.5*x);symsx;x=1:100;.专业.专注..word可编辑.f(x)=exp(-0.2*x).*sin(0.5*x);plot(x,f(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.1:2*pi;f(x)=exp(-0
7、.2*x).*sin(0.5*x);plot(x,f(x))symsx;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)^x')limit(f1)f1=sym('(1+a/x)^x'),x=inf;limit(f1)clcf1=sym('(1+1/x)^x');limit(f1,'inf')f1=sym('(1+1/x)^x');.专业.专注..wo
8、rd可编辑.limit(f1,'inf')clcf1=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'),holdonwhile(B-A)>0.2iff==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);endA=A+va*tk;B=B-vb*tk;plot(A,0,'R.',B,0,'g.'),pause
