非线性规划程序

dx999

#include<iostream.h> 3 x0 R$ P$ Z) b; Z#include<math.h> double fx(int n,double p[],double x) { * x$ r5 G: r2 Z; H9 Z double z=0,y=1; ' R4 h; v; E# Q int i=0,j=0; % a- b* z5 M! O& q0 y% A for(;i<n+1;i++) { for(;j<i;j++) ) M: ?# J. x$ s y*=x; 5 Y" t$ o2 E* U- s z+=*(p+i)*y; 8 K$ \- b @3 b4 i" d* J# j } 5 P3 w" C4 h; L1 R# w return z; }

int fi(int k) 5 s" V9 I: V5 b$ F! }. Y& z% z5 Z{ if(k==0) return 1; 9 f& y0 \" c2 C* {; v) d else { 4 \# j4 W+ N! `1 Z: r! t if(k==1) ' c. F- {" T6 C# [8 F, z. U+ K n return 1; * w- Y5 o6 v9 c/ C' S1 z else { / p+ j" _% v/ o( N; r" h int i=0,i1=1,i2=1,w=0; 5 k7 B D) Z; u7 W- O' Y for(;i<k-1;i++) 0 d/ @) L) m% R: E { w=i1+i2; i1=i2; ( L4 @# P0 m. V* `! H i2=w; + I! y; s0 l0 K3 u+ v { } return w; # E. ]& ?0 F# V% C } } - O; @9 P; O6 I# [, ^}

int r(double w,double a,double b) 3 d9 z% G: L7 ^/ t{ int i=0; + i- _3 b' Y1 o, P( b) s- [ for(;fi(i)<=(b-a)/w;i++) {} return i; } void main() 3 ]( Z6 r/ c9 M" y; _1 m{ cout<<"输入函数的最高次幂"; int n; cin>>n; double *p=new double[n+1]; double a,b,g=0,j=0,c=0; cout<<"依升降幂输入x项的系数"<<endl; 2 W; u5 R3 @* x. N, S for(int i=0;i<n+1;i++) - q1 M* k# | t" G: k, I cin>>*(p+i); # p# w! {% J' S6 t) e8 z cout<<"输入误差值"<<endl; double w; / @& I2 N1 K4 v+ [. G: w cin>>w; cout<<"输入初始区间"<<endl; 1 i& s2 X' A# B0 }, Y9 ?% I cin>>a>>b; cout<<"输入辨别常数"<<endl; cin>>c; n=r(w,a,b); " r7 D; ^% {$ ^2 ^8 v2 y i=0; cout<<fi(n-i-2)<<endl; ; V1 K& f& A. o M cout<<a<<" "<<b; 3 l( z9 E4 L9 d5 V- X g=a+fi(n-2-i)*(b-a)/fi(n-i); . k1 v1 \% [+ q/ F& F) `0 j j=a+fi(n-1-i)*(b-a)/fi(n-i); cout<<g<<" "<<j; 8 b" B9 Y' a' i0 { for(i=1;i<n-1;i++) { & |3 w+ F. R! c- Z5 H5 f( J 2 u- W- l2 J. T9 k* L if(fx(n,p,g)>fx(n,p,j)) { a=g; g=j; j=a+fi(n-1-i)*(b-a)/fi(n-i); ) t7 ?, v$ R$ X" e8 r" `5 w( ?+ r } ( S) N! ]$ ~1 o/ t$ d else 4 Q9 h. w4 [- h$ Z { b=j; j=g; % H4 v9 @, I$ [% O0 ` g=a+fi(n-2-i)*(b-a)/fi(n-i); 4 y _; ]5 K$ X4 {0 N } " C! }) s7 U5 i6 S. l% D5 v1 A; ?8 Z } 8 c. |: a( }2 X( U0 R% j4 j4 F; ^, q j=g+c; ; K" H/ Q i, u( r! }! t if(fx(n,p,g)>fx(n,p,j)) 9 z9 V9 F) K" n& }2 z5 a3 Z' V a=g; else b=g; cout<<"函数值在 ("<<g<<","<<j<<")之间；"; " M$ X2 n5 \) v$ D% p} % S' ?1 x) @% Q7 T; N, U c3 U

duanzengbin

C语言吗？

Allos

这个程序可读性太差

zhanfeng

如果加上注释会容易看一些！

ilikenba

用的是什么算法？

harleych

没什么用

leonsun2002

什么算法？

zhyi

有点像牛顿法求解非线性方程组

没有解时就找到某个区间

mengfanqi

可以慢慢研究!!!

garr

不是吧，用C语言求规划？汗。lingo不是能解吗？