資源描述:
《實(shí)驗(yàn)四 MATLAB數(shù)值計(jì)算》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、實(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
2、=sum(x>0.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ì)算0~900范圍內(nèi)整數(shù)點(diǎn)的正弦值和0~750范圍內(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
3、.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=
4、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次
5、多項(xiàng)式擬合的系數(shù)為:0.00540.0075-0.17450.00360.99940.00004.3197-10.20219.2877-3.25911.40270.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]))./DXdf=-0.4306-0.25750.25750.4306(2)用數(shù)值方法求積分g=inline('log(1+x)./(1+x.^2)');I=qua
6、dl(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.mfunctiondy=sys(t,y)dy=zeros(3,1);dy(1)=y(2)*y(3);dy(2)=-y(1)*y(
7、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é)習(xí)了多門數(shù)學(xué)課如數(shù)值分析、數(shù)學(xué)分析、最優(yōu)化和常微分等,各門課程里都有相關(guān)的問(wèn)題、方程、函數(shù)和解決這些的方法,實(shí)驗(yàn)四真正地讓我們學(xué)以致用,解決問(wèn)題的方法都可以通過(guò)MATLAB實(shí)現(xiàn)來(lái)求解問(wèn)題。