quad_exam.m

来自「matlab7.0从入门到精通(光盘内容)」· M 代码 · 共 28 行

%quad_exam.m
%求归一化高斯函数的在区间[-1 1]上的定积分,并求得到积分过程的中间结点
y=@(x)1/sqrt(pi)*exp(-x.^2);               %归一化高斯函数
quad(y,-1,1,2e-6,1)                        %求定积分，并显示中间迭代过程
fplot(y,[-1 1],'b');                       %画出函数
hold on;
%跟踪数据（运行完上面程序后，可以在命令行中复制这些数据）
trace = [ 9    -1.0000000000    5.43160000e-001     0.1804679399;
      11    -1.0000000000    2.71580000e-001     0.0728222057;
      13    -0.7284200000    2.71580000e-001     0.1076454255;
      15    -0.4568400000    9.13680000e-001     0.4817487615;
      17    -0.4568400000    4.56840000e-001     0.2408826755;
      19    -0.4568400000    2.28420000e-001     0.1142172651;
      21    -0.2284200000    2.28420000e-001     0.1266655031;
      23     0.0000000000    4.56840000e-001     0.2408826755;
      25     0.0000000000    2.28420000e-001     0.1266655031;
      27     0.2284200000    2.28420000e-001     0.1142172651;
      29     0.4568400000    5.43160000e-001     0.1804679399;
      31     0.4568400000    2.71580000e-001     0.1076454255;
      33     0.7284200000    2.71580000e-001     0.0728222057];
x1 = trace(:,2);                          %积分过程的中间结点
y1 = y(x1);                               %中间结点的函数值
plot(x1,y1,'ro');                         %画图
xlabel('x');
ylabel('y');
legend('高斯函数','求积分过程的中间节点');