PAGE
PAGE 1
数学实验作业三
题目:P93. 2).d;6);11)
日期:2003-
【实验目的】:
1、用MATLAB软件掌握梯形公式、辛普森公式和蒙特卡罗方法计算数值积分。
2、通过实例学习用数值积分和数值微分解决实际问题。
【实验内容】:
二:选择一些函数用梯形(trapz),辛普森(quad)和蒙特卡罗三种方法计算积分。改变步长(对梯形),改变精度要求(对辛普森),改变随机点数目(对蒙特卡罗),进行比较、分析。如下函数供选择参考:
a.
b.
c.
d.
【MATLAB源程序】
% 作业三2).b
clear;
a = 0; b = 2;
fun = inline(exp(3*x) .* sin(2*x));
format long e
% 以a1 为划分区间数做5次trapz积分
a1 = [10, 20, 40, 80, 160];
for i = 1:length(a1)
x1 = linspace(a,b,a1(i)+1);
answer1(i) = trapz(x1,feval(fun,x1));
end
disp( Trapz :)
disp( Steplength Answer)
disp([(b-a)./a1, answer1])
% 以a2为精度要求做5次quad积分
a2
原创力文档

文档评论(0)