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;
二重积分
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('就是闰年')else
disp('不是闰年')end
玫瑰线的绘制
相遇问题
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==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
数学实验
symsx;x=1:100;
f(x)=exp(-0.2*x).*sin(0.5*x);symsx;x=1:100;
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.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))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'),holdonwhile(b-a)>0.2iff==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;endk,tk,a,b 数学实验
f=[0913.5;0.100;00.20]x=[100;100;100]
y=x
forn=1:3x=f*xy=[y,x]endy
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*300y1=x1forn=1:3x1=f*x 1y1=[y1,x1]endy1
figure(4),bar(y1(1,:))figure(5),bar(y1(2,:))figure(6),bar(y1(3,:))
杨辉三角
functiony=yang(n)ifnargin==0,n=3;end
y=eye(n);y(:,1)=ones(n,1);fork=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按这样的运算进行下去:
