實(shí)驗(yàn)四 MATLAB數(shù)值計(jì)算.doc

實(shí)驗(yàn)四 MATLAB數(shù)值計(jì)算一、實(shí)驗(yàn)?zāi)康?.掌握數(shù)據(jù)統(tǒng)計(jì)和分析的方法2.掌握數(shù)據(jù)插值和曲線擬合的方法3.掌握求數(shù)值導(dǎo)數(shù)和數(shù)值積分的方法4.掌握代數(shù)方程數(shù)值求解的方法5.掌握常微分方程數(shù)值求解的方法二、實(shí)驗(yàn)內(nèi)容1利用matlab提供的rand函數(shù)生成30，000個(gè)符合均勻分布的隨機(jī)數(shù)，然后檢驗(yàn)隨機(jī)數(shù)的性質(zhì)：（1） 均值和標(biāo)準(zhǔn)方差；（2） 最大元素與最小元素；（3） 大于0.5的隨機(jī)數(shù)個(gè)數(shù)占總數(shù)額百分比。 x=rand(1,30000);Mean=mean(x) %均值Std=std(x,0,2)%標(biāo)準(zhǔn)方差MAX=max(x) %最大值MIN=min(x) %最小值n=sum(x0.5);percentage=n/30000 %大于0.5的隨機(jī)數(shù)個(gè)數(shù)占總數(shù)額百分比Mean = 0.5023Std = 0.2886MAX = 1.0000MIN = 2.4530e-005percentage = 0.50632按下表用3次樣條方法插值計(jì)算0900范圍內(nèi)整數(shù)點(diǎn)的正弦值和0750范圍內(nèi)整數(shù)點(diǎn)的正切值，然后用5次多項(xiàng)式擬合方法計(jì)算相同的函數(shù)值，并將兩種計(jì)算結(jié)果進(jìn)行比較。a度0153045607590Sina00.25880.50000.70710.86600.96591.0000Tana00.26790.57741.00001.73203.7320t1=0:15*pi/180:0.5*pi;t2=0:15*pi/180:75*pi/180;Sina=0,0.2588,0.5,0.7071,0.866,0.9659,1;Tana=0,0.2679,0.5774,1,1.732,3.732; x1=0:pi/180:0.5*pi; %3次樣條插值x2=0:pi/180:5/12*pi;T1=interp1(t1,Sina,x1,spline);T2=interp1(t2,Tana,x2,spline); p1=polyfit(t1,Sina,5); %多項(xiàng)式擬合disp(p1)L1=polyval(p1,t1);p2=polyfit(t2,Tana,5);disp(p2)L2=polyval(p2,t2); x=linspace(0,2*pi,100); subplot(2,1,1);plot(t1,L1,:o,x1,T1,-+,x,sin(x);title(sin(x);axis(0,0.5*pi,0,1);subplot(2,1,2);plot(t2,L2,-*,x2,T2,:.,x,tan(x);title(tan(x)axis(0,5/12*pi,0,4); sin（x）、tan（x）的5次多項(xiàng)式擬合的系數(shù)為： 0.0054 0.0075 -0.1745 0.0036 0.9994 0.0000 4.3197 -10.2021 9.2877 -3.2591 1.4027 0.00003（1）求函數(shù)在點(diǎn)的數(shù)值導(dǎo)數(shù)。x=pi*1/6,1/4,1/3,1/2;DX=pi*1/6,1/12,1/12,1/6;y=inline(sin(x).3+cos(x).3);df=diff(y(0,x)./DX df = -0.4306 -0.2575 0.2575 0.4306（2）用數(shù)值方法求積分g=inline(log(1+x)./(1+x.2);I=quadl(g,0,1)I = 0.27224. 求方程 在附近的根。f=inline(3*x+sin(x)-exp(x);fzero(fz,1.5)ans = 1.89005. 求函數(shù)在（0，1）內(nèi)的最小值。f=inline(x.3+cos(x)+x.*log(x)./exp(x);x,fval=fminbnd(f,0,1)x = 0.5223fval = 0.39746. 求微分方程的數(shù)值解，并繪制解的曲線先建立函數(shù)文件sys.mfunction dy=sys(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(3);dy(3)=-0.51*y(1)*y(2);取tf=20，再輸入命令：t0=0;tf=20;t,y=ode45(sys,t0,tf,0,1,1);plot(t,y(:,1),t,y(:,2),t,y(:,3)三、實(shí)驗(yàn)收獲與體會(huì)在學(xué)習(xí)MATLAB前，我們?cè)缫褜W(xué)