| 查看: 627 | 回复: 8 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】一个生物类的公式如何反编译（悬赏20）
|
|||
|
本人是学习生物类的，数学很糟糕，悬赏20个金币，求助一公式反编译。 X=363.1888+y* 23.25622-y的平方* 1.745972+y的立方*0.0552203， 要求输入x，得到y值。 如果有人能用vb实现同样的功能，同样有奖。 谢谢个位高手出手相助。 [ Last edited by csfn on 2008-11-30 at 11:53 ] |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
madonion
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 589.4
- 帖子: 307
- 在线: 26.3小时
- 虫号: 251164
- 注册: 2006-05-14
- 性别: GG
- 专业: 化学反应工程
★ ★ ★
suntao1982(金币+3,VIP+0):辛苦了，欢迎常来！
suntao1982(金币+3,VIP+0):辛苦了，欢迎常来！
|
找到一个现成的，不过是c语言的 用卡丹公式解一元三次方程 #include #include #include //////////////////函数--用卡丹公式解一元三次方程///////////////// void fun(double a,double b,double c,double d, double *real_y1,double *real_y2,double *real_y3, double *imag_y1,double *imag_y2,double *imag_y3) { double p,q,r,u,v,g,h,fai; p=(3.0*a*c-b*b)/(3*a*a); q=(2.0*pow(b,3.0)-9*a*b*c+27.0*a*a*d)/(27.0*pow(a,3.0)); r=b/(3.0*a); h=pow(q/2.0,2.0)+pow(p/3.0,3.0); g=sqrt(h); if(h>=0) { if(-q/2.0+g<0) u=-pow(fabs(-q/2.0+g),1.0/3.0); else u=pow((-q/2.0+g),1.0/3.0); if(-q/2.0-g<0) v=-pow(fabs(-q/2.0-g),1.0/3.0); else v=-pow((-q/2.0-g),1.0/3.0); if(h==0) { *real_y1=u+v-r; *imag_y1=0; *real_y2=-(u+v)/2-r; *imag_y2=0; *real_y3=-(u+v)/2-r; *imag_y3=0; } else { *real_y1=u+v-r; *imag_y1=0; *real_y2=-(u+v)/2; *imag_y2=sqrt(3.0)*(u-v)/2; *real_y3=-(u+v)/2; *imag_y3=-sqrt(3.0)*(u-v)/2; } } else { fai=acos((-q/2)/(sqrt(pow(fabs(p),3)/27))); *real_y1=2*sqrt(fabs(p)/3.0)*cos(fai/3.0)-r; *real_y2=-2*sqrt(fabs(p)/3.0)*cos((fai+3.1415926)/3.0)-r; *real_y3=-2*sqrt(fabs(p)/3.0)*cos((fai-3.1415926)/3.0)-r; *imag_y1=0; *imag_y2=0; *imag_y3=0; } } //////////////////////////////////主函数//////////////////////////////// void main() { double a,b,c,d; double real_x1,real_x2,real_x3; double *preal_x1=&real_x1; double *preal_x2=&real_x2; double *preal_x3=&real_x3; double imag_x1,imag_x2,imag_x3; double *pimag_x1=&imag_x1; double *pimag_x2=&imag_x2; double *pimag_x3=&imag_x3; cout<<"请输入方程的系数a,b,c,d:"<<"\n"< cout<<"系数a="; cin>>a; cout< cout<<"系数b="; cin>>b; cout< cout<<"系数c="; cin>>c; cout< cout<<"系数d="; cin>>d; cout< fun(a,b,c,d,preal_x1,preal_x2,preal_x3,pimag_x1,pimag_x2,pimag_x3); cout<<" "<<" "<<"根的实部"<<" "<<"根的虚部"<<"\n"< printf("x1 %.5f %.5f\n\n",real_x1,imag_x1); printf("x2 %.5f %.5f\n\n",real_x2,imag_x2); printf("x3 %.5f %.5f\n\n",real_x3,imag_x3); } |
6楼2008-07-19 20:05:46
595304082syl
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1975.5
- 散金: 1
- 帖子: 100
- 在线: 40.6小时
- 虫号: 374345
- 注册: 2007-05-17
- 性别: GG
- 专业: 药学
2楼2008-07-18 10:47:48
spkeey
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1244.9
- 散金: 3
- 帖子: 308
- 在线: 14小时
- 虫号: 320174
- 注册: 2007-03-10
- 性别: GG
- 专业: 理论和计算化学
★ ★
suntao1982(金币+2,VIP+0):xiexie
suntao1982(金币+2,VIP+0):xiexie
|
三次方程可能有一个值对应几个自变量的情况（不符合函数定义），只能说在某个区间上的反函数。 你的导数 x=23.25622-1.745972y1+0.0552203*y2 （二次可能经过X轴）算了下没有经过0。是单调的，不用分开讨论。。算你运气。。。。 可以求反函数。 数值上可以这样处理的： y=a+bx+cx2+dx3 等效于求 0=(a-y)+bx+cx2+dx3 的根。。 =========== 此有公式。 [ Last edited by spkeey on 2008-7-19 at 03:37 ] |
3楼2008-07-18 22:13:33
madonion
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 589.4
- 帖子: 307
- 在线: 26.3小时
- 虫号: 251164
- 注册: 2006-05-14
- 性别: GG
- 专业: 化学反应工程
★ ★ ★
suntao1982(金币+3,VIP+0):辛苦了，欢迎常来！
suntao1982(金币+3,VIP+0):辛苦了，欢迎常来！
|
卡丹公式 确定一般的三次方程的根的公式. 如果用现在的数学语言和符号,卡丹公式的结论可以借助于下面这样一种最基本的设想得出。 假如给我们一个一般的三次方程： ax3+3bx2+3cx+d=0 （1） 如果令 x=y-b/a 我们就把方程（1）推导成 y3+3py+2q=0 （2） 其中 p=c/a-b2/a2，2q=2b3/a3-3bc/a2+d/a 。 借助于等式 y=u-p/u 引入新变量u 。把这个表达式带入（2），得到： （u3）2+2qu3-p3=0 （3） 由此得 u3=-q±√（q2+p3）， 于是 y=3√（-q±√（q2+p3））-p/3√（-q±√（q2+p3）） 。 =3√（-q+√（q2+p3））+3√（-q-√（q2+p3）） 。 （最后这个等式里的两个立方根的积等于-p 。） 这就是著名的卡丹公式。 如果再由y转到x，那么，就能得到一个确定一般的三次方程的根的公式。 |
4楼2008-07-19 20:01:26
