1、數(shù)學(xué)實驗復(fù)習(xí)數(shù)據(jù)處理命令與常用函數(shù)數(shù)學(xué)實驗課典型問題實驗報告及思考題 max 求數(shù)組最大值 min 求數(shù)組最小值 mean 求數(shù)組的平均數(shù) median 求數(shù)組的中間數(shù) sort 數(shù)組按單增排序 sum 數(shù)組元素求和 prod 數(shù)組元素求積 hist 統(tǒng)計頻率 fix 向零方向取整zeros 全零矩陣.ones 全壹矩陣.eye 單位矩陣rand 均勻分布隨機數(shù)length 數(shù)組的長度det 求行列式eig 特征值和特征向量roots 求多項式的零點.polyval 計算多項式的值.數(shù)據(jù)處理命令與常用函數(shù)abs(x) 絕對值sqrt(x) 開平方exp(x) 自然指數(shù)log(x) e為底的對數(shù)

2、rat(x) 分?jǐn)?shù)表示sin 正弦函數(shù) asin 反正弦函數(shù) cos 余弦函數(shù) acos 反余弦函數(shù) tan 正切函數(shù) atan 反正切函數(shù)ezplot 函數(shù)簡單繪圖 fplot 函數(shù)繪圖.quad 數(shù)值積分.dblquad 數(shù)值重積分inline 構(gòu)造臨時函數(shù)fzero 求函數(shù)零點.fmin 求一元函數(shù)最小值fmins 求多元函數(shù)最小值 diff 符號求導(dǎo)數(shù) int 符號求積分 mod 求余數(shù) polyfit 數(shù)據(jù)多項式擬合.數(shù)據(jù)處理命令與常用函數(shù)v0=515;alpha=45*pi/180;T=2*v0*sin(alpha)/9.8;t= T*(0:16) /16;x=v0 *t *cos

3、(alpha);y=v0 *t *sin(alpha) -0.5*9.8*t.2;plot(x,y,x,y,r*)Xmax=x(17)例1. 發(fā)射角為450的拋射曲線繪制%計算飛行時間%飛行時刻%計算航點Xmax=2.7064e+004例5 生肖問題今年是鼠年還是狗年？n=input(input n: =);S=鼠?；⑼谬埳唏R羊猴雞狗豬;k=rem(n-4,12)+1; %取余，符號與n-4相同s=S(k);s=strcat(int2str(n),年是, s,年)input n:=2006 s= 2006年是狗年Input n:=2008 s=2008年是鼠年例6 天干地支問題？year=in

4、put(input year:=);S1=甲乙丙丁戊已庚辛壬癸;S2=子丑寅卯辰巳午未申酉戍亥;k1=mod(year-4,10)+1; %取余，符號與10相同 s1=S1(k1);k2=mod(year-4,12)+1; s2=S2(k2);calendar=strcat(int2str(year),年是,s1,s2,年)例7 程序功能如下：把邊長為2以原點為中心的正方形旋轉(zhuǎn)pi/24，將其壓縮(r=0.89),重復(fù)24次并繪圖 xy=-1 -1;1 -1;1 1;-1 1;-1 -1; A=cos(pi/24) -sin(pi/24); %創(chuàng)建正交矩陣 sin(pi/24) cos(pi/

5、24); x=xy(:,1);y=xy(:,2); %提取坐標(biāo)數(shù)據(jù)line(x,y),pause(1) %畫線并暫停一秒for k=1:24 %for循環(huán)，重復(fù)24次xy=.89*xy*A; %旋轉(zhuǎn)并壓縮 x=xy(:,1);y=xy(:,2); line(x,y),pause(1)endyear=input(input year:=);n1=year/4;n2=year/100;n3=year/400;if n1=fix(n1)&n2=fix(n2) disp(是潤年)elseif n1=fix(n1)&n3=fix(n3) disp(是潤年)else disp(不是潤年)end潤年條件有二

6、能被4整除，但不能被100整除；能被4整除，又能被400整除。例8 判潤年程序% fix（x）舍去小數(shù)取整； =不等于號，&邏輯與。while (B-A)0.2 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,BA=0; B=100; va=10;vb=8;vc=60;f=1;k=0;plot(A,0,ro,B,0,go),hold on例9 摩托車問題求解的函數(shù)文件例10. 計算曲線段 f(x)=

7、exp(a x)sin(b x), 繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面體積syms a b x %定義多個符號變量f=exp(a*x)*sin(b*x);f1=subs(f,a,-0.2);f2=subs(f1,b,0.5);V=pi*int(f2*f2,x,0,2*pi)double(V)V =pi*(-125/116*exp(-4/5*pi)+125/116)ans = 3.1111%subs為f，f1中參數(shù)a，b的代換；int（）為函數(shù)積分。f=inline(exp(-0.2*x).*sin(0.5*x);%構(gòu)造臨時函數(shù)t=(0:20)*pi/10;theta=t;r=f(t);x=t*ones(si

8、ze(t);y=r*cos(theta);z=r*sin(theta);mesh(x,y,z)colormap(0 0 0)axis offview(-17,54)例11. 繞X軸旋轉(zhuǎn)的旋轉(zhuǎn)曲面圖形繪制例12.牟合方蓋問題h=2*pi/100; t=0:h:pi; r=0:0.05:1;x=r*cos(t); y=r*sin(t); %極坐標(biāo)zz=sqrt(1-x.2);meshz(x,y,zz)colormap(0 0 1)axis offN=2000;P=rand(N,3); %N*3隨機矩陣x=P(:,1);y=P(:,2);z=P(:,3);II=find(x.2+y.2=1&x.2+z.2 z = 0.05。即求逆累積函數(shù)在 x=0.95 處的值 mu=170;sigam=6;z=norminv(0.95,mu,sigam)data=mu+sigam*randn(1000,1);II=find(data=z);F=length(II)/1