计算sinh的泰勒级数
该函数使用Taylor系列中的以下开发来计算sinh(x)
的值:
我想计算sinh(3)= 10.01787的值,但函数输出为9.我也得到这个警告:
1> main.c(24):警告C4244:'函数':从'double'转换为'int',可能丢失数据
这是我的代码:
int fattoriale(int n)
{
int risultato = 1;
if (n == 0)
{
return 1;
}
for (int i = 1; i < n + 1; i++)
{
risultato = risultato * i;
}
return risultato;
}
int esponenziale(int base, int esponente)
{
int risultato = 1;
for (int i = 0; i < esponente; i++)
{
risultato = risultato * base;
}
return risultato;
}
double seno_iperbolico(double x)
{
double risultato = 0, check = -1;
for (int n = 0; check != risultato; n++)
{
check = risultato;
risultato = risultato + (((esponenziale(x, ((2 * n) + 1))) / (fattoriale((2 * n) + 1))));
}
return risultato;
}
int main(void)
{
double numero = 1;
double risultato = seno_iperbolico(numero);
}
请帮我修复这个程序。
编译器警告你这种数据丢失实际上非常棒。
你看,当你调用这个:
esponenziale(x, ((2 * n) + 1))
由于您将double
(即x
)转换为int
,因此基本上会失去准确性。 这是因为esponenziale
的签名是int esponenziale(int base, int esponente)
。
将它改为double esponenziale(double base, int esponente)
, risultato
应该是double
,因为你正在从函数中返回它并使用/对其执行数学运算。
记住用int
分割一个double
会给你一个double
回馈。
编辑:根据ringø的评论,并看到它如何真正解决你的问题,你应该设置double fattoriale(int n)
并在double risultato = 1;
。
由于许多术语都是分数量,因此您正在失去精确度。 使用int
将打破小数部分。 根据需要用double
类型替换您的int
类型。
你的阶乘函数会溢出n
令人惊讶的小值。 对于16位int
, n
的最大值是7,对于32位是12,对于64位,它是19.溢出有signed
整型的行为是未定义的。 如果你的编译器支持它,你可以使用unsigned long long
或uint128_t
。 这会给你多一点时间。 但是,无论如何,如果你想转换成double
,那么你应该从一开始就使用double
。 请注意,一个IEEE764浮点加倍将在171处达到无穷大!
可以肯定的是,对于任何x
值, sinh
的Maclaurin展开的收敛半径是无限的。 所以x
任何值都可以工作,尽管收敛速度可能会很慢。 见http://math.cmu.edu/~bkell/21122-2011f/sinh-maclaurin.pdf。