MATLAB簡(jiǎn)單程序大全

1、MATLAB簡(jiǎn)單程序大全求特征值特征向量A=2 3 4;1 5 9;8 5 2det(A)A'rank(A)inv(A)rref(A)eig(A)%求特征值和特征向量衛(wèi)星運(yùn)行問題h=200,H=51000,R=6378;a=(h+H+2*R)/2;c=(H-h)/2;b=(a2-c2)(1/2);e=c/a;f=sqrt(1-exp(2).*cos(t)2);l=int(f,t,0,pi/2)L=4*a.*l動(dòng)態(tài)玫瑰線n=3;N=10000;theta=2*pi*(0:N)/N;r=cos(n*theta);x=r.*cos(theta);y=r.*sin(theta);comet(x

2、,y)二重積分syms x yf=x2*sin(y);int(int(f,x,0,1),y,0,pi)ezmesh(f,0,1,0,pi)函數(shù)畫圖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')求x2+y2=1和x2+z2=1所圍成的體積syms x y z Rr=1;Z=sqrt(1-x2);y0=Z;V=8*int(int(Z,y,0,y0),x,0,1)求導(dǎo)數(shù)及圖像f='1/(5+4*cos(x)'su

3、bplot(1,2,1);ezplot(f)f1=diff(f)subplot(1,2,2);ezplot(f1)繞x軸旋轉(zhuǎn)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)

4、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;endva=10;vb=8;vc=60;f=1;k=0; i

5、f 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數(shù)學(xué)實(shí)驗(yàn)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*

6、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)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');limit(f1,'inf')clcf1=s

7、ym('(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 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;endk,tk,A,B數(shù)學(xué)實(shí)驗(yàn)f=0 9 13.5;0.

8、1 0 0;0 0.2 0x=100;100;100y=xfor n=1:3 x=f*x y=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:3 x1=f*x1 y1=y1,x1endy1figure(4),bar(y1(1,:)figure(5),bar(y1(2,:)figure(6),bar(y1(3,:)楊輝三角function Y=yang(n)if nargin=0,n=3;endY=eye(n);Y

9、(:,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實(shí)例2 尋找水仙花數(shù)一、實(shí)驗(yàn)內(nèi)容數(shù)153有一個(gè)特殊的性質(zhì),即153=13+53+33,那么有沒有其他這樣的數(shù)呢？我們?cè)鯓永肕ATLAB程序設(shè)計(jì)快速尋求這些水仙花數(shù)呢？二、實(shí)驗(yàn)分析把任一數(shù)各位數(shù)字的立方和求出 ,再求所得的各位數(shù)字的立方和 ,依次進(jìn)行下去 ,總可以得到一個(gè)水仙花數(shù).如把76按這樣的運(yùn)算進(jìn)行下去: 結(jié)果得到一個(gè)水仙花數(shù)370。三、實(shí)驗(yàn)程序?qū)で笏苫〝?shù)經(jīng)典的程序for a=1:9 for b=0:9; for c=0:9; if a3+b3+c3=a*

10、100+b*10+c; sxh=a*100+b*10+c end end endend四、實(shí)驗(yàn)結(jié)果sxh =153，370, 371, 407。實(shí)例4 住房貸款的等額本息還款法計(jì)算一、實(shí)驗(yàn)內(nèi)容當(dāng)前在銀行放貸項(xiàng)目中，住房按揭貸款占相當(dāng)大的比重。這一貸款促進(jìn)了房地產(chǎn)事業(yè)的迅速發(fā)展，它不但改善了居民的居住條件，也促進(jìn)了建筑行業(yè)、鋼鐵行業(yè)和其他相關(guān)行業(yè)的發(fā)展，與此同時(shí)，也創(chuàng)造了大量的就業(yè)時(shí)機(jī)。在這中間，銀行利率起了巨大的杠桿作用，過高的利率會(huì)使按揭購房減少，房?jī)r(jià)下跌；過低的利率會(huì)使按揭貸款購房者數(shù)量增加，當(dāng)房源供不應(yīng)求時(shí)，受供求關(guān)系的影響，房?jī)r(jià)將上漲。假設(shè)某購房者向銀行貸款的金額為p0，銀行的月利率為a，貸款期限為 n個(gè)月，求每月的還款金額為r.二、實(shí)驗(yàn)分析 我們用等比級(jí)數(shù)求解： 第一月末貸款的本息之和為， 第二月末貸款的本息之和為：，第三月末貸款的本息之和為，第n月末貸款的本息之和為，考慮第n個(gè)月還清貸款，則p(n)=0,求得.假設(shè)銀行貸款20萬，月利率0.465%，借款期限為10年，即為120個(gè)月，問每月應(yīng)還金額和10年內(nèi)共計(jì)支付多少利息。三、實(shí)驗(yàn)程