1 函数的有界和无界
2 瑕点 (无界间断点)
3 无界函数的反常积分 (瑕积分)下面以两种图像相似的无界函数 (反比例函数) 为对象讲解瑕积分 . from sympy import latex, lambdify, exp ,integrate , oo ,sqrt from sympy.abc import x import numpy as np from matplotlib import pyplot as plt from IPython.display import display, Math
def show(left,right): display(Math(left+latex(right))) # 画图工具 def draw_invers_prop(a,b,f_x,y_label,int_label,fill_Flag=False):
f_x_func = lambdify(x,f_x) x_l_array = np.linspace(-4,-0.00001,200) x_r_array = np.linspace(0.00001,4,200) f_x_l_array = f_x_func(x_l_array) f_x_r_array = f_x_func(x_r_array)
fig, ax = plt.subplots() ax.plot(x_l_array, f_x_l_array, 'k', linewidth = 1) ax.plot(x_r_array, f_x_r_array, 'k', linewidth = 1) if((a<0) & fill_Flag) : ab_x_array = np.linspace(a,-0.00001,200) ab_y_array = f_x_func(ab_x_array) ax.fill_between(ab_x_array, ab_y_array, 0, color = '#EE9F20') if (int_label) : int_f_x = integrate(f_x, (x, a, 0)) if((b>0) & fill_Flag) : ab_x_array = np.linspace(0.00001,b,200) ab_y_array = f_x_func(ab_x_array) ax.fill_between(ab_x_array, ab_y_array, 0, color = '#EE9F20') if (int_label) : int_f_x = integrate(f_x, (x, 0, b))
if((a<0)&(b>0)& fill_Flag &(int_label!='')) : int_f_x = integrate(f_x, (x, a, b)) if(int_label): show(int_label,int_f_x)
ax.set_xlabel(r'$x$') ax.set_ylabel(y_label) ax.grid('--') ax.set_xlim(-4,4) ax.set_ylim(0,3)
y_label = r'$f(x)=\frac{1}{x^2}$' int_label='' f_x = x**(-2) draw_invers_prop(-4,4,f_x,y_label,int_label) 注意 :
3.1 积分区间为 且 为瑕点的反常积分
示例 3.1.1已知 是瑕点 , 求第一种无界函数在区间 上的反常积分 : 解 : 由图像可知 在区间 上连续 , 且 其存在原函数为 由微积分基本定理可得
示例 3.1.2已知 是瑕点 , 求第二种无界函数在区间 上的反常积分 : 解 : 由图像可知 在区间 上连续 其存在原函数为 由微积分基本定理可得 y_label = r'$f(x)=\frac{1}{x^2}$' int_label=' \int_{0}^{4}{\\frac{1}{x^2}}{\\rm d}x =' f_x = x**(-2) draw_invers_prop(0,4,f_x,y_label,int_label,True) 3.2 积分区间为 且 为瑕点的反常积分
示例 3.2.1已知 是瑕点 , 求第一种无界函数在区间 上的反常积分 : 解 : 由图像可知 在区间 上连续 , 且 其存在原函数为 由微积分基本定理可得
示例 3.2.2已知 是瑕点 , 求第二种无界函数在区间 上的反常积分 : 解 : 由图像可知 在区间 上连续 其存在原函数为 由微积分基本定理可得 y_label = r'$f(x)=\frac{1}{x^2}$' int_label=' \int_{-4}^{0}{\\frac{1}{x^2}}{\\rm d}x =' f_x = x**(-2) draw_invers_prop(-4,0,f_x,y_label,int_label,True) 3.3 积分区间为 和 且 为瑕点的反常积分
示例 3.3.1已知 是瑕点 , 求第一种无界函数在区间 上的反常积分 : 解 : 由定义可知
示例 3.3.2已知 是瑕点 , 求第二种无界函数在区间 上的反常积分 : 解 : 由定义可知 y_label = r'$f(x)=\frac{1}{x^2}$' int_label=' \int_{-4}^{4}{\\frac{1}{x^2}}{\\rm d}x =' f_x = x**(-2) draw_invers_prop(-4,4,f_x,y_label,int_label,True) |
|