您所在位置：主题：不知有没有对ＢＶＰ４Ｃ的ＣＯＤＥ 了解非常透彻的高手可以帮帮我？
qumarine
[专家分：0] 发布于 2009-03-05 07:13:00
大家好，不知有没有对ＢＶＰ４Ｃ的ＣＯＤＥ 了解非常透彻的高手可以帮帮我？
我发现ＣＯＤＥ中的方程和ＪＡＣＯＢＩＥＮ都不是很合理啊：
７１４行：
yip05 = (yi + yip1)/2 - (Fip1 - Fi)*(H/8);
这个对中点的近似是不是错了，应该是：yip05 = (yi + yip1)/2 - (Fip1 - Fi)*(H/４);
还有在计算ＪＡＣＯＢＩＥＮ的时候：
１０３４行：last_cols(rows,:) = -hi*dFdpar_ip05 + hi^2/12*Jip05*(dFdpar_ip1-dFdpar_i);
稍微按照它的方程计算下就知道这个不对啊！
偶非常非常地困惑，因为，我算了很多遍，怎么也想不通～可是貌似ＢＶＰ４Ｃ早已经过很多考验，大家还用得很起劲～～～可是为什么有这么明显的错误？？？而且它又是怎么行得通的捏？肯定是我哪里弄错了，所以我很期待有人能帮我解决这个问题，我又钻进死胡同了！谢谢大家！
我发现ＣＯＤＥ中的方程和ＪＡＣＯＢＩＥＮ都不是很合理啊：
７１４行：
yip05 = (yi + yip1)/2 - (Fip1 - Fi)*(H/8);
这个对中点的近似是不是错了，应该是：yip05 = (yi + yip1)/2 - (Fip1 - Fi)*(H/４);
还有在计算ＪＡＣＯＢＩＥＮ的时候：
１０３４行：last_cols(rows,:) = -hi*dFdpar_ip05 + hi^2/12*Jip05*(dFdpar_ip1-dFdpar_i);
稍微按照它的方程计算下就知道这个不对啊！
偶非常非常地困惑，因为，我算了很多遍，怎么也想不通～可是貌似ＢＶＰ４Ｃ早已经过很多考验，大家还用得很起劲～～～可是为什么有这么明显的错误？？？而且它又是怎么行得通的捏？肯定是我哪里弄错了，所以我很期待有人能帮我解决这个问题，我又钻进死胡同了！谢谢大家！