一、實驗名稱
newton-cotes型求積公式
二、實驗目的
學會newton-cotes型求積公式,並應用該演算法於實際問題。
三、實驗內容
求定積分
四、實驗要求
選擇等分份數,用復化simpson求積公式求上述定積分的誤差不超過的近似值,用matlab中的內部函式int求此定積分的準確值,與利用復化simpson求積公式計算的近似值進行比較。
五、實驗程式與輸出結果
在matalab的editor視窗中輸入以下程式:
function y=***simpson(fun,a,b,n)
z1=feval (fun,a)+ feval (fun,b);m=n/2;
h=(b-a)/(2*m); x=a;
z2=0; z3=0; x2=0; x3=0;
for k=2:2:2*m
x2=x+k*h; z2= z2+2*feval (fun,x2);
endfor k=3:2:2*m
x3=x+k*h; z3= z3+4*feval (fun,x3);
endy=(z1+z2+z3)*h/3;
然後儲存為然後儲存為***simpson.m的檔案
在matalab工作視窗命令視窗中輸入:
q2 =***simpson (@fun,0,pi,1000000000)
syms x
fi=int(exp(x).*cos(x),x,0,pi);
fs= double (fi)
wq2= double (abs(fi-q2) )
執行後結果:
q2=-12.0703,fs=-12.0703,
wq2=5.2654e-08
六、實驗結果分析
利用復化simpson求積公式計算執行後其結果為q2=-12.0703,利用內部函式求解的結果為fs=-12.0703,兩者的誤差為wq2=5.
2654e-08。從中可以看出誤差結果達到了1e-8級數,而相對應的n已經取到了10億次,再增大n對結果已經沒有太大變化。可見復化simpson要得到比較準確的結果需要運算的次數比較大。
第六章Spss實驗報告
第一題資料某公司專業生產電腦顯示器,通過使用3種不同的晶元,得出顯示器的某個關鍵引數如下表,每種晶元取10個觀測值 以 0.05的顯著性水平檢驗 不同晶元的引數相同 這一假設。試運用單因方差分析完成 分析結論 晶元b的均值最高,晶元a,c相近,這些可以在圖上看出。分析 分析 由於概率p值為0.305...
第六章學習小結數值分析
第6章數值積分 學習小結 1 本章學習體會 本章的最大收穫是學習了數值積分的一些常用的方法。例如牛頓 柯特斯公式是在等距節點的情形下的插值求積公式,其簡單的如梯形公式,辛普森公式等。再有事復化求積公式是改善求積公式精度的一種有效的方法,對於復化梯形公式 復化辛普森公式很常用。而高斯求積公式是一種高精...
第六章第六章財務計畫
6.1 資金 投資比例餅圖 希吉雅食品責任 成立初期,準備籌集資金100萬元。發起人自投60萬元,申請大學生創業貸款30萬元,10萬元尋求投資,企業固定資產作投資160萬元,向銀行貸款100萬元。共計註冊資本360萬元。投資比例如圖所示 圖8 1 投資比例 創業自籌資金由創業者個人以其個人名義籌集的...