带参数的分段的定积分---求帮助

赵明琳

本程序是从实验数据，对公式f（x，b）=2/pi*x*f(x)./(x.^2-b.^2)求定积分，得出g（b），最后画出g（b）的图形，f(x)已经通过二次多项式拟合分段得出，，所以积分变为分段的定积分了，b是参数，对x积分，写了一个程序，结果最后得出的图像为空白-------求帮忙---紧急，，，另x，b取值范围一样，怀疑问题出在奇异点，但程序没有报错。


x=0.3:0.01:0.8;
y1=[3.6093
3.7405
3.8674
3.9972
4.1174
4.2437
4.3728
4.4925
4.6145
4.7411
4.8606
4.98
5.1051
5.23
5.3473
5.4694
5.59
5.7162
5.8387
5.96
6.0778
6.198
6.3237
6.4475
6.5615
6.6869
6.8079
6.9245
7.0366
7.1456
7.2564
7.3687
7.4808
7.5887
7.6925
7.7938
7.8996
8.0079
8.1174
8.2216
8.3148
8.405
8.4991
8.5732
8.607
8.6195
8.6239
8.6127
8.5747
8.5148
8.4494
];
y=y1';
c=0;
for k=0:3:48
    x2=[0.3+k*0.01 0.31+k*0.01 0.32+k*0.01];
    y2=[y(k+1) y(k+2) y(k+3)];
    p=polyfit(x2,y2,2);
    syms a b
    f=b*(p(1)*b.^2+p(2)*b+p(3))./(b.^2-a.^2);
    n=int(f,b,0.3+k*0.01,0.32+k*0.01);
   c=c+n;
end
a=linspace(0.3,0.8);
ezplot(a,c)

赵明琳

自己顶一个，求解决啊，初学者伤不起啊