0

我正在编写一个程序，该程序从单独的信号和背景文件中读取波长和强度数据（因此每个文件都由多对波长和强度组成）。如您所见，我通过创建一个结构，然后在循环中使用 fscanf 将值分配给结构中的适当元素来做到这一点。一旦数据被读入，程序应该将其绘制在每个文件中记录的波长重叠的间隔上，即常见的波长范围。波长在存在这种重叠的地方完美对齐，并且已知以恒定的差异间隔开。因此，我辨别结构阵列的哪些元素适用的方法是确定两个文件中哪个最小波长更高，哪个最大波长更低。然后，对于具有较低最小值和较高最大值的文件，我会找到这与较高最小值/较低最大值之间的差异，然后将其除以恒定步长以确定要偏移多少元素。这是可行的，除非数学完成后，程序会返回一个完全无法解释的错误答案。

在下面的代码中，我通过计算一个元素的波长与其前一个元素的波长之间的差异，将恒定步长定义为 lambdastep。使用我的样本数据，它是 0.002，这是由 printf 确认的。但是，当我运行程序并除以 lambdastep 时，我得到了一个不正确的答案。当我运行除以 .002 的程序时，我得到了正确的答案。为什么会出现这种情况？我想不出任何解释。

#include<stdio.h>
#include<math.h>
#include<stdlib.h>
#include "plots.h"

struct spectrum{
    double lambda;
    double intensity;
};

main(){
double a=0,b=0,c=0,d=0,lambdastep,smin,smax,bmin,bmax,tmin,tmax,sintmin,bintmin,tintmin,sintmax,bintmax,tintmax,ymin,ymax;
int ns,nb,nt,i=0,sminel,smaxel,bminel,bmaxel,tminel,tmaxel;
double min(struct spectrum *a,int,int);
double max(struct spectrum *a,int,int);
FILE *Input;                                
Input = fopen("sig.dat","r");
FILE *InputII;                              
InputII = fopen("bck.dat","r");
fscanf(Input,"%d",&ns);
fscanf(InputII,"%d",&nb);
struct spectrum signal[ns];
struct spectrum background[nb];
struct spectrum *s = &signal[0];
struct spectrum *ba = &background[0];
s = malloc(ns*sizeof(struct spectrum));
ba = malloc(nb*sizeof(struct spectrum));
while( fscanf(Input,"%lf%lf",&a,&b) != EOF){
    signal[i].lambda = a;
    signal[i].intensity = b;
    i++;
}
i = 0;
while( fscanf(InputII,"%lf%lf",&c,&d) != EOF){
    background[i].lambda = c;
    background[i].intensity = d;
    i++;
}
for (i=0; i < ns ;i++){
    printf("%.7lf %.7lf\n", signal[i].lambda,signal[i].intensity);
}
printf("\n");
for (i=0; i < nb ;i++){
    printf("%.7lf %.7lf\n", background[i].lambda,background[i].intensity);
}
lambdastep = signal[1].lambda - signal[0].lambda;           //this is where I define lambdastep as the interval between two measurements
smin = signal[0].lambda;
smax = signal[ns-1].lambda;
bmin = background[0].lambda;
bmax = background[nb-1].lambda;
if (smin > bmin)
    tmin = smin;
else
    tmin = bmin;
if (smax > bmax)
    tmax = bmax;
else
    tmax = smax;
printf("%lf %lf %lf %lf %lf %lf %lf\n",lambdastep,smin,smax,bmin,bmax,tmin,tmax);   //here is where I confirm that it is .002, which is the expected value
sminel = (tmin-smin)/(lambdastep);  //sminel should be 27, but it returns 26 when lamdastep is used. it works right when .002 is directly entered , but not with lambdastep, even though i already confirmed they are exactly the same. why?
4

2 回答 2

1

sminel是一个整数，因此(tmin-smin)/lambdastep在计算结束时将转换为整数。

一个非常微小lambdastep的差异可能是获得例如 27.00001 和 26.99999 之间的差异；后者在转换为int.

尝试使用floor,ceil或round来更好地控制返回值的舍入。

于 2013-03-31T07:08:20.277 回答
0

它几乎肯定与浮点计算固有的不精确性有关。尝试打印出lambdastep许多有效数字——我打赌你会发现它的确切值比你想象的要大一些。

使用我的样本数据，它是.002，由 确认printf。

尝试打印出来(lambdastep == .002)。

于 2013-03-31T07:07:00.020 回答