Matlab上機(jī)實(shí)驗(yàn)報(bào)告(共4頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上Matlab上機(jī)實(shí)驗(yàn)報(bào)告姓名：王梅 學(xué)號(hào)： 上機(jī)時(shí)間：2013-5-141用符號(hào)工具箱計(jì)算：01(x4*(1-x)4)/(1+x2)dx這個(gè)積分能用數(shù)值積分計(jì)算嗎？如果能，應(yīng)該如何做？方法：符號(hào)計(jì)算和數(shù)值計(jì)算。 程序流程圖：新建一個(gè)指令窗口編輯指令將指令復(fù)制到操作窗口得出結(jié)果程序清單及運(yùn)行結(jié)果：>> syms xa=int(x4*(1-x)4)/(1+x2),0,1) a = 22/7 - pi >> f=inline('x.4.*(1-x).4./(1+x.*x)','x')f = Inline functio

2、n: f(x) = x.4.*(1-x).4./(1+x.*x)>> quad(f,0,1)ans =0.00132.誤差函數(shù)erf(x)的定義是： Erf(x)=2/(pi0.5)0xe-x2dx用指令quad計(jì)算當(dāng)x=0.1,0.21.0時(shí)相應(yīng)的erf函數(shù)值。經(jīng)結(jié)果與MATLAB的內(nèi)部函數(shù)erf(x)的結(jié)果進(jìn)行比較。方法：數(shù)值計(jì)算。程序流程圖：新建一個(gè)指令窗口編輯指令將指令復(fù)制到操作窗口得出結(jié)果程序清單及運(yùn)行結(jié)果：>> x=0.1:0.1:1.0f=inline('2./pi.0.5.*exp(-x.*x)','x')erf(x)fo

3、r x=0.1:0.1:1.0quad(f,0,x)endx = 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000f = Inline function: f(x) = 2./pi.0.5.*exp(-x.*x)ans = 0.1125 0.2227 0.3286 0.4284 0.5205 0.6039 0.6778 0.7421 0.7969 0.8427ans = 0.1125ans = 0.2227ans = 0.3286ans = 0.4284ans = 0.5205ans = 0.6039ans = 0.6778ans = 0.7421ans = 0.7969ans =0.84273.求微分方程xy+y-ex=0在初值條件y1=2e下的特接，并畫(huà)出解函數(shù)的圖形（解析解與數(shù)值解）程序流程圖：新建一個(gè)指令窗口編輯指令將指令復(fù)制到操作窗口得出結(jié)果程序按及運(yùn)行結(jié)果：y=dsolve('x*Dy+y-exp(x)=0','y(1)=2*exp(1