- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我想根据这个等式使用我自己的函数计算用户输入的正弦值:
sin(x) = sum_(i=0)^n (-1)^i * (x^(2 i + 1)/((2 i + 1)!))
我有这段代码,根据我的理解,我所做的与等式中所写的完全相同:
#include <stdio.h>
#include <math.h>
int faculty(int factor)
{
int result = 1;
if (factor > 0)
{
for (int i = 1; i <= factor; i++)
{
result = result * i;
}
}
else
{
result = 1;
}
return result;
}
double my_sin(double x, int n)
{
double my_sin = 0;
for (int j = 0; j < n ; j++)
{
double i = (double)j;
double faculty_j = (double)faculty(2*j+1);
my_sin = my_sin + (pow((-1.0), i) * (pow(x, (double)(2.0 * i + 1.0)) / faculty_j));
}
return my_sin;
}
int main()
{
int n = 0;
double x = 0;
printf("x=");
scanf("%lf", &x);
printf("n=");
scanf("%i", &n);
printf("sin(%i)=", (int)x);
printf("%lf\n", my_sin(x, n));
return 0;
}
但是,例如,当我使用 x = 8
和 n = 5
时,结果是 sin(8)=149
。我尝试调试代码一段时间了,但我不知道问题出在哪里,也不知道如何找出问题所在。
更新代码:
#include <stdio.h>
long int factorial(int factor)
{
long int result = 1;
if (factor > 0)
{
for (int i = 1; i <= factor; i++)
{
result = result * i;
}
}
else
{
result = 1;
}
return result;
}
double my_pow(double a, double b)
{
if (b == 0)
{
return 1;
}
double result = a;
double increment = a;
double i, j;
for (i = 1; i < b; i++)
{
for (j = 1; j < a; j++)
{
result += increment;
}
increment = result;
}
return result;
}
double my_sin(double x, int n)
{
double my_sin = 0;
for (int j = 0; j < n ; j++)
{
double i = (double)j;
double faculty_j = (double)factorial(2*i+1);
my_sin = my_sin + (my_pow((-1.0), i) * (my_pow(x, 2.0 * i + 1.0) / faculty_j));
}
return my_sin;
}
int main()
{
int n = 0;
double x = 0;
printf("x=");
scanf_s("%lf", &x);
printf("n=");
scanf_s("%i", &n);
printf("sin(%i)=", (int)x);
printf("%lf\n", my_sin(x, n));
return 0;
}
最佳答案
这里的问题不一定是你的代码。这只是您对等式的期望。我在 Wolfram Alpha 中测试了 n=5 的泰勒级数,并使用查询从该系列中减去 sin(x)
(sum (-1)^i * x^(2i+1)/(2i+1)!, i=0 to 5) - sin(x)
查看错误。
如您所见,当 x 约为 0 时,该级数的近似值非常好,但当 x 超过 3 时,误差开始变大,并且增长极快。在 x=8 时,错误的大小为 67,这显然是无用的。要获得合理的近似值,您需要使用 n=15 左右。只有五个太少了。
您在这里应该做的是利用 sin(x) = sin(x + k * 2 * PI) 的事实,其中 k 是任意整数。如果 x 是负数,您可以简单地使用 sin(x) = -sin(-x) 这一事实。我正在展示一个简单的例子来说明它是如何完成的:
double my_sin(double x, int n)
{
double my_sin = 0;
double sign = 1.0;
// sin(-x) = -sin(x)
// This does not improve accuracy. It's just to make the rest simpler
if(x<0) {
x=-x;
sign *= -1;
}
// sin(x) = sin(k*2*PI + x)
x -= 2*PI*floor(x/(2*PI));
// Continue as usual
return sign * my_sin;
}
您希望尽可能接近 x=0。为了进一步改进这一点,您可以使用 sin(x) = -sin(x+PI) 这一事实,如下所示。
以上足以将所有内容保持在 [0, PI] 区间内,但实际上我们可以通过使用一些更聪明的数学来做得更好。我们可以利用 sin(x) 关于 x=PI/2 + k*PI 对称的事实。所以 sin(PI/2 -x) = sin(PI/2 + x)。利用它将使您保持在 [0, PI/2] 区间内,当您处于该区间内时,n=2 足以使误差低于 0.005。对于 n=5,误差低于 10^-7。您可以在前面的步骤之后添加:
// sin(x) = -sin(x+PI)
if(x>PI) {
x-=PI;
sign *= -1;
}
// sin(PI/2 -x) = sin(PI/2 + x)
if(x>PI/2)
x = PI - x;
但是,您的代码中存在一个错误。请记住,j=n 的情况应该包括在内,所以更改
for (int j = 0; j < n ; j++)
到
for (int j = 0; j <= n ; j++)
另一种方法是始终在参数中添加 1 来调用函数。
只是为了清楚起见。请记住,泰勒级数要求 x
以弧度为单位。如果您使用度数,您会得到错误的结果。
只是为了完整起见,这里是完整的功能和一些其他的修复,这些修复不会影响正确性,但会影响性能和可读性:
double my_sin(double x, int n)
{
double ret = 0;
double sign = 1.0;
if(x<0) {
x=-x;
sign *= -1;
}
x -= 2*PI*floor(x/(2*PI));
if(x>PI) {
x-=PI;
sign *= -1;
}
if(x>PI/2)
x = PI - x;
size_t denominator = 1;
double numerator = x;
int s = -1;
for (int j = 0; j <= n ; j++) {
denominator *= 2*j + 1;
s *= -1;
ret += s * numerator / denominator;
numerator *= x*x;
}
return sign*ret;
}
假设您想要 3 个正确的十进制数字。然后我们可以利用分子比分母增长慢的事实并做这样的事情:
size_t denominator = 1;
double numerator = x;
double term;
int s = -1;
int j = 0;
double tolerance = 0.001;
do {
denominator *= 2*j + 1;
j++;
s *= -1;
term = numerator / denominator;
ret += s * term;
numerator *= x*x;
} while(abs(term) > tolerance);
关于c - 为什么我的 C 中的 sin 计算代码返回错误的值?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/65870518/
我已经使用 vue-cli 两个星期了,直到今天一切正常。我在本地建立这个项目。 https://drive.google.com/open?id=0BwGw1zyyKjW7S3RYWXRaX24tQ
您好,我正在尝试使用 python 库 pytesseract 从图像中提取文本。请找到代码: from PIL import Image from pytesseract import image_
我的错误 /usr/bin/ld: errno: TLS definition in /lib/libc.so.6 section .tbss mismatches non-TLS reference
我已经训练了一个模型,我正在尝试使用 predict函数但它返回以下错误。 Error in contrasts<-(*tmp*, value = contr.funs[1 + isOF[nn]])
根据Microsoft DataConnectors的信息我想通过 this ODBC driver 创建一个从 PowerBi 到 PostgreSQL 的连接器使用直接查询。我重用了 Micros
我已经为 SoundManagement 创建了一个包,其中有一个扩展 MediaPlayer 的类。我希望全局控制这个变量。这是我的代码: package soundmanagement; impo
我在Heroku上部署了一个应用程序。我正在使用免费服务。 我经常收到以下错误消息。 PG::Error: ERROR: out of memory 如果刷新浏览器,就可以了。但是随后,它又随机发生
我正在运行 LAMP 服务器,这个 .htaccess 给我一个 500 错误。其作用是过滤关键字并重定向到相应的域名。 Options +FollowSymLinks RewriteEngine
我有两个驱动器 A 和 B。使用 python 脚本,我在“A”驱动器中创建一些文件,并运行 powerscript,该脚本以 1 秒的间隔将驱动器 A 中的所有文件复制到驱动器 B。 我在 powe
下面的函数一直返回这个错误信息。我认为可能是 double_precision 字段类型导致了这种情况,我尝试使用 CAST,但要么不是这样,要么我没有做对...帮助? 这是错误: ERROR: i
这个问题已经有答案了: Syntax error due to using a reserved word as a table or column name in MySQL (1 个回答) 已关闭
我的数据库有这个小问题。 我创建了一个表“articoli”,其中包含商品的品牌、型号和价格。 每篇文章都由一个 id (ID_ARTICOLO)` 定义,它是一个自动递增字段。 好吧,现在当我尝试插
我是新来的。我目前正在 DeVry 在线学习中级 C++ 编程。我们正在使用 C++ Primer Plus 这本书,到目前为止我一直做得很好。我的老师最近向我们扔了一个曲线球。我目前的任务是这样的:
这个问题在这里已经有了答案: What is an undefined reference/unresolved external symbol error and how do I fix it?
我的网站中有一段代码有问题;此错误仅发生在 Internet Explorer 7 中。 我没有在这里发布我所有的 HTML/CSS 标记,而是发布了网站的一个版本 here . 如您所见,我在列中有
如果尝试在 USB 设备上构建 node.js 应用程序时在我的树莓派上使用 npm 时遇到一些问题。 package.json 看起来像这样: { "name" : "node-todo",
在 Python 中,您有 None单例,在某些情况下表现得很奇怪: >>> a = None >>> type(a) >>> isinstance(a,None) Traceback (most
这是我的 build.gradle (Module:app) 文件: apply plugin: 'com.android.application' android { compileSdkV
我是 android 的新手,我的项目刚才编译和运行正常,但在我尝试实现抽屉导航后,它给了我这个错误 FAILURE: Build failed with an exception. What wen
谁能解释一下?我想我正在做一些非常愚蠢的事情,并且急切地等待着启蒙。 我得到这个输出: phpversion() == 7.2.25-1+0~20191128.32+debian8~1.gbp108
我是一名优秀的程序员,十分优秀!