Sa函数在全时域积分怎么算
来源:学生作业帮 编辑:搜搜考试网作业帮 分类:综合作业 时间:2024/07/01 10:05:30
Sa函数在全时域积分怎么算
![Sa函数在全时域积分怎么算](/uploads/image/z/7409162-2-2.jpg?t=Sa%E5%87%BD%E6%95%B0%E5%9C%A8%E5%85%A8%E6%97%B6%E5%9F%9F%E7%A7%AF%E5%88%86%E6%80%8E%E4%B9%88%E7%AE%97)
理论解 3.141592653589793238.
计算机 只能算到有限区段,因为计算机有最大数值限制.
下面用辛普森公式积分.
#include
#include
#include
double fsimpf(double x)
{
if ( fabs(x) < 1e-10) return 1.0;
return sin(x) /x ;
}
double fsimp(double a,double b,double eps)
{
int n,k;
double h,t1,t2,s1,s2,ep,p,x;
n=1; h=b-a;
t1=h*(fsimpf(a)+fsimpf(b))/2.0;
s1=t1;
ep=eps+1.0;
while (ep>=eps)
{
p=0.0;
for (k=0;k
计算机 只能算到有限区段,因为计算机有最大数值限制.
下面用辛普森公式积分.
#include
#include
#include
double fsimpf(double x)
{
if ( fabs(x) < 1e-10) return 1.0;
return sin(x) /x ;
}
double fsimp(double a,double b,double eps)
{
int n,k;
double h,t1,t2,s1,s2,ep,p,x;
n=1; h=b-a;
t1=h*(fsimpf(a)+fsimpf(b))/2.0;
s1=t1;
ep=eps+1.0;
while (ep>=eps)
{
p=0.0;
for (k=0;k