应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Windows编程 » 【求助】求助：vb编程中用牛顿迭代解三次方程为什么只得到一个根？【已完成】
162/2返回列表上一页12
查看: 2808  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 2 个 程序强帖 ，点击这里进行查看

phychemlxd

金虫 (小有名气)
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
resonant(金币+5):欢迎交流，感谢回帖:-) 2010-04-21 18:34
liangzidou(金币+5):谢谢你的热情帮助，因为你回帖时金币已经赠送完了，就给你追加5个，币不多，略表谢意。 2010-04-22 18:00
wangen994(程序强帖+1):感谢你的参与，呵呵，欢迎常来 2010-05-09 21:49:52
看到你说热力学，来个C/C++版本的，如果你能知道压力对压缩因子/体积/密度的导数（二阶导数）信息的话，用迭代法求也不错，毕竟有时候三角函数计算也很耗时

#include
#include
#include
using namespace std;

typedef unsigned char Boolean;
#define False ((Boolean)'\0')
#define True  ((Boolean)'\1')

typedef int Index;
typedef double Real;

// The last parameter denotes complex roots is allowed or not
// X^3 + Coeff[0]*X^2 + Coeff[1]*X + Coeff[2] = 0
Index CubicRoot(Real *Coeff, Real *Root, Boolean Complex=True)
   {
   Real p=Coeff[0], q=Coeff[1], r=Coeff[2];
   Real Den_3, A, A3, B, B2, D, sqrt_D, M, N;

   Den_3 = 0.33333333333333333333333333;
   Real Den_27 = Den_3*Den_3*Den_3;

   A = Den_3 * (3.0*q - p*p);
   A3 = Den_27*A*A*A;
   B = 2.0*Den_27*p*p*p - Den_3*p*q + r;
   B2 = 0.25*B*B;

   D = A3 + B2;

   if (D >= 0.0)
      {
      sqrt_D = sqrt(D);
      M = pow((-0.5*B + sqrt_D), Den_3);
      N = pow((-0.5*B - sqrt_D), Den_3);
      if (D > 0.0) // 1 root
         {
         Root[0] = M + N;
         nRoot = 1;
         if (bComplex)
            {
            Root[1] = -0.5*(M + N);
            Root[2] = 0.5*sqrt(3.0)*(M - N);
            nRoot = 3;
            }
         }
      else // 2 roots
         {
         Root[0] = M + N;
         Root[1] = Root[2] = -0.5*Root[0];
         nRoot = 2;
         }
      }
   else // 3 real distinct roots
      {
      Real deg = sqrt(-B2/A3);
      Real PI_2 = 1.5707963267949; // 90C
      Real PI_3_2 = 2.0943951023932; // 120C
      Real Theta;
      if (B >= 0.0)
         {
         Theta = (PI_2 + atan(deg/sqrt(1.0 - deg*deg)))*Den_3;
         }
      else
         {
         Theta = atan(sqrt(1.0 - deg*deg)/deg)*Den_3;
         }
      Real factor = 2.0*sqrt(-Den_3*A);

      for (Index i=0; i<3; i++)
         {
         Root = factor*cos(Theta + PI_3_2*i);
         }
      nRoot = 3;
      }

   for (Index i=0; i<3; i++)
      {
      Root -= Den_3*p;
      }

   return nRoot;
   }

int main(int argc, char* argv[])
   {
   string Holding;

   Real a[3] = {-0.9964370251388, 0.0345649889308, -0.000135939438932387};
   Real root[3];
   Index nRoot = CubicRoot(a, root);
   cout<<"Root 1: "<    cout<<"Root 2: "<    cout<<"Root 3: "<    cin>>Holding;

   return 0;
   }
赞一下(3人) 回复此楼
修.齐.治.平
11楼2010-04-21 18:00:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangzidou

银虫 (小有名气)
是我用vb编程中要用到解立方型方程的这段程序，呵呵
赞一下 回复此楼
12楼2010-04-22 17:50:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangzidou

银虫 (小有名气)
引用回帖:
Originally posted by phychemlxd at 2010-04-21 18:00:29:
看到你说热力学，来个C/C++版本的，如果你能知道压力对压缩因子/体积/密度的导数（二阶导数）信息的话，用迭代法求也不错，毕竟有时候三角函数计算也很耗时

#include
#include

呵呵，谢谢你了，我没有系统的学过编程语言，不会用C语言，就是vb也是一知半解，不过还是很感谢你提供的解立方型方程的程序。
赞一下 回复此楼
13楼2010-04-22 17:55:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

phychemlxd

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
resonant(金币+1):感谢回帖交流：-） 2010-04-23 10:55
有热力学的问题我们可以探讨，如果是做热力学模型计算的话，我在方法上还是能提供一点点小小的帮忙的，呵呵。
赞一下(2人) 回复此楼
修.齐.治.平
14楼2010-04-23 09:16:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liangzidou

银虫 (小有名气)
★
jjdg(金币+1):欢迎常来 2010-04-25 01:40
引用回帖:
Originally posted by phychemlxd at 2010-04-23 09:16:17:
有热力学的问题我们可以探讨，如果是做热力学模型计算的话，我在方法上还是能提供一点点小小的帮忙的，呵呵。

嗯，好的，呵呵
回复此楼
15楼2010-04-24 23:23:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhimengfan

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
牛顿迭代法本身就是不完备的方法，不能保证找到所有根。
赞一下 回复此楼
16楼2012-07-11 22:01:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liangzidou 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 0XLacIJUOj8D 2026-08-24 6/300 2026-08-25 00:11 by BZKMTicpDhFj
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭