MATLAB簡單程序大全

1、百度文庫-讓每個人平等地提升自我1MATLAB 簡單程序大全求特征值特征向量A=2 3 4;1 5 9;8 5 2det(A)Aran k(A)inv(A)rref(A)eig(A)%求特征值和特征向量衛(wèi)星運行問題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=i nt(f,t,0,pi/2)L=4*a.*l動態(tài)玫瑰線n=3;N=10000;theta=2*pi*(0:N)/N;r=cos( n*theta);x=r.*cos(theta);y=r

2、.*s in (theta);comet(x,y)二重積分syms x yf=xA2*si n(y);in t(i nt(f,x,0,1),y,0,pi)ezmesh(f,0,1,0,pi)函數(shù)畫圖syms x; f=exp(-0.2*x)*si n(0.5*x);ezplot(f,0,8*pi)百度文庫-讓每個人平等地提升自我2玫瑰線theta=0:0.01:2*pi;r=cos(3*theta);polar(theta,r,r)求x2+yA2=1和xA2+zH=1所圍成的體積syms x y z Rr=1;Z=sqrt(1-xA2);y0=Z;V=8*i nt(i nt(Z,y,O,yO)

3、,x,O,1)求導(dǎo)數(shù)及圖像f=1/(5+4*cos(x);subplot(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).*si n(.5*t);theta=t;x=t* on es(size(t);y=r*cos(theta);z=r*s in (theta);mesh(x,y,z)colormap(0 0 0)某年是否閏年year=i nput(i nput year:=);n 1=year/4;n2=year/100;n3=year/400;if n1=fix (n 1

4、)&n 2=fix( n2)disp(是閏年)elseif n1=fix( n1)&n 3=fix( n3)disp(是閏年)elsedisp(不是閏年)End百度文庫-讓每個人平等地提升自我3玫瑰線的繪制theta=0:0.001:2*pi; rho=2*cos(3*theta);figure(1);polar(theta,rho) x=rho.*cos(theta); y=rho.*s in (theta);figure(2);comet(x,y)相遇問題function k=moto(A,B)if nargi n=0,A=0;B=100;e nd va=10;vb=8;v

5、c=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數(shù)學實驗syms x; x=1:100;f(x)=exp(-0.2*x).*si n(0.5*x); syms x; x=1:100;f(x)=exp(-0.2*x).*si n(0.5*x); plot(x,f(x)syms x; x=0:0.1:2*pi;f(x)=exp(-0.2*x)*si n(0.5*x); plot(x,f(x)syms x; x=0:0.1:2*p

6、i;f(x)=exp(-0.2*x).*si n(0.5*x);百度文庫-讓每個人平等地提升自我4plot(x,f(x)syms x; x=0:0.1:2*pi;f(x)=exp(-0.2*x).*si n(0.5*x); zplot(x,f(x) clc f=sym(si n( x)/x);limit(f)f1=sym(1+a/x)Ax)limit(fl) f1=sym(1+a/x)Ax),x=i nf;limit(fl)clcf1=sym(1+1/x)Ax);limit(f1,i nf)f1=sym(1+1/x)Ax);limit(f1,i nf)clcf1=sym(1+1/x)Ax);l

7、imit(f1,x,i nf)Clc相遇問題A=0;B=100; va=10;vb=8;vc=60;f=1;k=0;plot(A,0,ro,B,0,go),hold on while(B-A)0.2if f=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(1) f=-f;k=k+1;endk,tk,A,B數(shù)學實驗f=0 9 13.5;0.1 0 0;0 0.2 0 x=100;100;100y=xfor n=1:3x=f*xy=y,xendyfigure(1),bar

8、(y(1,:)百度文庫-讓每個人平等地提升自我5figure(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:3x仁f*x1y1= y1,x1endyifigure(4),bar(y1(1,:)figure(5),bar(y1(2,:)figure(6),bar(y1(3,:)楊輝三角function Y=ya ng(n)if nargin=0,n=3;e ndY=eye( n);Y(:,1)=on es( n,1);for k=3:nY(k,2:k-1)=Y(k-1,1:k-2

9、)+Y(k-1,2:k-1);End實例2尋找水仙花數(shù)一、實驗內(nèi)容數(shù)153有一個特殊的性質(zhì)，即153=1A3+5A3+3A3,那么有沒有其他這樣的數(shù) 呢？我們怎樣利用MATLAB程序設(shè)計快速尋求這些水仙花數(shù)呢？二、實驗分析把任一數(shù)各位數(shù)字的立方和求出，再求所得的各位數(shù)字的立方和，依次進行下去，總可以得到一個水仙花數(shù).如把76按這樣的運算進行下去:73635 59,5353939 79,9373931 801彳8303135 1 4,531343190,f9303730,733303370,3?7303370結(jié)果得到一個水仙花數(shù)370。三、實驗程序(尋求水仙花數(shù)經(jīng)典的程序)for a=1:9fo

10、r b=0:9;百度文庫-讓每個人平等地提升自我6for c=0:9;if aA3+bA3+cA3=a*100+b*10+c;sxh=a*100+b*10+cendendendend四、實驗結(jié)果sxh =153, 370, 371,407。實例4住房貸款的等額本息還款法計算一、實驗內(nèi)容當前在銀行放貸項目中，住房按揭貸款占相當大的比重。這一貸款促進了房 地產(chǎn)事業(yè)的迅速發(fā)展，它不但改善了居民的居住條件， 也促進了建筑行業(yè)、鋼鐵 行業(yè)和其他相關(guān)行業(yè)的發(fā)展，與此同時，也創(chuàng)造了大量的就業(yè)機會。在這中間， 銀行利率起了巨大的杠桿作用，過高的利率會使按揭購房減少，房價下跌；過低 的利率會使按揭貸款購房者數(shù)量

11、增加，當房源供不應(yīng)求時，受供求關(guān)系的影響， 房價將上漲。假設(shè)某購房者向銀行貸款的金額為p0,銀行的月利率為a,貸款期限為n個月，求每月的還款金額為r.二、實驗分析我們用等比級數(shù)求解：第一月末貸款的本息之和為p(1)p0(1 a) r，第二月末貸款的本息之和為：p(2)p(1)(1 a)rp0(1 a)2r(1 a)r，第三月末貸款的本息之和為p(3)p( 2)(1 a) r第n月末貸款的本息之和為p0(1、3a)r(1a)2r(1a) r，p(n) p(n 1)(1 a) rp0(1na)r(1a)n1(1 a) 1，考慮第n個月還清貸款，則p(n)=0，求得r pO*a(1 a)n/(1 a)n1.假設(shè)銀行貸款20萬，月利率0.465%,借款期限為10年，即為120個月，問每 月應(yīng)還金額和10年內(nèi)共計支付多少利息。三